Lattice Boltzmann Models for Flow and Transport in Saturated Karst

Flow and transport simulation in karst aquifers remains a significant challenge for the ground water modeling community. Darcy's law–based models cannot simulate the inertial flows characteristic of many karst aquifers. Eddies in these flows can strongly affect solute transport. The simple two-region conduit/matrix paradigm is inadequate for many purposes because it considers only a capacitance rather than a physical domain. Relatively new lattice Boltzmann methods (LBMs) are capable of solving inertial flows and associated solute transport in geometrically complex domains involving karst conduits and heterogeneous matrix rock. LBMs for flow and transport in heterogeneous porous media, which are needed to make the models applicable to large-scale problems, are still under development. Here we explore aspects of these future LBMs, present simple examples illustrating some of the processes that can be simulated, and compare the results with available analytical solutions. Simulations are contrived to mimic simple capacitance-based two-region models involving conduit (mobile) and matrix (immobile) regions and are compared against the analytical solution. There is a high correlation between LBM simulations and the analytical solution for two different mobile region fractions. In more realistic conduit/matrix simulation, the breakthrough curve showed classic features and the two-region model fit slightly better than the advection-dispersion equation (ADE). An LBM-based anisotropic dispersion solver is applied to simulate breakthrough curves from a heterogeneous porous medium, which fit the ADE solution. Finally, breakthrough from a karst-like system consisting of a conduit with inertial regime flow in a heterogeneous aquifer is compared with the advection-dispersion and two-region analytical solutions.     

Source Screening Module for Contaminant Transport Analysis Through Vadose and Saturated Zones

At complex sites there may be many potential sources of contaminants within the vadose zone. Screening‐level analyses are useful to identify which potential source areas should be the focus of detailed investigation and analysis. A source screening module (SSM) has been developed to support preliminary evaluation of the threat posed by vadose zone waste sites on groundwater quality. This tool implements analytical solutions to simulate contaminant transport through the unsaturated and saturated zones to predict time‐varying concentrations at potential groundwater receptors. The SSM integrates several transport processes in a single simulation that is implemented within a user‐friendly, Microsoft Excel? ‐ based interface.     

Numerical Simulation of Solute Transport in Saturated Porous Media with Bounded Domains

This work presents the first attempt to develop unconditionally stable, implicit finite difference solutions of one-sided spatial fractional advection-dispersion equation (s-FADE) by imposing the nonzero Dirichlet boundary condition (ND BC) or the nonzero fractional Robin boundary condition (NFR BC) at inlet boundary and the zero fractional Neumann boundary condition (ZFN BC) at outlet boundary. The results of the numerical studies performed using artificial solute transport parameters demonstrated that the numerical solution with the NFR BC as the inlet boundary produced much more realistic concentration values. The numerical solution with the NFR BC at the inlet boundary was capable of correctly describing the Fickian and non-Fickian behaviors of the solute transport at different α values, and it had the relatively same accuracy at different numbers of the spatial nodes. Also, the practical application of the numerical solution with the NFR BC as the inlet boundary was investigated by conducting tracer experiments in homogeneous and heterogeneous soil columns. According to the obtained results, this numerical solution described well solute transport in the homogenous and heterogeneous soils. The α values of the homogeneous and heterogeneous soils were obtained in the ranges of 1.849–1.999 and 1.248–1.570, respectively, which were in excellent agreement with the physical properties of the soils. In a nutshell, the numerical solution of the s-FADE with the NFR BC as the inlet boundary can be successfully applied to describe the solute transport in the homogeneous and heterogeneous soils with bounded spatial domains.     

饱和土体瞬态响应有限元分析

给出基于Biot多孔介质理论分析饱和土体在动载荷作用下瞬态响应的有限元公式,数值计算部分采用本文有限元法分别计算一维饱和土柱在两种不同类型动载荷作用下的瞬态响应,并将数值计算结果与文献中的解析解进行比较,二者结果十分吻合,从而验证本文方法的可行性。     

Simulation of Vapor Transport Through the Unsaturated Zone — Interpretation of Soil-Gas Surveys

Soil-gas surveys are becoming more widely accepted as a tool for the preliminary determination of the extent of soil and ground water contamination by volatile organic compounds (VOCs). The interpretation of the results of published soil-gas surveys has been necessarily limited and qualitative due to a lack of adequate models. There has been considerable effort in the recent past to characterize the transport and fate of pesticides in soil. However, the behavior of pesticides generally differ substantially from those of VOCs.
This paper presents a computer model developed to simulate the diffusive transport of VOC vapor through unsaturated soils using a two-dimensional, finite-difference, solution to Fick's second law of diffusion. An effective diffusion coefficient that incorporates the effects of tortuosity, moisture content, and soil organic carbon content is computed. Although the model has not been validated due to the unavailability of adequate field or laboratory data, nevertheless, sensitivity analyses demonstrate the importance of soil moisture and, secondarily, organic matter content in controlling the migration of VOC vapor through the unsaturated zone. The interpretation of soil-gas surveys can be complicated by unknown spatial heterogeneities in soil moisture and organic carbon content, temporal variations in moisture content, extent of contaminant migration as a non-aqueous phase liquid and by the unknown extent of VOC liquid and contaminated ground water.     

饱和土体-地下综合管廊结构地震响应分析

将土体视为固-液两相介质,基于饱和土体有效应力原理,建立饱和土体-地下综合管廊结构体系相互作用动力模型:在地应力平衡的静力状态下采用Duncan-Chang非线性弹性本构模型,在地震波作用的动力状态下采用Davidenkov非线性黏弹性本构模型;考虑饱和土体黏弹性动力人工边界条件,将地震动作用转化为作用在人工边界节点上的动力荷载。模型考察不同地震波时程、地震波加速度峰值、入射角度、孔隙率以及地应力场的影响,得出如下结论:(1)地震波的卓越周期与场地卓越周期相近时引起结构上的变形最大;随着地震波加速度峰值的增大结构变形增大;随着地震波入射角度的增加结构变形增大,地震波斜入射情况下产生的行波效应使得结构变形最大。(2)土体材料的孔隙水压力是影响地震中结构变形的主要因素之一。(3)将土体材料考虑为单相介质时结构上的变形要比考虑为固-液两相介质时大得多,直接将饱和土体场地中得到的地震波等效荷载施加到单相土介质-结构动力相互作用模型上,能够得到与完全基于有效应力法一致的结果。     

Sorption of ortho‐Phenylphenol to Soils

Conflicting sorption coefficients for ortho‐phenylphenol (OPP) have been reported in the literatures, which resulted in the conflicting assessments on OPP mobility in soil. To ascertain the sorption coefficient of OPP, batch experiments were performed based on OECD guideline 106, using three types of soils. Headspace solid‐phase microextraction (HS‐SPME) and GC‐MS were applied to the determination of OPP concentration in the liquid phase. The sorption isotherms obtained for all three soils under equilibrium conditions were described well, assuming linear sorption. The organic carbon normalized distribution coefficients (K_oc) ranged from 894 to 1703 L kg^?1, which suggested that OPP is moderately mobile in soil. The results also showed that the K_oc value of OPP can be predicted precisely from K_ow, whereas it was underestimated by one order of magnitude when water solubility is used.     

SV波在饱和土体自由表面的反射

基于饱和两相介质弹性波动方程分析SV波在饱和土体自由表面的反射问题,引入波动方程的势函数解答,求解出二维问题中SV波入射情况下饱和土体自由场的位移、速度、加速度和应力响应。在饱和土体自由场响应解析解基础上,建立SV波入射下饱和土体自由场静、动力有限元模型。建模中考虑了如下几方面因素:(1)在不同分析步,对土体单元赋予不同材料本构。通过*model change命令进行单元生死设定,从而实现在初始应力场平衡的静力状态下采用DuncanChang本构模型,而地震波动输入时采用Davidenkov动力本构模型;(2)采用多孔介质黏弹性人工边界条件,在人工边界上分别施加固相和液相介质的弹簧和阻尼来模拟饱和土体中能量的传播;(3)将地震波转化为作用在人工边界上的等效地震荷载,施加到人工边界节点上;(4)土体单元采用4结点平面应变孔压单元(CPE4P)。有限元计算与解析解比较结果表明:SV波在垂直入射和掠入射时,竖向位移响应为零;在45°左右入射时,水平位移响应最大;60°左右入射时,竖向位移响应最大。这些结论与解析解吻合较好,本文模型为建立土-结构动力相互作用模型打下良好的基础。     

Upstream Dispersion in Solute Transport Models: A Simple Evaluation and Reduction Methodology

This article outlines analytical solutions to quantify the length scale associated with “upstream dispersion,” the artificial movement of solutes in the opposite direction to groundwater flow, in solute transport models. Upstream dispersion is an unwanted artifact in common applications of the advection-dispersion equation (ADE) in problems involving groundwater flow in the direction of increasing solute concentrations. Simple formulae for estimating the one-dimensional distance of upstream dispersion are provided. These show that under idealized conditions (i.e., steady-state flow and transport, and a homogeneous aquifer), upstream dispersion may be a function of only longitudinal dispersivity. The scale of upstream dispersion in a selection of previously presented situations is approximated to highlight the utility of the presented formulae and the relevance of this ADE anomaly in common transport problems. Additionally, the analytical solution is applied in a hypothetical scenario to guide the modification of dispersion parameters to minimize upstream dispersion.     

关于孔隙相互孤立饱和土体弹性模量的讨论

对孔隙相互孤立饱和土体的弹性力学模型的部分观点进行了讨论。按Carroll的力学响应理论重新推导了此类土体的弹性本构关系和体积模量表达式。认为在准静态(很低频率)的波动问题中孔隙流体与固体骨架的相对运动为零,土体渗透系数为零或流体粘性为无穷大的条件是等价的。在此条件下没有必要区分孔隙是否连通,它们遵循相同的规律,而且其土体弹性体积模量的结果与Biot-Geertsma理论的无排水准静态的弹性体积模量相同。     

Transport and Biological Fate of Toluene in Low-Permeability Soils

The effect of simultaneous sorption, diffusion, and biodegradation on the fate and transport of toluene in low-permeability soil formations was examined. A transport model accounting for vapor and liquid sorption, vapor diffusions, and first-order biodegradation was developed to describe the movement of volatile solute in unsaturated soils. Modeling studies were followed with laboratory batch and column studies on fine-grained soil samples obtained from a gasoline-contaminated site. Batch experiments yielded the sorption and diffusion coefficients for generating theoretical solute transport profiles. Column studies were conducted to examine toluene sorption, diffusion, and biodegradation under aerobic and denitrifying conditions. Results from the column studies indicated that vapor sorption onto the soil was minimal due to the high moisture content of the soil. Comparison of model predictions with experimental results indicated that the SASK model, which is based on the resistivity theory, provided a more accurate prediction of the vapor phase tortuosity than the frequently used Millington-Quirk equation. Laboratory results of toluene concentration profiles matched well with the model predictions and yielded degradation rates comparable to those obtained in the field. Column studies, examining toluene biodegradation under aerobic and denitrifying conditions in low-permeability soils, indicated that the presence of excess nitrate in aerobic environments yielded higher solute degradation rates than those observed under exclusively aerobic systems.     

太湖地区丘陵旱地土壤磷的吸持解吸特征

旱地土壤中流失的磷是地表水体中磷的重要来源,本文通过模拟实验对比太湖地区丘陵旱地土壤和水稻土对磷的吸持解吸特征,结果表明虽然旱地土壤对磷的固定能力略高于水稻土,但由于旱地土壤的有效磷水平普遍高于水稻土,因而前者磷的吸持饱和度（DPS）要大大高于后者,这就决定了旱地土壤中的磷被淋溶或以溶解态随径流流失的风险和数量也同于水稻土,磷吸持饱和度是土壤磷素水平和土壤固磷能力的综合指标,更能反映土壤固相中的磷进入液相的难易程序,可以作为评价水－土界面磷迁移能力 重要指标。     

饱和黄土液化判别方法的两点发现

黄土液化实例多与现有可液化土地质年代规定不符。采用动三轴弯曲元试验设备对原状黄土饱和过程进行剪切波速跟踪测试,发现黄土浸水、结构先破坏再固结形成新的稳定结构的特点,证明经历过饱和的黄土已不能再视其为饱和前的地质年代。兰州马兰黄土中的黏粒有些是以黏土团块的形式存在,而黏土团块并不影响其他部分粉质土的液化,因此对于此类土应用黏粒含量进行液化判别时,应考虑将黏土团块不计入黏粒含量。     

Laboratory Study of Air Sparging of TCE-Contaminated Saturated Soils and Ground Water

Air sparging has proven to be an effective remediation technique for treating saturated soils and ground water contaminated by volatile organic compounds (VOCs). Since little is known about the system variables and mass transfer mechanisms important to air sparging, several researchers have recently performed laboratory investigations to study such issues. This paper presents the results of column experiments performed to investigate the behavior of dense nonaqueous phase liquids (DNAPFs). specifically trichloroethylene (TCE), during air sparging. The specific objectives of the study were (1) to compare the removal of dissolved TCE with the removal of dissolved light nonaqueous phase liquids (LNAPLs). such as benzene or toluene; (2) to determine the effect of injected air-flow rate on dissolved TCE removal; (3) to determine the effect of initial dissolved TCE concentration on removal efficiency; and (4) to determine the differences in removal between dissolved and pure-phase TCE. The test results showed that (1) the removal of dissolved TCE was similar to that of dissolved LNAPL: (2) increased air-injection rates led to increased TCE removal at lower ranges of air injection, but further increases at higher ranges of air injection did not increase the rate of removal, indicating a threshold removal rate had been reached; (3) increased initial concentration of dissolved TCE resulted in similar rates of removal: and (4) the removal of pure-phase TCE was difficult using a low air-injection rate, but higher air-injection rates led to easier removal.