Using fugacity to predict volatile emissions from layered materials with a clay/polymer diffusion barrier

Structural insulated panels (SIPs) have significant environmental and energy advantages. However, the tight structure that results may cause degraded indoor air quality and the potential release of volatile organic compounds (VOCs) from these layered materials must be considered. A physically based model for predicting VOC emissions from multi-layer materials is described. Fugacity is used to eliminate the concentration discontinuities at the interface between layers. This avoids an obstacle associated with numerically simulating mass transfer in composite materials. The numerical model is verified for a double-layer system by comparing predicted concentrations to those obtained with a previously published analytical model. In addition, hexanal emissions from multi-layer SIPs are simulated to demonstrate the usefulness of the fugacity approach. Finally, the multi-layer model is used to investigate the impact that clay/polyurethane nanocomposite diffusion barriers can have on VOC emissions. Indoor gas-phase concentrations can be greatly reduced with a barrier layer on the surface, thereby minimizing the environmental impact of SIPs.     

Numerical modeling of coupled variably saturated fluid flow and reactive transport with fast and slow chemical reactions

The couplings among chemical reaction rates, advective and diffusive transport in fractured media or soils, and changes in hydraulic properties due to precipitation and dissolution within fractures and in rock matrix are important for both nuclear waste disposal and remediation of contaminated sites. This paper describes the development and application of LEHGC2.0, a mechanistically based numerical model for simulation of coupled fluid flow and reactive chemical transport, including both fast and slow reactions in variably saturated media. Theoretical bases and numerical implementations are summarized, and two example problems are demonstrated. The first example deals with the effect of precipitation/dissolution on fluid flow and matrix diffusion in a two-dimensional fractured media. Because of the precipitation and decreased diffusion of solute from the fracture into the matrix, retardation in the fractured medium is not as large as the case wherein interactions between chemical reactions and transport are not considered. The second example focuses on a complicated but realistic advective-dispersive-reactive transport problem. This example exemplifies the need for innovative numerical algorithms to solve problems involving stiff geochemical reactions.     

Simulation of solute transport across low-permeability barrier walls

Low-permeability, non-reactive barrier walls are often used to contain contaminants in an aquifer. Rates of solute transport through such barriers are typically many orders of magnitude slower than rates through the aquifer. Nevertheless, the success of remedial actions may be sensitive to these low rates of transport. Two numerical simulation methods for representing low-permeability barriers in a finite-difference groundwater-flow and transport model were tested. In the first method, the hydraulic properties of the barrier were represented directly on grid cells and in the second method, the intercell hydraulic-conductance values were adjusted to approximate the reduction in horizontal flow, allowing use of a coarser and computationally efficient grid. The alternative methods were tested and evaluated on the basis of hypothetical test problems and a field case involving tetrachloroethylene (PCE) contamination at a Superfund site in New Hampshire. For all cases, advective transport across the barrier was negligible, but preexisting numerical approaches to calculate dispersion yielded dispersive fluxes that were greater than expected. A transport model (MODFLOW-GWT) was modified to (1) allow different dispersive and diffusive properties to be assigned to the barrier than the adjacent aquifer and (2) more accurately calculate dispersion from concentration gradients and solute fluxes near barriers. The new approach yields reasonable and accurate concentrations for the test cases.     

A model for contaminant mass flux in capped sediment under consolidation

The paper presents a model for contaminant transport and flux through a consolidating subaqueous sediment and overlying cap. The formulation is based on the effect of consolidation and excess pore pressure dissipation on transient, nonlinear advective component of transport through sediment and the cap. The consolidation is induced by the buoyant weight of the cap when it is placed on the contaminated sediments. One equation is presented for advective-diffusive transport through the sediment that is dependent upon soil/contaminant properties and transient advective velocity, which is calculated from a second equation based on the Terzaghi consolidation theory. A third equation is provided to describe the transport of contaminants in the cap. The parameters, including advective velocity, and boundary conditions used for contaminant transport through the cap are derived from the solution of the first two equations. The finite difference method is used to solve the system of equations for consolidation and contaminant transport. A hypothetical case is analyzed to demonstrate the formulation, and the results show that advection due to consolidation can accelerate breakthrough of contaminant through the cap by orders of magnitude. The derivation and results show that consolidation should be included for cap design, and that reactive caps are essential for delaying and reducing dissolved contaminant flux.     

An analytical solution for vertical transport of volatile chemicals in the vadose zone

An analytical solution is presented for one-dimensional vertical transport of volatile chemicals through the vadose zone to groundwater. The solution accounts for the important transport mechanisms of the steady advection of water and gas, diffusion and dispersion in water and gas, as well as adsorption, and first-order degradation. By assuming a linear, equilibrium partitioning between water, gas and the adsorbed chemical phases, the dependent variable in the mathematical model becomes the total resident concentration. The general solution was derived for cases having a constant initial total concentration over a discrete depth interval and a zero initial total concentration elsewhere. A zero concentration gradient is assumed at the groundwater table. Examples are given to demonstrate the application of the new solution for calculating the case of a non-uniform initial source concentration, and estimating the transport of chemicals to the groundwater and the atmosphere. The solution was also used to verify a numerical code called VLEACH. We discovered an error in VLEACH, and found that the new solution agreed very well with the numerical results by corrected VLEACH. A simplified solution to predict the migration of volatile organic chemical due to the gas density effect has shown that a high source concentration may lead to significant downward advective gas-phase transport in a soil with a high air-permeability.     

A numerical model for transient-hysteretic flow and solute transport in unsaturated porous media

A two-dimensional flow and transport model was developed for simulating transient water flow and nonreactive solute transport in heterogeneous, unsaturated porous media containing air and water. The model is composed of a unique combination of robust and accurate numerical algorithms for solving the Richards', Darcy flux, and advection-dispersion equations. The mixed form of Richards' equation is solved using a finite-element formulation and a modified Picard iteration scheme. Mass lumping is employed to improve solution convergence and stability behavior. The flow algorithm accounts for hysteresis in the pressure head-water content relationship. Darcy fluxes are approximated with a Galerkin and Petrov-Galerkin finite-element method developed for random heterogeneous porous media. The transport equation is solved using an Eulerian-Lagrangian method. A multi-step, fourth-order Runge-Kutta, reverse particle tracking technique and a quadratic-linear interpolation scheme are shown to be superior for determining the advective concentration. A Galerkin finite-element method is used for approximating the dispersive flux. The unsaturated flow and transport model was applied to a variety of rigorous problems and was found to produce accurate, mass-conserving solutions when compared to analytical solutions and published numerical results.     

Determination of spatially-resolved porosity, tracer distributions and diffusion coefficients in porous media using MRI measurements and numerical simulations

This work is focused on measuring the concentration distribution of a conservative tracer in a homogeneous synthetic porous material and in heterogeneous natural sandstone using MRI techniques, and on the use of spatially resolved porosity data to define spatially variable diffusion coefficients in heterogeneous media. The measurements are made by employing SPRITE, a fast MRI method that yields quantitative, spatially-resolved tracer concentrations in porous media. Diffusion experiments involving the migration of H(2)O into D(2)O-saturated porous media are conducted. One-dimensional spatial distributions of H(2)O-tracer concentrations acquired from experiments with the homogeneous synthetic calcium silicate are fitted with the one-dimensional analytical solution of Fick's second law to confirm that the experimental method provides results that are consistent with expectations for Fickian diffusion in porous media. The MRI-measured concentration profiles match well with the solution for Fick's second law and provide a pore-water diffusion coefficient of 1.75×10(-9)m(2)s(-1). The experimental approach was then extended to evaluate diffusion in a heterogeneous natural sandstone in three dimensions. The relatively high hydraulic conductivity of the sandstone, and the contrast in fluid density between the H(2)O tracer and the D(2)O pore fluid, lead to solute transport by a combination of diffusion and density-driven advection. The MRI measurements of spatially distributed tracer concentration, combined with numerical simulations allow for the identification of the respective influences of advection and diffusion. The experimental data are interpreted with the aid of MIN3P-D - a multicomponent reactive transport code that includes the coupled processes of diffusion and density-driven advection. The model defines local diffusion coefficients as a function of spatially resolved porosity measurements. The D(e) values calculated for the heterogeneous sandstone and used to simulate diffusive and advective transport range from 5.4×10(-12) to 1.0×10(-10)m(2)s(-1). These methods have broad applicability to studies of contaminant migration in geological materials.     

Design and laboratory testing of a chamber device to measure total flux of volatile organic compounds from the unsaturated zone under natural conditions

To determine if an aquifer contaminated with volatile organic compounds (VOCs) has potential for natural remediation, all natural processes affecting the fate and transport of VOCs in the subsurface must be identified and quantified. This research addresses the quantification of air-phase volatile organic compounds (VOCs) leaving the unsaturated zone soil gas and entering the atmosphere-including the additional flux provided by advective soil-gas movement induced by barometric pumping. A simple and easy-to-use device for measuring VOC flux under natural conditions is presented. The vertical flux chamber (VFC) was designed using numerical simulations and evaluated in the laboratory. Mass-balance numerical simulations based on continuously stirred tank reactor equations (CSTR) provided information on flux measurement performance of several sampling configurations with the final chamber configuration measuring greater than 96% of model-simulated fluxes. A laboratory device was constructed to evaluate the flux chamber under both diffusion-only and advection-plus-diffusion transport conditions. The flux chamber measured an average of 82% of 15 diffusion-only fluxes and an average of 95% of 15 additional advection-plus-diffusion flux experiments. The vertical flux chamber has the capability of providing reliable measurement of VOC flux from the unsaturated zone under both diffusion and advection transport conditions.     

Analytical model of solute transport by unsteady unsaturated gravitational infiltration

Penetration of reactive solute into a soil during a cycle of water infiltration and redistribution is investigated by deriving analytical closed form solutions for fluid flux, moisture content and contaminant concentration. The solution is developed for gravitational flow and advective transport and is applied to two scenarios of solute applications encountered in the applications: a finite pulse of solute dissolved in irrigation water and an instantaneous pulse broadcasted onto the soil surface. Through comparison to simulations of Richards' flow, capillary suction is shown to have contrasting effects on the upper and lower boundaries of the fluid pulse, speeding penetration of the wetting front and reducing the rate of drying. This leads to agreement between the analytical and numerical solutions for typical field and experimental conditions. The analytical solution is further incorporated into a stochastic column model of flow and transport to compute mean solute concentration in a heterogeneous field. An unusual phenomenon of plume contraction is observed at long times of solute propagation during the drying stage. The mean concentration profiles match those of the Monte-Carlo simulations for capillary length scales typical of sandy soils.     

Numerical simulation of air pollution in urban areas: Model development

A three-dimensional, grid-based numerical air pollution model for the estimation of air pollutant concentrations in an urban area is developed. Based on the continuity equation, the modeling system incorporates the combined influences of advective transport, turbulent diffusion, chemical transformation, source emissions and surface removal of air contaminants. Recent developments in plume rise and plume penetration processes, objective wind field analysis procedures and numerical solution techniques incorporated into the model are described.     

Dynamics of sulfur dioxide absorption in a raindrop falling at terminal velocity

Sulfur dioxide absorption dynamics in a raindrop are studied numerically by means of a fully numerical simulation method (FNSM) in which a composite orthogonal grid system consisting of both gas- and liquid-phase is adopted. When a raindrop with fixed radius falls in association with terminal velocity, a recirculation bubble always accompanies the gas-phase flow field in the aft region of the drop. With regard to the drop internal flow structure it has a drastic variation with drop size. When the drop radius is small (e.g. r_s=200 μm), only a single vortex motion is seen inside the drop. Under such a situation, sulfur dioxide absorbed from the interface is mainly transported from the area in front of the aft stagnation point. In contrast, as the drop is relatively large, say, r_s=500 μm, it is of interest to find that a double-vortex motion, composed of a primary and a secondary vortexes, is clearly observed. As a result, the onset of SO₂ transport process occurs at the drop’s surface near where the two vortexes meet. By defining a drop mass transport number, it indicates that the mass transported via internal circulation is always much faster than that by mass diffusion and the latter is highly relevant to the drop radius. Accordingly, the SO₂ transport dynamics in a raindrop is essentially determined by drop size.     

Diffusive partitioning tracer test for the quantification of nonaqueous phase liquid (NAPL) in the vadose zone: Performance evaluation for heterogeneous NAPL distribution

A partitioning tracer test based on gas-phase diffusion in the vadose zone yields estimates of the residual nonaqueous phase liquid (NAPL) saturation. The present paper investigates this technique further by studying diffusive tracer breakthrough curves in the vadose zone for a heterogeneous NAPL distribution. Tracer experiments were performed in a lysimeter with a horizontal layer of artificial kerosene embedded in unsaturated sand. Tracer disappearance curves at the injection point and tracer breakthrough curves at some distance from the injection point were measured inside and outside of the NAPL layer. A numerical code was used to generate independent model predictions based on the physicochemical sand, NAPL, and tracer properties. The measured and modeled tracer breakthrough curves were in good agreement confirming the validity of important modeling assumptions such as negligible sorption of chlorofluorocarbon (CFC) tracers to the uncontaminated sand and their fast reversible partitioning between the soil air and the NAPL phase. Subsequently, the model was used to investigate different configurations of NAPL contamination. The experimental and model results show that the tracer disappearance curves of a single-well diffusive partitioning tracer test (DPTT) are dominated by the near-field presence of NAPL around the tip of the soil gas probe. In contrast, breakthrough curves of inter-well tracer tests reflect the NAPL saturation in between the probes, although there is no unique interpretation of the tracer signals if the NAPL distribution is heterogeneous. Numerical modeling is useful for the planning of a DPTT application. Simulations suggest that several cubic meters of soil can be investigated with a single inter-well partitioning tracer test of 24-hour duration by placing the injection point in the center of the investigated soil volume and probes at up to 1 m distance for the monitoring of gaseous tracers.     

Oxygen Stratification in Industrial Boiler Stacks

The DSS/2 software system was used to numerically simulate oxygen stratification effects in a typical industrial boiler stack. The model developed for the simulation involves turbulent mass transfer in a tube with a constant, high mass flux from the walls of the tube. The concentration profile is considered to be developing, but the velocity profile is assumed to be fully developed. The model is for a gaseous system at constant temperature and pressure; the transport properties are shown to be constant. The universal velocity profile for turbulent flow and von Karman’s expressions for the turbulent eddy diffusivity are incorporated in the model. An expression for a wall concentration profile corresponding to conditions of high mass flux at the wall is developed. The model is then expressed as a system of partial differential equations which are solved using the numerical method of lines together with the DSS/2 numerical integration system. The resulting concentration profiles are presented for various flux rates which might be encountered in a typical stack. These profiles indicate that a potential for oxygen stratification in stacks does exist; they also suggest that careful placement of the oxygen cell probe can help minimize such errors.     

Formulation of a soil-pesticide transport model based on a compartmental approach

A semianalytical soil-pesticide transport model is formulated based on a compartmental approach to determine spatial and temporal variations of pesticide residues across a soil profile. The compartmental model is implemented by drawing an analogy between a series of continuous-flow stirred tank reactors and a soil horizon that consists of multiple perfectly mixed compartments. The analogy is strengthened by exploiting a relation between the compartment series and the conventional convective-dispersive equation (CDE) for vertical transport in the soil. Consequently, the number of compartments in the model formulation is not free, but dictated as a function of transport parameters. The model formulation allows consideration of arbitrary boundary value specifications and also, for some cases, spatially varying initial concentration profiles. Sorption kinetics is represented via a two-site model that involves a linear sorption isotherm and a first-order irreversible sorption or a radial diffusive penetrating model. For these three cases, analysis of the compartmental model allows the resultant concentration profiles to be expressed in terms of the Poisson distribution. When a nonlinear kinetic sorption model is used to simulate the sorption processes, an analytical solution is not found and a numerical approach is required.     

Multidimensional, multicomponent, subsurface reactive transport in nonuniform velocity fields: code verification using an advective reactive streamtube approach

High performance computing has made possible the development of high resolution, multidimensional, multicomponent reactive transport models that can be used to analyze complex geochemical environments. However, as increasingly complex processes are included in these models, the accuracy of the numerical formulation coupling the nonlinear processes becomes difficult to verify. Analytical solutions are not available for realistically complex problems and benchmark solutions are not generally available for specific problems. We present an advective reactive streamtube (ARS) transport technique that efficiently provides accurate solutions of nonlinear multicomponent reactive transport in nonuniform multidimensional velocity fields. These solutions can be compared with results from Eulerian-based advection-dispersion-reaction models to evaluate the accuracy of the numerical formulation used. The ARS technique includes mixed equilibrium and kinetic complexation and precipitation-dissolution reactions subject to the following assumptions: (1) transport is purely advective (i.e., no explicit diffusion or dispersion), and (2) chemistry is described by a canonical system of reactions that evolves with time and is unaffected by position in space. Results from the ARS technique are compared with results from the massively parallel, multicomponent reactive transport model MCTRACKER on a test problem involving irreversible oxidation of organic carbon and reaction of the oxidation products with two immobile mineral phases, gypsum and calcite, and fifteen aqueous complexes. Truncation error, operator splitting error, and the nonlinear transformation of these errors in the high-resolution reactive transport model are identified for this problem.     

Anomalous transport of colloids and solutes in a shear zone

Transport experiments with colloids and radionuclides in a shear zone were conducted during the Colloid and Radionuclide Retardation experiment (CRR) at Nagra's Grimsel Test Site. Breakthrough curves of bentonite colloids and uranine, a non-sorbing solute, were measured in an asymmetric dipole flow field. The colloid breakthrough is earlier than that of uranine. Both breakthrough curves show anomalously long late time tails and the slope of the late time tails for the colloids is slightly higher. Anomalous late time tails are commonly associated with matrix diffusion processes; the diffusive interaction of solutes transported in open channels with the adjacent porous rock matrix or zones of stagnant water. The breakthrough curves for different colloid size classes are very similar and show no signs of fractionation due to their (size-dependent) diffusivity. It is proposed that tailing of the colloids is mainly caused by the structure of the flow field and that for the colloid transport, matrix diffusion is of minor importance. This has consequences for the interpretation of the uranine breakthrough. Comparisons of experimental results with numerical studies and with the evaluation of the colloid breakthrough with continuous time random theory imply that the tailing in the conservative solute breakthrough in this shear zone is not only caused by matrix diffusion. Part of the tailing can be attributed to advective transport in fracture networks and advection in low velocity regions. Models based on the advection-dispersion equation and matrix diffusion do not properly describe the temporal and spatial evolution of colloid and solute transport in such systems with a consistent set of parameters.     

Investigation of hydrogeologic processes in a dipping layer structure: 2. Transport and biodegradation of organics

Numerical simulation tools have been used to study the dominating processes during transport of aromatic hydrocarbons in the unsaturated soil zone. Simulations were based on field observations at an experimental site located on a glacial delta plain with pronounced layered sedimentary structures. A numerical model for transport in the unsaturated zone, SWMS-3D, has been extended to incorporate coupled multispecies transport, microbial degradation following Monod kinetics and gas diffusive transport of oxygen and hydrocarbons. The flow field parameters were derived from previous work using nonreactive tracers. Breakthrough curves (BTC) from the hydrocarbon field experiment were used to determine sorption parameters and Monod kinetic parameters using a fitting procedure. The numerical simulations revealed that the assumption of homogeneous layers resulted in deviations from the field observations. The deviations were more pronounced with incorporation of reactive transport, compared with earlier work on nonreactive transport. To be able to model reasonable BTC, sorption had to be reduced compared to laboratory experiments. The initial biomass and the maximum utilisation rate could be adjusted to capture both the initial lag phase and the overall degradation rate. Nevertheless, local oxygen limitation is predicted by the model, which was not observed in the field experiment. Incorporation of evaporation and diffusive gas transport of the hydrocarbons did not significantly change the local oxygen demand. The main cause of the observed discrepancies between model and field are attributed to channelling as a result of small-scale heterogeneities such as biopores.     

A mountain-scale thermal-hydrologic model for simulating fluid flow and heat transfer in unsaturated fractured rock

A multidimensional, mountain-scale, thermal-hydrologic (TH) numerical model is presented for investigating unsaturated flow behavior in response to decay heat from the proposed radioactive waste repository in the Yucca Mountain unsaturated zone (UZ), The model, consisting of both two-dimensional (2-D) and three-dimensional (3-D) representations of the UZ repository system, is based on the current repository design, drift layout, thermal loading scenario, and estimated current and future climate conditions. This mountain-scale TH model evaluates the coupled TH processes related to mountain-scale UZ flow. It also simulates the impact of radioactive waste heat release on the natural hydrogeological system, including heat-driven processes occurring near and far away from the emplacement tunnels or drifts. The model simulates predict thermally perturbed liquid saturation, gas- and liquid-phase fluxes, and water and rock temperature elevations, as well as the changes in water flux driven by evaporation/condensation processes and drainage between drifts. These simulations provide insights into mountain-scale thermally perturbed flow fields under thermal loading conditions.     

Biodegradation of hydrocarbons vapors: Comparison of laboratory studies and field investigations in the vadose zone at the emplaced fuel source experiment, Airbase Vaerløse, Denmark

The natural attenuation of volatile organic compounds (VOCs) in the unsaturated zone can only be predicted when information about microbial biodegradation rates and kinetics are known. This study aimed at determining first-order rate coefficients for the aerobic biodegradation of 13 volatile petroleum hydrocarbons which were artificially emplaced as a liquid mixture during a field experiment in an unsaturated sandy soil. Apparent first-order biodegradation rate coefficients were estimated by comparing the spatial evolution of the resulting vapor plumes to an analytical reactive transport model. Two independent reactive numerical model approaches have been used to simulate the diffusive migration of VOC vapors and to estimate degradation rate coefficients. Supplementary laboratory column and microcosm experiments were performed with the sandy soil at room temperature under aerobic conditions. First-order kinetics adequately matched the lab column profiles for most of the compounds. Consistent compound-specific apparent first-order rate coefficients were obtained by the three models and the lab column experiment, except for benzene. Laboratory microcosm experiments lacked of sensitivity for slowly degrading compounds and underestimated degradation rates by up to a factor of 5. Addition of NH3 vapor was shown to increase the degradation rates for some VOCs in the laboratory microcosms. All field models suggested a significantly higher degradation rate for benzene than the rates measured in the lab, suggesting that the field microbial community was superior in developing benzene degrading activity.