A positivity preserving method for simulation of steam injection for NAPL site remediation

Numerical modelling of steam injection methods for cleanup of non-aqueous phase liquid (NAPL) contamination of groundwater requires consideration of multiphase, multicomponent convective and dispersive transport. Standard techniques do not ensure that the solution of the discrete equations has positive mole fractions, for finite mesh sizes. Negative mole fractions may cause the simulation to abort due to failure of the Newton iteration. A method for alleviating this problem is described. This method ensures that the mole fractions are positive, and results in an error that is the same size as the usual finite element discretization error. Example computations are presented for cartesian and axisymmetric two-dimensional geometries.     

Estimation of Scalar Source/Sink Distributions in Plant Canopies Using Lagrangian Dispersion Analysis: Corrections for Atmospheric Stability and Comparison with a Multilayer Canopy Model

Source/sink distributions of heat, water vapour andCO₂ within a rice canopy were inferred using aninverse Lagrangian dispersion analysis and measuredmean profiles of temperature, specific humidity andCO₂ mixing ratio. Monin–Obukhov similarity theorywas used to account for the effects of atmosphericstability on _w(z), the standard deviation ofvertical velocity and _L(z), the Lagrangian timescale of the turbulence. Classical surface layer scaling was applied in the inertial sublayer (z > z_ruf)using the similarity parameter = (z - d)/L, where z is height above ground, d is the zero plane displacementheight for momentum, L is the Obukhov length,and z_ruf 2.3h_c, where h_c iscanopy height. A single length scale h_c, was usedfor the stability parameter 3 = h_c/L in the height range 0.25 < z/h_c < 2.5. This choice is justified by mixing layer theory, which shows that within the roughness sublayer there is one dominant turbulence length scaledetermined by the degree of inflection in the windprofile at the canopy top. In the absence of theoretical or experimental evidence for guidance,standard Monin–Obukhov similarity functions, with = h_c/L, were used to calculate the stabilitydependence of _w(z) and _L(z) in the roughness sublayer. For z/h_c < 0.25 the turbulence length and time scales are influenced by the presence of the lowersurface, and stability effects are minimal. With theseassumptions there was excellent agreement between eddycovariance flux measurements and deductions from theinverse Lagrangian analysis. Stability correctionswere particularly necessary for night time fluxes whenthe atmosphere was stably stratified.The inverse Lagrangian analysis provides a useful toolfor testing and refining multilayer canopy models usedto predict radiation absorption, energy partitioningand CO₂ exchanges within the canopy and at thesoil surface. Comparison of model predictions withsource strengths deduced from the inverse analysisgave good results. Observed discrepancies may be dueto incorrect specification of the turbulent timescales and vertical velocity fluctuations close to theground. Further investigation of turbulencecharacteristics within plant canopies is required toresolve these issues.     

半拉格朗日、半隐式欧拉方程组大气数值模式研究

介绍了一个具有较高时间积分效率的三维弹性大气数值模式。其中,控制方程为能够描述大气非静力及可压缩性的欧拉方程组,时间积分采用了高效率的“半拉格朗日、半隐式”方案。通过引用一个简化的“云物理过程”参数化方案,针对大气中深厚湿对流过程的若干问题进行了数值试验,结果表明,模式能够对该过程中的一些基本和复杂现象进行有效地模拟。试验结果还表明,动力学框架具有很好的稳定性,能够实现高效率的时间积分,它还具有较好的频散特性,保证了空间计算的精度,从而完好地描述了风暴的形态。     

Experimental determination of saltating glass particle dispersion in a turbulent boundary layer

A horizontal saltation layer of glass particles in air is investigated experimentally over a flat bed and also over a triangular ridge in a wind tunnel. Particle concentrations are measured by light scattering diffusion (LSD) and digital image processing, and velocities using particle image velocimetry (PIV). All the statistical moments of the particle concentration are determined such as mean concentration, root mean square concentration fluctuations, skewness and flatness coefficients. Over the flat bed, it is confirmed that the mean concentration decreases exponentially with height, the mean dispersion height being a significant length scale. It is shown that the concentration distribution follows quite well a lognormal distribution. Over the ridge, measurements were made at the top of the ridge and in the cavity region and are compared with measurements without the ridge. On the hill crest, particles are retarded, the saltation layer decreases in thickness and concentration is increased. Downwind of the ridge, particle flow behaves like a jet, in particular no particle return flow is observed. Copyright © 2006 John Wiley & Sons, Ltd.     

东沟特大型斑岩钼矿床的物化探找矿效果

对河南东沟钼矿床的区域物化探成果进行整理、研究,运用物化探资料循序渐进地揭示出钼矿床的空间分布规律。东沟特大型钼矿床的发现,不仅为我国的钼矿产业今后几十年的稳步发展将产生深远的影响,而且为在东秦岭—大别山钼成矿带进一步运用物化探作为充分必要手段,为寻找更多的钼矿资源提供了难得的成功经验。     

Determination of Porosity and Permeability out of Dispersion Analyses of Borehole Wall Displacements

Waves carried in the borehole can be measured by the use of controlled sources inside the borehole and receivers located in the same borehole as well. This article gives the theoretical background and develops a new method to determine rock-physical parameters out of such controlled measurements of waves carried inside the borehole. Theoretical dispersion curves of the group velocities are matched to dispersion analyses of registrations obtained in the borehole. This inversion process intends to determine the in situ porosity and permeability of the rock.     

Turbulence Scale Dependence of the Richardson Constant in Lagrangian Stochastic Models

We investigate the relative dispersion properties of the well-mixed class of Lagrangian stochastic models. Dimensional analysis shows that, given a model in the class, its properties depend solely on a non-dimensional parameter, which measures the relative weight of Lagrangian-to-Eulerian scales. This parameter is formulated in terms of Kolmogorov constants, and model properties are then studied by modifying its value in a range that contains the experimental variability. Large variations are found for the quantity, g* = 2gC₀^{− 1}, where g is the Richardson constant.     

CFD simulation of airflow over a regular array of cubes. Part II: analysis of spatial average properties

In the first part of this study, results of a computational fluid dynamics simulation over an array of cubes have been validated against a set of wind-tunnel measurements. In Part II, such numerical results are used to investigate spatially-averaged properties of the flow and passive tracer dispersion that are of interest for high resolution urban mesoscale modelling (e.g. non resolved obstacle approaches). The results show that vertical profiles of mean horizontal wind are linear within the canopy and logarithmic above. The drag coefficient, derived from the numerical results using the classical formula for the drag force, is height dependent (it decreases with height). However, a modification of the formula is proposed (accounting for subgrid velocity scales) that makes the drag coefficient constant with height. Results also show that the dispersive fluxes are similar in magnitude to the turbulent fluxes, and that they play a very important role within the canopy. Vertical profiles of turbulent length scales (to be used in k–l closure schemes, where k is the turbulent kinetic energy and l a turbulent length scale) are also derived. Finally the distribution of the values around the mean over the reference volumes are analysed for wind and tracer concentrations.     

北京市朝阳区沙子营垃圾污染物在潜水含水层中迁移规律的现场模拟实验

选择北京平原区水文地质、环境地质等方面都比较典型的沙子营垃圾堆放场,建立了其水文地质模型。在充分收集资料、分析得出部分计算所需参数后,采用现场弥散实验、勘查取样测试等方法,求得了该含水层的弥散系数等参数;用二维非稳定流溶质运移方程对污染物在此含水层中的迁移扩散规律、速度和污染范围等进行了模拟计算;采用现场钻探、取样测试分析等方法,评价了该含水层的实际污染状况。实验模拟计算和现场调查结果表明:污染物在潜水含水层中的运移规律遵循二维非稳定流场中的溶质运移方程,污染物在潜水含水层中的运移速度约为86.25m/a,迁移扩散主要发生在地下水流向上,侧向扩散宽度极小,是地下水流向上的1/17。     

Migration and dispersion of trace elements in the rock–soil–plant system in areas underlain by black shales and slates of the Okchon Zone, Korea

The Okchon black shale in Korea provides a typical example of natural geological materials enriched with potentially toxic elements. The Chung-Joo, Duk-Pyung, Geum-Kwan, I-Won, Bo-Eun and Chu-Bu areas are underlain by these black shales and slates of the Guryongsan Formation or the Changri Formation, which are parts of the Okchon Group in the central part of the southern Korean Peninsula. In order to investigate the enrichment levels and dispersion patterns of potentially toxic elements in the rock–soil–plant system, environmental geochemical surveys were undertaken in the above six study areas in the Okchon Zone. After appropriate preparation, rock and soil samples were analyzed for potentially toxic elements by instrumental neutron activation analysis (INAA) and inductively coupled plasma atomic emission spectrometry (ICP-AES), and plant samples by atomic absorption spectrometry (AAS). In particular, Ba, Cd, Mo, V and U in Okchon black shales are highly enriched, and their mean concentrations are significantly higher than those in black slates. These elements are geochemically associated, and might be enriched simultaneously. The highest mean concentrations of 42.0 μg g⁻¹ As, 2100 μg g⁻¹ Ba, 10.9 μg g⁻¹ Cd, 213 μg g⁻¹ Mo, 83 μg g⁻¹ U, 938 μg g⁻¹ V and 394 μg g⁻¹ Zn are found in black shales from the Duk-Pyung area. Mean concentrations of As, Mo and U in soils overlying black shales occurring in the Duk-Pyung area (30 μg g⁻¹ As, 24 μg g⁻¹ Mo and 50 μg g⁻¹ U) and Chu-Bu area (39 μg g⁻¹ As, 15 μg g⁻¹ Mo and 27 μg g⁻¹ U) are higher than the permissible level. Enrichment index values of the six study areas decrease in the order of Duk-Pyung > Chu-Bu > Bo-Eun > Chung-Joo > Geum-Kwan = I-Won areas. Relationships between trace element concentrations in soils and plants are significantly correlated, and the biological absorption coefficients (BAC) in plants are in the order of Cd > Zn = Cu > Pb, which suggests that Cd is more bioavailable to plants than the other elements. Cadmium concentrations in plant species decrease in the order of chinese cabbage > red pepper > soybean = sesame > rice stalk > corn > rice grain. From the result of sequential extraction analysis of soils, relatively high proportions of Cu, Pb and Zn are present as residual fractions, and that of Cd as non-residual fractions. Cadmium occurs predominantly as exchangeable/water-acid soluble phase in soils, and this is in agreement with the findings of high Cd concentrations in plants.