Macrodispersivity and Large-scale Hydrogeologic Variability

Although groundwater velocities vary over a wide range of spatial scales it is generally only feasible to model the largest variations explicitly. Smaller-scale velocity variability must be accounted for indirectly, usually by increasing the magnitude of the dispersivity tensor (i.e. by introducing a so-called macrodispersivity). Most macrodispersion theories tacitly assume that a macrodispersivity tensor which works well when there is only small-scale velocity variability will also work well when there is larger-scale variability. We analyze this assumption in a high resolution numerical experiment which simulates solute transport through a two-scale velocity field. Our results confirm that a transport model which uses an appropriately adjusted macrodispersivity can reproduce the large-scale features of a solute plume when the velocity varies only over small scales. However, if the velocity field includes both small and large-scale components, the macrodispersivity term does not appear to be able to capture all of the effects of small-scale variability. In this case the predicted plume is more well mixed and consistently underestimates peak solute concentrations at all times. We believe that this result can be best explained by scale interactions resulting from the nonlinear transformation from velocity to concentration. However, additional analysis will be required to test this hypothesis.     

Role of Molecular Diffusion in Contaminant Migration and Recovery in an Alluvial Aquifer System

Highly-resolved simulations and flow and transport in an alluvial system at the Lawrence Livermore National Laboratory (LLNL) site explore the role of diffusion in the migration and recovery of a conservative solute. Heterogeneity is resolved to the hydrofacies scale with a discretization of 10.0, 5.0 and 0.5m in the strike, dip and vertical directions of the alluvial-fan system. Transport simulations rely on recently developed random-walk techniques that accurately account for local dispersion processes at interfaces between materials with contrasting hydraulic and transport properties. Solute migration and recovery by pump and treat are shown to be highly sensitive to magnitude of effective diffusion coefficient. Further, transport appears significantly more sensitive to the diffusion coefficient than to local-scale dispersion processes represented by a dispersivity coefficient. Predicted hold back of solute mass near source locations during ambient migration and pump-and-treat remediation is consistent with observations at LLNL, and reminiscent of observations at the MADE site of Columbus Air Force Base, Mississippi. Results confirm the important role of diffusion in low-conductivity materials and, consequently, its impact on efficacy of pump-and-treat and other remedial technologies. In a typical alluvial system on a decadal time scale this process is, in part, fundamentally nonreversible because the average thickness of low-K hydrofacies is considerably greater than the mean-square length of penetration of the solute plume.     

稳定分层湍流中逆梯度输运的关联分析研究

本文采用关联分析方法研究了稳定温度分层湍流中的结构特性、输运特性,以及热量、动量逆梯度输运现象的尺度效应及其参数演化．首先采用大涡模拟方法对稳定分层湍流中的结构特性和输运特性进行了分析,将逆梯度输运发生的时间尺度作为已知条件,结合关联量分析方法在波数空间中的解析解,对逆梯度输运现象的尺度效应进行了分析研究．结果发现,稳定分层强度较大的流动中发生垂向热量及动量逆梯度输运现象,发生的结构尺度与关联分析所发现垂向热量、动量逆梯度输运的波数形成了呼应．随着分层强度增加,热量、动量的输运强度均受抑制,与逆梯度输运关联的流场结构尺度减小,同样的效应也发生在流场结构向下游演化的过程中．     

Experiments on Solute Transport in Aggregated Porous Media: Are Diffusions Within Aggregates and Hydrodynamic Dispersion Independent?

Two region models for solute transport in porous media assume that hydrodynamic dispersion in mobile water and solute diffusion within immobile water regions are independent. Experimental and theoretical results for transport through a macropore indicate that hydrodynamic dispersion and solute exchange are interdependent. Experiments were carried out to investigate this problem for a column packed with spherical porous aggregates. The effective diffusion coefficient of a tracer within the agreggates was determined from specific experiments. The dispersivity of the bed was determined from experiments carried out with a column filled with nonporous beads. We took advantage of the dependence of hydrodynamic dispersion on density ratios between the invading and displaced solutions to obtain a set of breakthrough curves corresponding to situations where the diffusion coefficient remains constant, whereas the dispersivity varies. Simulations reproduce correctly the experiments. Small discrepancies are noted that can be corrected either by increasing the dispersion coefficient or by fitting the external mass transfer coefficient. Increased dispersion coefficients probably reveal a modification of Taylor dispersion due to solute exchange. The fitted external mass transfer coefficients are close to the values obtained with classical correlations of the chemical engineering literature.     

Linear and non-linear double diffusive convection in a rotating porous layer using a thermal non-equilibrium model

Double diffusive convection in a fluid-saturated rotating porous layer heated from below and cooled from above is studied when the fluid and solid phases are not in local thermal equilibrium, using both linear and non-linear stability analyses. The Darcy model that includes the time derivative and Coriolis terms is employed as momentum equation. A two-field model that represents the fluid and solid phase temperature fields separately is used for energy equation. The onset criterion for stationary, oscillatory and finite amplitude convection is derived analytically. It is found that small inter-phase heat transfer coefficient has significant effect on the stability of the system. There is a competition between the processes of thermal and solute diffusions that causes the convection to set in through either oscillatory or finite amplitude mode rather than stationary. The effect of solute Rayleigh number, porosity modified conductivity ratio, Lewis number, diffusivity ratio, Vadasz number and Taylor number on the stability of the system is investigated. The non-linear theory based on the truncated representation of Fourier series method predicts the occurrence of subcritical instability in the form of finite amplitude motions. The effect of thermal non-equilibrium on heat and mass transfer is also brought out.     

Scale properties of turbulent transport and coherent structure in stably stratified flows

The empirical mode decomposition (EMD) is used to study the scale properties of turbulent transport and coherent structures based on velocity and temperature time series in stably stratified turbulence. The analysis is focused on the scale properties of intermittency and coherent structures in different modes and the contributions of energy-contained coherent structures to turbulent scalar counter-gradient transport (CGT). It is inferred that the velocity intermittency is scattered to more modes with the development of the stratified flow, and the intermittency is enhanced by the vertical stratification, especially in small scales. The anisotropy of the field is presented due to different time scales of coherent structures of streamwise and vertical velocities. There is global counter-gradient heat transport close to the turbulence-generated grid, and there is local counter-gradient heat transport at certain modes in different positions. Coherent structures play a principal role in the turbulent vertical transport of temperature.     

Analytical Solutions for Macrodispersion in a 3D Heterogeneous Porous Medium with Random Hydraulic Conductivity and Dispersivity

A stochastic analysis of macrodispersion for conservative solute transport in three-dimension (3D) heterogeneous statistically isotropic and anisotropic porous media when both hydraulic conductivity and local dispersivity are random is presented. Analytical expressions of macrodispersivity are derived using Laplace and Fourier transforms. The effects of various parameters such as ratio of transverse to longitudinal local dispersivity, correlation length ratio, correlation coefficient and direction of flow on asymptotic macrodispersion are studied. The behaviour of growth of macrodispersivity in preasymptotic stage is also shown in this paper. The variation in local dispersion coefficient causes change in transverse macrodispersivity. The consideration of random dispersivity along with random hydraulic conductivity indicates that the total dispersion is affected and important in the case when the hydraulic conductivity and dispersivity are correlated. It is observed that the pre-asymptotic behavior of the macrodispersivity is not sensitive to the choice of spectral density functions.     

结构层热扩散系数对超薄涂层热力性能的影响

将非傅立叶热传导模型(用于超薄热涂层)与傅立叶热传导模型(用于结构层)相结合 求解温度场,运用有限元法求解热涂层热应力和裂纹驱动力,并分析结构层材料热扩散系数 的变化对热涂层的热力学性能(温度场、应力场和断裂性能)的影响. 研究表明,结构层材 料性能变化对温度场的影响主要表现在热冲击后期,对热应力和裂纹尖端驱动力后期的变化 也有一定的影响.     

Transient Stochastic Analysis of Biodegradable Contaminant Transport: First-Order Decay

A stochastic methodology was used to analyze the field-scale transport of solutes in heterogeneous aquifers with first-order biodegradation. Spectral methods and perturbation techniques were utilized to develop expressions for the field-scale effective parameters in the mean transport model. Expressions were obtained for the longitudinal and transverse macrodispersivity coefficients, and effective velocity and an effective decay parameter for statistically anisotropic, and isotropic, heterogeneous porous medium, respectively. The behavior of these parameters was described as function of time and log K correlation scale. The expressions for asymptotic values of the dispersivity coefficients and effective velocity and decay parameters were also derived for the isotropic case.     

Generalized Taylor-Aris moment analysis of the transport of sorbing solutes through porous media with spatially-periodic retardation factor

Taylor-Aris dispersion theory, as generalized by Brenner, is employed to investigate the macroscopic behavior of sorbing solute transport in a three-dimensional, hydraulically homogeneous porous medium under steady, unidirectional flow. The porous medium is considered to possess spatially periodic geochemical characteristics in all three directions, where the spatial periods define a rectangular parallelepiped or a unit-element. The spatially-variable geochemical parameters of the solid matrix are incorporated into the transport equation by a spatially-periodic distribution coefficient and consequently a spatially-periodic retardation factor. Expressions for the effective or large-time coefficients governing the macroscopic solute transport are derived for solute sorbing according to a linear equilibrium isotherm as well as for the case of a first-order kinetic sorption relationship. The results indicate that for the case of a chemical equilibrium sorption isotherm the longitudinal macrodispersion incorporates a second term that accounts for the eflect of averaging the distribution coefficient over the volume of a unit element. Furthermore, for the case of a kinetic sorption relation, the longitudinal macrodispersion expression includes a third term that accounts for the effect of the first-order sorption rate. Therefore, increased solute spreading is expected if the local chemical equilibrium assumption is not valid. The derived expressions of the apparent parameters governing the macroscopic solute transport under local equilibrium conditions agreed reasonably with the results of numerical computations using particle tracking techniques.     

Linear and nonlinear double diffusive convection in a rotating sparsely packed porous layer using a thermal non-equilibrium model

Double diffusive convection in a fluid-saturated rotating porous layer is studied when the fluid and solid phases are not in local thermal equilibrium, using both linear and nonlinear stability analyses. The Brinkman model that includes the Coriolis term is employed as the momentum equation. A two-field model that represents the fluid and solid phase temperature fields separately is used for the energy equation. The onset criterion for stationary, oscillatory, and finite amplitude convection is derived analytically. It is found that small inter-phase heat transfer coefficient has significant effect on the stability of the system. There is a competition between the processes of thermal diffusion, solute diffusion, and rotation that causes the convection to set in through either oscillatory or finite amplitude mode rather than stationary. The effect of solute Rayleigh number, porosity modified conductivity ratio, Lewis number, diffusivity ratio, Vadasz number, and Taylor number on the stability of the system is investigated. The nonlinear theory based on the truncated representation of Fourier series method predicts the occurrence of subcritical instability in the form of finite amplitude motions. The effect of thermal non-equilibrium on heat and mass transfer is also brought out. Some of the convection systems previously reported in the literature is shown to be special cases of the system presented in this study.     

Thermal Dispersion Effects on Combined Convection in Non-Newtonian Fluids Along a Nonisothermal Vertical Plate in a Porous Medium

The method of non-similarity solution is used to study the influence of thermal dispersion on combined convection from vertical surfaces in a porous medium saturated with a power-law type non-Newtonian fluid. The coefficient of thermal diffusivity has been assumed to be the sum of molecular diffusivity and the dispersion thermal diffusivity due to mechanical dispersion. The transformed conservation laws are solved numerically for the case of variable surface heat flux conditions. Results for the details of the velocity and temperature fields as well as the Nusselt number have been presented.     

Effects of Material Symmetry on the Coefficients of Transport in Anisotropic Porous Media

The objective of this article is to highlight certain features of a number of coefficients that appear in models of phenomena of transport in anisotropic porous media, especially the coefficient of dispersion the second-rank tensor D _ij , and the dispersivity coefficient, the fourth-rank tensor a _ijkl , that appear in models of solute transport. Although we shall focus on the transport of mass of a dissolved chemical species in a fluid phase that occupies the void space, or part of it, the same discussion is also applicable to transport coefficients that appear in models that describe the advective mass flux of a fluid and the diffusive transport of other extensive quantities, like heat. The case of coupled processes, e.g. the simultaneous transport of heat and mass of a chemical species, are also considered. The entire discussion will be at the macroscopic level, at which a porous medium domain is visualized as a homogenized continuum.     

Experimental investigation of solute transport in large,homogeneous and heterogeneous,saturated soil columns

Laboratory tracer experiments were conducted to investigate solute transport in 12.5-m long, horizontally placed soil columns during steady saturated water flow. Two columns having cross-sectional areas of 10×10cm² were used: a uniformly packed homogeneous sandy column and a heterogeneous column containing layered, mixed, and lenticular formations of various shapes and sizes. The heterogeneous soil column gradually changed, on average, from coarse-textured at one end to fine-textured at the other end. NaCl breakthrough curves (BTC's) in the columns were measured with electrical conductivity probes inserted at 50- or 100-cm intervals. Observed BTC's in the homogeneous sandy column were relatively smooth and sigmoidal (S-shaped), while those in the heterogeneous column were very irregular, nonsigmoidal, and exhibited extensive tailing. Effective average pore-water velocities (v _eff) and dispersion coefficients (D _eff) were estimated simultaneously by fitting an analytical solution of the convection-dispersion equation to the observed BTC's. Velocity variations in the heterogeneous medium were found to be much larger than those in the homogeneous sand. Values of the dispersivity,=D _eff/v _eff, for the homogeneous sandy column ranged from 0.1 to 5.0 cm, while those for the heterogeneous column were as high as 200cm. The dispersivity for transport in both columns increased with travel distance or travel time, thus exhibiting scale-dependency. The heterogeneous soil column also showed the effects of preferential flow, i.e., some locations in the column showed earlier solute breakthrough than several locations closer to the inlet boundary. Spatial fluctuations in the dispersivity could be explained qualitatively by the particular makeup of the heterogeneities in the column.     

Non-passive Transport of Volatile Organic Compounds Infiltrated from a Surface Disk Source

An unsaturated flow and non-passive transport model for water-soluble organic compounds has been implemented in cylindrical coordinates, with a top boundary condition that accounts for different zones of the surface that can be under infiltration or volatilization independently. We simulated two-dimensional infiltration of aqueous mixtures of methanol from a disk source, its redistribution and volatilization in both homogeneous and heterogeneous soils. Simulations showed that the incoming composition significantly affects volumetric liquid content and concentration profiles, as well as the fraction of infiltrated mass of methanol that is released to the atmosphere. Concentration-dependent viscosity had the major impact on the liquid flow. The differences in volumetric liquid content and normalized concentration of methanol became more pronounced during transport through a soil composed of a clay lens embedded within a main matrix of sandy clay loam texture. Dispersion in the liquid-phase was the predominant transport mechanism when dispersivity at saturation was set to 7.8 cm. However, for dispersivity of 1.0 cm, changes in composition led to changes in surface tension inducing significantly higher liquid flow. In this case, liquid-phase advection was the most active transport mechanism for homogeneous soils and highly concentrated infiltrating mixtures.     

Migration of salts in the unsaturated zone caused by heating

Heat-transfer phenomena as well as moisture movement in unsaturated soils due to thermal gradients, have been extensively studied during the last four decades. Less attention has been devoted to the transport and redistribution of solutes caused by heating.Solar radiation, radioactive waste repositories, underground energy storage, buried electric cables and steam pipes, disposal of waste heat from power plants are examples of heat sources in the soil.Soil-water properties, such as surface tension, viscosity, density, as well as the equilibrium composition of phases, depend on temperature. Hence, nonuniform heating of a soil partially saturated by saline water has an effect on such processes as water flow under capillary and gravitational forces, evaporation, condensation and diffusion of vapor and transport and precipitation of salts.A mathematical model is presented for the migration of salts in the vadoze zone in the soil under nonisothermal conditions, taking into account the above-mentioned phenomena. The physical assumptions underlying the model are briefly discussed.The study of a particular case shows that under certain conditions, a heat source may attract dissolved salts, and cause their precipitation in the hot area.     

A Solute Transport in Stratified Media

Two-dimensional and steady solute transport in a stratified porous formation is analysed under assumption that the effect of pore-scale dispersion is negligible. The longitudinal dispersion produced as a result of the vertical variation of hydraulic conductivity is analysed by averaging the variability of a solute flux concentration and conductivity. The evolution of the solute flux concentration is expressed with respect to the correlated variable, that is the travel (arrival) time at a fixed location and the averaging procedure is constructed to satisfy the boundary condition where the inlet concentration is a known function of time. In such a statement, a velocity-averaged solute flux concentration is described by a conventional dispersion model (CDM) with a dispersion coefficient which is a function of the arrival time. It is demonstrated that such CDM satisfies the assumption that hydraulic conductivity of the layers is gamma distributed with the parameter of distribution which is chosen to represent a reasonable value of the field scale solute dispersion. The overall behaviour of the model is illustrated by several examples of two-dimensional mass transport.     

Effect of thermal dispersion on free convection in a fluid saturated porous medium

The present article considers a numerical study of thermal dispersion effect on the non-Darcy natural convection over a vertical flat plate in a fluid saturated porous medium. Forchheimer extension is considered in the flow equations. The coefficient of thermal diffusivity has been assumed to be the sum of molecular diffusivity and the dispersion thermal diffusivity due to mechanical dispersion. The non-dimensional governing equations are solved by the finite element method (FEM) with a Newton–Raphson solver. Numerical results for the details of the stream function, velocity and temperature contours and profiles as well as heat transfer rates in terms of Nusselt number are obtained. The study shows that the increase in thermal dispersion coefficient of the porous medium results in more heat energy to disperse away in the normal direction to the wall. This induces more fluid to flow along the wall, enhancing the heat transfer coefficient particularly near the wall.     

Analysis of hydrodynamic and thermal dispersion in porous media by means of a local approach

A pore scale analysis is implemented in this numerical study to investigate the behavior of microscopic inertia and thermal dispersion in a porous medium with a periodic structure. The macroscopic characteristics of the transport phenomena are evaluated with an averaging technique of the controlling variables at a pore scale level in an elementary cell of the porous structure. The Darcy–Forchheimer model describes the fluid motion through the porous medium while the continuity and Navier–Stokes equations are applied within the unit cell. An average energy equation is employed for the thermal part of the porous medium. The macroscopic pressure loss is computed in order to evaluate the dominant microscopic inertial effects. Local fluctuations of velocity and temperature at the pore scale are instrumental in the quantification of the thermal dispersion through the total effective thermal diffusivity. The numerical results demonstrate that microscopic inertia contributes significantly to the magnitude of the macroscopic pressure loss, in some instances with as much as 70%. Depending on the nature of the porous medium, the thermal dispersion may have a marked bearing on the heat transfer, particularly in the streamwise direction for a highly conducting fluid and certain values of the Peclet number.