任意阶次磁多极子场力线方程的解析表示

本文证明任意阶次磁多极子场力线方程总可以用有限形式解析地表示,并对偶极族和四极族磁多极子场力线的位形给予注记。     

Kepler Equation solver

Kepler's Equation is solved over the entire range of elliptic motion by a fifth-order refinement of the solution of a cubic equation. This method is not iterative, and requires only four transcendental function evaluations: a square root, a cube root, and two trigonometric functions. The maximum relative error of the algorithm is less than one part in 10¹⁸, exceeding the capability of double-precision computer arithmetic. Roundoff errors in double-precision implementation of the algorithm are addressed, and procedures to avoid them are developed.     

流体阻力对软管空间平衡形态的影响

深海采矿系统中的一个关键环节是连接中继舱和集矿机的软管。软管空间形态的非线性分析是一个复杂的问题。本文采用荷载增量法与牛顿一拉夫逊相结合的方法研究了软管在稳定水流作用下位移的反应,利用非线性有限元理论与Morison流体阻力计算公式分析了流体阻力对软管空间形态的影响。该计算模型可为软管输送性能的设计提供理论参考依据。     

北京市房价与地价的动态关系——基于结构方程模型的实证分析

城市房价与地价之间的关系错综复杂,不仅受多种因素的交织影响,相互之间也存在动态关系。研究房价与地价关系的传统方法,如Granger因果检验和回归分析等,无法刻画房价与地价之间多维的网络状关系,相比之下,结构方程模型能同时处理多个内生潜变量,且不受观测指标共线性的影响,为刻画地价与房价的交互作用提供了新的工具。从住房与土地市场的供需传导机制出发,推导出房价与地价的结构模型,以北京市为例,运用2003-2013年居住用地价格和2014年在售楼盘价格,与北京市GIS电子地图相匹配,提取商服中心可达性、公共交通可达性、道路可达性、商服繁华度、设施便利性等解释变量,构建地价与房价结构方程模型,分析二者之间的结构关系。     

Comparison of two different equations of state for application of carbon dioxide sequestration

We suggest that different equations of state (EOS) algorithms can and frequently will provide very different predictions of CO₂ migration following injection for sequestration. Rather than carry out an exhaustive examination of all EOS algorithms available, we elected to evaluate this general hypothesis by making detailed comparisons of simulation results of two very common EOS algorithms. We simulated and compared CO₂ migration patterns using two fundamentally different EOS algorithms – Modified Redlich-Kwong EOS (MRKEOS) and Span and Wagner EOS (SWEOS). In general, the predictions of thermophysical properties for both algorithms are close, except for a contrast in the predicted fugacity coefficient of CO₂, which subsequently propagates to a contrast in predicted solubility in water/brine. Typically, MRKEOS underestimates solubility of CO₂ compared to both SWEOS and experimental solubility data. In simulations of CO₂ migration, dissolution rates of separate-phase CO₂ predicted from the two EOS algorithms were significantly different, even for small contrasts in predicted fluid properties from EOS algorithms, resulting in markedly different migration patterns.     

An easily implemented agro-hydrological procedure with dynamic root simulation for water transfer in the crop–soil system: Validation and application

Models for water transfer in the crop–soil system are key components of agro-hydrological models for irrigation, fertilizer and pesticide practices. Many of the hydrological models for water transfer in the crop–soil system are either too approximate due to oversimplified algorithms or employ complex numerical schemes. In this paper we developed a simple and sufficiently accurate algorithm which can be easily adopted in agro-hydrological models for the simulation of water dynamics. We used a dual crop coefficient approach proposed by the FAO for estimating potential evaporation and transpiration, and a dynamic model for calculating relative root length distribution on a daily basis. In a small time step of 0.001 d, we implemented algorithms separately for actual evaporation, root water uptake and soil water content redistribution by decoupling these processes. The Richards equation describing soil water movement was solved using an integration strategy over the soil layers instead of complex numerical schemes. This drastically simplified the procedures of modeling soil water and led to much shorter computer codes. The validity of the proposed model was tested against data from field experiments on two contrasting soils cropped with wheat. Good agreement was achieved between measurement and simulation of soil water content in various depths collected at intervals during crop growth. This indicates that the model is satisfactory in simulating water transfer in the crop–soil system, and therefore can reliably be adopted in agro-hydrological models. Finally we demonstrated how the developed model could be used to study the effect of changes in the environment such as lowering the groundwater table caused by the construction of a motorway on crop transpiration.     

MATHEMATICAL MODEL OF OVERLAND FLOW AND MECHANISM OF SOIL CONSERVATION FOR FORESTED STEEP HILLSLOPE

ＭＡＴＨＥＭＡＴＩＣＡＬＭＯＤＥＬＯＦＯＶＥＲＬＡＮＤＦＬＯＷＡＮＤＭＥＣＨＡＮＩＳＭＯＦＳＯＩＬＣＯＮＳＥＲＶＡＴＩＯＮＦＯＲＦＯＲＥＳＴＥＤＳＴＥＥＰＨＩＬＬＳＬＯＰＥ（１）ＡｎａｌｙｔｉｃａｌＳｏｌｕｔｉｏｎｔｏｔｈｅＯｖｅｒｌａｎｄＦｌｏｗｏ...     

任意阶次磁多极子场力线方程的解析表示

本文证明任意阶次磁多极子场力线方程总可以用有限形式解析地表示，并对偶极族和四极族磁多极子场力线的位形给予注记。     

用辐射传输方程改进算法反演LandSAT8海表温度的效果检验和适用条件分析

为检验辐射传输方程改进算法推广用于反演LandSAT8 海表温度(SST)的可行性及适用条件,本研究在修订算法部分参数的基础上,分别反演出算法改进前后的LandSAT8 SST并进行比较。MODIS SST验证表明改进算法对大气透过率偏差和SST偏差均有明显改善,除2019-08-21图像外,其他3景图像SST偏差在0.5℃左右;浮标SST验证表明改进算法在近岸海域对SST偏差改善效果同样显著,其平均偏差bias和均方根误差rmse分别减少到0.14、0.18℃,基本可以忽略不计;进一步研究发现,改进算法SST反演精度与研究海域大气透过率分布均匀程度呈显著的正相关关系,R²[bias(SST)]=0.920 7,R²[rmse(SST)]=0.934 0。使用改进算法遥感监测电厂温排水：嵩屿电厂表层海水温升羽流不明显;而后石电厂表层温升现象最显著,温升幅度和扩散影响范围最大,在其排水口附近存在稳定的高温水体(温升＞3℃);晋江电厂表层温升羽流呈扇形分布,流轴短,主要集中在排水口附近。     

P–V Equations of State and the relative stabilities of serpentine varieties

Serpentines are hydrous phyllosilicates which form by hydration of Mg–Fe minerals. The reasons for the occurrence of the structural varieties lizardite and chrysotile, with respect to the variety antigorite, stable at high pressure, are not yet fully elucidated, and their relative stability fields are not quantitatively defined. In order to increase the database of thermodynamic properties of serpentines, the P–V Equations of State (EoS) of lizardite and chrysotile were determined at ambient temperature up to 10 GPa, by in situ synchrotron X-ray diffraction in a diamond-anvil cell. Neither amorphization nor hysteresis was observed during compression and decompression, and no phase transition was resolved in lizardite. In chrysotile, a reversible change in compression mechanism, possibly due to an unresolved phase transition, occurs above 5 GPa. Both varieties exhibit strong anisotropic compression, with the c axis three times more compressible than the others. Fits to ambient temperature Birch–Murnaghan EoS gave for lizardite V ₀=180.92(3) Å³, K ₀ = 71.0(19) GPa and K′ ₀=3.2(6), and for chrysotile up to 5 GPa, V ₀ = 730.57(31) Å³ and K ₀ = 62.8(24) GPa (K′ ₀ fixed to 4). Compared to the structural variety antigorite is stable at high pressure (HP) (Hilairet et al. 2006), the c axis is more compressible in these varieties, whereas the a and b axes are less compressible. These differences are attributed to the less anisotropic distribution of stiff covalent bonds in the corrugated structure of antigorite. The three varieties have almost identical bulk compressibility curves. Thus the compressibility has negligible influence on the relative stability fields of the serpentine varieties. They are dominated by first-order thermodynamic properties, which stabilizes antigorite at high temperature with respect to lizardite, and by out-of-equilibrium phenomena for metastable chrysotile (Evans 2004).