Modeling Variable Density Flow and Solute Transport in Porous Medium: 2. Re‐Evaluation of the Salt Dome Flow Problem

Case 5, Level 1 of the international HYDROCOIN groundwater flow modeling project is an example of idealized flow over a salt dome. The groundwater flow is strongly coupled to solute transport since density variations in this example are large (20%).Several independent teams simulated this problem using different models. Results obtained by different codes can be contradictory. We develop a new numerical model based on the mixed hybrid finite elements approximation for flow, which provides a good approximation of the velocity, and the discontinuous finite elements approximation to solve the advection equation, which gives a good approximation of concentration even when the dispersion tensor is very small. We use the new numerical model to simulate the salt dome flow problem.In this paper we study the effect of molecular diffusion and we compare linear and nonlinear dispersion equations. We show the importance of the discretization of the boundary condition on the extent of recirculation and the final salt distribution. We study also the salt dome flow problem with a more realistic dispersion (very small dispersion tensor). Our results are different to prior works with regard to the magnitude of recirculation and the final concentration distribution. In all cases, we obtain recirculation in the lower part of the domain, even for only dispersive fluxes at the boundary. When the dispersion tensor becomes very small, the magnitude of recirculation is small. Swept forward displacement could be reproduced by using finite difference method to compute the dispersive fluxes instead of mixed hybrid finite elements.     

Modeling Variable Density Flow and Solute Transport in Porous Medium: 1. Numerical Model and Verification

A new numerical model for the resolution of density coupled flow and transport in porous media is presented. The model is based on the mixed hybrid finite elements (MHFE) and discontinuous finite elements (DFE) methods. MHFE is used to solve the flow equation and the dispersive part of the transport equation. This method is more accurate in the calculation of velocities and ensures continuity of fluxes from one element to the adjacent one. DFE is used to solve the convective part of the transport equation. Combined with a slope limiting procedure, it avoids numerical instabilities and creates a very limited numerical dispersion, even for high grid Peclet number.Flow and transport equations are coupled by a standard iterative scheme. Residual based criterion is used to stop the iterations. Simulations of an unstable equilibrium show the effects of the criteria used to stop the iterations and the stopping criterion in the solver. The effects are more important for finer grids than for coarser grids.The numerical model is verified by the simulation of standard benchmarks: the Henry and the Elder test cases. A good agreement is found between the revised semianalytical Henry solution and the numerical solution. The Elder test case was also studied. The simulations were similar to those presented in previous works but with significantly less unknowns (i.e. coarser grids). These results show the efficiency of the used numerical schemes.     

气相色谱-质谱法测定地下水中多氯萘

建立了地下水中1-氯萘、2-氯萘、1,4-二氯萘、1,2,3,4-四氯萘、1,3,5,7-四氯萘、1,2,3,5,7-五氯萘、1,2,3,5,6,7-六氯萘、1,2,3,4,5,6,7-七氯萘和八氯萘9种多氯萘(PCNs)的气相色谱-质谱(GC-MS)分析方法。对比研究了液液萃取(LLE)和固相萃取(SPE)萃取地下水中PCNs的提取效率,优选二氯甲烷-液液萃取为PCNs检测的前处理方法。在优化条件下,9种PCNs的线性范围为5~100μg/L,各组分的相关系数(r)大于0.995,方法检出限(S/N=3)为4.21~7.41 ng/L,地下水的平均加标回收率为70.7%~112%,相对标准偏差(RSD,n=5)均小于9.9%。该方法已用于地下水样中多氯萘的检测。     

Analysis of Concentration as Sampled in

The aim of this paper is to analyze the statistical properties of solute concentration in natural aquifers as sampled in observation wells, having a small diameter in comparison with the characteristic size of the heterogeneity in hydraulic properties. The analysis, in Langragian framework, takes advantage of the reverse formulation, where, instead of considering the destination of the injected particles, the origin of the particle being sampled is sought. In the case of small values of the log-conductivity variance _Y², it allows the derivation of an analytical expression for concentration mean, variance and pdf, while for aquifer characterized by high value in _Y², a numerical analysis based on a Monte Carlo approach using a reverse scheme is developed and applied for values of _Y² up to 2. In this case, the use of a Beta function to fit the concentration pdf proves valid for practical applications. The comparison between the numerical and the analytical results defines the range of validity of the analytical ones. The relative role of large-scale dispersion processes and pore-scale effects is analyzed in terms of global variance in order to point out limits and accuracy of the Eulerian scheme in comparison with the Lagrangian one.     

Stochastic Analysis of Macrodispersion in a Semi‐Confined Aquifer

In this paper, the macroscopic dispersion resulting from one and twodimensional flows through a semiconfined aquifer with spatially variable hydraulic conductivity K which is represented by a stationary (statistically homogeneous) random process is analyzed using the spectral representation technique. Stochastic fluctuation equations of the steady flow and solute transport are solved to construct the macroscopic dispersive flux and evaluate the resulting macrodispersivity tensor in terms of the leakage factor and input covariances describing the hydraulic conductivity in a semiconfined aquifer bounded by a leaky layer above and an impervious stratum below. The macrodispersivity tensor is studied using some convenient forms of the log hydraulic conductivity process. The sensitivity of the resulting macrodispersivity to the input covariances is discussed along with the influence of the leakage factor for both one and twodimensional flows. It is found that the longitudinal macrodispersivities are increased due to the presence of leakage, while the transverse macrodispersivities are reduced due to leakage.     

The Kinematic Flow Structure for the Gvirtzman–Gorelick In Situ VOC Remediation System

Potential theory and Stokes' stream function techniques are used to investigate the flow structure around the recirculation system developed by Gvirtzman and Gorelick (1992, 1993), which consists of an extraction well and a gallery (trench) for the recharge of treated water to the aquifer. Analytical formulas are derived for the drawdown, velocity, and stream function for a model in which the extraction well is modeled as a uniformly distributed line sink and the gallery is modeled as a uniformly distributed ring source. Travel times are reported for water particles traveling along the streamlines containing 50 and 90% of the flow for various degrees of well penetration and various radii of the ring source. The travel times along the streamline resulting in the shortest travel time (not necessarily the shortest path) are also reported for various degrees of well penetration and various radii of the ring source. The method completely eliminates the use of numerical finite-difference or finite-element methods and can be used for optimization of technological parameters of this remediation system.     

基于复杂网络的并联式项目群结构脆弱性评估

运用复杂网络研究并联式项目群的结构脆弱性。将项目群中的任务用节点表示、任务之间的依赖关系用边表示,依据任务的工期确定点权、任务之间的依赖强度确定边权,则并联式项目群被抽象为一个有向加权网络。在分析并联式项目群网络拓扑特征的基础上,构建并联式项目群脆弱性评估模型,并对一个并联式项目群进行脆弱性评估。仿真结果表明该方法可以对并联式项目群脆弱性进行定量化评估、能够有效识别结构中的脆弱环节。本研究为保障项目群顺利实施、抵御风险提供了新的思路。     

Modelling groundwater changes due to fluctuating dam discharge

In this contribution, two numerical methods are used to predict the free surface changes in a sand bar due to fluctuations in river stage. One is a fixed-mesh, finite-element seepage formulation including Biot's consolidation theory, and the other is a boundary element method solution of the Laplace equation. Both models give overall predictions that are in good agreement with field data recorded at an instrumented sand bar in the Colorado River subjected to stage fluctuations from operation of the Glen Canyon Dam. The boundary element method appears to offer significant advantage in data preparation and computational times over the finite-element method for the problem studied in this paper.     

低频电磁脉冲对复杂电大目标耦合效应

针对复杂电大目标电磁脉冲耦合效应缺乏定量研究的问题,基于时域有限差分方法(FDTD)建立了一套电磁脉冲耦合效应建模与评估分析流程,开展了低频电磁脉冲激励下的仿真试验,获取和分析了特性数据,进而提出了一种表面加载透明导电材料的低频电磁脉冲激励防护方法。试验结果表明:基于功率密度的易损性分析具有良好的工程适用性,而且防护方法对低频电磁脉冲激励具有较好的屏蔽效能。     

Detection and analysis of cyclotrimethylenetrinitramine (RDX) in environmental samples by surface‐enhanced Raman spectroscopy

Techniques for rapid and sensitive detection of energetics such as cyclotrimethylenetrinitramine (RDX) are needed both for environmental and security screening applications. Here we report the use of surface‐enhanced Raman scattering (SERS) spectroscopy to detect traces of RDX with good sensitivity and reproducibility. Using gold (Au) nanoparticles (∼90–100 nm in diameter) as SERS substrates, RDX was detectable at concentrations as low as 0.15 mg/l in a contaminated groundwater sample. This detection limit is about two orders of magnitude lower than those reported previously using SERS techniques. A surface enhancement factor of ∼6 × 10⁴ was obtained. This research further demonstrates the potential for using SERS as a rapid, in situ field screening tool for energetics detection when coupled with a portable Raman spectrometer. Copyright © 2010 John Wiley & Sons, Ltd.