A Finite Volume Scheme for the Transport of Radionucleides in Porous Media

The COUPLEX1 Test case (Bourgeat et al., 2003) is devoted to the comparison of numerical schemes on a convection–diffusion–reaction problem. We first show that the results of the simulation can be mainly predicted by a simple analysis of the data. A finite volume scheme, with three different treatments of the convective term, is then shown to deliver accurate and stable results under a low computational cost.     

川西亚高山三个森林群落的湿林冠蒸发速率

于3个假说和林冠上方2m处的气象变量，采用Penman-Monteith组合模型估算了一个生长季节内川西亚高山林区分别以云杉(SF)、冷杉(FF)和白桦(BF)为优势树种的3个林分的湿林冠蒸发速率(Er)。研究结果表明，SF、FF和BF的湿林冠蒸发量(E)分别为44．51mm、88．51mm和57．8mm，分别占总降雨量的9．2％、16．6％和10．2％。与SF和BF相比，FF具有最高的月平均Er和蒸发比例。SF、FF和BF的平均Er分别为0．097mm／h(变化范围：0．028-0．487mm／h)、0．242mm／h(变化范围：0．068～0．711mm／h)和0．149mm／h(0．060～0．576mm／h)。最高和最低的月平均Er分别在6月(SF、FF和BF分别为0．120mm／h、0．317mm／h和0．169mm／h)和10月(SF、FF和BF分别为0．083mm／h、0．187mm／h和0．101mm／h)。8：00至16：00期间的平均点Er显著高于0：00至8：00以及16：00至0：00期间的平均Er。Er显著的日变化和月变化主要归因于林冠上方的太阳辐射、空气温度和相对湿度的变化。     

信息时代地图技术美的研究与实践

通过分析现代地图表现形式及存在的问题，探讨了信息时代地图制图技术美的含义和地图技术美的表现形式。并在技术美创新理念指导下，对利用计算机技术和图形图像处理技术来改进地图符号设计，提高地图表现力进行了一些理论探讨和技术尝试．同时也得到了一些有益的启示。     

用遗传算法反演地壳的变密度模型

用青藏高原布格重力异常资料 ,结合亚东 格尔木地学断面得到的Moho面结构 ,反演了青藏高原中部地壳内密度随深度变化的指数分布函数 ,试验得出地壳表层与地幔的密度差约为 - 0 .92 6 g·cm-3 ,地壳的平均密度值为 2 .75 0 g·cm-3 。     

基于分层的面状地图符号设计

根据面状符号的结构特征和实现机制 ,提出了基于分层的面状地图符号设计方法 ,把面状符号分为若干层 ,每层采用一种填充方式 ,并在基于VC 的Cartosymbol软件中实现。结果表明 ,该方法既便于地理现象的表达 ,又便于计算机实现。     

时间动态地图模型

就时间动态地图模型进行了研究，认为动态地图由于表达空间维数的增加，使得地图作为客观现实世界的概念模型更加完善，并且作为物质模型，动态地图的可视化应与空间认知的计算机表达紧密融合。在制作时间动态地图时，不同的视觉变量具有不同的动态表现能力，根据时态现象的特征进行选取。     

Nonlinear free thermal convection in a spherical shell:A variable viscosity model

Introduction The velocity field of surface plate motion can be split into a poloidal and a toroidal parts.At the Earth′s surface,the toroidal component is manifested by the existence of transform faults,and the poloidal component by the presence of convergence and divergence,i.e.spreading and subduc-tion zones.They have coupled each other and completely depicted the characteristics of plate tec-tonic motions.The mechanism of poloidal field has been studied fairly clearly which is related to …     

Mathematical model for hydraulically aided electrokinetic remediation of aquifer and removal of nonanionic copper

One of the most cost-effective in situ technologies for soil and groundwater (i.e., aquifer) remediation is electrokinetic remediation. In electrokinetic remediation, electromigration due to electric field is combined with hydromigration due to hydraulic flow by purge water to remove pollutants from aquifers through the pore water. This study aims at investigating theoretically the role of electromigration (as active movement) of pollutants and the role of hydromigration (as passive movement) of pollutants in electrokinetic remediation, and making it clear that the control variables for electrokinetic remediation are the applied voltage and the hydraulic flow rate. These aims are pursued by construction of a mathematical model based on physico-chemical considerations and by model simulations of the electrokinetic remediation applied to the virtual aquifer polluted by heavy metals of copper sulfate. According to numerical simulations with the model: (1) heavy metal (nonanionic copper) is removed from the upstream anode region and accumulated in the downstream cathode region; (2) to carry away the heavy metal outside the aquifer (global removal), hydromigration by purge water flow is essential; and (3) electromigration contributes mainly to the redistribution of heavy metals within the aquifer (local removal and local accumulation).     

Refined modeling and movement characteristics analyses of irregularly shaped particles

Irregularly shaped (IRS) particles widely exist in many engineering and industrial fields. The macro physical and mechanical properties of the particle system are governed by the interaction between the particles in the system. The interaction between IRS particles is more complicated because of their complex geometric shape with extremely irregular and co‐existed concave and convex surfaces. These particles may interlock each other, making the sliding and friction of IRS particles more complex than that of particles with regular shape. In order to study the interaction of IRS particles more efficiently, a refined method of constructing discrete element model based on computed tomography scanning of IRS particles is proposed. Three parameters were introduced to control the accuracy and the number of packing spheres. Subsequently, the inertia tensor of the IRS particle model was optimized. Finally, laboratory and numerical open bottom cylinder tests were carried out to verify the refined modeling method. The influence of particle shape, particle position, and mesoscopic friction coefficient on the interaction of particles was also simulated. It is noteworthy that with the increase of mesoscopic friction coefficient, the fluidity of IRS particle assembly decreases, and intermittent limit equilibrium state may appear. Copyright © 2014 John Wiley & Sons, Ltd.     

Discrete element method modeling of inherently anisotropic rocks under uniaxial compression loading

A new numerical approach is proposed in this study to model the mechanical behaviors of inherently anisotropic rocks in which the rock matrix is represented as bonded particle model, and the intrinsic anisotropy is imposed by replacing any parallel bonds dipping within a certain angle range with smooth‐joint contacts. A series of numerical models with β = 0°, 15°, 30°, 45°, 60°, 75°, and 90° are constructed and tested (β is defined as the angle between the normal of weak layers and the maximum principal stress direction). The effect of smooth‐joint parameters on the uniaxial compression strength and Young's modulus is investigated systematically. The simulation results reveal that the normal strength of smooth‐joint mainly affects the behaviors at high anisotropy angles (β > 45°), while the shear strength plays an important role at medium anisotropy angles (30°–75°). The normal stiffness controls the mechanical behaviors at low anisotropy angles. The angle range of parallel bonds being replaced plays an important role on defining the degree of anisotropy. Step‐by‐step procedures for the calibration of micro parameters are recommended. The numerical model is calibrated to reproduce the behaviors of different anisotropic rocks. Detailed analyses are conducted to investigate the brittle failure process by looking at stress‐strain behaviors, increment of micro cracks, initiation and propagation of fractures. Most of these responses agree well with previous experimental findings and can provide new insights into the micro mechanisms related to the anisotropic deformation and failure behaviors. The numerical approach is then applied to simulate the stress‐induced borehole breakouts in anisotropic rock formations at reduced scale. The effect of rock anisotropy and stress anisotropy can be captured. Copyright © 2015 John Wiley & Sons, Ltd.