Accuracy of numerical simulations of contaminant transport in heterogeneous aquifers: A comparative study

This work deals with a comparison of different numerical schemes for the simulation of contaminant transport in heterogeneous porous media. The numerical methods under consideration are Galerkin finite element (GFE), finite volume (FV), and mixed hybrid finite element (MHFE). Concerning the GFE we use linear and quadratic finite elements with and without upwind stabilization. Besides the classical MHFE a new and an upwind scheme are tested. We consider higher order finite volume schemes as well as two time discretization methods: backward Euler (BE) and the second order backward differentiation formula BDF (2). It is well known that numerical (or artificial) diffusion may cause large errors. Moreover, when the Péclet number is large, a numerical code without some stabilising techniques produces oscillating solutions. Upwind schemes increase the stability but show more numerical diffusion. In this paper we quantify the numerical diffusion for the different discretization schemes and its dependency on the Péclet number. We consider an academic example and a realistic simulation of solute transport in heterogeneous aquifer. In the latter case, the stochastic estimates used as reference were obtained with global random walk (GRW) simulations, free of numerical diffusion. The results presented can be used by researchers to test their numerical schemes and stabilization techniques for simulation of contaminant transport in groundwater.     

黑龙江省暴雪时空分布特征与发生风险研究

暴雪是黑龙江省常见的灾害性天气之一,基于1970～2006年黑龙江省78个市（县）的逐日降雪资料,对暴雪初、终日分布规律,暴雪的时间、空间分布特征进行分析,采用信息扩散理论计算了各市（县）发生不同暴雪日数的概率风险估计值,并结合GIS技术进行风险区划．研究表明：黑龙江省暴雪初始日期和终止日期37a来变化不大,初始日出现...     

Noble Gas Diffusion Mechanism in Lunar Soil Simulant Grains:Results from ~4He~+ Implantation and Extraction Experiments

Experiments on ion implantation were performed in order to better characterize diffusion of noble gases in lunar soil.4He+ at 50 keV with 5×1016 ions/cm2 was implanted into lunar simulants and crystal ilmenite.Helium in the samples was released by stepwise heating experiments.Based on the data,we calculated the helium diffusion coefficient and activation energy.Lunar simulants display similar 4He release patterns in curve shape as lunar soil,but release temperatures are a little lower.This is probably a con...     

深海环境中溢油输移扩散的初步数值模拟

通过对溢油在深海环境中的输移过程及行为特点的分析,初步建立基于拉格朗日积分法的深海溢油模型.该模型除了能够模拟油气混合物在真实深海环境中的共同输移与分离输移扩散过程,还考虑了石油溶解、气体溶解、天然气水合物形成与分解等行为变化对溢油运动轨迹的影响.应用该模型初步数值模拟了一次实际深海溢油试验,结果表明溢油在水下的空间分...     

太湖水污染扩散插件式模拟系统

插件式太湖水污染扩散模拟系统,基于三维水体污染扩散模型,结合GIS技术和三维可视化技术,再现了水体污染物运移扩散的时空分布特征与动态演进过程.为了保证系统的可扩展性,需要一种松散的耦合机制.插件技术可以有效地降低系统功能模块间的耦合,把插件技术引入GIS的二次开发可以赋予系统动态可扩展性.本文首先,介绍插件式系统的一般...     

双重介质中核素迁移模型的解析解及其应用

给出了具有单裂隙的双重介质中，裂隙一端有分时注入源情形下的迁移模型的解析解及一些简化条件下的解析解，并利用解析解研究了核素在双重介质中的迁移问题。     

橄榄石中H的结合机制及扩散行为

名义上无水矿物(nominally anhydrous minerals,NAMs)中以点缺陷形式存在的结构水,因其对于矿物物理化学性质的显著影响而受到越来越广泛的关注。橄榄石是上地幔中含量最丰富的矿物,水在橄榄石中的存在形式、扩散机制和速率对上地幔的流变学和电导性质有着重要的影响。因此,对于橄榄石中H的结合机制及其扩散机制和速率的了解有助于理解地球深部的水循环,也有助于构建上地幔流变结构以及解释电导测量结果。本文总结了近几十年来学界对于H在橄榄石中存在的缺陷类型与OH红外光谱吸收峰之间的对应关系,以及H在矿物晶格中的扩散机制、扩散速率等重要问题的研究成果,并探讨了现有研究中依然存在的问题。     

Flux and dispersion of gases from the “Drachenschlund” hydrothermal vent at 8°18' S, 13°30' W on the Mid-Atlantic Ridge

Hydrothermal emission of mantle helium appears to be directly related to magma production rate, but other processes can generate methane and hydrogen on mid-ocean ridges. In an on-going effort to characterize these processes in the South Atlantic, the flux and distribution of these gases were investigated in the vicinity of a powerful black smoker recently discovered at 8°17.9' S, 13°30.4' W. The vent lies on the shoulder of an oblique offset in the Mid-Atlantic Ridge and discharges high concentrations of methane and hydrogen. Measurements during expeditions in 2004 and 2006 show that the ratio of CH₄ to ³He in the neutrally buoyant plume is quite high, 4 × 10⁸. The CTD stations were accompanied by velocity measurements with lowered acoustic Doppler current profilers (LADCP), and from these data we estimate the methane transport to have been 0.5 mol s^− 1 in a WSW-trending plume that seems to develop during the ebb tidal phase. This transport is an order of magnitude greater than the source of CH₄ calculated from its concentration in the vent fluid and the rise height of the plume. From this range of methane fluxes, the source of ³He is estimated to be between 0.14 and 1.2 nmol s^− 1. In either case, the ³He source is significantly lower than expected from the spreading rate of the Mid-Atlantic Ridge. From the inventory of methane in the rift valley adjacent to the vent, it appears that the average specific rate of oxidation is 2.6 to 23 yr^− 1, corresponding to a turnover time between 140 and 16 days. Vertical profiles of methane in the surrounding region often exhibited Gaussian-like distributions, and the variances appear to increase with distance from the vent. Using a Gaussian plume model, we obtained a range of vertical eddy diffusivities between 0.009 and 0.08 m²m² s^− 1. These high values may be due to tidally driven internal waves across the promontory on which the vent is located.     

井水温度微动态及其形成机制

介绍了我国地震地下水温度观测概况及已获得的微动态现象, 即水温固体潮、 水温同震效应、 水温前兆异常等, 并对其形成机理做了理论分析。 提出水温微动态形成机理研究应从二个地温梯度、 二个对流、 一个传导、 一个扩散等热系统与热过程去考虑。     

Water and nutrient fluxes off Northwest Africa

A historical data set is used to describe the coastal transition zone off Northwest Africa during spring 1973 and fall 1975, from 17° to 26°N, with special emphasis on the interaction between subtropical (North Atlantic Central Waters) and tropical (South Atlantic Central Waters) gyres. The near-surface geostrophic circulation, relative to 300 m, is quite complex. Major features are a large cyclonic pattern north of Cape Blanc (21°N) and offshore flow at the Cape Verde front. The large cyclone occurs in the region of most intense winds, and resembles a large meander of the baroclinic southward upwelling jet. The Cape Verde frontal system displays substantial interleaving that may partly originate as mesoscale features at the coastal upwelling front. Property–property diagrams show that the front is an effective barrier to all properties except temperature. The analysis of the Turner angle suggests that the frontal system is characterized by large heat horizontal diffusion as a result of intense double diffusion, which results in the smoothing of the temperature horizontal gradients. Nine cross-shore sections are used to calculate along-shore geostrophic water-mass and nutrient transports and to infer exchanges between the coastal transition zone and the deep ocean (import: deep ocean to transition zone; export: transition zone to deep ocean). These exchanges compare well with mean wind-induced transports and actual geostrophic cross-shore transport estimates. The region is divided into three areas: southern (18–21°N), central (21–23.5°N), and northern (23.5–26°N). In the northern area geostrophic import is roughly compensated with wind-induced export during both seasons. In the central area geostrophic import is greater than wind-induced export during spring, resulting in net import of both water (0.8 Sv) and nitrate (14 kmol s⁻¹), but during fall both factors again roughly cancel. In the southern area geostrophy and wind join to export water and nutrients during both seasons, they increase from 0.6 Sv and 3 kmol s⁻¹ during fall to 2.9 Sv and 53 kmol s⁻¹ during spring.