CONTAMINANT TRANSPORT IN SEDIMENT UNDER THE INFLUENCE OF ADVECTIVE FLUX
Chemical flux across the sediment/water interface is controlled by a combination of diffusive, dispersive and advective processes. The advective process is a function of submarine groundwater discharge and tidal effects. In areas where surface water interacts with groundwater, ...
Technology Transfer Automated Retrieval System (TEKTRAN)
Analytical solutions of the advection-dispersion equation and related models are indispensable for predicting or analyzing contaminant transport processes in streams and rivers, as well as in other surface water bodies. Many useful analytical solutions originated in disciplines other than surface-w...
It is well known that the fate and transport of contaminants in the subsurface are controlled by complex processes including advection, dispersion-diffusion, and chemical reactions. However, the interplay between the physical transport processes and chemical reactions, and their...
Efficient mass transport by optical advection
NASA Astrophysics Data System (ADS)
Kajorndejnukul, Veerachart; Sukhov, Sergey; Dogariu, Aristide
2015-10-01
Advection is critical for efficient mass transport. For instance, bare diffusion cannot explain the spatial and temporal scales of some of the cellular processes. The regulation of intracellular functions is strongly influenced by the transport of mass at low Reynolds numbers where viscous drag dominates inertia. Mimicking the efficacy and specificity of the cellular machinery has been a long time pursuit and, due to inherent flexibility, optical manipulation is of particular interest. However, optical forces are relatively small and cannot significantly modify diffusion properties. Here we show that the effectiveness of microparticle transport can be dramatically enhanced by recycling the optical energy through an effective optical advection process. We demonstrate theoretically and experimentally that this new advection mechanism permits an efficient control of collective and directional mass transport in colloidal systems. The cooperative long-range interaction between large numbers of particles can be optically manipulated to create complex flow patterns, enabling efficient and tunable transport in microfluidic lab-on-chip platforms.
This manual describes the next generation of the modular three-dimensional transport model, MT3D, with significantly expanded capabilities, including the addition of (a) a third-order total-variation-diminishing (TVD) scheme for solving the advection term that is mass conservativ...
Efficient mass transport by optical advection
Kajorndejnukul, Veerachart; Sukhov, Sergey; Dogariu, Aristide
2015-01-01
Advection is critical for efficient mass transport. For instance, bare diffusion cannot explain the spatial and temporal scales of some of the cellular processes. The regulation of intracellular functions is strongly influenced by the transport of mass at low Reynolds numbers where viscous drag dominates inertia. Mimicking the efficacy and specificity of the cellular machinery has been a long time pursuit and, due to inherent flexibility, optical manipulation is of particular interest. However, optical forces are relatively small and cannot significantly modify diffusion properties. Here we show that the effectiveness of microparticle transport can be dramatically enhanced by recycling the optical energy through an effective optical advection process. We demonstrate theoretically and experimentally that this new advection mechanism permits an efficient control of collective and directional mass transport in colloidal systems. The cooperative long-range interaction between large numbers of particles can be optically manipulated to create complex flow patterns, enabling efficient and tunable transport in microfluidic lab-on-chip platforms. PMID:26440069
Advective turbulent transport in the fluid plasma
NASA Astrophysics Data System (ADS)
Min, Byung-Hoon; An, Chan-Yong; Kim, Chang-Bae
2013-10-01
The Hasegawa-Wakatani model (HWM) has been employed in pedagogical analyses of the physics behind the behavior of the tokamak plasmas. In addition to the geometric simplicity HWM has an appealing feature of sustaining autonomous quasi-steady state, unstable modes providing the power that is being transported by the nonlinear interactions and is eventually dissipated by the collisional damping at small scales. Emergence of the zonal flow out of the turbulence is a main candidate to cause the transition from the low plasma confinement to the high mode. In the study of such LH transition with the HWM, the adiabaticity parameter has been shown to play an important role in forcing the zonal flow that results in the regulation of the drift-wave turbulence. Instead of concentrating on the physics of the feedback loop between the turbulence and the zonal flow the present study focuses on the presence of the advective transport of the energy. Numerical simulations of HWM are performed and the connections between the advective transport and the zonal flow will be presented. This work was supported by the Supercpmputing Center/Korea Institute of Science and Technology Information with supercomputing resources including technical support (KSC-2013-C1-009).
Chaotic Advection, Fluid Spreading, and Groundwater Contaminant Plumes
NASA Astrophysics Data System (ADS)
Neupauer, R. M.; Mays, D. C.
2011-12-01
In situ remediation of contaminated groundwater requires degradation reactions at the interface between the contaminant plume and an injected treatment solution containing chemical or biological amendments. Therefore a promising approach to accelerate in situ remediation is to elongate the interface between the contaminant plume and treatment solution through fluid spreading. The literature on chaotic advection describes how to accomplish spreading in laminar flows, which lack the turbulent eddies that provide spreading in streams or engineered reactors. A key result from the literature on chaotic advection is that spreading is inherent in the vicinity of certain periodic points, which are points to which fluid particles return in successive iterations of chaotic flows. Specifically, spreading is enhanced near the stable and unstable manifolds associated with hyperbolic periodic points. We investigate the transient flow created with a four-well system in which wells are operated sequentially as either injection wells or extraction wells. In particular, we identify the periodic points and demonstrate that fluid spreading occurs nearby. For appropriately designed injection and extraction sequences, the periodic points are located near the interface between the contaminant plume and treatment solution, leading to elongation of the interface, with expected benefits of enhanced reaction and accelerated remediation.
Garges, J.A.; Baehr, A.L.
1998-01-01
The relative importance of advection and dispersion for both solute and vapor transport can be determined from type curves or concentration, flux, or cumulative flux. The dimensionless form of the type curves provides a means to directly evaluate the importance of mass transport by advection relative to that of mass transport by diffusion and dispersion. Type curves based on an analytical solution to the advection-dispersion equation are plotted in terms of dimensionless time and Peclet number. Flux and cumulative flux type curves provide additional rationale for transport regime determination in addition to the traditional concentration type curves. The extension of type curves to include vapor transport with phase partitioning in the unsaturated zone is a new development. Type curves for negative Peclet numbers also are presented. A negative Peclet number characterizes a problem in which one direction of flow is toward the contamination source, and thereby diffusion and advection can act in opposite directions. Examples are the diffusion of solutes away from the downgradient edge of a pump-and-treat capture zone, the upward diffusion of vapors through the unsaturated zone with recharge, and the diffusion of solutes through a low hydraulic conductivity cutoff wall with an inward advective gradient.
Backward fractional advection dispersion model for contaminant source prediction
NASA Astrophysics Data System (ADS)
Zhang, Yong; Meerschaert, Mark M.; Neupauer, Roseanna M.
2016-04-01
The forward Fractional Advection Dispersion Equation (FADE) provides a useful model for non-Fickian transport in heterogeneous porous media. The space FADE captures the long leading tail, skewness, and fast spreading typically seen in concentration profiles from field data. This paper develops the corresponding backward FADE model, to identify source location and release time. The backward method is developed from the theory of inverse problems, and then explained from a stochastic point of view. The resultant backward FADE differs significantly from the traditional backward Advection Dispersion Equation (ADE) because the fractional derivative is not self-adjoint and the probability density function for backward locations is highly skewed. Finally, the method is validated using tracer data from a well-known field experiment, where the peak of the backward FADE curve predicts source release time, while the median or a range of percentiles can be used to determine the most likely source location for the observed plume. The backward ADE cannot reliably identify the source in this application, since the forward ADE does not provide an adequate fit to the concentration data.
Advective and diffusive cosmic ray transport in galactic haloes
NASA Astrophysics Data System (ADS)
Heesen, Volker; Dettmar, Ralf-Jürgen; Krause, Marita; Beck, Rainer; Stein, Yelena
2016-05-01
We present 1D cosmic ray transport models, numerically solving equations of pure advection and diffusion for the electrons and calculating synchrotron emission spectra. We find that for exponential halo magnetic field distributions advection leads to approximately exponential radio continuum intensity profiles, whereas diffusion leads to profiles that can be better approximated by a Gaussian function. Accordingly, the vertical radio spectral profiles for advection are approximately linear, whereas for diffusion they are of `parabolic' shape. We compare our models with deep Australia Telescope Compact Array observations of two edge-on galaxies, NGC 7090 and 7462, at λλ 22 and 6 cm. Our result is that the cosmic ray transport in NGC 7090 is advection dominated with V=150^{+80}_{-30} km s^{-1}, and that the one in NGC 7462 is diffusion dominated with D=3.0± 1.0 × 10^{28}E_GeV^{0.5} cm^2 s^{-1}. NGC 7090 has both a thin and thick radio disc with respective magnetic field scale heights of hB1 = 0.8 ± 0.1 kpc and hB2 = 4.7 ± 1.0 kpc. NGC 7462 has only a thick radio disc with hB2 = 3.8 ± 1.0 kpc. In both galaxies, the magnetic field scale heights are significantly smaller than what estimates from energy equipartition would suggest. A non-negligible fraction of cosmic ray electrons can escape from NGC 7090, so that this galaxy is not an electron calorimeter.
Advective Removal of Intraparticle Uranium from Contaminated Vadose Zone Sediments, Hanford, USA
Ilton, Eugene S.; Qafoku, Nikolla; Liu, Chongxuan; Moore, D. A.; Zachara, John M.
2008-03-01
A column study on U contaminated vadose zone sediments from the Hanford Site, WA, was performed in order to aid the development of a model for predicting U(VI) release rates under a dynamic flow regime and for variable geochemical conditions. The sediments of interest are adjacent to and below tank BX-102, part of the BX tank farm that contained high level liquid radioactive waste. Two sediments, with different U(VI) loadings and intraparticle large fracture vs. smaller fracture ratios, were reacted with three different solutions. The primary reservoir for U(VI) appears to be a micron-sized nanocrystalline Na-U-Si phase, possibly Na-boltwoodite, that nucleated and grew on plagioclase grains that line fractures within sand-sized granitic clasts. The solutions were all calcite saturated and in equilibrium with atmospheric CO2, where one solution was simply DI-water, the second was a synthetic ground water (SGW) with elevated Na, and the third was the same SGW but with both elevated Na and Si. The latter two solutions were employed, in part, to test the effect of saturation state on U(VI) release. For both sediments and all three electrolytes, there was an initial rapid release of U(VI) to the advecting solution followed by a plateau of low U(VI) concentration. U(VI) effluent concentration increased during subsequent stop flow (SF) events. The electrolytes with elevated Na and Si appreciably depressed U(VI) concentrations relative to DI water. The effluent data for both sediments and all three electrolytes was simulated reasonably well by a three domain model (the advecting fluid, fractures, and matrix) that coupled U(VI) dissolution rates, intraparticle U(VI) diffusion, and interparticle advective transport of U(VI); where key transport and dissolution processes had been parameterized in previous batch studies. For the calcite-saturated DI-water, U(VI) concentrations in the effluent remained far below saturation with respect to Na-boltwoodite and release of U(VI) to
Lewis-Brown, Jean C.; Carleton, Glen B.; Imbrigiotta, Thomas E.
2006-01-01
Volatile organic compounds, predominantly trichloroethylene and its degradation products, have been detected in ground water at the Naval Air Warfare Center (NAWC), West Trenton, New Jersey. An air-stripping pump-and-treat system has been in operation at the NAWC since 1998. An existing ground-water-flow model was used to evaluate the effect of a change in the configuration of the network of recovery wells in the pump-and-treat system on flow paths of contaminated ground water. The NAWC is underlain by a fractured-rock aquifer composed of dipping layers of sedimentary rocks of the Lockatong and Stockton Formations. Hydraulic and solute-transport properties of the part of the aquifer composed of the Lockatong Formation were measured using aquifer tests and tracer tests. The heterogeneity of the rocks causes a wide range of values of each parameter measured. Transmissivity ranges from 95 to 1,300 feet squared per day; the storage coefficient ranges from 9 x 10-5 to 5 x 10-3; and the effective porosity ranges from 0.0003 to 0.002. The average linear velocity of contaminated ground water was determined for ambient conditions (when no wells at the site are pumped) using an existing ground-water-flow model, particle-tracking techniques, and the porosity values determined in this study. The average linear velocity of flow paths beginning at each contaminated well and ending at the streams where the flow paths terminate ranges from 0.08 to 130 feet per day. As a result of a change in the pump-and-treat system (adding a 165-foot-deep well pumped at 5 gallons per minute and reducing the pumping rate at a nearby 41-foot-deep well by the same amount), water in the vicinity of three 100- to 165-foot-deep wells flows to the deep well rather than the shallower well.
Technology Transfer Automated Retrieval System (TEKTRAN)
Understanding and modeling transport of solutes in porous media is a critical issue in the environmental protection. Contaminants from various industrial and agricultural sources can travel in soil and ground water and eventually affect human and animal health. The parabolic advective-dispersive equ...
Sediment transport in a surface-advected estuarine plume
NASA Astrophysics Data System (ADS)
Yao, H. Y.; Leonardi, N.; Li, J. F.; Fagherazzi, S.
2016-03-01
The interplay between suspended-sediment transport and plume hydrodynamics in a surface-advected estuarine plume is studied using a three-dimensional numerical model. Our analysis focuses on the formation of a sediment-rich alongshore current and on the effect of sediments on the structure of the recirculating freshwater bulge. We introduce the ratio Y between the traveling time of sediment along the bulge edge and the settling timescale. When Y <1, suspended sediments enter the alongshore coastal current. When Y >1 the sediments are deposited within the bulge. We find that a critical range of settling velocities exist above which no transport in the costal current is allowed. Critical settling-velocity values increase with river discharge. Therefore, low magnitude and long-lasting floods promote sediment sorting in the continental shelf. We further find that, for a given flood duration, intermediate flood magnitudes at the limit between subcritical and supercritical flow maximize the alongshore sediment transport. Similarly, for a fixed input of water and sediments, intermediate discharge durations maximize alongshore sediment transport.
Xu, Bruce S; Lollar, Barbara Sherwood; Passeport, Elodie; Sleep, Brent E
2016-04-15
Aqueous phase diffusion-related isotope fractionation (DRIF) for carbon isotopes was investigated for common groundwater contaminants in systems in which transport could be considered to be one-dimensional. This paper focuses not only on theoretically observable DRIF effects in these systems but introduces the important concept of constraining "observable" DRIF based on constraints imposed by the scale of measurements in the field, and on standard limits of detection and analytical uncertainty. Specifically, constraints for the detection of DRIF were determined in terms of the diffusive fractionation factor, the initial concentration of contaminants (C0), the method detection limit (MDL) for isotopic analysis, the transport time, and the ratio of the longitudinal mechanical dispersion coefficient to effective molecular diffusion coefficient (Dmech/Deff). The results allow a determination of field conditions under which DRIF may be an important factor in the use of stable carbon isotope measurements for evaluation of contaminant transport and transformation for one-dimensional advective-dispersive transport. This study demonstrates that for diffusion-dominated transport of BTEX, MTBE, and chlorinated ethenes, DRIF effects are only detectable for the smaller molar mass compounds such as vinyl chloride for C0/MDL ratios of 50 or higher. Much larger C0/MDL ratios, corresponding to higher source concentrations or lower detection limits, are necessary for DRIF to be detectable for the higher molar mass compounds. The distance over which DRIF is observable for VC is small (less than 1m) for a relatively young diffusive plume (<100years), and DRIF will not easily be detected by using the conventional sampling approach with "typical" well spacing (at least several meters). With contaminant transport by advection, mechanical dispersion, and molecular diffusion this study suggests that in field sites where Dmech/Deff is larger than 10, DRIF effects will likely not be
Fractional Advective-Dispersive Equation as a Model of Solute Transport in Porous Media
Technology Transfer Automated Retrieval System (TEKTRAN)
Understanding and modeling transport of solutes in porous media is a critical issue in the environmental protection. The common model is the advective-dispersive equation (ADE) describing the superposition of the advective transport and the Brownian motion in water-filled pore space. Deviations from...
NASA Astrophysics Data System (ADS)
Ancey, Christophe; Bohorquez, Patricio; Heyman, Joris
2016-04-01
The advection-diffusion equation arises quite often in the context of sediment transport, e.g., for describing time and space variations in the particle activity (the solid volume of particles in motion per unit streambed area). Stochastic models can also be used to derive this equation, with the significant advantage that they provide information on the statistical properties of particle activity. Stochastic models are quite useful when sediment transport exhibits large fluctuations (typically at low transport rates), making the measurement of mean values difficult. We develop an approach based on birth-death Markov processes, which involves monitoring the evolution of the number of particles moving within an array of cells of finite length. While the topic has been explored in detail for diffusion-reaction systems, the treatment of advection has received little attention. We show that particle advection produces nonlocal effects, which are more or less significant depending on the cell size and particle velocity. Albeit nonlocal, these effects look like (local) diffusion and add to the intrinsic particle diffusion (dispersal due to velocity fluctuations), with the important consequence that local measurements depend on both the intrinsic properties of particle displacement and the dimensions of the measurement system.
Thermally driven advection for radioxenon transport from an underground nuclear explosion
NASA Astrophysics Data System (ADS)
Sun, Yunwei; Carrigan, Charles R.
2016-05-01
Barometric pumping is a ubiquitous process resulting in migration of gases in the subsurface that has been studied as the primary mechanism for noble gas transport from an underground nuclear explosion (UNE). However, at early times following a UNE, advection driven by explosion residual heat is relevant to noble gas transport. A rigorous measure is needed for demonstrating how, when, and where advection is important. In this paper three physical processes of uncertain magnitude (oscillatory advection, matrix diffusion, and thermally driven advection) are parameterized by using boundary conditions, system properties, and source term strength. Sobol' sensitivity analysis is conducted to evaluate the importance of all physical processes influencing the xenon signals. This study indicates that thermally driven advection plays a more important role in producing xenon signals than oscillatory advection and matrix diffusion at early times following a UNE, and xenon isotopic ratios are observed to have both time and spatial dependence.
Magnetic helicity transport in the advective gauge family
Candelaresi, Simon; Brandenburg, Axel; Hubbard, Alexander; Mitra, Dhrubaditya
2011-01-15
Magnetic helicity fluxes are investigated in a family of gauges in which the contribution from ideal magnetohydrodynamics takes the form of a purely advective flux. Numerical simulations of magnetohydrodynamic turbulence in this advective gauge family exhibit instabilities triggered by the build-up of unphysical irrotational contributions to the magnetic vector potential. As a remedy, the vector potential is evolved in a numerically well behaved gauge, from which the advective vector potential is obtained by a gauge transformation. In the kinematic regime, the magnetic helicity density evolves similarly to a passive scalar when resistivity is small and turbulent mixing is mild, i.e., when the fluid Reynolds number is not too large. In the dynamical regime, resistive contributions to the magnetic helicity flux in the advective gauge are found to be significant owing to the development of small length scales in the irrotational part of the magnetic vector potential.
Analytical solution for the advection-dispersion transport equation in layered media
Technology Transfer Automated Retrieval System (TEKTRAN)
The advection-dispersion transport equation with first-order decay was solved analytically for multi-layered media using the classic integral transform technique (CITT). The solution procedure used an associated non-self-adjoint advection-diffusion eigenvalue problem that had the same form and coef...
Theory of advection-driven long range biotic transport
Technology Transfer Automated Retrieval System (TEKTRAN)
We propose a simple mechanistic model to examine the effects of advective flow on the spread of fungal diseases spread by wind-blown spores. The model is defined by a set of two coupled non-linear partial differential equations for spore densities. One equation describes the long-distance advectiv...
An exact peak capturing and essentially oscillation-free (EPCOF) algorithm, consisting of advection-dispersion decoupling, backward method of characteristics, forward node tracking, and adaptive local grid refinement, is developed to solve transport equations. This algorithm repr...
Kim, Hong; Benson, Craig H
2004-07-01
The relative contributions of four mechanisms of oxygen transport in multilayer composite (MLC) caps placed over oxygen-consuming mine waste were evaluated using numerical and analytical methods. MLC caps are defined here as caps consisting of earthen and geosynthetic (polymeric) components where a composite barrier layer consisting of a geomembrane (1-2 mm thick polymeric sheet) overlying a clay layer is the primary barrier to transport. The transport mechanisms that were considered are gas-phase advective transport, gas-phase diffusive transport, liquid-phase advective transport via infiltrating precipitation and liquid-phase diffusive transport. A numerical model was developed to simulate gas-phase advective-diffusive transport of oxygen through a multilayer cap containing seven layers. This model was also used to simulate oxygen diffusion in the liquid phase. An approximate analytical method was used to compute the advective flux of oxygen in the liquid phase. The numerical model was verified for limiting cases using an analytical solution. Comparisons were also made between model predictions and field data for earthen caps reported by others. Results of the analysis show that the dominant mechanism for oxygen transport through MLC caps is gas-phase diffusion. For the cases that were considered, the gas-phase diffusive flux typically comprises at least 99% of the total oxygen flux. Thus, designers of MLC caps should focus on design elements and features that will limit diffusion of gas-phase oxygen. PMID:15145567
NASA Astrophysics Data System (ADS)
Maghrebi, M.; Jankovic, I.; Rabideau, A. J.; Allen-King, R. M.; Weissmann, G. S.
2011-12-01
Effects of three key transport mechanisms (advection, diffusion and sorption) on transport and contaminant tailing of chlorinated solvents have been investigated using a numerical model. Thousands of model simulations have been conducted for various combinations of transport parameters that govern three key mechanisms in order to quantify tailing and relative importance of each mechanism. Hydraulic conductivity model contains a single inclusion of constant conductivity K embedded in a homogeneous anisotropic background of conductivity Kh,Kv. The inclusion is shaped as an oblate ellipsoid and subject to uniform flow. The background represents "average" conductivity of a heterogeneous formation while inclusion is used to represent geologic units that are more or less conductive than the background. The ratio of long to short semi-axis of the inclusion (a/b) models the ratio of horizontal to vertical integral scales (Ih/Iv) of different geologic units, where integral scales can be obtained, for example, using indicator variograms. The flow solution for present problem is obtained analytically as a closed form solution with exact expressions for Darcy velocity valid both inside and outside the inclusion. Sorption is modeled as an equilibrium process governed by a linear isotherm. The effects on transport and tailing are accounted for using retardation factors. Sorption heterogeneity is created by allowing different values of retardation factor for the interior (Ri) and the exterior of the inclusion (Rb). Diffusive displacements have been added to retarded advective displacements using random walk method. Peclet number, defined as Pe=U Ih/D (U is the groundwater velocity, D is the molecular diffusion coefficient for chlorinated solvents), is used to quantify the diffusion process. Very large numbers of particles (hundreds of thousands) have been tracked using very small time steps (as small as a/10,000) to provide sufficient resolution to breakthrough curves and to
General solution of a fractional diffusion-advection equation for solar cosmic-ray transport
NASA Astrophysics Data System (ADS)
Rocca, M. C.; Plastino, A. R.; Plastino, A.; Ferri, G. L.; de Paoli, A.
2016-04-01
In this effort we exactly solve the fractional diffusion-advection equation for solar cosmic-ray transport and give its general solution in terms of hypergeometric distributions. Numerical analysis of this equation shows that its solutions resemble power-laws.
Scaling of geochemical reaction rates via advective solute transport.
Hunt, A G; Ghanbarian, B; Skinner, T E; Ewing, R P
2015-07-01
Transport in porous media is quite complex, and still yields occasional surprises. In geological porous media, the rate at which chemical reactions (e.g., weathering and dissolution) occur is found to diminish by orders of magnitude with increasing time or distance. The temporal rates of laboratory experiments and field observations differ, and extrapolating from laboratory experiments (in months) to field rates (in millions of years) can lead to order-of-magnitude errors. The reactions are transport-limited, but characterizing them using standard solute transport expressions can yield results in agreement with experiment only if spurious assumptions and parameters are introduced. We previously developed a theory of non-reactive solute transport based on applying critical path analysis to the cluster statistics of percolation. The fractal structure of the clusters can be used to generate solute distributions in both time and space. Solute velocities calculated from the temporal evolution of that distribution have the same time dependence as reaction-rate scaling in a wide range of field studies and laboratory experiments, covering some 10 decades in time. The present theory thus both explains a wide range of experiments, and also predicts changes in the scaling behavior in individual systems with increasing time and/or length scales. No other theory captures these variations in scaling by invoking a single physical mechanism. Because the successfully predicted chemical reactions include known results for silicate weathering rates, our theory provides a framework for understanding changes in the global carbon cycle, including its effects on extinctions, climate change, soil production, and denudation rates. It further provides a basis for understanding the fundamental time scales of hydrology and shallow geochemistry, as well as the basis of industrial agriculture. PMID:26232976
Scaling of geochemical reaction rates via advective solute transport
NASA Astrophysics Data System (ADS)
Hunt, A. G.; Ghanbarian, B.; Skinner, T. E.; Ewing, R. P.
2015-07-01
Transport in porous media is quite complex, and still yields occasional surprises. In geological porous media, the rate at which chemical reactions (e.g., weathering and dissolution) occur is found to diminish by orders of magnitude with increasing time or distance. The temporal rates of laboratory experiments and field observations differ, and extrapolating from laboratory experiments (in months) to field rates (in millions of years) can lead to order-of-magnitude errors. The reactions are transport-limited, but characterizing them using standard solute transport expressions can yield results in agreement with experiment only if spurious assumptions and parameters are introduced. We previously developed a theory of non-reactive solute transport based on applying critical path analysis to the cluster statistics of percolation. The fractal structure of the clusters can be used to generate solute distributions in both time and space. Solute velocities calculated from the temporal evolution of that distribution have the same time dependence as reaction-rate scaling in a wide range of field studies and laboratory experiments, covering some 10 decades in time. The present theory thus both explains a wide range of experiments, and also predicts changes in the scaling behavior in individual systems with increasing time and/or length scales. No other theory captures these variations in scaling by invoking a single physical mechanism. Because the successfully predicted chemical reactions include known results for silicate weathering rates, our theory provides a framework for understanding changes in the global carbon cycle, including its effects on extinctions, climate change, soil production, and denudation rates. It further provides a basis for understanding the fundamental time scales of hydrology and shallow geochemistry, as well as the basis of industrial agriculture.
Technology Transfer Automated Retrieval System (TEKTRAN)
Analytical solutions of the advection-dispersion solute transport equation remain useful for a large number of applications in science and engineering. In this paper we extend the Duhamel theorem, originally established for diffusion type problems, to the case of advective-dispersive transport subj...
Wellman, Dawn M.; Zachara, John M.; Liu, Chongxuan; Qafoku, Nikolla; Smith, Steven C.; Forrester, Steven W.
2008-11-03
Sedimentary, hydrologic, and geochemical variations in the Hanford subsurface environment, as well as compositional differences in contaminating waste streams, have created vast differences in the migration and mobility of uranium within the subsurface environment. A series of hydraulically-saturated and -unsaturated column experiments were performed to i.) assess the effect of water content on the advective desorption and migration of uranium from contaminated sediments, and ii.) evaluate the uranium concentration that can develop in porewater and/or groundwater as a result of desorption/dissolution reactions. Flow rate and moisture content were varied to evaluate the influence of contact time, pore water velocity, and macropore desaturation on aqueous uranium concentrations. Sediments were collected from the T-TX-TY tank farm complex and the 300 Area Process Ponds located on the Hanford Site, southeastern Washington State. The sediments vary in depth, mineralogy, and in contamination events. Experiments were conducted under mildly alkaline/calcareous conditions representative of conditions commonly encountered at repository sites across the arid western United States and, in particular, the Hanford site. Results illustrate the release of uranium from these sediments is kinetically controlled and low water contents encountered within the Hanford vadose zone result in the formation of mobile-immobile water regimes, which isolate a fraction of the reactive sites within the sediments, effectively reducing the concentration of uranium released into migrating porewaters.
NASA Astrophysics Data System (ADS)
Dvoretskaya, Olga A.; Kondratenko, Peter S.
2009-04-01
We study the transport of impurity particles on a comb structure in the presence of advection. The main body concentration and asymptotic concentration distributions are obtained. Seven different transport regimes occur on the comb structure with finite teeth: classical diffusion, advection, quasidiffusion, subdiffusion, slow classical diffusion, and two kinds of slow advection. Quasidiffusion deserves special attention. It is characterized by a linear growth of the mean-square displacement. However, quasidiffusion is an anomalous transport regime. We established that a change in transport regimes in time leads to a change in regimes in space. Concentration tails have a cascade structure, namely, consisting of several parts.
NASA Astrophysics Data System (ADS)
Chauhan, R. P.; Kumar, Amit
The present work is aimed that out of diffusive and advective transport which is dominant process for indoor radon entry under normal room conditions. For this purpose the radon diffusion coefficient and permeability of concrete were measured by specially designed experimental set up. The radon diffusion coefficient of concrete was measured by continuous radon monitor. The measured value was (3.78 ± 0.39)×10-8 m2/s and found independent of the radon gas concentration in source chamber. The radon permeability of concrete varied between 1.85×10-17 to 1.36×10-15 m2 for the bulk pressure difference fewer than 20 Pa to 73.3 kPa. From the measured diffusion coefficient and absolute permeability, the radon flux from the concrete surface having concentrations gradient 12-40 kBq/m3 and typical floor thickness 0.1 m was calculated by the application of Fick and Darcy laws. Using the measured flux attributable to diffusive and advective transport, the indoor radon concentration for a typical Indian model room having dimension (5×6×7) m3 was calculated under average room ventilation (0.63 h-1). The results showed that the contribution of diffusive transport through intact concrete is dominant over the advective transport, as expected from the low values of concrete permeability.
NASA Astrophysics Data System (ADS)
Guihéneuf, N.; Boisson, A.; Bour, O.; Le Borgne, T.; Marechal, J.; Nigon, B.; Wajiddudin, M.; Ahmed, S.
2013-12-01
The prediction of transport in weathered and fractured rocks is critical as it represents the primary control of contaminant transfer from the subsurface in many parts of the world. This is the case in Southern India, where the subsurface is composed mainly of weathered and fractured granite and where the overexploitation of the groundwater resource since the 70's has led to high water table depletion and strong groundwater quality deterioration. One key issue for modelling transport in such systems is to quantify the respective role of advective heterogeneities and matrix diffusion, which can both lead to strongly non Fickian transport properties. We investigate this question by analysing tracer test experiments performed under different flow configurations at a fractured granite experimental site located in Andhra Pradesh (India). We performed both convergent and push-pull tracer tests within the same fracture and at different scales. Three convergent tracer tests were performed with a solution of fluorescein for different pumping rate and for different distances between injection and pumping boreholes: 6, 30 and 41 meters. To evaluate diffusive process, we performed two long-duration push-pull tests (push time of 3 hours) with a solution of two conservative tracers of different diffusion coefficient (fluorescein and sodium chloride). We performed also six others push-pull tests with only fluorescein but for a variable push times of 14 min and 55 min with or without resting time of about 60 min. The late-time behaviour on the breakthrough curves (BTCs) obtained for all convergent tracer tests showed a power-law slope of -2. Two of them showed an inflexion in the BTCs suggesting the existence of two independent flow paths and thus a highly channelized flow. The long-duration push-pull tests showed similar late-time behaviour with a power-law slope of -2.2 for both tracers. The six others push-pull tests showed a variation of power-law exponent from -3 to -2
Technology Transfer Automated Retrieval System (TEKTRAN)
This paper presents a formal exact solution of the linear advection-diffusion transport equation with constant coefficients for both transient and steady-state regimes. A classical mathematical substitution transforms the original advection-diffusion equation into an exclusively diffusive equation. ...
NASA Astrophysics Data System (ADS)
Künzli, Pierre; Tsunematsu, Kae; Albuquerque, Paul; Falcone, Jean-Luc; Chopard, Bastien; Bonadonna, Costanza
2016-04-01
Volcanic ash transport and dispersal models typically describe particle motion via a turbulent velocity field. Particles are advected inside this field from the moment they leave the vent of the volcano until they deposit on the ground. Several techniques exist to simulate particles in an advection field such as finite difference Eulerian, Lagrangian-puff or pure Lagrangian techniques. In this paper, we present a new flexible simulation tool called TETRAS (TEphra TRAnsport Simulator) based on a hybrid Eulerian-Lagrangian model. This scheme offers the advantages of being numerically stable with no numerical diffusion and easily parallelizable. It also allows us to output particle atmospheric concentration or ground mass load at any given time. The model is validated using the advection-diffusion analytical equation. We also obtained a good agreement with field observations of the tephra deposit associated with the 2450 BP Pululagua (Ecuador) and the 1996 Ruapehu (New Zealand) eruptions. As this kind of model can lead to computationally intensive simulations, a parallelization on a distributed memory architecture was developed. A related performance model, taking into account load imbalance, is proposed and its accuracy tested.
Numerical advection algorithms and their role in atmospheric transport and chemistry models
NASA Technical Reports Server (NTRS)
Rood, Richard B.
1987-01-01
During the last 35 years, well over 100 algorithms for modeling advection processes have been described and tested. This review summarizes the development and improvements that have taken place. The nature of the errors caused by numerical approximation to the advection equation are highlighted. Then the particular devices that have been proposed to remedy these errors are discussed. The extensive literature comparing transport algorithms is reviewed. Although there is no clear cut 'best' algorithm, several conclusions can be made. Spectral and pseudospectral techniques consistently provide the highest degree of accuracy, but expense and difficulties assuring positive mixing ratios are serious drawbacks. Schemes which consider fluid slabs bounded by grid points (volume schemes), rather than the simple specification of constituent values at the grid points, provide accurate positive definite results.
NASA Astrophysics Data System (ADS)
Ancey, C.; Bohorquez, P.; Heyman, J.
2015-12-01
The advection-diffusion equation is one of the most widespread equations in physics. It arises quite often in the context of sediment transport, e.g., for describing time and space variations in the particle activity (the solid volume of particles in motion per unit streambed area). Phenomenological laws are usually sufficient to derive this equation and interpret its terms. Stochastic models can also be used to derive it, with the significant advantage that they provide information on the statistical properties of particle activity. These models are quite useful when sediment transport exhibits large fluctuations (typically at low transport rates), making the measurement of mean values difficult. Among these stochastic models, the most common approach consists of random walk models. For instance, they have been used to model the random displacement of tracers in rivers. Here we explore an alternative approach, which involves monitoring the evolution of the number of particles moving within an array of cells of finite length. Birth-death Markov processes are well suited to this objective. While the topic has been explored in detail for diffusion-reaction systems, the treatment of advection has received no attention. We therefore look into the possibility of deriving the advection-diffusion equation (with a source term) within the framework of birth-death Markov processes. We show that in the continuum limit (when the cell size becomes vanishingly small), we can derive an advection-diffusion equation for particle activity. Yet while this derivation is formally valid in the continuum limit, it runs into difficulty in practical applications involving cells or meshes of finite length. Indeed, within our stochastic framework, particle advection produces nonlocal effects, which are more or less significant depending on the cell size and particle velocity. Albeit nonlocal, these effects look like (local) diffusion and add to the intrinsic particle diffusion (dispersal due
Benchmarking of a Markov multizone model of contaminant transport.
Jones, Rachael M; Nicas, Mark
2014-10-01
A Markov chain model previously applied to the simulation of advection and diffusion process of gaseous contaminants is extended to three-dimensional transport of particulates in indoor environments. The model framework and assumptions are described. The performance of the Markov model is benchmarked against simple conventional models of contaminant transport. The Markov model is able to replicate elutriation predictions of particle deposition with distance from a point source, and the stirred settling of respirable particles. Comparisons with turbulent eddy diffusion models indicate that the Markov model exhibits numerical diffusion in the first seconds after release, but over time accurately predicts mean lateral dispersion. The Markov model exhibits some instability with grid length aspect when turbulence is incorporated by way of the turbulent diffusion coefficient, and advection is present. However, the magnitude of prediction error may be tolerable for some applications and can be avoided by incorporating turbulence by way of fluctuating velocity (e.g. turbulence intensity). PMID:25143517
NASA Technical Reports Server (NTRS)
Leonard, B. P.
1988-01-01
A fresh approach is taken to the embarrassingly difficult problem of adequately modeling simple pure advection. An explicit conservative control-volume formation makes use of a universal limiter for transient interpolation modeling of the advective transport equations. This ULTIMATE conservative difference scheme is applied to unsteady, one-dimensional scalar pure advection at constant velocity, using three critical test profiles: an isolated sine-squared wave, a discontinuous step, and a semi-ellipse. The goal, of course, is to devise a single robust scheme which achieves sharp monotonic resolution of the step without corrupting the other profiles. The semi-ellipse is particularly challenging because of its combination of sudden and gradual changes in gradient. The ULTIMATE strategy can be applied to explicit conservation schemes of any order of accuracy. Second-order schemes are unsatisfactory, showing steepening and clipping typical of currently popular so-called high resolution shock-capturing of TVD schemes. The ULTIMATE third-order upwind scheme is highly satisfactory for most flows of practical importance. Higher order methods give predictably better step resolution, although even-order schemes generate a (monotonic) waviness in the difficult semi-ellipse simulation. Little is to be gained above ULTIMATE fifth-order upwinding which gives results close to the ultimate for which one might hope.
Advective transport and decomposition of chain-forming planktonic diatoms in permeable sediments
NASA Astrophysics Data System (ADS)
Ehrenhauss, Sandra; Huettel, Markus
2004-09-01
In laboratory chamber experiments we demonstrate that permeable sediments (>7×10 -12 m 2) exposed to boundary flows filter chain-forming coastal bloom diatoms ( Skeletonema costatum and Thalassiosira rotula) from the water column, causing rapid transfer of fresh organic particulate matter into sediment layers as deep as 5 cm within 72 h. The penetration depth of the diatoms depends on the permeability of the bed and the length of the chains. Long chains were not transported as deep into the sediment as short chains or single cells. The fast advective transfer of phytoplankton cells into sandy sediments may be an important process facilitating organic matter uptake and preventing resuspension of deposited organic material in high-energy coastal environments. High advective flushing rates in medium- and coarse-grained sandy sediments enhanced the mineralisation of the trapped diatoms (2300 to 3200 μmol C m -2 d -1), stimulated benthic oxygen consumption (2300 to 3000 μmol O 2 m -2 d -1), as well as nitrification (up to 20 μmol NO 3- m -2 d -1), relative to sediment where diffusion dominated the solute exchange. Advective solute exchange rates that increase with increasing permeability prevent the accumulation of Si(OH) 4 near the dissolving frustules and in the pore water, leading to an effective recycling of dissolved silica to the production process in the water column (95 to 101 μmol Si(OH) 4 m -2 d -1). This process may also enhance dissolution rates of the deposited opal in coarse-grained sands by maintaining higher degrees of undersaturation than in fine-grained sediments. Our results suggest that advective filtration of planktonic diatoms into permeable sediments increases mineralisation and recycling of Si(OH) 4 and organic matter in high energetic shelf areas.
Analytical solutions of a fractional diffusion-advection equation for solar cosmic-ray transport
Litvinenko, Yuri E.; Effenberger, Frederic
2014-12-01
Motivated by recent applications of superdiffusive transport models to shock-accelerated particle distributions in the heliosphere, we analytically solve a one-dimensional fractional diffusion-advection equation for the particle density. We derive an exact Fourier transform solution, simplify it in a weak diffusion approximation, and compare the new solution with previously available analytical results and with a semi-numerical solution based on a Fourier series expansion. We apply the results to the problem of describing the transport of energetic particles, accelerated at a traveling heliospheric shock. Our analysis shows that significant errors may result from assuming an infinite initial distance between the shock and the observer. We argue that the shock travel time should be a parameter of a realistic superdiffusive transport model.
Analytical Solutions of a Fractional Diffusion-advection Equation for Solar Cosmic-Ray Transport
NASA Astrophysics Data System (ADS)
Litvinenko, Yuri E.; Effenberger, Frederic
2014-12-01
Motivated by recent applications of superdiffusive transport models to shock-accelerated particle distributions in the heliosphere, we analytically solve a one-dimensional fractional diffusion-advection equation for the particle density. We derive an exact Fourier transform solution, simplify it in a weak diffusion approximation, and compare the new solution with previously available analytical results and with a semi-numerical solution based on a Fourier series expansion. We apply the results to the problem of describing the transport of energetic particles, accelerated at a traveling heliospheric shock. Our analysis shows that significant errors may result from assuming an infinite initial distance between the shock and the observer. We argue that the shock travel time should be a parameter of a realistic superdiffusive transport model.
Modeling of advection-diffusion-reaction processes using transport dissipative particle dynamics
NASA Astrophysics Data System (ADS)
Li, Zhen; Yazdani, Alireza; Tartakovsky, Alexandre; Karniadakis, George Em
2015-11-01
We present a transport dissipative particle dynamics (tDPD) model for simulating mesoscopic problems involving advection-diffusion-reaction (ADR) processes, along with a methodology for implementation of the correct Dirichlet and Neumann boundary conditions in tDPD simulations. In particular, the transport of concentration is modeled by a Fickian flux and a random flux between tDPD particles, and the advection is implicitly considered by the movements of Lagrangian particles. To validate the proposed tDPD model and the boundary conditions, three benchmark simulations of one-dimensional diffusion with different boundary conditions are performed, and the results show excellent agreement with the theoretical solutions. Also, two-dimensional simulations of ADR systems are performed and the tDPD simulations agree well with the results obtained by the spectral element method. Finally, an application of tDPD to the spatio-temporal dynamics of blood coagulation involving twenty-five reacting species is performed to demonstrate the promising biological applications of the tDPD model. Supported by the DOE Center on Mathematics for Mesoscopic Modeling of Materials (CM4) and an INCITE grant.
Simulation of the advective methane transport and AOM in Shenhu area, the Northern South China Sea
NASA Astrophysics Data System (ADS)
Liu, L.; Wu, N.
2012-04-01
Anaerobic Oxidation of Methane (AOM) occurs in the transition zone between the presence of sulfate and methane. This reaction is an important process for methane and the global carbon cycle. Methane gas hydrates bearing sediments were recovered in Shenhu Area, the Northern South China Sea, and methane advective transport was detected in this area as well. A one dimension numerical simulation tool was implemented to study the AOM process combined with the advective methane transport in Shenhu Area according to the local drilling data and geochemical information. The modeled results suggest that local methane flux will be consumed in the sediment column via dissolution, sorption and AOM reaction. A portion of methane will enter water column and possibly atmosphere if the methane flux was one order of magnitude higher than current level. Furthermore, the calculated rates of AOM in Shenhu area range similar to that of gas hydrate mounds in Mexico Golf. However, AOM is ability to consume more methane than that in Golf of Mexico due to the lower permeable sediment associated with a deeper sulfate methane transition layer.
NASA Astrophysics Data System (ADS)
Gillibrand, P. A.; Herzfeld, M.
2016-05-01
We present a flux-form semi-Lagrangian (FFSL) advection scheme designed for offline scalar transport simulation with coastal ocean models using curvilinear horizontal coordinates. The scheme conserves mass, overcoming problems of mass conservation typically experienced with offline transport models, and permits long time steps (relative to the Courant number) to be used by the offline model. These attributes make the method attractive for offline simulation of tracers in biogeochemical or sediment transport models using archived flow fields from hydrodynamic models. We describe the FFSL scheme, and test it on two idealised domains and one real domain, the Great Barrier Reef in Australia. For comparison, we also include simulations using a traditional semi-Lagrangian advection scheme for the offline simulations. We compare tracer distributions predicted by the offline FFSL transport scheme with those predicted by the original hydrodynamic model, assess the conservation of mass in all cases and contrast the computational efficiency of the schemes. We find that the FFSL scheme produced very good agreement with the distributions of tracer predicted by the hydrodynamic model, and conserved mass with an error of a fraction of one percent. In terms of computational speed, the FFSL scheme was comparable with the semi-Lagrangian method and an order of magnitude faster than the full hydrodynamic model, even when the latter ran in parallel on multiple cores. The FFSL scheme presented here therefore offers a viable mass-conserving and computationally-efficient alternative to traditional semi-Lagrangian schemes for offline scalar transport simulation in coastal models.
Exact PDF equations and closure approximations for advective-reactive transport
Venturi, D.; Tartakovsky, Daniel M.; Tartakovsky, Alexandre M.; Karniadakis, George E.
2013-06-01
Mathematical models of advection–reaction phenomena rely on advective flow velocity and (bio) chemical reaction rates that are notoriously random. By using functional integral methods, we derive exact evolution equations for the probability density function (PDF) of the state variables of the advection–reaction system in the presence of random transport velocity and random reaction rates with rather arbitrary distributions. These PDF equations are solved analytically for transport with deterministic flow velocity and a linear reaction rate represented mathematically by a heterog eneous and strongly-correlated random field. Our analytical solution is then used to investigate the accuracy and robustness of the recently proposed large-eddy diffusivity (LED) closure approximation [1]. We find that the solution to the LED-based PDF equation, which is exact for uncorrelated reaction rates, is accurate even in the presence of strong correlations and it provides an upper bound of predictive uncertainty.
NASA Technical Reports Server (NTRS)
Zhou, YE; Vahala, George
1993-01-01
The advection of a passive scalar by incompressible turbulence is considered using recursive renormalization group procedures in the differential sub grid shell thickness limit. It is shown explicitly that the higher order nonlinearities induced by the recursive renormalization group procedure preserve Galilean invariance. Differential equations, valid for the entire resolvable wave number k range, are determined for the eddy viscosity and eddy diffusivity coefficients, and it is shown that higher order nonlinearities do not contribute as k goes to 0, but have an essential role as k goes to k(sub c) the cutoff wave number separating the resolvable scales from the sub grid scales. The recursive renormalization transport coefficients and the associated eddy Prandtl number are in good agreement with the k-dependent transport coefficients derived from closure theories and experiments.
Space shuttle exhaust plumes in the lower thermosphere: Advective transport and diffusive spreading
NASA Astrophysics Data System (ADS)
Stevens, Michael H.; Lossow, Stefan; Siskind, David E.; Meier, R. R.; Randall, Cora E.; Russell, James M.; Urban, Jo; Murtagh, Donal
2014-02-01
The space shuttle main engine plume deposited between 100 and 115 km altitude is a valuable tracer for global-scale dynamical processes. Several studies have shown that this plume can reach the Arctic or Antarctic to form bursts of polar mesospheric clouds (PMCs) within a few days. The rapid transport of the shuttle plume is currently not reproduced by general circulation models and is not well understood. To help delineate the issues, we present the complete satellite datasets of shuttle plume observations by the Sounding of the Atmosphere using Broadband Emission Radiometry instrument and the Sub-Millimeter Radiometer instrument. From 2002 to 2011 these two instruments observed 27 shuttle plumes in over 600 limb scans of water vapor emission, from which we derive both advective meridional transport and diffusive spreading. Each plume is deposited at virtually the same place off the United States east coast so our results are relevant to northern mid-latitudes. We find that the advective transport for the first 6-18 h following deposition depends on the local time (LT) of launch: shuttle plumes deposited later in the day (~13-22 LT) typically move south whereas they otherwise typically move north. For these younger plumes rapid transport is most favorable for launches at 6 and 18 LT, when the displacement is 10° in latitude corresponding to an average wind speed of 30 m/s. For plumes between 18 and 30 h old some show average sustained meridional speeds of 30 m/s. For plumes between 30 and 54 h old the observations suggest a seasonal dependence to the meridional transport, peaking near the beginning of year at 24 m/s. The diffusive spreading of the plume superimposed on the transport is on average 23 m/s in 24 h. The plume observations show large variations in both meridional transport and diffusive spreading so that accurate modeling requires knowledge of the winds specific to each case. The combination of transport and spreading from the STS-118 plume in August
NASA Technical Reports Server (NTRS)
Waugh, Darryn W.; Plumb, R. Alan
1994-01-01
We present a trajectory technique, contour advection with surgery (CAS), for tracing the evolution of material contours in a specified (including observed) evolving flow. CAS uses the algorithms developed by Dritschel for contour dynamics/surgery to trace the evolution of specified contours. The contours are represented by a series of particles, which are advected by a specified, gridded, wind distribution. The resolution of the contours is preserved by continually adjusting the number of particles, and finescale features are produced that are not present in the input data (and cannot easily be generated using standard trajectory techniques). The reliability, and dependence on the spatial and temporal resolution of the wind field, of the CAS procedure is examined by comparisons with high-resolution numerical data (from contour dynamics calculations and from a general circulation model), and with routine stratospheric analyses. These comparisons show that the large-scale motions dominate the deformation field and that CAS can accurately reproduce small scales from low-resolution wind fields. The CAS technique therefore enables examination of atmospheric tracer transport at previously unattainable resolution.
Technology Transfer Automated Retrieval System (TEKTRAN)
Contaminant transport processes in streams, rivers, and other surface water bodies can be analyzed or predicted using the advection-dispersion equation and related transport models. In part 1 of this two-part series we presented a large number of one- and multi-dimensional analytical solutions of t...
Mass transport contamination study
NASA Technical Reports Server (NTRS)
Robertson, S. J.
1972-01-01
A theoretical analysis was performed to determine the effects of outgassing and waste dumping on the contamination field around an orbiting spacecraft. The spacecraft was assumed to be spherical in shape with the mass flow emitting uniformly from the spherical surface at a constant rate and in a D'Lambertian spatial distribution. The outflow of gases were assumed to be neutrally charged and of a single species with a molecular weight characteristic of a composite of the actual species involved in the mass flow. The theoretical analysis showed that, for outgassing only, less than 1.5 percent of the outgas products will return to the Skylab spacecraft as a result of intermolecular collisions. When the total mass flow from the spacecraft, including waste dumps and reaction control motor firings, was considered, it was estimated that about 30 percent will return to the spacecraft.
Subsurface Flow and Contaminant Transport
Energy Science and Technology Software Center (ESTSC)
2000-09-19
FACT is a transient three-dimensional, finite element code for simulating isothermal groundwater flow, moisture movement, and solute transport in variably and/or fully saturated subsurface porous media. Both single and dual-domain transport formulations are available. Transport mechanisms considered include advection, hydrodynamic dispersion, linear adsorption, mobile/immobile mass transfer and first-order degradation. A wide range of acquifier conditions and remediation systems commonly encountered in the field can be simulated. Notable boundary condition (BC) options include, a combined rechargemore » and drain BC for simulating recirculation wells, and a head dependent well BC that computes flow based on specified drawdown. The code is designed to handle highly heterogenous, multi-layer, acquifer systems in a numerically efficient manner. Subsurface structure is represented with vertically distorted rectangular brick elements in a Cartesian system. The groundwater flow equation is approximated using the Bubnov-Galerkin finite element method in conjunction with an efficient symmetric Preconditioned Conjugate Gradient (PCG) ICCG matrix solver. The solute transport equation is approximated using an upstream weighted residual finite element method designed to alleviate numerical oscillation. An efficient asymmetric PCG (ORTHOMIN) matrix solver is employed for transport. For both the flow and transport equations, element matrices are computed from either influence coefficient formulas for speed, or two point Gauss-Legendre quadrature for accuracy. Non-linear flow problems can be solved using either Newton-Ralphson linearization or Picard iteration, with under-relaxation formulas to further enhance convergence. Dynamic memory allocation is implemented using Fortran 90 constructs. FACT coding is clean and modular.« less
Transport dissipative particle dynamics model for mesoscopic advection-diffusion-reaction problems.
Li, Zhen; Yazdani, Alireza; Tartakovsky, Alexandre; Karniadakis, George Em
2015-07-01
We present a transport dissipative particle dynamics (tDPD) model for simulating mesoscopic problems involving advection-diffusion-reaction (ADR) processes, along with a methodology for implementation of the correct Dirichlet and Neumann boundary conditions in tDPD simulations. tDPD is an extension of the classic dissipative particle dynamics (DPD) framework with extra variables for describing the evolution of concentration fields. The transport of concentration is modeled by a Fickian flux and a random flux between tDPD particles, and the advection is implicitly considered by the movements of these Lagrangian particles. An analytical formula is proposed to relate the tDPD parameters to the effective diffusion coefficient. To validate the present tDPD model and the boundary conditions, we perform three tDPD simulations of one-dimensional diffusion with different boundary conditions, and the results show excellent agreement with the theoretical solutions. We also performed two-dimensional simulations of ADR systems and the tDPD simulations agree well with the results obtained by the spectral element method. Finally, we present an application of the tDPD model to the dynamic process of blood coagulation involving 25 reacting species in order to demonstrate the potential of tDPD in simulating biological dynamics at the mesoscale. We find that the tDPD solution of this comprehensive 25-species coagulation model is only twice as computationally expensive as the conventional DPD simulation of the hydrodynamics only, which is a significant advantage over available continuum solvers. PMID:26156459
Transport dissipative particle dynamics model for mesoscopic advection-diffusion-reaction problems
NASA Astrophysics Data System (ADS)
Li, Zhen; Yazdani, Alireza; Tartakovsky, Alexandre; Karniadakis, George Em
2015-07-01
We present a transport dissipative particle dynamics (tDPD) model for simulating mesoscopic problems involving advection-diffusion-reaction (ADR) processes, along with a methodology for implementation of the correct Dirichlet and Neumann boundary conditions in tDPD simulations. tDPD is an extension of the classic dissipative particle dynamics (DPD) framework with extra variables for describing the evolution of concentration fields. The transport of concentration is modeled by a Fickian flux and a random flux between tDPD particles, and the advection is implicitly considered by the movements of these Lagrangian particles. An analytical formula is proposed to relate the tDPD parameters to the effective diffusion coefficient. To validate the present tDPD model and the boundary conditions, we perform three tDPD simulations of one-dimensional diffusion with different boundary conditions, and the results show excellent agreement with the theoretical solutions. We also performed two-dimensional simulations of ADR systems and the tDPD simulations agree well with the results obtained by the spectral element method. Finally, we present an application of the tDPD model to the dynamic process of blood coagulation involving 25 reacting species in order to demonstrate the potential of tDPD in simulating biological dynamics at the mesoscale. We find that the tDPD solution of this comprehensive 25-species coagulation model is only twice as computationally expensive as the conventional DPD simulation of the hydrodynamics only, which is a significant advantage over available continuum solvers.
Transport dissipative particle dynamics model for mesoscopic advection-diffusion-reaction problems
Yazdani, Alireza; Tartakovsky, Alexandre; Karniadakis, George Em
2015-01-01
We present a transport dissipative particle dynamics (tDPD) model for simulating mesoscopic problems involving advection-diffusion-reaction (ADR) processes, along with a methodology for implementation of the correct Dirichlet and Neumann boundary conditions in tDPD simulations. tDPD is an extension of the classic dissipative particle dynamics (DPD) framework with extra variables for describing the evolution of concentration fields. The transport of concentration is modeled by a Fickian flux and a random flux between tDPD particles, and the advection is implicitly considered by the movements of these Lagrangian particles. An analytical formula is proposed to relate the tDPD parameters to the effective diffusion coefficient. To validate the present tDPD model and the boundary conditions, we perform three tDPD simulations of one-dimensional diffusion with different boundary conditions, and the results show excellent agreement with the theoretical solutions. We also performed two-dimensional simulations of ADR systems and the tDPD simulations agree well with the results obtained by the spectral element method. Finally, we present an application of the tDPD model to the dynamic process of blood coagulation involving 25 reacting species in order to demonstrate the potential of tDPD in simulating biological dynamics at the mesoscale. We find that the tDPD solution of this comprehensive 25-species coagulation model is only twice as computationally expensive as the conventional DPD simulation of the hydrodynamics only, which is a significant advantage over available continuum solvers. PMID:26156459
Transport and Recruitment of Blue Crab Larvae:a Model with Advection and Mortality
NASA Astrophysics Data System (ADS)
Garvine, R. W.; Epifanio, C. E.; Epifanio, C. C.; Wong, K.-C.
1997-07-01
The present paper develops a mathematical model for the transport and recruitment of blue crab (Callinectes sapidus) larvae, and applies it to the inner continental shelf of the Middle Atlantic Bight near Delaware Bay, U.S.A. Blue crab larvae develop through seven or eight planktonic zoeal stages to a megalopa stage suitable for recruitment to adult populations of east coast estuaries. The larvae are concentrated near the surface, and the currents are primarily forced by alongshelf winds and river discharge through major estuaries. Model currents are prescribed based on a realistic synthesis of their observed relationship to wind and river discharge. Besides the resulting advection, particle diffusion and biological mortality are added to determine the fate of larvae released from their parent estuary. Groups of particles were released across the source region of the outflowing buoyancy-driven current in the model estuary mouth. Most larvae were swept alongshelf to the south with the buoyancy-driven coastal current, and thus were lost as recruits to the population of their parent estuary. However, some larvae released close to the seaward edge of the emerging coastal current were able to cross the coastal current front and move seaward into inner shelf water during upwelling-favorable (northward) wind events. Some of these, in turn, were suitably placed near the parent estuary mouth so that they could be advected landward as megalopae into the estuary during a subsequent downwelling-favorable (southward) wind event and thus join the adult population. The model results for megalopae returns were computed from consecutive daily release of 1000 particles, and were compared with 4 years of blue crab megalopa settlement data for Delaware Bay. The model results for 1989 and 1990 matched the observed data remarkably well, with both years showing dominance by a single return event of a few days duration. For 1991 and 1992, the observed results showed multiple return events
NASA Astrophysics Data System (ADS)
Fan, Niannian; Singh, Arvind; Guala, Michele; Foufoula-Georgiou, Efi; Wu, Baosheng
2016-04-01
Bed load transport is a highly stochastic, multiscale process, where particle advection and diffusion regimes are governed by the dynamics of each sediment grain during its motion and resting states. Having a quantitative understanding of the macroscale behavior emerging from the microscale interactions is important for proper model selection in the absence of individual grain-scale observations. Here we develop a semimechanistic sediment transport model based on individual particle dynamics, which incorporates the episodic movement (steps separated by rests) of sediment particles and study their macroscale behavior. By incorporating different types of probability distribution functions (PDFs) of particle resting times Tr, under the assumption of thin-tailed PDF of particle velocities, we study the emergent behavior of particle advection and diffusion regimes across a wide range of spatial and temporal scales. For exponential PDFs of resting times Tr, we observe normal advection and diffusion at long time scales. For a power-law PDF of resting times (i.e., f>(Tr>)˜Tr-ν), the tail thickness parameter ν is observed to affect the advection regimes (both sub and normal advective), and the diffusion regimes (both subdiffusive and superdiffusive). By comparing our semimechanistic model with two random walk models in the literature, we further suggest that in order to reproduce accurately the emerging diffusive regimes, the resting time model has to be coupled with a particle motion model able to produce finite particle velocities during steps, as the episodic model discussed here.
Technology Transfer Automated Retrieval System (TEKTRAN)
The classical model to describe solute transport in soil is based on the advective-dispersive equation where Fick’s law is used to explain dispersion. From the microscopic point of view this is equivalent to consider that the motion of the particles of solute may be simulated by the Brownian motion....
The predictability of advection-dominated flux-transport solar dynamo models
Sanchez, Sabrina; Fournier, Alexandre; Aubert, Julien
2014-01-20
Space weather is a matter of practical importance in our modern society. Predictions of forecoming solar cycles mean amplitude and duration are currently being made based on flux-transport numerical models of the solar dynamo. Interested in the forecast horizon of such studies, we quantify the predictability window of a representative, advection-dominated, flux-transport dynamo model by investigating its sensitivity to initial conditions and control parameters through a perturbation analysis. We measure the rate associated with the exponential growth of an initial perturbation of the model trajectory, which yields a characteristic timescale known as the e-folding time τ {sub e}. The e-folding time is shown to decrease with the strength of the α-effect, and to increase with the magnitude of the imposed meridional circulation. Comparing the e-folding time with the solar cycle periodicity, we obtain an average estimate for τ {sub e} equal to 2.76 solar cycle durations. From a practical point of view, the perturbations analyzed in this work can be interpreted as uncertainties affecting either the observations or the physical model itself. After reviewing these, we discuss their implications for solar cycle prediction.
Transport dissipative particle dynamics model for mesoscopic advection- diffusion-reaction problems
Zhen, Li; Yazdani, Alireza; Tartakovsky, Alexandre M.; Karniadakis, George E.
2015-07-07
We present a transport dissipative particle dynamics (tDPD) model for simulating mesoscopic problems involving advection-diffusion-reaction (ADR) processes, along with a methodology for implementation of the correct Dirichlet and Neumann boundary conditions in tDPD simulations. tDPD is an extension of the classic DPD framework with extra variables for describing the evolution of concentration fields. The transport of concentration is modeled by a Fickian flux and a random flux between particles, and an analytical formula is proposed to relate the mesoscopic concentration friction to the effective diffusion coefficient. To validate the present tDPD model and the boundary conditions, we perform three tDPD simulations of one-dimensional diffusion with different boundary conditions, and the results show excellent agreement with the theoretical solutions. We also performed two-dimensional simulations of ADR systems and the tDPD simulations agree well with the results obtained by the spectral element method. Finally, we present an application of the tDPD model to the dynamic process of blood coagulation involving 25 reacting species in order to demonstrate the potential of tDPD in simulating biological dynamics at the mesoscale. We find that the tDPD solution of this comprehensive 25-species coagulation model is only twice as computationally expensive as the DPD simulation of the hydrodynamics only, which is a significant advantage over available continuum solvers.
Lewis, F.M.; Voss, C.I.; Rubin, Jacob
1986-01-01
A model was developed that can simulate the effect of certain chemical and sorption reactions simultaneously among solutes involved in advective-dispersive transport through porous media. The model is based on a methodology that utilizes physical-chemical relationships in the development of the basic solute mass-balance equations; however, the form of these equations allows their solution to be obtained by methods that do not depend on the chemical processes. The chemical environment is governed by the condition of local chemical equilibrium, and may be defined either by the linear sorption of a single species and two soluble complexation reactions which also involve that species, or binary ion exchange and one complexation reaction involving a common ion. Partial differential equations that describe solute mass balance entirely in the liquid phase are developed for each tenad (a chemical entity whose total mass is independent of the reaction process) in terms of their total dissolved concentration. These equations are solved numerically in two dimensions through the modification of an existing groundwater flow/transport computer code. (Author 's abstract)
Experimental study of advective-diffusive gaseous CO2 transport through porous media
NASA Astrophysics Data System (ADS)
Basirat, Farzad; Sharma, Prabhakar; Niemi, Auli; Fagerlund, Fritjof
2014-05-01
Leakage of gaseous CO2 into the shallow subsurface system is one of the main concerns associated with geologic storage resources. A better understanding of CO2 leakage in the shallow subsurface plays an important role for developing leakage monitoring programs. CO2 may reach the unsaturated zone by different leak mechanisms such as exsolution from CO2 supersaturated water and continuous bubbling or gas flow along a leakage path. In the unsaturated zone, the CO2 is heavier than air and may accumulate below the ground surface and move laterally. We developed a small-scale experiment setup to study the possible gaseous CO2 transport mechanisms with different controlled conditions. In this study, the experiment setup was applied to measure CO2 distributions in time and space through homogenous dry sand in which the CO2 concentrations through the domain were measured by sensitive gas sensors. The preliminary analysis of the result suggests that the transport and distribution of gaseous CO2 is spatially and temporally sensitive for the selected experimental conditions of gas flow rate and porous media. To better understand the advection and diffusion processes through the unsaturated zone, the experimental results are coupled with the dusty gas model (DGM) of Mason et al. (1967). The dusty gas model's constitutive relationships are integrated into a numerical model for multicomponent gas mixture flow and transport in porous media. The DGM considers interactions between all gaseous species and Knudsen diffusion which is important in fine grained soils. Results from the applied model were consistent with the experimental breakthrough curves obtained in this study.
McGraw R.
2012-03-01
Moment methods are finding increasing usage for simulations of particle population balance in box models and in more complex flows including two-phase flows. These highly efficient methods have nevertheless had little impact to date for multi-moment representation of aerosols and clouds in atmospheric models. There are evidently two reasons for this: First, atmospheric models, especially if the goal is to simulate climate, tend to be extremely complex and take many man-years to develop. Thus there is considerable inertia to the implementation of novel approaches. Second, and more fundamental, the nonlinear transport algorithms designed to reduce numerical diffusion during advection of various species (tracers) from cell to cell, in the typically coarse grid arrays of these models, can and occasionally do fail to preserve correlations between the moments. Other correlated tracers such as isotopic abundances, composition of aerosol mixtures, hydrometeor phase, etc., are subject to this same fate. In the case of moments, this loss of correlation can and occasionally does give rise to unphysical moment sets. When this happens the simulation can come to a halt. Following a brief description and review of moment methods, the goal of this paper is to present two new approaches that both test moment sequences for validity and correct them when they fail. The new approaches work on individual grid cells without requiring stored information from previous time-steps or neighboring cells.
Reliability analysis of contaminant transport in saturated porous media
Jang, Yeon-Soo; Sitar, N.; Der Kiureghian, A. )
1994-08-01
An approach to probabilistic analysis of contaminant transport based on first- and second-order reliability methods (FORM and SORM) is presented. In addition, system reliability methodology is introduced for the analysis of problems with more than one limit state function. Conventional one-dimensional finite difference and two-dimensional finite element models are coupled with the FORM and SORM algorithms to perform reliability analyses of advection-dominated contaminant transport. A comparison of the results of FORM and SORM analyses with the results of Monte Carlo simulations shows that FORM tends to overestimate the probability of exceedence in spatially variable domains. However, SORM accounts for the nonlinearity of the limit state surface and remains accurate, giving results consistent with Monte Carlo simulation. Finally, while the analyses presented here considered relatively simple problems, the methodology is shown to have the necessary flexibility for application to problems of practical interest.
Anderman, Evan R.; Hill, Mary Catherine
2001-01-01
Observations of the advective component of contaminant transport in steady-state flow fields can provide important information for the calibration of ground-water flow models. This report documents the Advective-Transport Observation (ADV2) Package, version 2, which allows advective-transport observations to be used in the three-dimensional ground-water flow parameter-estimation model MODFLOW-2000. The ADV2 Package is compatible with some of the features in the Layer-Property Flow and Hydrogeologic-Unit Flow Packages, but is not compatible with the Block-Centered Flow or Generalized Finite-Difference Packages. The particle-tracking routine used in the ADV2 Package duplicates the semi-analytical method of MODPATH, as shown in a sample problem. Particles can be tracked in a forward or backward direction, and effects such as retardation can be simulated through manipulation of the effective-porosity value used to calculate velocity. Particles can be discharged at cells that are considered to be weak sinks, in which the sink applied does not capture all the water flowing into the cell, using one of two criteria: (1) if there is any outflow to a boundary condition such as a well or surface-water feature, or (2) if the outflow exceeds a user specified fraction of the cell budget. Although effective porosity could be included as a parameter in the regression, this capability is not included in this package. The weighted sum-of-squares objective function, which is minimized in the Parameter-Estimation Process, was augmented to include the square of the weighted x-, y-, and z-components of the differences between the simulated and observed advective-front locations at defined times, thereby including the direction of travel as well as the overall travel distance in the calibration process. The sensitivities of the particle movement to the parameters needed to minimize the objective function are calculated for any particle location using the exact sensitivity
NASA Astrophysics Data System (ADS)
Comolli, Alessandro; Moussey, Charlie; Dentz, Marco
2016-04-01
Transport processes in groundwater systems are strongly affected by the presence of heterogeneity. The heterogeneity leads to non-Fickian features, that manifest themselves in the heavy-tailed breakthrough curves, as well as in the non-linear growth of the mean squared displacement and in the non-Gaussian plumes of solute particles. The causes of non-Fickian transport can be the heterogeneity in the flow fields and the processes of mass exchange between mobile and immobile phases, such as sorption/desorption reactions and diffusive mass transfer. Here, we present a Continuous Time Random Walk (CTRW) model that describes the transport of solutes in d-dimensional systems by taking into account both heterogeneous advection and mobile-immobile mass transfer. In order to account for these processes in the CTRW, the heterogeneities are mapped onto a distribution of transition times, which can be decomposed into advective transition times and trapping times, the latter being treated as a compound Poisson process. While advective transition times are related to the Eulerian flow velocities and, thus, to the conductivity distribution, trapping times depend on the sorption/desorption time scale, in case of reactive problems, or on the distribution of diffusion times in the immobile zones. Since the trapping time scale is typically much larger than the advective time scale, we observe the existence of two temporal regimes. The pre-asymptotic regime is defined by a characteristic time scale at which the properties of transport are fully determined by the heterogeneity of the advective field. On the other hand, in the asymptotic regime both the heterogeneity and the mass exchange processes play a role in conditioning the behaviour of transport. We consider different scenarios to discuss the relative importance of the advective heterogeneity and the mass transfer for the occurrence of non-Fickian transport. For each case we calculate analytically the scalings of the breakthrough
NASA Astrophysics Data System (ADS)
Bullen, T.; Izbicki, J.
2007-12-01
groundwater samples having <4 ppb Cr(VI), taken to be representative of regional groundwater, and the contaminated water do not pass through the remainder of the data, discounting a simple advective mixing scenario. We hypothesize a more likely scenario that involves both Cr(VI) reduction and advective mixing. As the plume initially expands downgradient, Cr(VI) in water at the leading edge encounters reductant in the aquifer resulting in limited Cr(VI) reduction. As a result of reduction, δ53Cr of Cr(VI) remaining in solution at the leading edge increases along the "reduction" trend from 0 to ~+2‰. Inevitable mixing of this water at the leading edge with regional groundwater results in a suitable mixing end-member to combine with Cr(VI) within the plume in order to explain the bulk of the remaining data. Neither Cr(VI) reduction nor advective mixing of plume and regional groundwaters can explain the data on their own, implying an interplay of at least these two processes during plume evolution. Ellis, A.S., Johnson, T.M. and Bullen, T.D. 2002, Science, 295, 2060-2062.
Advective transport observations with MODPATH-OBS--documentation of the MODPATH observation process
Hanson, R.T.; Kauffman, L.K.; Hill, M.C.; Dickinson, J.E.; Mehl, S.W.
2013-01-01
The MODPATH-OBS computer program described in this report is designed to calculate simulated equivalents for observations related to advective groundwater transport that can be represented in a quantitative way by using simulated particle-tracking data. The simulated equivalents supported by MODPATH-OBS are (1) distance from a source location at a defined time, or proximity to an observed location; (2) time of travel from an initial location to defined locations, areas, or volumes of the simulated system; (3) concentrations used to simulate groundwater age; and (4) percentages of water derived from contributing source areas. Although particle tracking only simulates the advective component of conservative transport, effects of non-conservative processes such as retardation can be approximated through manipulation of the effective-porosity value used to calculate velocity based on the properties of selected conservative tracers. This program can also account for simple decay or production, but it cannot account for diffusion. Dispersion can be represented through direct simulation of subsurface heterogeneity and the use of many particles. MODPATH-OBS acts as a postprocessor to MODPATH, so that the sequence of model runs generally required is MODFLOW, MODPATH, and MODPATH-OBS. The version of MODFLOW and MODPATH that support the version of MODPATH-OBS presented in this report are MODFLOW-2005 or MODFLOW-LGR, and MODPATH-LGR. MODFLOW-LGR is derived from MODFLOW-2005, MODPATH 5, and MODPATH 6 and supports local grid refinement. MODPATH-LGR is derived from MODPATH 5. It supports the forward and backward tracking of particles through locally refined grids and provides the output needed for MODPATH_OBS. For a single grid and no observations, MODPATH-LGR results are equivalent to MODPATH 5. MODPATH-LGR and MODPATH-OBS simulations can use nearly all of the capabilities of MODFLOW-2005 and MODFLOW-LGR; for example, simulations may be steady-state, transient, or a combination
Rigorous upper bounds for transport due to passive advection by inhomogeneous turbulence
Krommes, J.A.; Smith, R.A.
1987-05-01
A variational procedure, due originally to Howard and explored by Busse and others for self-consistent turbulence problems, is employed to determine rigorous upper bounds for the advection of a passive scalar through an inhomogeneous turbulent slab with arbitrary generalized Reynolds number R and Kubo number K. In the basic version of the method, the steady-state energy balance is used as a constraint; the resulting bound, though rigorous, is independent of K. A pedagogical reference model (one dimension, K = infinity) is described in detail; the bound compares favorably with the exact solution. The direct-interaction approximation is also worked out for this model; it is somewhat more accurate than the bound, but requires considerably more labor to solve. For the basic bound, a general formalism is presented for several dimensions, finite correlation length, and reasonably general boundary conditions. Part of the general method, in which a Green's function technique is employed, applies to self-consistent as well as to passive problems, and thereby generalizes previous results in the fluid literature. The formalism is extended for the first time to include time-dependent constraints, and a bound is deduced which explicitly depends on K and has the correct physical scalings in all regimes of R and K. Two applications from the theory of turbulent plasmas ae described: flux in velocity space, and test particle transport in stochastic magnetic fields. For the velocity space problem the simplest bound reproduces Dupree's original scaling for the strong turbulence diffusion coefficient. For the case of stochastic magnetic fields, the scaling of the bounds is described for the magnetic diffusion coefficient as well as for the particle diffusion coefficient in the so-called collisionless, fluid, and double-streaming regimes.
Ackerman, D.J.
1995-01-01
Quantitative estimates of ground-water flow directions and traveltimes for advective flow were developed for the regional aquifer system of the eastern Snake River Plain, Idaho. The work included: (1) descriptions of compartments in the aquifer that function as intermediate and regional flow systems, (2) descriptions of pathlines for flow originating at or near the water table, and (3) quantitative estimates of traveltimes for advective transport originating at or near the water table. A particle-tracking postprocessing program was used to compute pathlines on the basis of output from an existing three-dimensional steady-state flow model. The flow model uses 1980 conditions to approximate average annual conditions for 1950-80. The advective transport model required additional information about the nature of flow across model boundaries, aquifer thickness, and porosity. Porosity of two types of basalt strata has been reported for more than 1,500 individual cores from test holes, wells, and outcrops near the south side of the Idaho National Engineering Laboratory. The central 80 percent of samples had porosities of 0.08 to 0.25, the central 50 percent of samples, O. 11 to 0.21. Calibration of the model involved choosing a value for porosity that yielded the best solution. Two radiologic contaminants, iodine-129 and tritium, both introduced to the flow system about 40 years ago, are relatively conservative tracers. Iodine- 129 was considered to be more useful because of a lower analytical detection limit, longer half-life, and longer flow path. The calibration value for porosity was 0.21. Most flow in the aquifer is contained within a regional-scale compartment and follows paths that discharge to the Snake River downstream from Milner Dam. Two intermediate-scale compartments exist along the southeast side of the aquifer and near Mud Lake.One intermediate-scale compartment along the southeast side of the aquifer discharges to the Snake River near American Fails
Monger, Gregg R.; Duncan, Candice Morrison; Brusseau, Mark L.
2015-01-01
A gas-phase tracer test (GTT) was conducted at a landfill in Tucson, AZ, to help elucidate the impact of landfill gas generation on the transport and fate of chlorinated aliphatic volatile organic contaminants (VOCs). Sulfur hexafluoride (SF6) was used as the non-reactive gas tracer. Gas samples were collected from a multiport monitoring well located 15.2 m from the injection well, and analyzed for SF6, CH4, CO2, and VOCs. The travel times determined for SF6 from the tracer test are approximately two to ten times smaller than estimated travel times that incorporate transport by only gas-phase diffusion. In addition, significant concentrations of CH4 and CO2 were measured, indicating production of landfill gas. Based on these results, it is hypothesized that the enhanced rates of transport observed for SF6 are caused by advective transport associated with landfill gas generation. The rates of transport varied vertically, which is attributed to multiple factors including spatial variability of water content, refuse mass, refuse permeability, and gas generation. PMID:26380532
NASA Technical Reports Server (NTRS)
Shia, Run-Lie; Ha, Yuk Lung; Wen, Jun-Shan; Yung, Yuk L.
1990-01-01
Extensive testing of the advective scheme proposed by Prather (1986) has been carried out in support of the California Institute of Technology-Jet Propulsion Laboratory two-dimensional model of the middle atmosphere. The original scheme is generalized to include higher-order moments. In addition, it is shown how well the scheme works in the presence of chemistry as well as eddy diffusion. Six types of numerical experiments including simple clock motion and pure advection in two dimensions have been investigated in detail. By comparison with analytic solutions, it is shown that the new algorithm can faithfully preserve concentration profiles, has essentially no numerical diffusion, and is superior to a typical fourth-order finite difference scheme.
A study of turbulent transport of an advective nature in a fluid plasma
NASA Astrophysics Data System (ADS)
Min, Byunghoon; An, Chan-Yong; Kim, Chang-Bae
2014-08-01
The advective nature of the electrostatic turbulent flux of plasma energy in Fourier space is studied numerically in a nearly adiabatic state. Such a state is represented by the Hasegawa-Mima equation, which is driven by a noise that may model the destabilization due to the phase mismatch of the plasma density and the electric potential. The noise is assumed to be Gaussian and not to be invariant under reflection along a direction ŝ. The flux density induced by such noise is found to be anisotropic: While it is random along ŝ, it is not along the perpendicular direction ŝ ⊥, and the flux is not diffusive. The renormalized response may be approximated as advective, with the velocity being proportional to ( kρ s )2, in the Fourier space.
An efficient horizontal advection scheme for the modeling of global transport of constituents
Hundsdorfer, W.; Spee, E.J.
1995-12-01
In this paper the authors consider a dimensional-splitting scheme for horizontal advection on a sphere with a uniform longitude-latitude grid. The 1D subprocesses that arise within the splitting are solved with an explicit finite-volume type scheme, which is made unconditionally stable by allowing the stencil to vary with the Courant numbers. The scheme is made positive by flux limiting. For the inaccuracies at the poles some special measures are discussed. Numerical tests show that the scheme is almost shape preserving and conservative, and it gives accurate results at low computational costs. 23 refs., 7 figs., 1 tab.
Richon, Patrick; Perrier, Frédéric; Koirala, Bharat Prasad; Girault, Frédéric; Bhattarai, Mukunda; Sapkota, Soma Nath
2011-02-01
Temporal variation of radon-222 concentration was studied at the Syabru-Bensi hot springs, located on the Main Central Thrust zone in Central Nepal. This site is characterized by several carbon dioxide discharges having maximum fluxes larger than 10 kg m(-2) d(-1). Radon concentration was monitored with autonomous Barasol™ probes between January 2008 and November 2009 in two small natural cavities with high CO(2) concentration and at six locations in the soil: four points having a high flux, and two background reference points. At the reference points, dominated by radon diffusion, radon concentration was stable from January to May, with mean values of 22 ± 6.9 and 37 ± 5.5 kBq m(-3), but was affected by a large increase, of about a factor of 2 and 1.6, respectively, during the monsoon season from June to September. At the points dominated by CO(2) advection, by contrast, radon concentration showed higher mean values 39.0 ± 2.6 to 78 ± 1.4 kBq m(-3), remarkably stable throughout the year with small long-term variation, including a possible modulation of period around 6 months. A significant difference between the diffusion dominated reference points and the advection-dominated points also emerged when studying the diurnal S(1) and semi-diurnal S(2) periodic components. At the advection-dominated points, radon concentration did not exhibit S(1) or S(2) components. At the reference points, however, the S(2) component, associated with barometric tide, could be identified during the dry season, but only when the probe was installed at shallow depth. The S(1) component, associated with thermal and possibly barometric diurnal forcing, was systematically observed, especially during monsoon season. The remarkable short-term and long-term temporal stability of the radon concentration at the advection-dominated points, which suggests a strong pressure source at depth, may be an important asset to detect possible temporal variations associated with the
Modeling Facilitated Contaminant Transport by Mobile Bacteria
NASA Astrophysics Data System (ADS)
Corapcioglu, M. Yavuz; Kim, Seunghyun
1995-01-01
Introduction of exogenous biocolloids such as genetically engineered bacteria in a bioremediation operation can enhance the transport of contaminants in groundwater by reducing the retardation effects. Because of their colloidal size and favorable surface conditions, bacteria are efficient contaminant carriers. In cases where contaminants have a low mobility in porous media because of their high partition with solid matrix, facilitated contaminant transport by mobile bacteria can create high contaminant fluxes. When metabolically active mobile bacteria are present in a subsurface environment, the system can be treated as consisting of three phases: water phase, bacterial phase, and stationary solid matrix phase. In this work a mathematical model based on mass balance equations is developed to describe the facilitated transport and fate of a contaminant and bacteria in a porous medium. Bacterial partition between the bulk solution and the stationary solid matrix and contaminant partition among three phases are represented by expressions in terms of measurable quantities. Solutions were obtained to provide estimates of contaminant and bacterial concentrations. A dimensional analysis of the transport model was utilized to estimate model parameters from the experimental data and to assess the effect of several parameters on model behavior. The model results matched favorably with experimental data of Jenkins and Lion (1993). The presence of mobile bacteria enhances the contaminant transport. However, bacterial consumption of the contaminant, which serves as a bacterial nutrient, can attenuate the contaminant mobility. The work presented in this paper is the first three-phase model to include the effects of substrate metabolism on the fate of groundwater contaminants.
Advective excess Ba transport as shown from sediment and trap geochemical signatures
Fagel, N.; Andre, L.; Dehairs, F.
1999-08-01
The authors report the results of a geochemical study of sediment and trap material. Major and trace elements (Zr, Ba, rare earth elements, and Th) were analyzed on bulk sedimentary material collected along the NE Atlantic margin. The aim is to test the widespread use of Ba-barite as a proxy for paleoproductivity in a continental margin area. This environment is of great interest because atmospheric-oceanic exchanges are important. In sediments, the geochemical signatures remain close to an upper crust reference, with flat shale-normalized rare earth elements patterns and constant elementary ratios. The calculated biogenic fraction of Ba or excess Ba (20--45%) remains lower than the excess Ba record in trap material (80--99%). The evolution of the geochemical signature along the margin reflects variable dilution of a detrital Post Archean Australian Shale-like component by a biogenic carbonaceous seawater-derived component. The trap material displays a wide range of variation in its trace element content (e.g., Ba {approximately}150--3,000 ppm, Zr {approximately}2--100 ppm), except for the abyssal site, which is characterized by constant signature. In the two other sites, all of the trace element contents increase with water depth and present pronounced seasonal changes at each sampled water depth. The amount of excess Ba also increases in the deepest traps, and its evolution throughout the year mimics the change of the other analyzed trace elements. In contrast, its relationships with particulate organic carbon are not obvious. In terms of fluxes, two periods of enhanced excess Ba fluxes are observed: (1) excess Ba flux increases with the detrital-like elements like Th especially during winter, and (2) excess Ba flux is enhanced without any change for the other trace elements during spring. To explain the first case, a supply through lateral advection is proposed. Such transient input of significant excess Ba flux will have a great impact on the yearly averaged
Code System to Solve for Release and Transport of Contaminants through Saturated/Unsaturated Media.
Energy Science and Technology Software Center (ESTSC)
1996-03-07
The BLT code solves for release and transport of contaminants from containerized wastes. Each container may have unique properties (i.e., time to failure or localized failure, e.g. pitting) and each waste form may have unique release properties. Release from the waste form is limited by one of four physical or chemical constraints: solubility, diffusion, dissolution, and surface wash-off with partitioning. The release from the waste form acts as a source for transport in the advection/dispersionmore » equation. Transport is modeled in two-dimensions through the groundwater pathway from subsurface disposal.« less
Lichtner, P.C.; Helgeson, H.C.
1986-06-20
A general formulation of multi-phase fluid flow coupled to chemical reactions was developed based on a continuum description of porous media. A preliminary version of the computer code MCCTM was constructed which implemented the general equations for a single phase fluid. The computer code MCCTM incorporates mass transport by advection-diffusion/dispersion in a one-dimensional porous medium coupled to reversible and irreversible, homogeneous and heterogeneous chemical reactions. These reactions include aqueous complexing, oxidation/reduction reactions, ion exchange, and hydrolysis reactions of stoichiometric minerals. The code MCCTM uses a fully implicit finite difference algorithm. The code was tested against analytical calculations. Applications of the code included investigation of the propagation of sharp chemical reaction fronts, metasomatic alteration of microcline at elevated temperatures and pressures, and ion-exchange in a porous column. Finally numerical calculations describing fluid flow in crystalline rock in the presence of a temperature gradient were compared with experimental results for quartzite.
Mohanram, Arvind; Ray, Chittaranjan; Harvey, Ronald W; Metge, David W; Ryan, Joseph N; Chorover, Jon; Eberl, D D
2010-10-01
In order to gain more information about the fate of Cryptosporidium parvum oocysts in tropical volcanic soils, the transport and attachment behaviors of oocysts and oocyst-sized polystyrene microspheres were studied in the presence of two soils. These soils were chosen because of their differing chemical and physical properties, i.e., an organic-rich (43-46% by mass) volcanic ash-derived soil from the island of Hawaii, and a red, iron (22-29% by mass), aluminum (29-45% by mass), and clay-rich (68-76% by mass) volcanic soil from the island of Oahu. A third agricultural soil, an organic- (13% by mass) and quartz-rich (40% by mass) soil from Illinois, was included for reference. In 10-cm long flow-through columns, oocysts and microspheres advecting through the red volcanic soil were almost completely (98% and 99%) immobilized. The modest breakthrough resulted from preferential flow-path structure inadvertently created by soil-particle aggregation during the re-wetting process. Although a high (99%) removal of oocysts and microsphere within the volcanic ash soil occurred initially, further examination revealed that transport was merely retarded because of highly reversible interactions with grain surfaces. Judging from the slope of the substantive and protracted tail of the breakthrough curve for the 1.8-μm microspheres, almost all (>99%) predictably would be recovered within ∼4000 pore volumes. This suggests that once contaminated, the volcanic ash soil could serve as a reservoir for subsequent contamination of groundwater, at least for pathogens of similar size or smaller. Because of the highly reversible nature of organic colloid immobilization in this soil type, C. parvum could contaminate surface water should overland flow during heavy precipitation events pick up near-surface grains to which they are attached. Surprisingly, oocyst and microsphere attachment to the reference soil from Illinois appeared to be at least as sensitive to changes in pH as was
Mohanram, A.; Ray, C.; Harvey, R.W.; Metge, D.W.; Ryan, J.N.; Chorover, J.; Eberl, D.D.
2010-01-01
In order to gain more information about the fate of Cryptosporidium parvum oocysts in tropical volcanic soils, the transport and attachment behaviors of oocysts and oocyst-sized polystyrene microspheres were studied in the presence of two soils. These soils were chosen because of their differing chemical and physical properties, i.e., an organic-rich (43-46% by mass) volcanic ash-derived soil from the island of Hawaii, and a red, iron (22-29% by mass), aluminum (29-45% by mass), and clay-rich (68-76% by mass) volcanic soil from the island of Oahu. A third agricultural soil, an organic- (13% by mass) and quartz-rich (40% by mass) soil from Illinois, was included for reference. In 10-cm long flow-through columns, oocysts and microspheres advecting through the red volcanic soil were almost completely (98% and 99%) immobilized. The modest breakthrough resulted from preferential flow-path structure inadvertently created by soil-particle aggregation during the re-wetting process. Although a high (99%) removal of oocysts and microsphere within the volcanic ash soil occurred initially, further examination revealed that transport was merely retarded because of highly reversible interactions with grain surfaces. Judging from the slope of the substantive and protracted tail of the breakthrough curve for the 1.8-??m microspheres, almost all (>99%) predictably would be recovered within ~4000 pore volumes. This suggests that once contaminated, the volcanic ash soil could serve as a reservoir for subsequent contamination of groundwater, at least for pathogens of similar size or smaller. Because of the highly reversible nature of organic colloid immobilization in this soil type, C. parvum could contaminate surface water should overland flow during heavy precipitation events pick up near-surface grains to which they are attached. Surprisingly, oocyst and microsphere attachment to the reference soil from Illinois appeared to be at least as sensitive to changes in pH as was observed
Iwasaki, T; Nabi, M; Shimizu, Y; Kimura, I
2015-01-01
A numerical model capable of simulating the transfer of (137)Cs in rivers associated with transport of fine sediment is presented. The accident at Fukushima Dai-ichi Nuclear Power Plant (FDNPP) released radionuclides into the atmosphere, and after fallout several radionuclides in them, such as radiocesium ((134)Cs, (137)Cs) and radioiodine ((131)I) were adsorbed on surface soil particles around FDNPP and transported by surface water. To understand the transport and deposition of the radioactive contaminant along with surface soil particles and its flux to the ocean, we modeled the transport of the (137)Cs contaminant by computing the water flow and the associated washload and suspended load transport. We have developed a two-dimensional model to simulate the plane flow structure, sediment transport and associated (137)Cs contaminant transport in rivers by combining a shallow water flow model and an advection-diffusion equation for the transport of sediment. The proposed model has been applied to the lower reach of Abukuma River, which is the main river in the highly contaminated area around FDNPP. The numerical results indicate that most (137)Cs supplied from the upstream river reach with washload would directly reach to Pacific Ocean. In contrast, washload-oriented (137)Cs supplied from the upstream river basin has a limited role in the radioactive contamination in the river. The results also suggest that the proposed framework of computational model can be a potential tool for understanding the sediment-oriented (137)Cs behavior in rivers. PMID:24909793
Arctic seabirds transport marine-derived contaminants.
Blais, Jules M; Kimpe, Lynda E; McMahon, Dominique; Keatley, Bronwyn E; Mallory, Mark L; Douglas, Marianne S V; Smol, John P
2005-07-15
Long-range atmospheric transport of pollutants is generally assumed to be the main vector for arctic contamination, because local pollution sources are rare. We show that arctic seabirds, which occupy high trophic levels in marine food webs, are the dominant vectors for the transport of marine-derived contaminants to coastal ponds. The sediments of ponds most affected by seabirds had 60 times higher DDT, 25 times higher mercury, and 10 times higher hexachlorobenzene concentrations than nearby control sites. Bird guano greatly stimulates biological productivity in these extreme environments but also serves as a major source of industrial and agricultural pollutants in these remote ecosystems. PMID:16020729
Technology Transfer Automated Retrieval System (TEKTRAN)
It has been reported that this model cannot take into account several important features of solute movement through soil. Recently, a new model has been suggested that results in a solute transport equation with fractional spatial derivatives, or FADE. We have assembled a database on published solu...
Parkhurst, D.L.
1995-01-01
PHREEQC is a computer program written in the C pwgranuning language that is designed to perform a wide variety of aqueous geochemical calculations. PHREEQC is based on an ion-association aqueous model and has capabilities for (1) speciation and saturation-index calculations, (2) reaction-path and advective-transport calculations involving specified irreversible reactions, mixing of solutions, mineral and gas equilibria surface-complex-ation reactions, and ion-exchange reactions, and (3) inverse modeling, which finds sets of mineral and gas mole transfers that account for composition differences between waters, within specified compositional uncertainties. PHREEQC is derived from the Fortran program PHREEQE, but it has been completely rewritten in C with the addition many new capabilities. New features include the capabilities to use redox couples to distribute redox elements among their valence states in speciation calculations; to model ion-exchange and surface-compiexation reactions; to model reactions with a fixed-pressure, multicomponent gas phase (that is, a gas bubble); to calculate the mass of water in the aqueous phase during reaction and transport calculations; to keep track of the moles of minerals present in the solid phases and determine antomaticaHy the thermodynamically stable phase assemblage; to simulate advective transport in combination with PHREEQC's reaction-modeling capability; and to make inverse modeling calculations that allow for uncertainties in the analytical data. The user interface is improved through the use of a simplified approach to redox reactions, which includes explicit mole-balance equations for hydrogen and oxygen; the use of a revised input that is modular and completely free format; and the use of mineral names and standard chemical symbolism rather than index numbers. The use of (2 eliminates nearly all limitations on army sizes, including numbers of elements, aqueous species, solutions, phases, and lengths of character
Kim, Chang-Bae; Krommes, J.A.
1988-08-01
The work of Krommes and Smith on rigorous upper bounds for the turbulent transport of a passively advected scalar (/ital Ann. Phys./ 177:246 (1987)) is extended in two directions: (1) For their ''reference model,'' improved upper bounds are obtained by utilizing more sophisticated two-time constraints which include the effects of cross-correlations up to fourth order. Numerical solutions of the model stochastic differential equation are also obtained; they show that the new bounds compare quite favorably with the exact results, even at large Reynolds and Kubo numbers. (2) The theory is extended to take account of a finite spatial autocorrelation length L/sub c/. As a reasonably generic example, the problem of particle transport due to statistically specified stochastic magnetic fields in a collisionless turbulent plasma is revisited. A bound is obtained which reduces for small L/sub c/ to the quasilinear limit and for large L/sub c/ to the strong turbulence limit, and which provides a reasonable and rigorous interpolation for intermediate values of L/sub c/. 18 refs., 6 figs.
ADVECTIVE TRANSPORT OF INTERSTELLAR PLASMA INTO THE HELIOSPHERE ACROSS THE RECONNECTING HELIOPAUSE
Strumik, M.; Grzedzielski, S.; Czechowski, A.; Macek, W. M.; Ratkiewicz, R.
2014-02-10
We discuss results of magnetohydrodynamical model simulations of plasma dynamics in the proximity of the heliopause (HP). The model is shown to fit details of the magnetic field variations observed by the Voyager 1 spacecraft during the transition from the heliosphere to the local interstellar medium (LISM). We propose an interpretation of magnetic field structures observed by Voyager 1 in terms of fine-scale physical processes. Our simulations reveal an effective transport mechanism of relatively dense LISM plasma across the reconnecting HP into the heliosphere. The mechanism is associated with annihilation of magnetic sectors in the heliospheric plasma near the HP.
Advective heat transport associated to regional Earth degassing in central Apennine (Italy)
NASA Astrophysics Data System (ADS)
Cardellini, Carlo; Chiodini, Giovanni; Caliro, Stefano; Chiarabba, Claudio; Frondini, Francesco
2013-04-01
The main springs of central Italy Apennines were investigated, in order to compute the amount of heat transported by groundwaters and to compute the fraction of heat due to the geothermal heat flux. The 46 investigated springs represent a significant portion of the permeable structures of the Apennine being characterised by a cumulative flow rate of 130 m3/s, i.e. ~ 50% of the water discharged in this sector of the Apennines. The groundwaters are characterised by relatively low temperatures, but the occurrence of an heat anomaly is evident when the differences between the temperatures of springs and recharge waters are compared with the corresponding altitude difference. A total amount of heat of ~ 2.1 × 109 J/s has been estimated to be transported by these groundwaters. Most of this heat (57%) is given by geothermal warming while the remaining 43% is due to gravitational potential energy dissipation. The computed geothermal warming implies very high heat flux, with values higher than 300 mW/m2, in a large sector of the Apennines which was considered to date be characterised by normal to low conductive heat flux. The same area is affected by high fluxes of CO2 from a deep source and the strict correlation between the geothermal warming and the input of deep CO2-rich fluids is testified by the fact that all the thermally anomalous groundwaters are also affected by the input of deeply derived CO2 contrary to those not thermally anomalous which display any input of deeply derived CO2. This correspondence reasonably suggest the geothermal heat is transported from depth by CO2 rich fluids, which enter the aquifers and mix with infiltrating waters. The amount of geothermal heat transported by central Apennine cold groundwaters is in absolute very high. It results the double than the hydrothermal heat discharge of the US Cascade Range (~1×103 MW) and is about the half of the total heat discharged at Yellowstone, one of the largest hydrothermal system of the world (5-6
McCaulou, D.R.
1993-10-01
Replicate column experiments were done to quantify the effects of temperature and bacterial motility on advective transport through repacked, but otherwise unaltered, natural aquifer sediment. The bacteria used in this study, A0500, was a flagellated, spore-forming rod isolated from the deep subsurface at DOE`s Savannah River Laboratory. Motility was controlled by turning on flagellar metabolism at 18{degrees}C but off at 40{degrees}C. Microspheres were used to independently quantify the effects of temperature on the sticking efficiency ({alpha}), estimated using a steady-state filtration model. The observed greater microsphere removal at the higher temperature agreed with the physical-chemical model, but bacteria removal at 18{degrees}C was only half that at 4{degrees}C. The sticking efficiency for non-motile A0500 (4{degrees}C) was over three times that of the motile A0500 (18{degrees}C), 0.073 versus 0.022 respectively. Analysis of complete breakthrough curves using a non-steady, kinetically limited, transport model to estimate the time scales of attachment and detachment suggested that motile A 0500 bacteria traveled twice as far as non-motile A 0500 bacteria before becoming attached. Once attached, non-motile colloids detached on the time scale of 9 to 17 days. The time scale for detachment of motile A0500 bacteria was shorter, 4 to 5 days. Results indicate that bacterial attachment was reversible and detachment was enhanced by bacterial motifity. The kinetic energy of bacterial motility changed the attachment-detachment kinetics in favor of the detached state. The chemical factors responsible for the enhanced transport are not known. However, motility may have caused weakly held bacteria to detach from the secondary minimum, and possibly from the primary minimum, as described by DLVO theory.
Advective heat transport associated with regional Earth degassing in central Apennine (Italy)
NASA Astrophysics Data System (ADS)
Chiodini, G.; Cardellini, C.; Caliro, S.; Chiarabba, C.; Frondini, F.
2013-07-01
In this work we show that the main springs of the central Apennine transport a total amount of heat of ˜2.2×109 J s-1. Most of this heat (57%) is the result of geothermal warming while the remaining 43% is due to gravitational potential energy dissipation. This result indicates that a large area of the central Apennines is very hot with heat flux values >300 mW m-2. These values are higher than those measured in the magmatic and famously geothermal provinces of Tuscany and Latium and about 1/3 of the total heat discharged at Yellowstone. This finding is surprising because the central Apennines have been thought to be a relatively cold area. Translated by CO2 rich fluids, this heat anomaly suggests the existence of a thermal source such as a large magmatic intrusion at depth. Recent tomographic images of the area support the presence of such an intrusion visible as a broad negative velocity anomaly in seismic waves. Our results indicate that the thermal regime of tectonically active areas of the Earth, where meteoric waters infiltrate and deeply circulate, should be revised on the basis of mass and energy balances of the groundwater systems.
MULTIMEDIA CONTAMINANT FATE, TRANSPORT, AND EXPOSURE MODEL (MMSOILS)
The Multimedia Contaminant Fate, Transport, and Exposure Model (MMSOILS) estimates the human exposure and health risk associated with releases of contamination from hazardous waste sites. The methodology consists of a multimedia model that addresses the transport of a chemical in...
Kurylyk, Barret L.; McKenzie, Jeffrey M; MacQuarrie, Kerry T. B.; Voss, Clifford I.
2014-01-01
Numerous cold regions water flow and energy transport models have emerged in recent years. Dissimilarities often exist in their mathematical formulations and/or numerical solution techniques, but few analytical solutions exist for benchmarking flow and energy transport models that include pore water phase change. This paper presents a detailed derivation of the Lunardini solution, an approximate analytical solution for predicting soil thawing subject to conduction, advection, and phase change. Fifteen thawing scenarios are examined by considering differences in porosity, surface temperature, Darcy velocity, and initial temperature. The accuracy of the Lunardini solution is shown to be proportional to the Stefan number. The analytical solution results obtained for soil thawing scenarios with water flow and advection are compared to those obtained from the finite element model SUTRA. Three problems, two involving the Lunardini solution and one involving the classic Neumann solution, are recommended as standard benchmarks for future model development and testing.
Alcaraz, Mar; García-Gil, Alejandro; Vázquez-Suñé, Enric; Velasco, Violeta
2016-02-01
Borehole Heat Exchangers (BHEs) are increasingly being used to exploit shallow geothermal energy. This paper presents a new methodology to provide a response to the need for a regional quantification of the geothermal potential that can be extracted by BHEs and the associated environmental impacts. A set of analytical solutions facilitates accurate calculation of the heat exchange of BHEs with the ground and its environmental impacts. For the first time, advection and dispersion heat transport mechanisms and the temporal evolution from the start of operation of the BHE are taken into account in the regional estimation of shallow geothermal resources. This methodology is integrated in a GIS environment, which facilitates the management of input and output data at a regional scale. An example of the methodology's application is presented for Barcelona, in Spain. As a result of the application, it is possible to show the strengths and improvements of this methodology in the development of potential maps of low temperature geothermal energy as well as maps of environmental impacts. The minimum and maximum energy potential values for the study site are 50 and 1800 W/m(2) for a drilled depth of 100 m, proportionally to Darcy velocity. Regarding to thermal impacts, the higher the groundwater velocity and the energy potential, the higher the size of the thermal plume after 6 months of exploitation, whose length ranges from 10 to 27 m long. A sensitivity analysis was carried out in the calculation of heat exchange rate and its impacts for different scenarios and for a wide range of Darcy velocities. The results of this analysis lead to the conclusion that the consideration of dispersion effects and temporal evolution of the exploitation prevent significant differences up to a factor 2.5 in the heat exchange rate accuracy and up to several orders of magnitude in the impacts generated. PMID:26605833
TNT transport and fate in contaminated soil
Comfort, S.D.; Shea, P.J.; Hundal, L.S.
1995-11-01
Past disposal practices at munitions production plants have contaminated terrestrial and aquatk ecosystems with 2,4,6-trinitrotoluene (TNT). We determined TNT transport, degradation, and long-term sorption characteristics in soil. Transport experiments were conducted with repacked, unsaturated soil columns containing uncontaminated soil or layers of contaminated and uncontaminated soil. Uncontaminated soil columns received multiple pore volumes (22-50) of a TNT-{sup 3}H{sub 2}O pulse, containing 70 or 6.3 mg TNT L{sup -1} at a constant pore water velocity. TNT breakthrough curves (BTCs) never reached initial solute pulse concentrations. Apex concentrations (C/C{sub o}) were between 0.6 and 0.8 for an initial pulse of 70 mg TNT L{sup -1} and 0.2 to 0.3 for the 6.3 mg TNT L{sup -1} pulse. Earlier TNT breakthrough was observed at the higher pulse concentration. This mobility difference was predicted from the nonlinear adsorption isotherm determined for TNT sorption. In all experiments, a significant fraction of added TNT was recovered as amino degradates of TNT. Mass balance estimates indicated 81% of the added TNT was recovered (as TNT and amino degradates) from columns receiving the 70 mg TNT L{sup -1} pulse compared to 35% from columns receiving the 6.3 mg TNT L{sup -1} pulse. Most of the unaccountable TNT was hypothesized to be unextractable. This was supported by a 168-d sorption experiment, which found that within 14d, 80% of {sup 14}C activity (added as {sup 14}C-TNT) was adsorbed and roughly 40% unextractable. Our observations illustrate that TNT sorption and degradation are concentration-dependent and the assumptions of linear adsorption and adsorption-desorption singularity commonly used in transport modeling, may not be valid for predicting TNT transport in munitions-contaminated soils. 29 refs., 6 figs., 7 tabs.
Wagner, Brian J.; Gorelick, Steven M.
1986-01-01
A simulation nonlinear multiple-regression methodology for estimating parameters that characterize the transport of contaminants is developed and demonstrated. Finite difference containment transport simulation is combined with a nonlinear weighted least squares multiple-regression procedure. The technique provides optimal parameter estimates and gives statistics for assessing the reliability of these estimates under certain general assumptions about the distributions of the random measurement errors. Monte Carlo analysis is used to estimate parameter reliability for a hypothetical homogeneous soil column for which concentration data contain large random measurement errors. The value of data collected spatially versus data collected temporally was investigated for estimation of velocity, dispersion coefficient, effective porosity, first-order decay rate, and zero-order production. The use of spatial data gave estimates that were 2-3 times more reliable than estimates based on temporal data for all parameters except velocity. (Estimated author abstract) Refs.
Can contaminant transport models predict breakthrough?
Peng, Wei-Shyuan; Hampton, Duane R.; Konikow, Leonard F.; Kambham, Kiran; Benegar, Jeffery J.
2000-01-01
A solute breakthrough curve measured during a two-well tracer test was successfully predicted in 1986 using specialized contaminant transport models. Water was injected into a confined, unconsolidated sand aquifer and pumped out 125 feet (38.3 m) away at the same steady rate. The injected water was spiked with bromide for over three days; the outflow concentration was monitored for a month. Based on previous tests, the horizontal hydraulic conductivity of the thick aquifer varied by a factor of seven among 12 layers. Assuming stratified flow with small dispersivities, two research groups accurately predicted breakthrough with three-dimensional (12-layer) models using curvilinear elements following the arc-shaped flowlines in this test. Can contaminant transport models commonly used in industry, that use rectangular blocks, also reproduce this breakthrough curve? The two-well test was simulated with four MODFLOW-based models, MT3D (FD and HMOC options), MODFLOWT, MOC3D, and MODFLOW-SURFACT. Using the same 12 layers and small dispersivity used in the successful 1986 simulations, these models fit almost as accurately as the models using curvilinear blocks. Subtle variations in the curves illustrate differences among the codes. Sensitivities of the results to number and size of grid blocks, number of layers, boundary conditions, and values of dispersivity and porosity are briefly presented. The fit between calculated and measured breakthrough curves degenerated as the number of layers and/or grid blocks decreased, reflecting a loss of model predictive power as the level of characterization lessened. Therefore, the breakthrough curve for most field sites can be predicted only qualitatively due to limited characterization of the hydrogeology and contaminant source strength.
Packaging and transportation of radioactively contaminated lead
Gleason, Eugene; Holden, Gerard
2007-07-01
Under the management of the Nuclear Decommissioning Authority (NDA) the government of the United Kingdom has launched an ambitious program to remediate the nation's nuclear waste legacy. Over a twenty-five year period NDA plans to decommission several first generation nuclear power plants and other radioactive facilities. The use innovative, safe 'fit for purpose' technologies will be a major part of this complex program. This paper will present a case study of a recently completed project undertaken in support of the nuclear decommissioning activities at the Sellafield site in the United Kingdom. The focus is on an innovative application of new packaging technology developed for the safe transportation of radioactively contaminated lead objects. Several companies collaborated on the project and contributed to its safe and successful conclusion. These companies include British Nuclear Group, Gravatom Engineering, W. F. Bowker Transport, Atlantic Container Lines, MHF Logistical Solutions and Energy Solutions. New containers and a new innovative inter-modal packaging system to transport the radioactive lead were developed and demonstrated during the project. The project also demonstrated the potential contribution of international nuclear recycling activities as a safe, economic and feasible technical option for nuclear decommissioning in the United Kingdom. (authors)
Shapiro, A.M.; Renken, R.A.; Harvey, R.W.; Zygnerski, M.R.; Metge, D.W.
2008-01-01
A tracer experiment, using a nonreactive tracer, was conducted as part of an investigation of the potential for chemical and pathogen migration to public supply wells that draw groundwater from the highly transmissive karst limestone of the Biscayne aquifer in southeastern Florida. The tracer was injected into the formation over approximately 1 h, and its recovery was monitored at a pumping well approximately 100 m from the injection well. The first detection of the tracer occurred after approximately 5 h, and the peak concentration occurred at about 8 h after the injection. The tracer was still detected in the production well more than 6 days after injection, and only 42% of the tracer mass was recovered. It is hypothesized that a combination of chemical diffusion and slow advection resulted in significant retention of the tracer in the formation, despite the high transmissivity of the karst limestone. The tail of the breakthrough curve exhibited a straight-line behavior with a slope of -2 on a log-log plot of concentration versus time. The -2 slope is hypothesized to be a function of slow advection, where the velocities of flow paths are hypothesized to range over several orders of magnitude. The flow paths having the slowest velocities result in a response similar to chemical diffusion. Chemical diffusion, due to chemical gradients, is still ongoing during the declining limb of the breakthrough curve, but this process is dwarfed by the magnitude of the mass flux by slow advection.
COLLOIDAL CONSIDERATIONS IN GROUNDWATER SAMPLING AND CONTAMINANT TRANSPORT PREDICTIONS
The association of contaminants with suspended colloidal material in groundwater is a possible transport mechanism and a complicating factor for accurate estimations of the aqueous geochemistry of subsurface systems. esearch to date indicates colloidal facilitated transport of co...
NASA Astrophysics Data System (ADS)
Zhou, R.; Zhan, H.
2015-12-01
With the consideration of advection, dispersion, adsorption and first order decay in the fracture and rock matrix in a single fracture model, a new semi-analytical solution is derived using the Mobile-Immobile Method. It can be used to estimate the concentration at any location at any time precisely within the fracture and rock matrix. Most fractures found underground are filled with the conglomerate, sand, clay and other kinds of possible porous media. The existence of those filling ingredients leads to the isolated pore space within the fracture, which is also called immobile zone. Certain assumptions have be made: the diffusion is the only way that the contamination travels from the fracture to the matrix as the large permeability difference between them; the diffusive transport is dominant in the rock matrix while the advective-dispersive transport plays the major role in the fracture. Experimental data have been collected from literatures to compare the performance of this semi-analytical solution from the classical analytical solution. The comparison shows that the semi-analytical solution simulates it better when the mobile zone percentage is limited. Also, the effects of matrix diffusion, dispersivity and Darcy velocity in the fracture, fracture aperture, first order mass transfer rate and mobile zone percentage on solute transport are demonstrated through the sensitivity analysis, concentration profiles and breakthrough curves. By modifying the boundary conditions and adding an advection term in the rock matrix governing equation, this model can be extended to a two-layer solute transport model.
Feed gas contaminant removal in ion transport membrane systems
Underwood, Richard Paul; Makitka, III, Alexander; Carolan, Michael Francis
2012-04-03
An oxygen ion transport membrane process wherein a heated oxygen-containing gas having one or more contaminants is contacted with a reactive solid material to remove the one or more contaminants. The reactive solid material is provided as a deposit on a support. The one or more contaminant compounds in the heated oxygen-containing gas react with the reactive solid material. The contaminant-depleted oxygen-containing gas is contacted with a membrane, and oxygen is transported through the membrane to provide transported oxygen.
Hill, M.C.; Ely, D.M.; Tiedeman, C.R.; O'Brien, G.M.; D'Agnese, F.A.; Faunt, C.C.
2001-08-01
When a model is calibrated by nonlinear regression, calculated diagnostic statistics and measures of uncertainty provide a wealth of information about many aspects of the system. This report presents a method of ranking the likely importance of existing observation locations using measures of prediction uncertainty. It is suggested that continued monitoring is warranted at more important locations, and unwarranted or less warranted at less important locations. The report develops the methodology and then demonstrates it using the hydraulic-head observation locations of a three-layer model of the Death Valley regional flow system (DVRFS). The predictions of interest are subsurface transport from beneath Yucca Mountain and 14 underground Test Area (UGTA) sites. The advective component of transport is considered because it is the component most affected by the system dynamics represented by the regional-scale model being used. The problem is addressed using the capabilities of the U.S. Geological Survey computer program MODFLOW-2000, with its ADVective-Travel Observation (ADV) Package, and an additional computer program developed for this work.
NASA Astrophysics Data System (ADS)
Zhu, Yan; Shi, Liangsheng; Yang, Jinzhong; Wu, Jingwei; Mao, Deqiang
2013-09-01
An efficient integrated modeling approach is developed to simulate the contaminant transport in the subsurface system. The unsaturated zone is divided into a number of horizontal sub-areas according to the atmospheric boundary conditions, land use types and hydrological conditions. Solute migration through the unsaturated zone of each sub-area is assumed to be vertical and can be represented by the one-dimensional advection-dispersion equation, which is then coupled to the three-dimensional advection-dispersion equation representing the subsequent groundwater transport. The finite element method is adopted to discretize the vertical solute equation, while the hybrid finite element and finite difference method is used to discretize the three-dimensional saturated solute transport equation, which is split into the horizontal and vertical equations based on the concept of the horizontal/vertical splitting. The unsaturated and saturated solute transport equations are combined into a unified matrix by the mass balance analysis for the adjacent nodes located at the one-dimensional soil column and at the water table. Two hypothetical cases and two field cases are simulated to test the validity of the model with the results compared with those from HYDRUS-1D, SWMS2D and the measured data. The limitations of the model are discussed as well. The analysis of the four cases demonstrates that the proposed model can calculate the water flow and solute transport reasonably even with complex boundary and variable topography conditions. It also shows that the model is efficient to simulate the water flow and solute transport in regional-scale areas with small computational costs. However, the model will lose accuracy when the lateral dispersion effect is dominant in the unsaturated zone.
Ackerman, Daniel J.; Rousseau, Joseph P.; Rattray, Gordon W.; Fisher, Jason C.
2010-01-01
Three-dimensional steady-state and transient models of groundwater flow and advective transport in the eastern Snake River Plain aquifer were developed by the U.S. Geological Survey in cooperation with the U.S. Department of Energy. The steady-state and transient flow models cover an area of 1,940 square miles that includes most of the 890 square miles of the Idaho National Laboratory (INL). A 50-year history of waste disposal at the INL has resulted in measurable concentrations of waste contaminants in the eastern Snake River Plain aquifer. Model results can be used in numerical simulations to evaluate the movement of contaminants in the aquifer. Saturated flow in the eastern Snake River Plain aquifer was simulated using the MODFLOW-2000 groundwater flow model. Steady-state flow was simulated to represent conditions in 1980 with average streamflow infiltration from 1966-80 for the Big Lost River, the major variable inflow to the system. The transient flow model simulates groundwater flow between 1980 and 1995, a period that included a 5-year wet cycle (1982-86) followed by an 8-year dry cycle (1987-94). Specified flows into or out of the active model grid define the conditions on all boundaries except the southwest (outflow) boundary, which is simulated with head-dependent flow. In the transient flow model, streamflow infiltration was the major stress, and was variable in time and location. The models were calibrated by adjusting aquifer hydraulic properties to match simulated and observed heads or head differences using the parameter-estimation program incorporated in MODFLOW-2000. Various summary, regression, and inferential statistics, in addition to comparisons of model properties and simulated head to measured properties and head, were used to evaluate the model calibration. Model parameters estimated for the steady-state calibration included hydraulic conductivity for seven of nine hydrogeologic zones and a global value of vertical anisotropy. Parameters
Koch, M.
1995-12-31
A new mesh-adaptive 1D collocation technique has been developed to efficiently solve transient advection-dominated transport problems in porous media that are governed by a hyperbolic/parabolic (singularly perturbed) PDE. After spatial discretization a singularly perturbed ODE is obtained which is solved by a modification of the COLNEW ODE-collocation code. The latter also contains an adaptive mesh procedure that has been enhanced here to resolve linear and nonlinear transport flow problems with steep fronts where regular FD and FE methods often fail. An implicit first-order backward Euler and a third-order Taylor-Donea technique are employed for the time integration. Numerical simulations on a variety of high Peclet-number transport phenomena as they occur in realistic porous media flow situations are presented. Examples include classical linear advection-diffusion, nonlinear adsorption, two-phase Buckley-Leverett flow without and with capillary forces (Rapoport-Leas equation) and Burgers` equation for inviscid fluid flow. In most of these examples sharp fronts and/or shocks develop which are resolved in an oscillation-free manner by the present adaptive collocation method. The backward Euler method has some amount of numerical dissipation is observed when the time-steps are too large. The third-order Taylor-Donea technique is less dissipative but is more prone to numerical oscillations. The simulations show that for the efficient solution of nonlinear singularly perturbed PDE`s governing flow transport a careful balance must be struck between the optimal mesh adaptation, the nonlinear iteration method and the time-stepping procedure. More theoretical research is needed with this regard.
NASA Astrophysics Data System (ADS)
Périard, Yann; José Gumiere, Silvio; Rousseau, Alain N.; Caron, Jean
2013-04-01
Certain contaminants may travel faster through soils when they are sorbed to subsurface colloidal particles. Indeed, subsurface colloids may act as carriers of some contaminants accelerating their translocation through the soil into the water table. This phenomenon is known as colloid-facilitated contaminant transport. It plays a significant role in contaminant transport in soils and has been recognized as a source of groundwater contamination. From a mechanistic point of view, the attachment/detachment of the colloidal particles from the soil matrix or from the air-water interface and the straining process may modify the hydraulic properties of the porous media. Šimůnek et al. (2006) developed a model that can simulate the colloid-facilitated contaminant transport in variably saturated porous media. The model is based on the solution of a modified advection-dispersion equation that accounts for several processes, namely: straining, exclusion and attachement/detachement kinetics of colloids through the soil matrix. The solutions of these governing, partial differential equations are obtained using a standard Galerkin-type, linear finite element scheme, implemented in the HYDRUS-2D/3D software (Šimůnek et al., 2012). Modeling colloid transport through the soil and the interaction of colloids with the soil matrix and other contaminants is complex and requires the characterization of many model parameters. In practice, it is very difficult to assess actual transport parameter values, so they are often calibrated. However, before calibration, one needs to know which parameters have the greatest impact on output variables. This kind of information can be obtained through a sensitivity analysis of the model. The main objective of this work is to perform local and global sensitivity analyses of the colloid-facilitated contaminant transport module of HYDRUS. Sensitivity analysis was performed in two steps: (i) we applied a screening method based on Morris' elementary
A Generalized Model for Transport of Contaminants in Soil by Electric Fields
Paz-Garcia, Juan M.; Baek, Kitae; Alshawabkeh, Iyad D.; Alshawabkeh, Akram N.
2012-01-01
A generalized model applicable to soils contaminated with multiple species under enhanced boundary conditions during treatment by electric fields is presented. The partial differential equations describing species transport are developed by applying the law of mass conservation to their fluxes. Transport, due to migration, advection and diffusion, of each aqueous component and complex species are combined to produce one partial differential equation hat describes transport of the total analytical concentrations of component species which are the primary dependent variables. This transport couples with geochemical reactions such as aqueous equilibrium, sorption, precipitation and dissolution. The enhanced model is used to simulate electrokinetic cleanup of lead and copper contaminants at an Army Firing Range. Acid enhancement is achieved by the use of adipic acid to neutralize the basic front produced for the cathode electrochemical reaction. The model is able to simulate enhanced application of the process by modifying the boundary conditions. The model showed that kinetics of geochemical reactions, such as metals dissolution/leaching and redox reactions might be significant for realistic prediction of enhanced electrokinetic extraction of metals in real world applications. PMID:22242884
A generalized model for transport of contaminants in soil by electric fields.
Paz-Garcia, Juan Manuel; Baek, Kitae; Alshawabkeh, Iyad D; Alshawabkeh, Akram N
2012-01-01
A generalized model applicable to soils contaminated with multiple species under enhanced boundary conditions during treatment by electric fields is presented. The partial differential equations describing species transport are developed by applying the law of mass conservation to their fluxes. Transport, due to migration, advection and diffusion, of each aqueous component and complex species are combined to produce one partial differential equation that describes transport of the total analytical concentrations of component species which are the primary dependent variables. This transport couples with geochemical reactions such as aqueous equilibrium, sorption, precipitation and dissolution. The enhanced model is used to simulate electrokinetic cleanup of lead and copper contaminants at an Army Firing Range. Acid enhancement is achieved by the use of adipic acid to neutralize the basic front produced for the cathode electrochemical reaction. The model is able to simulate enhanced application of the process by modifying the boundary conditions. The model showed that kinetics of geochemical reactions, such as metals dissolution/leaching and redox reactions, may be significant for realistic prediction of enhanced electrokinetic extraction of metals in real-world applications. PMID:22242884
Update on Advection-Diffusion Purge Flow Model
NASA Technical Reports Server (NTRS)
Brieda, Lubos
2015-01-01
Gaseous purge is commonly used in sensitive spacecraft optical or electronic instruments to prevent infiltration of contaminants and/or water vapor. Typically, purge is sized using simplistic zero-dimensional models that do not take into account instrument geometry, surface effects, and the dependence of diffusive flux on the concentration gradient. For this reason, an axisymmetric computational fluid dynamics (CFD) simulation was recently developed to model contaminant infiltration and removal by purge. The solver uses a combined Navier-Stokes and Advection-Diffusion approach. In this talk, we report on updates in the model, namely inclusion of a particulate transport model.
Experimental investigation of contaminant transport in porous media. Research report
Wang, J.C.; Booker, J.R.; Carter, J.P.
1998-10-01
When numerical methods are applied to simulate a real contaminant transport problem, the values of a number of key parameters such as porosity, hydrodynamic dispersion coefficient or dispersivity and Darcy velocity or seepage velocity or seepage velocity are needed. In this paper, two different experimental programs, involving two types of column test and a well simulation test, were designed to demonstrate that the theory developed to explain contaminant transport in porous media is capable of representing the actual phenomenon of contaminant migration in soil. It is demonstrated that experiments can also be carried out to determine the properties necessary to model a real case of contaminant migration in porous media.
Smith, J.A.; Sahoo, D.; Mclellan, H.M.; Imbrigiotta, T.E.
1997-01-01
Transport of a nonionic surfactant (Triton X-100) at aqueous concentrations less than 400 mg/L through a trichloroethene-contaminated sand-and-gravel aquifer at Picatinny Arsenal, NJ, has been studied through a series of laboratory and field experiments. In the laboratory, batch and column experiments were conducted to quantify the rate and amount of Triton X-100 sorption to the aquifer sediments. In the field, a 400 mg/L aqueous Triton X-100 solution was injected into the aquifer at a rate of 26.5 L/min for a 35-d period. The transport of Triton X-100 was monitored by sampling and analysis of groundwater at six locations surrounding the injection well. Equilibrium batch sorption experiments showed that Triton X-100 sorbs strongly and nonlinearly to the field soil with the sharpest inflection point of the isotherm occurring at an equilibrium aqueous Triton X-100 concentration close to critical micelle concentration. Batch, soil column, and field experimental data were analyzed with zero-, one-, and two- dimensional (respectively) transient solute transport models with either equilibrium or rate-limited sorption. These analyses reveal that Triton X- 100 sorption to the aquifer solids is slow relative to advective and dispersive transport and that an equilibrium sorption model cannot simulate accurately the observed soil column and field data. Comparison of kinetic sorption parameters from batch, column, and field transport data indicate that both physical heterogeneities and Triton X-100 mass transfer between water and soil contribute to the kinetic transport effects.Transport of a nonionic surfactant (Triton X-100) at aqueous concentrations less than 400 mg/L through a trichloroethene-contaminated sand-and-gravel aquifer was studied. Equilibrium batch sorption experiments showed that Triton X-100 sorbs strongly and nonlinearly to the field soil with the sharpest inflection point of the isotherm occurring at an equilibrium aqueous Triton X-100 concentration close to
FACILITATED TRANSPORT OF INORGANIC CONTAMINANTS IN GROUNDWATER: PART II. COLLOIDAL TRANSPORT
This project consisted of both field and laboratory components. Field studies evaluated routine sampling procedures for determination of aqueous inorganicgeochemistry and assessment of contaminant transport by colloidal mobility. Research at three different metal-contaminated sit...
CONTAMINATED SEDIMENT TRANSPORT AND FATE MODELING
Many Superfund sites include rivers, reservoirs, and other surface bodies of water (and the adjacent floodplains) that are highly contaminated with PCBs, metals, and other toxic chemicals. Examples of contaminated waters are the Hudson, Fox, Housatonic and Clark Fork Rivers, Lake...
Woodman, N.D. Rees-White, T.C.; Stringfellow, A.M.; Beaven, R.P.; Hudson, A.P.
2015-04-15
Highlights: • Multiple tracers were applied to saturated MSW to test dual-porosity properties. • Lithium demonstrated to be non-conservative as a tracer. • 260 mm diameter column too small to test transport properties of MSW. • The classical advection-dispersion mode was rejected due to high dispersivity. • Characteristic diffusion times did not vary with the tracer. - Abstract: Two column tests were performed in conditions emulating vertical flow beneath the leachate table in a biologically active landfill to determine dominant transport mechanisms occurring in landfills. An improved understanding of contaminant transport process in wastes is required for developing better predictions about potential length of the long term aftercare of landfills, currently measured in timescales of centuries. Three tracers (lithium, bromide and deuterium) were used. Lithium did not behave conservatively. Given that lithium has been used extensively for tracing in landfill wastes, the tracer itself and the findings of previous tests which assume that it has behaved conservatively may need revisiting. The smaller column test could not be fitted with continuum models, probably because the volume of waste was below a representative elemental volume. Modelling compared advection-dispersion (AD), dual porosity (DP) and hybrid AD–DP models. Of these models, the DP model was found to be the most suitable. Although there is good evidence to suggest that diffusion is an important transport mechanism, the breakthrough curves of the different tracers did not differ from each other as would be predicted based on the free-water diffusion coefficients. This suggested that solute diffusion in wastes requires further study.
Kim, Song-Bae; Corapcioglu, M Yavuz
2002-12-01
In riverbank filtration, contaminant transport is affected by colloidal particles such as dissolved organic matter (DOM) and bacterial particles. In addition, the subsurface heterogeneity influences the behavior of contaminant transport in riverbank filtration. A mathematical model is developed to describe the contaminant transport in dual-porosity media in the presence of DOM and bacteria as mobile colloids. In the model development, a porous medium is divided into the mobile and immobile regions to consider the presence of ineffective micropores in physically heterogeneous riverbanks. We assume that the contaminant transport in the mobile region is controlled by the advection and dispersion while the contaminant transport in the immobile region occurs due to the molecular diffusion. The contaminant transfer between the mobile and immobile regions takes place by diffusive mass transfer. The mobile region is conceptualized as a four-phase system: two mobile colloidal phases, an aqueous phase, and a solid matrix. The complete set of governing equations is solved numerically with a fully implicit finite difference method. The model results show that in riverbank filtration, the contaminant can migrate further than expected due to the presence of DOM and bacteria. In addition, the contaminant mobility increases further in the presence of the immobile region in aquifers. A sensitivity analysis shows that in dual-porosity media, earlier breakthrough of the contaminant takes place as the volumetric fraction of the mobile region decreases. It is also demonstrated that as the contaminant mass transfer rate coefficient between the mobile and immobile regions increases, the contaminant concentration gradient between the two regions reverses at earlier pore volumes. The contaminant mass transfer coefficient between the mobile and immobile regions mainly controls the tailing effect of the contaminant breakthrough. The contaminant breakthrough curves are sensitive to changes in
Preliminary study of niobium alloy contamination by transport through helium
NASA Technical Reports Server (NTRS)
Scheuermann, C. M.; Moore, T. J.; Wheeler, D. R.
1987-01-01
Transport of gaseous contaminants through the working fluid to or from sensitive refractory alloys is theoretically possible during long time operation of Brayton and Stirling space power generation systems which use a gas as the working fluid. A test was designed which could give an answer to whether transport of contaminants through the working fluid was a potential major problem. The findings of that preliminary study are summarized.
Contaminant transport from an array of sources
Kim, C.L.; Chambre, P.L.; Lee, W.W.L.; Pigford, T.H.
1987-04-01
This document shows analytic solutions to the problem of contaminant dispersion from an array of point sources in a waste disposal site. These solutions are for waste sources in a fluid-saturated porous medium, and may be for isotropic or anisotropic dispersion. The solutions are illustrated through isopleths of contaminants for a planar array of point sources perpendicular to ground-water flow. The concentration fields several meters away from this plane can be approximated by equivalent plane sources. 2 refs., 4 figs.
Investigation of contaminant transport from the saginaw confined disposal facility
Velleux, M.L.; Rathbun, J.E.; Kreis, R.G.; Martin, J.L.; Mac, M.J.
1993-01-01
Pilot biomonitoring and modeling studies were conducted at the Saginaw Confined Disposal Facility (CDF), Saginaw Bay, Lake Huron, during 1987 to develop methods to assess the potential for or magnitude of (1) contaminant transport from the dike interior to the outside environment, (2) impacts of CDF disposal on the water column and sediments, and (3) impacts of CDF disposal on aquatic biota living in the outdike zone. Polychlorinated biphenyls (PCBs) were selected for study due to their presence in the sediments of the Saginaw River/Bay ecosystem. A mathematical model of near-field contaminant transport through the dike walls was constructed. Model predictions indicate that the rate of contaminant transport through the dike is expected to be small, amounting to less than 0.25 kg of PCBs after 5,000 days of simulation. A mathematical model of the farfield impacts of CDF transport was also constructed. Model predictions indicate that the incremental increase in steady-state, water column PCB concentrations in Saginaw Bay is expected to be approximately 0.05 ng/L per kg of PCB transported from the CDF. A biomonitoring program was developed to assess contaminant transport through dike walls and its impact on contaminant concentrations in biological tissues.
Nilson, R.H.; Lie, K.H. )
1987-12-01
A double-porosity model is used to describe the oscillatory gas motion and associated contaminant transport induced by cyclical variations in the barometric pressure at the surface of a fractured porous medium. Flow along the fractures and within the permeable matrix blocks is locally one-dimensional. The interaction between fractures and blocks includes the Darcian seepage of fluid as well as the Fickian diffusion of contaminant. To guard against artificial numerical diffusion, the FRAM filtering remedy and methodology of Chapman is used in calculating the advective fluxes along fractures and within blocks. The entire system of equations, including the fracture/matrix interaction terms, is solved by a largely implicit non-computational time step is large compared to the cross-block transit time of Darcian pressure waves. The numerical accuracy is tested by comparison with exact solutions for oscillatory and unidirectional flows, some of which include Darcian seepage or Fickian diffusion interaction between the fracture and the matrix. The method is used to estimate the rate of transport of radioactive gases through the rubblized chimney produced by an underground nuclear explosion.
NASA Astrophysics Data System (ADS)
Ngo-Cong, D.; Mohammed, F. J.; Strunin, D. V.; Skvortsov, A. T.; Mai-Duy, N.; Tran-Cong, T.
2015-06-01
The contaminant transport process governed by the advection-diffusion equation plays an important role in modelling industrial and environmental flows. In this article, our aim is to accurately reduce the 2-D advection-diffusion equation governing the dispersion of a contaminant in a turbulent open channel flow to its 1-D approximation. The 1-D model helps to quickly estimate the horizontal size of contaminant clouds based on the values of the model coefficients. We derive these coefficients analytically and investigate numerically the model convergence. The derivation is based on the centre manifold theory to obtain successively more accurate approximations in a consistent manner. Two types of the average velocity profile are considered: the classical logarithmic profile and the power profile. We further develop the one-dimensional integrated radial basis function network method as a numerical approach to obtain the numerical solutions to both the original 2-D equation and the approximate 1-D equations. We compare the solutions of the original models with their centre-manifold approximations at very large Reynolds numbers. The numerical results obtained from the approximate 1-D models are in good agreement with those of the original 2-D model for both the logarithmic and power velocity profiles.
Merritt, M.L.
1993-01-01
The simulation of the transport of injected freshwater in a thin brackish aquifer, overlain and underlain by confining layers containing more saline water, is shown to be influenced by the choice of the finite-difference approximation method, the algorithm for representing vertical advective and dispersive fluxes, and the values assigned to parametric coefficients that specify the degree of vertical dispersion and molecular diffusion that occurs. Computed potable water recovery efficiencies will differ depending upon the choice of algorithm and approximation method, as will dispersion coefficients estimated based on the calibration of simulations to match measured data. A comparison of centered and backward finite-difference approximation methods shows that substantially different transition zones between injected and native waters are depicted by the different methods, and computed recovery efficiencies vary greatly. Standard and experimental algorithms and a variety of values for molecular diffusivity, transverse dispersivity, and vertical scaling factor were compared in simulations of freshwater storage in a thin brackish aquifer. Computed recovery efficiencies vary considerably, and appreciable differences are observed in the distribution of injected freshwater in the various cases tested. The results demonstrate both a qualitatively different description of transport using the experimental algorithms and the interrelated influences of molecular diffusion and transverse dispersion on simulated recovery efficiency. When simulating natural aquifer flow in cross-section, flushing of the aquifer occurred for all tested coefficient choices using both standard and experimental algorithms. ?? 1993.
Transport of trace contaminants through porous media
NASA Technical Reports Server (NTRS)
Madey, R.
1975-01-01
Research accomplishments in the following areas are discussed: (1) the calibration of the gas chromatograph for acetaldehyde and ethanol; (2) the development of data reduction and analysis methods; (3) the generation and analysis of experimental data for the transport of 100 ppm acetaldehyde through a cylindrical bed packed with activated carbon granules; (4) the generation and analysis of experimental data for the transport of 100 ppm ethanol through a cylindrical bed packed with activated carbon granules; and (5) a comparison of the volume adsorption capacity of activated carbon for 100 ppm concentrations of acetaldehyde, ethanol, and acetone. Activities in progress and planned activities are reviewed.
Effects of vegetation on contaminant transport in surface flows
Green, R.; Govindaraju, R.S.; Erickson, L.E.; Roig, L.
1996-12-31
It is well known that vegetation reduces off-site contamination that would result from surface flows. A significant portion of heavy metal contamination occurs at abandoned mine sites due to sediment movement. The effects of vegetation on sediment transport and surface runoff are reviewed, with an emphasis on factors that can reduce or prevent the movement of such metals in mine tailings. Several mathematical models for sediment transport in surface flows are briefly discussed, including advantages and limitations of the Universal Soil-Loss Equation and CREAMS model. Reported experimental and field data on contaminant transport in surface flows are reviewed and evaluated, as well as studies in treating the bioavailability of heavy metals in attempts to reduce metal phytotoxicity or decreasing the potential for entrance of the metals into the food chain via vegetation. Pollutants of concern include lead, zinc, and cadmium. 55 refs.
Preliminary study of niobium alloy contamination by transport through helium
NASA Technical Reports Server (NTRS)
Scheuermann, Coulson M.; Moore, Thomas J.; Wheeler, Donald R.
1987-01-01
Preliminary tests were conducted to determine if interstitial element transport through a circulating helium working fluid was a potential problem in Brayton and Stirling space power systems. Test specimens exposed to a thermal gradient for up to 3000-hr included Nb-1%Zr, a Sm-Co alloy (referred to as SmCo in this paper), Hiperco 50 steel, and alumina to simulate various engine components of the Brayton and Stirling systems. Results indicate that helium transport of interstitial contaminants can be minimized over a 7-yr life with a monometallic Nb-1%Zr design. Exposure with other materials indicated a potential for interstitial contaminant transport. Determination of contamination kinetics and the effects on structural integrity will require additional testing.
Contaminant transport from Elliott and Commencement Bays. Technical memo
Curl, H.C.; Baker, E.T.; Bates, T.S.; Cannon, G.A.; Feely, R.A.
1988-04-01
Contaminant transport from Elliott and Commencement Bays to the main basin of the Puget Sound was investigated by the Pacific Marine Environmental Laboratory during the spring and summer of 1985 and January of 1986. Distributions of water properties (salinity, suspended particulate matter, and toxic trace metals and organics) were mapped during periods of high river runoff and during combined sewer outflow after heavy rainfall. Results indicate that: (1) dissolved contaminants remain in the very thin, fresh-water plume and are transported through the bays into the main basin quite rapidly; (2) PCB and DDT isomers were undetectable in either bay; and (3) there was no evidence that resuspension and transport of contaminated bottom sediments was taking place in Elliot Bay. These results must be qualified due to the short periods during which measurements were taken.
Wagner, B.J.
1992-01-01
Parameter estimation and contaminant source characterization are key steps in the development of a coupled groundwater flow and contaminant transport simulation model. Here a methodologyfor simultaneous model parameter estimation and source characterization is presented. The parameter estimation/source characterization inverse model combines groundwater flow and contaminant transport simulation with non-linear maximum likelihood estimation to determine optimal estimates of the unknown model parameters and source characteristics based on measurements of hydraulic head and contaminant concentration. First-order uncertainty analysis provides a means for assessing the reliability of the maximum likelihood estimates and evaluating the accuracy and reliability of the flow and transport model predictions. A series of hypothetical examples is presented to demonstrate the ability of the inverse model to solve the combined parameter estimation/source characterization inverse problem. Hydraulic conductivities, effective porosity, longitudinal and transverse dispersivities, boundary flux, and contaminant flux at the source are estimated for a two-dimensional groundwater system. In addition, characterization of the history of contaminant disposal or location of the contaminant source is demonstrated. Finally, the problem of estimating the statistical parameters that describe the errors associated with the head and concentration data is addressed. A stage-wise estimation procedure is used to jointly estimate these statistical parameters along with the unknown model parameters and source characteristics. ?? 1992.
Subsurface barrier design alternatives for confinement and controlled advection flow
Phillips, S.J.; Stewart, W.E.; Alexander, R.G.; Cantrell, K.J.; McLaughlin, T.J.
1994-02-01
Various technologies and designs are being considered to serve as subsurface barriers to confine or control contaminant migration from underground waste storage or disposal structures containing radioactive and hazardous wastes. Alternatives including direct-coupled flood and controlled advection designs are described as preconceptual examples. Prototype geotechnical equipment for testing and demonstration of these alternative designs tested at the Hanford Geotechnical Development and Test Facility and the Hanford Small-Tube Lysimeter Facility include mobile high-pressure injectors and pumps, mobile transport and pumping units, vibratory and impact pile drivers, and mobile batching systems. Preliminary laboratory testing of barrier materials and additive sequestering agents have been completed and are described.
NASA Astrophysics Data System (ADS)
Baechler, S.; Croisé, J.; Altmann, S.
2012-12-01
Chemico-osmosis is a recognized phenomenon taking place in clay mineral-rich sedimentary formations and a number of questions have been raised concerning its potential effects on pressure fields in and around underground radioactive waste repositories installed in such formations. Certain radioactive waste packages contain large quantities of nitrate salts whose release might result in the presence of highly concentrated salt solutions in the disposal cells, during their resaturation after closure of the facility. This would lead to large solute concentration gradients within the formation's porewater which could then potentially induce significant chemico-osmotic fluxes. In this paper, we assess the impact of chemico-osmotic fluxes on the water pressure during the post-closure period of a typical disposal cell for intermediate-level, long-lived bituminised radioactive waste in the Callovo-Oxfordian Clay formation. A numerical model of chemico-osmotic water flow and solute transport has been developed based on the work of Bader and Kooi (2005) [5], and including Bresler's dependence of osmotic efficiency on concentration and compaction state [9]. Model validity has been extended to highly concentrated solutions by incorporating a concentration-dependent activity coefficient, based on the Pitzer's equations. Results show that due to the strong dependence of the osmotic coefficient on concentration, the impact of chemico-osmosis on water flow and on the pressure field around the disposal cell is relatively low. A maximum overpressure of the order of 1 MPa was obtained. No difference in the simulation results were noticed for disposal cell solutions having concentrations higher than 1 M NaNO3. Differences between simulations were found to be almost entirely due to Bresler's relationship i.e., the model of the dependence between osmotic efficiency and concentration, and only slightly on the activity coefficient correction. Questions remain regarding the appropriate
Zhong, L; Szecsody, J E; Truex, M J; Williams, M D; Liu, Y
2015-05-30
Use of gas-phase amendments for in situ remediation of inorganic contaminants in unsaturated sediments of the vadose zone may be advantageous, but there has been limited development and testing of gas remediation technologies. Treatment with ammonia gas has a potential for use in treating inorganic contaminants (such as uranium) because it induces a high pore-water pH, causing mineral dissolution and subsequent formation of stable precipitates that decrease the mobility of some contaminants. For field application of this treatment, further knowledge of ammonia transport in porous media and the geochemical reactions induced by ammonia treatment is needed. Laboratory studies were conducted to support calculations needed for field treatment design, to quantify advective and diffusive ammonia transport in unsaturated sediments, to evaluate inter-phase (gas/sediment/pore water) reactions, and to study reaction-induced pore-water chemistry changes as a function of ammonia delivery conditions, such as flow rate, gas concentration, and water content. Uranium-contaminated sediment was treated with ammonia gas to demonstrate U immobilization. Ammonia gas quickly partitions into sediment pore water and increases the pH up to 13.2. Injected ammonia gas advection front movement can be reasonably predicted by gas flow rate and equilibrium partitioning. The ammonia gas diffusion rate is a function of the water content in the sediment. Sodium, aluminum, and silica pore-water concentrations increase upon exposure to ammonia and then decline as aluminosilicates precipitate when the pH declines due to buffering. Up to 85% of the water-leachable U was immobilized by ammonia treatment. PMID:25723886
The effects of a perturbed source on contaminant transport near the Weldon Spring quarry
Tomasko, D.
1989-03-01
The effects of a perturbed contamination source at the Weldon Spring quarry in St. Charles County, Missouri, on downstream solute concentrations were investigated using one-dimensional analytical solutions to an advection-dispersion equation developed for both constant-strength and multiple-stepped source functions. A sensitivity study using parameter base-case values and ranges consistent with the geologic conceptualization of the quarry area indicates that the parameters having the greatest effect on predicted concentrations are the distance from the quarry to the point of interest, the average linear groundwater velocity, the contaminant retardation coefficient, and the amplitude and duration of the source perturbation caused by response action activities. Use of base-case parameter value and realistic values for the amplitude and duration of the source perturbation produced a small effect on solute concentrations near the western extremity of the nearby municipal well field, as well as small uncertainties in the predicted results for the assumed model. The effect of simplifying assumptions made in deriving the analytic solution is unknown: use of a multidimensional flow and transport model and additional field work are needed to validate the model. 13 refs., 18 figs.
NASA Astrophysics Data System (ADS)
Stremler, Mark A.; Haselton, F. R.; Aref, Hassan
2004-05-01
Chaotic advection can play an important role in efficient microfluidic mixers. We discuss a design paradigm that exploits chaotic advection and illustrate by two recent examples, namely enhancing gene expression profiling and constructing an in-line microfluidic mixing channel, how application of this paradigm has led to successful micromixers. We suggest that 'designing for chaos', that is, basing practical mixer design on chaotic advection analysis, is a promising approach to adopt in this developing field which otherwise has little to guide it and is constrained by issues of scale and manufacturability.
Rockhold, Mark L.; White, Mark D.; Freeman, Eugene J.
2004-10-12
This letter report documents initial numerical analyses conducted by PNNL to provide support for a feasibility study on decommissioning of the canyon buildings at Hanford. The 221-U facility is the first of the major canyon buildings to be decommissioned. The specific objective of this modeling effort was to provide estimates of potential rates of migration of residual contaminants out of the 221-U facility during the first 40 years after decommissioning. If minimal contaminant migration is predicted to occur from the facility during this time period, then the structure may be deemed to provide a level of groundwater protection that is essentially equivalent to the liner and leachate collection systems that are required at conventional landfills. The STOMP code was used to simulate transport of selected radionuclides out of a canyon building, representative of the 221-U facility after decommissioning, for a period of 40 years. Simulation results indicate that none of the selected radionuclides that were modeled migrated beyond the concrete structure of the facility during the 40-year period of interest. Jacques (2001) identified other potential contaminants in the 221-U facility that were not modeled, however, including kerosene, phenol, and various metals. Modeling of these contaminants was beyond the scope of this preliminary effort due to increased complexity. Simulation results indicate that contaminant release from the canyon buildings will be diffusion controlled at early times. Advection is expected to become much more important at later times, after contaminants have diffused out of the facility and into the surrounding soil environment. After contaminants have diffused out of the facility, surface infiltration covers will become very important for mitigating further transport of contaminants in the underlying vadose zone and groundwater.
CONTAMINANT TRANSPORT IN PARALLEL FRACTURED MEDIA: SUDICKY AND FRIND REVISITED
This paper is concerned with a modified, nondimensional form of the parallel fracture, contaminant transport model of Sudicky and Frind (1982). The modifications include the boundary condition at the fracture wall, expressed by a parameter, and the power-law relationship between...
STABILITY AND TRANSPORT OF INORGANIC COLLOIDS THROUGH CONTAMINATED AQUIFER MATERIAL
Laboratory columns using contaminated natural aquifer material from Globe, Arizona, were used to investigate the transport of inorganic colloids under saturated flow conditions. e2O3 radio-labeled spherical colloids of various diameters were synthesized and introduced into the co...
INVESTIGATION OF CONTAMINANT TRANSPORT FROM THE SAGINAW CONFINED DISPOSAL FACILITY
Pilot biomonitoring and monitoring studies were conducted at the Saginaw Confined Disposal Facility (CDF), Saginaw Bay, Lake Huron, during 1987 to develop methods to assess the potential for magnitude of 1) contaminant transport from the dike interior to the outside environment, ...
CONTAMINANT TRANSPORT IN PARALLEL FRACTURED MEDIA: SUDICKY AND FRIND REVISITED
This paper is concerned with a modified, nondimensional form of the parallel fracture, contaminant transport model of Sudicky and Frind (1982). The modifications include the boundary condition at the fracture wall, expressed by a parameter , and the power-law relationship betwe...
Energy Science and Technology Software Center (ESTSC)
1995-07-07
DUST solves for release and transport of contaminants from containerized wastes. Each container may have unique properties (i.e., time to failure or localized failure, e.g., pitting) and each waste form may have unique release properties. Release from the waste form is limited by one of four physical or chemical restraints: solubility, diffusion, dissolution, and surface wash-off with partitioning. The release from the waste form acts as a source for transport in the advection/dispersion equation. Transportmore » is modeled in one-dimension through the groundwater pathway from subsurface disposal. RNUCL.DAT, database of half-lives, solubility limits, and atomic mass for selected radionuclides, is included in this package.« less
A contaminant transport model for wetlands accounting for distinct residence time bimodality
NASA Astrophysics Data System (ADS)
Musner, T.; Bottacin-Busolin, A.; Zaramella, M.; Marion, A.
2014-07-01
Vegetation plays a major role in controlling the fate of contaminants in natural and constructed wetlands. Estimating the efficiency of contaminant removal of a wetland requires separate knowledge of the residence time statistics in the main flow channels, where the flow velocity is relatively higher, and in the more densely vegetated zones, where the velocity is smaller and most of the biochemical transformations occur. A conceptual wetland characterized by a main flow channel (MFC) and lateral vegetated zones (LVZs) is modeled here using a two-dimensional depth-averaged hydrodynamic and advection-dispersion model. The effect of vegetation is described as a flow resistance represented in the hydrodynamic model as a function of the stem density. Simulations are performed for a given flow discharge and for increasing values of the ratio between the vegetation density in the LVZs and in the MFC. Residence time distributions (RTDs) of a nonreactive tracer are derived from numerical simulations of the solute breakthrough curves (BTCs) resulting from a continuous concentration input. Results show that increasing vegetation densities produce an increasingly pronounced bimodality of the RTDs. At longer times, the RTDs decrease exponentially, with different timescales depending on the stem density ratio and other system parameters. The overall residence time distribution can be decomposed into a first component associated with the relatively fast transport in the MFC, and a second component associated with the slower transport in the LVZs. The weight of each temporal component is related to the exchange flux at the MFC-LVZ interface. A one-dimensional transport model is proposed that is capable to reproduce the RTDs predicted by the depth-averaged model, and the relationship between model and system parameters is investigated using a combination of direct and inverse modeling approaches.
Contaminant Transport Through Subsurface Material from the DOE Hanford Reservation
Pace, M.N.; Mayes, M.A.; Jardine, P.M.; Fendorf, S.E.; Nehlhorn, T.L.; Yin, X.P.; Ladd, J.; Teerlink, J.; Zachara, J.M.
2003-03-26
Accelerated migration of contaminants in the vadose zone has been observed beneath tank farms at the U.S. Department of Energy's Hanford Reservation. This paper focuses on the geochemical processes controlling the fate and transport of contaminants in the sediments beneath the Hanford tank farms. Laboratory scale batch sorption experiments and saturated transport experiments were conducted using reactive tracers U(VI), Sr, Cs, Co and Cr(VI) to investigate geochemical processes controlling the rates and mechanisms of sorption to Hanford subsurface material. Results indicate that the rate of sorption is influenced by changes in solution chemistry such as ionic strength, pH and presence of competing cations. Sediment characteristics such as mineralogy, iron content and cation/anion exchange capacity coupled with the dynamics of flow impact the number of sites available for sorption. Investigative approaches using a combination of batch and transport experiments will contribute to the conceptual and Hanford vadose zone.
Modelling contaminant transport for pumping wells in riverbank filtration systems.
Mustafa, Shaymaa; Bahar, Arifah; Aziz, Zainal Abdul; Suratman, Saim
2016-01-01
Analytical study of the influence of both the pumping well discharge rate and pumping time on contaminant transport and attenuation is significant for hydrological and environmental science applications. This article provides an analytical solution for investigating the influence of both pumping time and travelling time together for one-dimensional contaminant transport in riverbank filtration systems by using the Green's function approach. The basic aim of the model is to understand how the pumping time and pumping rate, which control the travelling time, can affect the contaminant concentration in riverbank filtration systems. Results of analytical solutions are compared with the results obtained using a MODFLOW numerical model. Graphically, it is found that both analytical and numerical solutions have almost the same behaviour. Additionally, the graphs indicate that any increase in the pumping rate or simulation pumping time should increase the contamination in groundwater. The results from the proposed analytical model are well matched with the data collected from a riverbank filtration site in France. After this validation, the model is then applied to the first pilot project of a riverbank filtration system conducted in Malaysia. Sensitivity analysis results highlight the importance of degradation rates of contaminants on groundwater quality, for which higher utilization rates lead to the faster consumption of pollutants. PMID:26433356
Concentration distribution of contaminant transport in wetland flows
NASA Astrophysics Data System (ADS)
Wu, Zi; Fu, Xudong; Wang, Guangqian
2015-06-01
Study on contaminant transport in wetland flows is of fundamental importance. Recent investigation on scalar transport in laminar tube flows (Wu and Chen, 2014. J. Fluid Mech., 740: 196-213.) indicates that the vertical concentration difference in wetland flows may be remarkable for a very long time, which cannot be captured by the extensively applied one-dimensional Taylor dispersion model. To understand detailed information for the vertical distribution of contaminant in wetland flows, for the first time, the present paper deduces an analytical solution for the multi-dimensional concentration distribution by the method of mean concentration expansion. The solution is verified by both our analytical and numerical results. Representing the effects of vegetation in wetlands, the unique dimensionless parameter α can cause the longitudinal contraction of the contaminant cloud and the change of the shape of the concentration contours. By these complicated effects, it is shown unexpectedly that the maximum vertical concentration difference remains nearly unaffected, although its longitudinal position may change. Thus the slow-decaying transient effect (Wu and Chen, 2014. J. Hydrol., 519: 1974-1984.) is shown also apply to the process of contaminant transport in wetland flows.
NASA Astrophysics Data System (ADS)
Makedonska, N.; Karra, S.; Painter, S. L.; Viswanathan, H.; Gable, C. W.
2014-12-01
's Manual: A Massively Parallel Reactive Flow andTransport Model for describing Surface and Subsurface Processes, 2014. [3] Makedonska N., Painter S.L., Karra S., and Gable C.W., NumericalExperiments on Advective Transport in Large Three-Dimensional DFNs,Abstract H53A-1398 ,2013, AGU, San-Francisco, CA, 9-13 Dec.
Monitoring Potential Transport of Radioactive Contaminants in Shallow Ephemeral Channels
Miller, Julianne J.; Mizell, Steve A.; Nikolich, George; Campbell, Scott A.
2012-02-01
The U.S. Department of Energy (DOE) National Nuclear Security Administration (NNSA), Nevada Site Office (NSO), Environmental Restoration Soils Activity has authorized the Desert Research Institute (DRI) to conduct field assessments of potential sediment transport of contaminated soil from Corrective Action Unit (CAU) 550, Area 8 Smoky Contamination Area (CA), during precipitation runoff events. CAU 550 includes Corrective Action Sites (CASs) 08-23-03, 08-23-04, 08-23-06, and 08-23-07; these CASs are associated with tests designated Ceres, Smoky, Oberon, and Titania, respectively.
Analysis Of Multispectral Imagery And Modeling Contaminant Transport
NASA Astrophysics Data System (ADS)
Irvine, J. M.; Becker, N. M.; Brumby, S.; David, N. A.
2003-12-01
A significant concern in the monitoring of hazardous waste is the potential for contaminants to migrate into locations where their presence poses greater environmental risks. The transport modeling performed in this study demonstrates the joint use of remotely sensed multispectral imagery and mathematical modeling to assess the surface migration of contaminants. KINEROS, an event-driven model of surface runoff and sediment transport, was used to assess uranium transport for various rain events. While our specific application was uranium transport, the methods apply to surface transport of any substance of concern. The model inputs include parameters related to the size and slope of watershed components, vegetation, and soil conditions. One distinct set of model inputs was derived from remotely sensed imagery data and another from site-specific knowledge. To derive the parameters of the KINEROS model from remotely sensed data, classification analysis was performed on IKONOS four-band multispectral imagery of the watershed. A system known as GENIE, developed by Los Alamos National Laboratory, employs genetics algorithms to evolve classifiers based on small, user-selected training samples. The classification analysis derived by employing GENIE provided insight into the correct KINEROS parameters for various sub-elements of the watershed. The model results offer valuable information about portions of the watershed that contributed the most to contaminant transport. These methods are applicable to numerous sites where possible transport of waste materials or other hazardous substances poses an environmental risk. Consequently, the approach presented here is relevant to homeland security and emergency response scenarios, as well as long-term environmental monitoring applications. Because the approach rests on the analysis of remote sensing data, the techniques can be used to monitor a range of sites and can reduce costs of data collection for model calibration.
Multimedia transport of organic contaminants and exposure modeling
Layton, D.W.; McKone, T.E.
1988-01-01
Human exposures to organic contaminants in the environment are a complex function of human factors, physicochemical properties of the contaminants, and characteristics of the environmental media in which the contaminants reside. One subject of interest in the screening of organic chemicals for the purpose of identifying exposure pathways of potential concern is the relationship between exposures and contaminant properties. To study such relationships, a multimedia environmental model termed GEOTOX is used to predict the equilibrium partitioning and transport of ''reference'' organic chemicals between compartments representing different media (i.e., soil layers, ground water, air, biota, etc.) of a contaminated landscape. Reference chemicals, which are added to the surface soil of a landscape, are defined by properties consisting of the Henry's law constant, soil water-soil organic carbon partition coefficient, and bioconcentration factors. The steady-state concentrations of the chemical in the GEOTOX compartments are then used to estimate lifetime exposures (in mg/kg-d) to the compartments for individuals living in the contaminated landscape. Exposure pathways addressed include ingestion, inhalation, and dermal absorption. Local sensitivity analyses are performed to determine which chemical and landscape properties have the greatest effect on the exposure estimates. 9 refs., 4 figs., 3 tabs.
NASA Astrophysics Data System (ADS)
Wilkins, M. J.; Pearce, C.; Zhang, C.; Heald, S.; Fredrickson, J. K.; Zachara, J. M.
2011-12-01
Microenvironments and transition zones dominate the subsurface biogeochemical cycling of key contaminants, with strong effects resulting from the coupling of chemical reactions, physical transport and microbiological processes. Understanding the impact of pore-scale environments (e.g. spatial heterogeneity, chemical gradients, and redox potential) is essential for modeling contaminant fate and transport in the subsurface. The driver for biogeochemical processes at the pore scale changes from macroscopic advection to microscale diffusion, and this has a significant effect on the retention of soluble, highly mobile contaminants such as U(VI). Here, etched-silicon microfluidic models with defined chemistry, mineralogy, microbiology, and flow regimes are used for the incremental development of complex microenvironments that approach real-world systems. We demonstrate the colonization of such pore spaces by an anaerobic Fe(III)-reducing bacterium, the enzymatic reduction of a bioavailable Fe(III) phase within this environment, and the subsequent effects of both oxidized and reduced Fe phases on uranium biogeochemistry under flow conditions using both X-ray Microprobe (XMP) and X-ray Absorption Spectroscopy (XAS). Precipitated Fe(III) phases within the microfluidic model were most effectively reduced in the presence of an electron shuttle (e.g. AQDS), with Fe(II) ions adsorbing onto mineral precipitates and surfaces. In the absence of Fe, U(VI) was effectively reduced by the microbial population to insoluble U(IV), which was precipitated in discrete regions associated within biomass. In the presence of both oxidized and reduced Fe phases however, differing effects were observed with regards to U behavior; oxidized U(VI) was frequently adsorbed to poorly crystalline Fe(III), and reduced U(IV) associated with more reduced regions of the microscale flow cell. In the future, the flexibility in the design of the microfluidic models, in combination with advanced
IMPACT OF REDOX DISEQUILIBRIA ON CONTAMINANT TRANSPORT AND REMEDIATION IN SUBSURFACE SYSTEMS
Partitioning to mineral surfaces exerts significant control on inorganic contaminant transport in subsurface systems. Remedial technologies for in-situ treatment of subsurface contamination are frequently designed to optimize the efficiency of contaminant partitioning to solid s...
Transport mechanisms of contaminants released from fine sediment in rivers
NASA Astrophysics Data System (ADS)
Cheng, Pengda; Zhu, Hongwei; Zhong, Baochang; Wang, Daozeng
2015-12-01
Contaminants released from sediment into rivers are one of the main problems to study in environmental hydrodynamics. For contaminants released into the overlying water under different hydrodynamic conditions, the mechanical mechanisms involved can be roughly divided into convective diffusion, molecular diffusion, and adsorption/desorption. Because of the obvious environmental influence of fine sediment (D_{90}= 0.06 mm), non-cohesive fine sediment, and cohesive fine sediment are researched in this paper, and phosphorus is chosen for a typical adsorption of a contaminant. Through theoretical analysis of the contaminant release process, according to different hydraulic conditions, the contaminant release coupling mathematical model can be established by the N-S equation, the Darcy equation, the solute transport equation, and the adsorption/desorption equation. Then, the experiments are completed in an open water flume. The simulation results and experimental results show that convective diffusion dominates the contaminant release both in non-cohesive and cohesive fine sediment after their suspension, and that they contribute more than 90 % of the total release. Molecular diffusion and desorption have more of a contribution for contaminant release from unsuspended sediment. In unsuspension sediment, convective diffusion is about 10-50 times larger than molecular diffusion during the initial stages under high velocity; it is close to molecular diffusion in the later stages. Convective diffusion is about 6 times larger than molecular diffusion during the initial stages under low velocity, it is about a quarter of molecular diffusion in later stages, and has a similar level with desorption/adsorption. In unsuspended sediment, a seepage boundary layer exists below the water-sediment interface, and various release mechanisms in that layer mostly dominate the contaminant release process. In non-cohesive fine sediment, the depth of that layer increases linearly with shear
Amanzi and Akuna: Two New Community Codes for Subsurface Contaminant Flow and Transport
NASA Astrophysics Data System (ADS)
Dixon, P. R.; Moulton, J. D.; Gorton, I.; Meza, J.; Freshley, M.
2011-12-01
The Advanced Simulation Capability for Environmental Management (ASCEM) program is developing a modular and extensible open-source set of tools for understanding the fate and transport of contaminants in natural and engineered systems. These tools not only support a fundamental shift toward standardized assessments of performance and risk for the Department of Energy Office of Environmental Management (DOE-EM) cleanup and closure decisions, but establish a modern high-quality code base for a growing interdisciplinary community. Specifically, ASCEM is leveraging advances and expertise from applied mathematics, computer and computational sciences, and the geosciences, in this new development. A toolset named Akuna will provide capabilities for data management, visualization, conceptual model development, uncertainty quantification, parameter estimation, risk analysis, and decision support. Akuna will integrate with Amanzi, a flexible high performance computing simulator, which is designed to leverage the growing parallelism in modern systems. This talk will describe the approach that we have taken to develop this new open-source capability, including issues of intellectual property, licensing, the developers tool chain, and the users tool chain. The modular and extensible design will be discussed, highlighting the potential for collaboration and inclusion of recent modeling and algorithmic advances. In addition, it will discuss the advantages and challenges of relying on an open-source model that leverages a wide variety of open-source efforts from other programs. Results from early prototype development will be presented to highlight the potential of these new tools to contaminated subsurface environments, including calculations for variably saturated flow, advection of non-reactive species and the reactive-transport of 17 different chemical species on both structured and unstructured meshes.
Mass transfer model of nanoparticle-facilitated contaminant transport in saturated porous media.
Johari, Wan Lutfi Wan; Diamessis, Peter J; Lion, Leonard W
2010-02-01
A one-dimensional model has been evaluated for transport of hydrophobic contaminants, such as polycyclic aromatic hydrocarbon (PAH) compounds, facilitated by synthetic amphiphilic polyurethane (APU) nanoparticles in porous media. APU particles synthesized from poly(ethylene glycol)-modified urethane acrylate (PMUA) precursor chains have been shown to enhance the desorption rate and mobility of phenanthrene (PHEN) in soil. A reversible process governed by attachment and detachment rates was considered to describe the PMUA binding in soil in addition to PMUA transport through advection and dispersion. Ultimately, an irreversible second-order PMUA attachment rate in which the fractional soil saturation capacity with PMUA was a rate control was found to be adequate to describe the retention of PMUA particles. A gamma-distributed site model (GS) was used to describe the spectrum of physical/chemical constraints for PHEN transfer from solid to aqueous phases. Instantaneous equilibrium was assumed for PMUA-PHEN interactions. The coupled model for PMUA and PHEN behavior successfully described the enhanced elution profile of PHEN by PMUA. Sensitivity analysis was performed to analyze the significance of model parameters on model predictions. The adjustable parameter alpha in the gamma-distribution shapes the contaminant desorption distribution profile as well as elution and breakthrough curves. Model simulations show the use of PMUA can be also expected to improve the release rate of PHEN in soils with higher organic carbon content. The percentage removal of PHEN mass over time is shown to be influenced by the concentration of PMUA added and this information can be used to optimize cost and time require to accomplish a desired remediation goal. PMID:19406449
NASA Astrophysics Data System (ADS)
Parker, Jack C.; Kim, Ungtae
2015-11-01
The mono-continuum advection-dispersion equation (mADE) is commonly regarded as unsuitable for application to media that exhibit rapid breakthrough and extended tailing associated with diffusion between high and low permeability regions. This paper demonstrates that the mADE can be successfully used to model such conditions if certain issues are addressed. First, since hydrodynamic dispersion, unlike molecular diffusion, cannot occur upstream of the contaminant source, models must be formulated to prevent "back-dispersion." Second, large variations in aquifer permeability will result in differences between volume-weighted average concentration (resident concentration) and flow-weighted average concentration (flux concentration). Water samples taken from wells may be regarded as flux concentrations, while soil samples may be analyzed to determine resident concentrations. While the mADE is usually derived in terms of resident concentration, it is known that a mADE of the same mathematical form may be written in terms of flux concentration. However, when solving the latter, the mathematical transformation of a flux boundary condition applied to the resident mADE becomes a concentration type boundary condition for the flux mADE. Initial conditions must also be consistent with the form of the mADE that is to be solved. Thus, careful attention must be given to the type of concentration data that is available, whether resident or flux concentrations are to be simulated, and to boundary and initial conditions. We present 3-D analytical solutions for resident and flux concentrations, discuss methods of solving numerical models to obtain resident and flux concentrations, and compare results for hypothetical problems. We also present an upscaling method for computing "effective" dispersivities and other mADE model parameters in terms of physically meaningful parameters in a diffusion-limited mobile-immobile model. Application of the latter to previously published studies of
Parker, Jack C; Kim, Ungtae
2015-11-01
The mono-continuum advection-dispersion equation (mADE) is commonly regarded as unsuitable for application to media that exhibit rapid breakthrough and extended tailing associated with diffusion between high and low permeability regions. This paper demonstrates that the mADE can be successfully used to model such conditions if certain issues are addressed. First, since hydrodynamic dispersion, unlike molecular diffusion, cannot occur upstream of the contaminant source, models must be formulated to prevent "back-dispersion." Second, large variations in aquifer permeability will result in differences between volume-weighted average concentration (resident concentration) and flow-weighted average concentration (flux concentration). Water samples taken from wells may be regarded as flux concentrations, while soil samples may be analyzed to determine resident concentrations. While the mADE is usually derived in terms of resident concentration, it is known that a mADE of the same mathematical form may be written in terms of flux concentration. However, when solving the latter, the mathematical transformation of a flux boundary condition applied to the resident mADE becomes a concentration type boundary condition for the flux mADE. Initial conditions must also be consistent with the form of the mADE that is to be solved. Thus, careful attention must be given to the type of concentration data that is available, whether resident or flux concentrations are to be simulated, and to boundary and initial conditions. We present 3-D analytical solutions for resident and flux concentrations, discuss methods of solving numerical models to obtain resident and flux concentrations, and compare results for hypothetical problems. We also present an upscaling method for computing "effective" dispersivities and other mADE model parameters in terms of physically meaningful parameters in a diffusion-limited mobile-immobile model. Application of the latter to previously published studies of
Simulation of contaminated sediment transport in White Oak Creek basin
Bao, Y.; Clapp, R.B.; Brenkert, A.L.; Moore, T.D.; Fontaine, T.A.
1995-12-31
This paper presents a systematic approach to management of the contaminated sediments in the White Oak Creek watershed at Oak Ridge National Laboratory near Oak Ridge, Tennessee. The primary contaminant of concern is radioactive cesium-137 ({sup 137}Cs), which binds to soil and sediment particles. The key components in the approach include an intensive sampling and monitoring system for flood events; modeling of hydrological processes, sediment transport, and contaminant flux movement; and a decision framework with a detailed human health risk analysis. Emphasis is placed on modeling of watershed rainfall-runoff and contaminated sediment transport during flooding periods using the Hydrologic Simulation Program- Fortran (HSPF) model. Because a large number of parameters are required in HSPF modeling, the major effort in the modeling process is the calibration of model parameters to make simulation results and measured values agree as closely as possible. An optimization model incorporating the concepts of an expert system was developed to improve calibration results and efficiency. Over a five-year simulation period, the simulated flows match the observed values well. Simulated total amount of sediment loads at various locations during storms match with the observed values within a factor of 1.5. Simulated annual releases of {sup 137}Cs off-site locations match the data within a factor of 2 for the five-year period. The comprehensive modeling approach can provide a valuable tool for decision makers to quantitatively analyze sediment erosion, deposition, and transport; exposure risk related to radionuclides in contaminated sediment; and various management strategies.
Prediction methodology for contaminant transport from rangeland watersheds
Devaurs, M.A.; Springer, E.P.; Lane, L.J.; Langhorst, G.J.
1988-01-01
Weather on arid and semiarid lands can be extremely variable. Runoff is generally emphermeral, and high intensity, short-duration rainfall events are the major stimulus for runoff events. Transport of sediment and associated contaminants occurs with these infrequent events. Incorporation of variability in weather into any prediction technology is essential to provide accurate representations of climate-induced uncertainty in predictions of hydrologic response. The objective of this study is to investigate a method for including short-term climatic variations in analyses for contaminant transport from rangeland watersheds in arid/semiarid regions. Short term is defined here as a twenty to fifty time frame and it is assumed that lone term climatic fluctuations are not observed during this time. Also, most weather records are available for this time period; predictions of greater length are extrapolations of existing records unless corroborative data for longer term trends are collected. Predictions are being made with condensable uncertainty in the weather inputs even if the models for water, sediment, and contaminant transport are perfectly unknown. This study will incorporate uncertainty in weather inputs into the prediction process and address the ramifications of this uncertainty. Uncertainty introduced by improper model or parameter specification is only briefly addressed.
Modeling of contaminant transport in underground coal gasification
Lanhe Yang; Xing Zhang
2009-01-15
In order to study and discuss the impact of contaminants produced from underground coal gasification on groundwater, a coupled seepage-thermodynamics-transport model for underground gasification was developed on the basis of mass and energy conservation and pollutant-transport mechanisms, the mathematical model was solved by the upstream weighted multisell balance method, and the model was calibrated and verified against the experimental site data. The experiment showed that because of the effects of temperature on the surrounding rock of the gasification panel the measured pore-water-pressure was higher than the simulated one; except for in the high temperature zone where the simulation errors of temperature, pore water pressure, and contaminant concentration were relatively high, the simulation values of the overall gasification panel were well fitted with the measured values. As the gasification experiment progressed, the influence range of temperature field expanded, the gradient of groundwater pressure decreased, and the migration velocity of pollutant increased. Eleven months and twenty months after the test, the differences between maximum and minimum water pressure were 2.4 and 1.8 MPa, respectively, and the migration velocities of contaminants were 0.24-0.38 m/d and 0.27-0.46 m/d, respectively. It was concluded that the numerical simulation of the transport process for pollutants from underground coal gasification was valid. 42 refs., 13 figs., 1 tab.
Attenuation of Selected Emerging Contaminants During River Transport
NASA Astrophysics Data System (ADS)
Reinhard, M.; Gross, B.; Hadeler, A.
2002-12-01
The ubiquitous occurrence of emerging (non-regulated) contaminants in the aquatic environment is of concern because some of these chemicals are biologically active at low concentrations and a potential threat to wildlife and human health.. Emerging contaminants include a diverse range of chemicals, including pharmaceuticals, natural and synthetic hormones and industrial surfactants, such as alkylphenol ethoxylates (APEO) and their metabolites. To address the ecotoxicological impact of these chemicals, it is necessary to know their sources, removal efficiencies during wastewater treatment, and their behavior in the environment. In this study, the fate of selected emerging contaminants in the Santa Ana River (SAR) in Southern California was investigated. The SAR originates in the San Bernardino Mountains and flows 80 miles into the Pacific Ocean. The SAR flow stems mainly from storm runoff, wastewater treatment effluents and several other minor sources. During the dry season, SAR flow is dominated by effluent from public wastewater treatment plants. Input into the SAR was studied by analyzing samples from four major treatment plants that employ different tertiary treatment processes. To assess the fate during river water transport and during wetland treatment, samples from six sites along the river were analyzed. Effluent samples were analyzed every two months, river water every four months. River samples were taken considering the flow velocity, which is approximately 1 mile per hour. The analytical method involves solid-phase extraction using C-18 cartridges and extraction of three fractions. Samples were analyzed with and without further derivatization using GC/MS and GC/MS/MS. Results indicate significant contaminant removal during river transport, presumably by photochemical oxidation. Within a distance of nine miles, pharmaceuticals, plasticizers, flame retardants, APEOs and metabolites were attenuated with removal rates ranging from 76% for a flame retardant
Varank, Gamze; Demir, Ahmet; Yetilmezsoy, Kaan; Bilgili, M. Sinan; Top, Selin; Sekman, Elif
2011-11-15
Highlights: > We conduct 1D advection-dispersion modeling to estimate transport parameters. > We examine fourteen phenolic compounds and three inorganic contaminants. > 2-MP, 2,4-DCP, 2,6-DCP, 2,4,5-TCP, 2,3,4,6-TeCP have the highest coefficients. > Dispersion coefficients of Cu are determined to be higher than Zn and Fe. > Transport of phenolics can be prevented by zeolite and bentonite in landfill liners. - Abstract: One-dimensional (1D) advection-dispersion transport modeling was conducted as a conceptual approach for the estimation of the transport parameters of fourteen different phenolic compounds (phenol, 2-CP, 2-MP, 3-MP, 4-MP, 2-NP, 4-NP, 2,4-DNP, 2,4-DCP, 2,6-DCP, 2,4,5-TCP, 2,4,6-TCP, 2,3,4,6-TeCP, PCP) and three different inorganic contaminants (Cu, Zn, Fe) migrating downward through the several liner systems. Four identical pilot-scale landfill reactors (0.25 m{sup 3}) with different composite liners (R1: 0.10 + 0.10 m of compacted clay liner (CCL), L{sub e} = 0.20 m, k{sub e} = 1 x 10{sup -8} m/s, R2: 0.002-m-thick damaged high-density polyethylene (HDPE) geomembrane overlying 0.10 + 0.10 m of CCL, L{sub e} = 0.20 m, k{sub e} = 1 x 10{sup -8} m/s, R3: 0.002-m-thick damaged HDPE geomembrane overlying a 0.02-m-thick bentonite layer encapsulated between 0.10 + 0.10 m CCL, L{sub e} = 0.22 m, k{sub e} = 1 x 10{sup -8} m/s, R4: 0.002-m-thick damaged HDPE geomembrane overlying a 0.02-m-thick zeolite layer encapsulated between 0.10 + 0.10 m CCL, L{sub e} = 0.22 m, k{sub e} = 4.24 x 10{sup -7} m/s) were simultaneously run for a period of about 540 days to investigate the nature of diffusive and advective transport of the selected organic and inorganic contaminants. The results of 1D transport model showed that the highest molecular diffusion coefficients, ranging from 4.77 x 10{sup -10} to 10.67 x 10{sup -10} m{sup 2}/s, were estimated for phenol (R4), 2-MP (R1), 2,4-DNP (R2), 2,4-DCP (R1), 2,6-DCP (R2), 2,4,5-TCP (R2) and 2,3,4,6-TeCP (R1). For all reactors
NASA Astrophysics Data System (ADS)
Tóth, Gábor; Keppens, Rony
2012-07-01
The Versatile Advection Code (VAC) is a freely available general hydrodynamic and magnetohydrodynamic simulation software that works in 1, 2 or 3 dimensions on Cartesian and logically Cartesian grids. VAC runs on any Unix/Linux system with a Fortran 90 (or 77) compiler and Perl interpreter. VAC can run on parallel machines using either the Message Passing Interface (MPI) library or a High Performance Fortran (HPF) compiler.
Characterization of Anomalous Contaminant Transport via Push-Pull Tracer Tests
NASA Astrophysics Data System (ADS)
Hansen, S. K.; Vesselinov, V. V.; Berkowitz, B.
2015-12-01
Push-pull (single-well-injection-withdrawal) tracer tests are widely used as an economical means of characterizing field-scale solute transport properties such as sorption and dispersion. Typically, these are analyzed by means of analytic solutions that assume transport obeys the radial advection-dispersion equation. We revisit this approach as: (1) Recognition of the ubiquity of anomalous transport and its impact on contaminant remediation necessitates the use of new methods to characterize it, and (2) Improved computational power and numerical methods have rendered reliance on analytical solutions obsolete. Here, we present a technique for characterizing diffusion-driven anomalous transport (i.e., anomalous transport driven by a "trapping" process whose trapping and release statistics are independent of the groundwater flow velocity). Examples include diffusion into low permeability zones, kinetic sorption, and matrix diffusion. Using field observations, we simultaneously calibrate an exponential probability distribution for time spent on a single sojourn in the mobile domain and a truncated power law probability distribution for time spent on a single sojourn in the immobile domain via a stochastic global optimization technique. The calibrated distributions, being independent of the flow regime, are applicable to the same domain under any flow conditions, including linear flow. In the context of the continuous time random walk (CTRW), one may simply define a transition to represent a single trap-and-release cycle, and directly compute the spatiotemporal transition distribution that defines the CTRW from the two calibrated distributions and the local seepage velocity (so that existing CTRW transport theory applies). A test of our methodology against a push-pull test from the MADE site demonstrated fitting performance comparable to that of a 3-D MODFLOW/MT3DMS model with a variety of hydraulic conductivity zones and explicit treatment of mobile-immobile mass
Transport of contaminants in the planetary boundary layer
NASA Technical Reports Server (NTRS)
Lee, I. Y.; Swan, P. R.
1978-01-01
A planetary boundary layer model is described and used to simulate PBL phenomena including cloud formation and pollution transport in the San Francisco Bay Area. The effect of events in the PBL on air pollution is considered, and governing equations for the average momentum, potential temperature, water vapor mixing ratio, and air contaminants are presented. These equations are derived by integrating the basic equations vertically through the mixed layer. Characteristics of the day selected for simulation are reported, and the results suggest that the diurnally cyclic features of the mesoscale motion, including clouds and air pollution, can be simulated in a readily interpretable way with the model.
Karniadakis, George Em
2014-03-11
The main objective of this project is to develop new computational tools for uncertainty quantifica- tion (UQ) of systems governed by stochastic partial differential equations (SPDEs) with applications to advection-diffusion-reaction systems. We pursue two complementary approaches: (1) generalized polynomial chaos and its extensions and (2) a new theory on deriving PDF equations for systems subject to color noise. The focus of the current work is on high-dimensional systems involving tens or hundreds of uncertain parameters.
Subsurface Flow and Contaminant Transport Documentation and User's Guide
Aleman, S.E.
1999-07-28
This report documents a finite element code designed to model subsurface flow and contaminant transport, named FACT. FACT is a transient three-dimensional, finite element code designed to simulate isothermal groundwater flow, moisture movement, and solute transport in variably saturated and fully saturated subsurface porous media. The code is designed specifically to handle complex multi-layer and/or heterogeneous aquifer systems in an efficient manner and accommodates a wide range of boundary conditions. Additionally, 1-D and 2-D (in Cartesian coordinates) problems are handled in FACT by simply limiting the number of elements in a particular direction(s) to one. The governing equations in FACT are formulated only in Cartesian coordinates.
Analysis of Contaminant Transport through the Vadose and Saturated Zones for Source Screening
NASA Astrophysics Data System (ADS)
Bedekar, V.; Neville, C. J.; Tonkin, M. J.
2010-12-01
At complex sites there may be many potential source areas. Screening level analyses are useful to identify which of the source areas should be the focus of detailed investigation and analysis. A screening tool has been developed to evaluate the threat posed by waste sites on groundwater quality. This tool implements analytical solutions to simulate contaminant transport through the vadose and saturated zones and predict time-varying concentrations at potential groundwater receptors. The screening tool is developed within a user friendly, Microsoft ExcelTM based interface; however, care has been taken to implement rigorous solutions. The screening tool considers the following mechanisms: (a) Partitioning of soil contamination in to an equivalent dissolved concentration. For a time-invariant source, the solution is generalized from [3] for sorption and decay. For a time-varying source, the solution represents a special, degenerate, case of a solution implemented in ATRANS [2]; (b) One-dimensional (1D) transport of the dissolved contamination through the vadose zone considering 1D dispersion, equilibrium sorption, and first order transformation reactions. Steady state infiltration and moisture content are assumed; (c) Blending (mixing) of ambient water quality in the saturated zone with the contaminated water leaching from the vadose zone; and (d) Three-dimensional (3D) transport through the saturated zone using the formulation provided in [2], considering advection, dispersion, sorption, and first-order transformation reactions. The solution is derived using integral transform methods, following approaches adopted in [1] and [4]. Independent verification showed that the analytical techniques implemented in this study generate solutions that closely approximate those obtained using sophisticated numerical approaches, with a systematic over-estimate of the likely impact to groundwater that (predictably) stems from the use of a 1D approximation in the vadose zone. As a
Antidiffusive velocities for multipass donor cell advection
Margolin, L.; Smolarkiewicz, P.K.
1999-01-01
Multidimensional positive definite advection transport algorithm (MPDATA) is an iterative process for approximating the advection equation, which uses a donor cell approximation to compensate for the truncation error of the originally specified donor cell scheme. This step may be repeated an arbitrary number of times, leading to successfully more accurate solutions to the advection equation. In this paper, the authors show how to sum the successive approximations analytically to find a single antidiffusive velocity that represents the effects of an arbitrary number of passes. The analysis is first done in one dimension to illustrate the method and then is repeated in two dimensions. The existence of cross terms in the truncation analysis of the two-dimensional equations introduces an extra complication into the calculation. The authors discuss the implementation of the antidiffusive velocities and provide some examples of applications, including a third-order accurate scheme.
CONTAMINANT TRANSPORT IN SEDIMENT UNDER THE INFLUENCE OF SUBMARINE GROUNDWATER DISCHARGE
Theoretical estimations and laboratory studies suggest that capping can effectively retard contaminant transport from sediments under undisturbed conditions. However, contaminated near-shore areas, commonly selected as capping sites, are frequently subjected to submarine groundwa...
NASA Astrophysics Data System (ADS)
Göransson, G.; Larson, M.; Bendz, D.; Åkesson, M.
2011-12-01
Landslides of contaminated soil into surface water represent an overlooked exposure pathway that has not been addressed properly in existing risk analysis for landslide hazard, contaminated land, or river basin management. A landslide of contaminated soil into surface water implies an instantaneous exposure of the water to the contaminated soil, dramatically changing the prerequisites for the mobilisation and transport of pollutants. In this study, an analytical approach is taken to simulate the transport of suspended matter released in connection with landslides into rivers. Different analytical solutions to the advection-dispersion equation (ADE) were tested against the measured data from the shallow rotational, retrogressive landslide in clayey sediments that took place in 1993 on the Göta River, SW Sweden. The landslide encompassed three distinct events, namely an initial submerged slide, followed by a main slide, and a retrogressive slide. These slides generated three distinct and non-Gaussian peaks in the online turbidity recordings at the freshwater intake downstream the slide area. To our knowledge, this registration of the impact in a river of the sediment release from a landslide is one of the few of its kind in the world, and unique for Sweden considering the low frequency of landslide events, making it highly useful for evaluating how appropriate the ADE is to describe a landslide into surface water. The results yielded realistic predictions of the measured concentration variation, after proper calibration. For the three individual slides it was estimated that a total of about 0.6% (515 000 kg) of the total landslide mass went into suspension/was suspended and was transported downstream. This release corresponds to about 1 to 2% of the annual suspended sediment delivery for that river stretch. The studied landslide partly involved an industrial area and by applying the analytical solution for the transport of metals in the sediments it was found that
NASA Astrophysics Data System (ADS)
Göransson, G.; Larson, M.; Bendz, D.; Åkesson, M.
2012-07-01
Landslides of contaminated soil into surface water represent an overlooked exposure pathway that has not been addressed properly in existing risk analysis for landslide hazard, contaminated land, or river basin management. A landslide of contaminated soil into surface water implies an instantaneous exposure of the water to the soil, dramatically changing the prerequisites for the mobilisation and transport of pollutants. In this study, an analytical approach is taken to simulate the transport of suspended matter released in connection with landslides into rivers. Different analytical solutions to the advection-dispersion equation (ADE) were tested against the measured data from the shallow rotational, retrogressive landslide in clayey sediments that took place in 1993 on the Göta River, SW Sweden. The landslide encompassed three distinct events, namely an initial submerged slide, followed by a main slide, and a retrogressive slide. These slides generated three distinct and non-Gaussian peaks in the online turbidity recordings at the freshwater intake downstream the slide area. To our knowledge, this registration of the impact on a river of the sediment release from a landslide is one of few of its kind in the world and unique for Sweden. Considering the low frequency of such events, the data from this landslide are highly useful for evaluating how appropriate the ADE is to describe the effects of landslides into surface water. The results yielded realistic predictions of the measured variation in suspended particle matter (SPM) concentration, after proper calibration. For the three individual slides it was estimated that a total of about 0.6% of the total landslide mass went into suspension and was transported downstream. This release corresponds to about 1 to 2% of the annual suspended sediment transport for that river stretch. The studied landslide partly involved an industrial area, and by applying the analytical solution to estimate the transport of metals in
Xie, Haijian; Chen, Yunmin; Thomas, Hywel R; Sedighi, Majid; Masum, Shakil A; Ran, Qihua
2016-02-01
A field investigation of contaminant transport beneath and around an uncontrolled landfill site in Huainan in China is presented in this paper. The research aimed at studying the migration of some chemicals present in the landfill leachate into the surrounding clayey soils after 17 years of landfill operation. The concentrations of chloride and sodium ions in the pore water of soil samples collected at depths up to 15 m were obtained through an extensive site investigation. The contents of organic matter in the soil samples were also determined. A two-dimensional numerical study of the reactive transport of sodium and chloride ion in the soil strata beneath and outside the landfill is also presented. The numerical modelling approach adopted is based on finite element/finite difference techniques. The domain size of approximately 300 × 30 m has been analysed and major chemical transport parameters/mechanisms are established via a series of calibration exercises. Numerical simulations were then performed to predict the long-term behaviour of the landfill in relation to the chemicals studied. The lateral migration distance of the chloride ions was more than 40 m which indicates that the advection and mechanical dispersion are the dominant mechanism controlling the contaminant transport at this site. The results obtained from the analysis of chloride and sodium migration also indicated a non-uniform advective flow regime of ions with depth, which were localised in the first few metres of the soil beneath the disposal site. The results of long-term simulations of contaminant transport indicated that the concentrations of ions can be 10 to 30 times larger than that related to the allowable limit of concentration values. The results of this study may be of application and interest in the assessment of potential groundwater and soil contamination at this site with a late Pleistocene clayey soil. The obtained transport properties of the soils and the contaminant transport
A Linear Systems Approach to Segmented Watershed Contaminant Transport
NASA Astrophysics Data System (ADS)
Carleton, J. N.
2013-12-01
The U.S. Environmental Protection Agency (USEPA) employs simulation models to estimate concentrations of pesticide residues in surface waters for risk assessment. These models have historically been used to simulate runoff loadings from homogeneous landscapes to isolated, well-mixed lentic systems that generically represent vulnerable waters. Recent efforts to refine this approach in terms of realism and geographic specificity have focused on enhancing the level of detail of the landscape representation, rather than that of receiving water hydrology. Linear systems theory and transfer function based approaches have been applied by various investigators to the representation of contaminant leaching through soils, and to surface water hydrology (e.g., unit hydrographs), but rarely to contaminant transport either within surface waters, or through multi-compartment systems such as stream networks. This poster describes a straightforward approach to simulating watersheds as segmented into collections of linked water bodies. The approach employs convolution integrals, impulse response functions, and the Discrete Fourier Transform to propagate concentration time series from upstream to downstream locations. Given knowledge only of estimated mean stream residence times, with appropriately-scaled segmentations of catchments, realistic representations of concentration dynamics are shown to be achievable. These representations are based upon high-frequency atrazine monitoring data sets collected over common time periods from upstream and downstream locations within the same small watersheds. Simulated concentrations are shown to match measured concentrations well in both the temporal and spectral domains without the need for calibration, and despite inherent simplifying assumptions such as steady flow. The approach may have utility for enhancing surface water hydrologic representation in contaminant modeling used for regulatory purposes.
Coliform contamination of a coastal embayment: Sources and transport pathways
Weiskel, P.K.; Howes, B.L.; Heufelder, G.R.
1996-01-01
Fecal bacterial contamination of nearshore waters has direct economic impacts to coastal communities through the loss of shellfisheries and restrictions of recreational uses. We conducted seasonal measurements of fecal coliform (FC) sources and transport pathways contributing to FC contamination of Buttermilk Bay, a shallow embayment adjacent to Buzzards Bay, MA. Typical of most coastal embayments, there were no direct sewage discharges (i.e., outfalls), and fecal bacteria from human, domestic animal, and wildlife pools entered open waters primarily through direct deposition or after transport through surface waters or groundwaters. Direct fecal coliform inputs to bay waters occurred primarily in winter (December-March) from waterfowl, ~33 x 1012 FC yr-1 or ~67% of the total annual loading. Effects of waterfowl inputs on bay FC densities were mitigated by their seasonality, wide distribution across the bay surface, and the apparent limited dispersal from fecal pellets. On-site disposal of sewage by septic systems was the single largest FC source in the watershed-embayment system, 460 x 1012 FC yr-1, but due to attenuation during subsurface transport only a minute fraction, < 0.006 x 1012 FC yr-1, reached bay waters (<0.01% of annual input to bay). Instead, surface water flows, via storm drains and natural streams under both wet- and dry-weather conditions, contributed the major terrestrial input, 12 x 1012 FC yr-1 (24% of annual input), all from animal sources. Since most of the surface water FC inputs were associated with periodic, short-duration rain events with discharge concentrated in nearshore zones, wet-weather flows were found to have a disproportionately high impact on nearshore FC levels. Elution of FC from shoreline deposits of decaying vegetation (wrack) comprised an additional coliform source. Both laboratory and field experiments suggest significant elution of bacteria from wrack, ~3 x 1012 FC yr-1 on a bay-wide basis (6% of annual input), primarily
LES validation for contaminant transport in urban areas
NASA Astrophysics Data System (ADS)
Hertwig, D.; Leitl, B.; Schatzmann, M.; Patnaik, G.
2010-09-01
Contaminant transport in urban areas poses a major challenge with respect to its simulation with computational fluid dynamics (CFD) models. The use of time-resolved approaches like large-eddy simulation (LES) can provide insight into transient flow and dispersion regimes, which are strongly influenced by the urban geometry. LES models have the potential to resolve the characteristic unsteady flow features and their impact on plume dynamics, whereas standard industrial codes based on Reynolds-averaged Navier-Stokes (RANS) equations can only yield steady state solutions. However, the potential to simulate the energetically dominating part of an inherently unsteady turbulent flow with LES also sets higher requirements for validation strategies. This includes that the evaluation of the model performance must go beyond comparisons of first and second order statistics which were adequate for RANS models and currently provide the basis for most of the validation metrics used as a standard. With regard to an a posteriori validation of model results for atmospheric boundary layer (ABL) flow and dispersion in complex geometry, laboratory data from boundary-layer wind tunnels are of special value. Since inflow and boundary conditions are well-defined, systematic laboratory studies provide high statistical confidence levels of measured quantities. The potential of field measurements - in this regard - is limited due to the natural atmospheric variability. In order to verify the realistic simulation of the spatio-temporal behavior of turbulent eddies, transient flow phenomena have to be characterized in experimental validation data sets as well. This topic is closely linked to structure identification and the characterization of organized motions in ABL flows, for which advanced analysis strategies like wavelet transforms, orthogonal decomposition, or stochastic estimation can be employed. Systematic comparisons of wind-tunnel measurements and LES simulation results are planned
Monitoring and modeling contaminated sediment transport in the White Oak Creek watershed
Fontaine, T.A.
1991-11-01
Over the past 47 years, operations and waste disposal activities at Oak Ridge National Laboratory have resulted in the contamination of the White Oak Creek drainage system. The containments presenting the highest risk to human health and the environment are particle reactive and are associated with the soils and sediments in White Oak Creek. During floods, the erosion of these sediments results in the transport of contaminants out of the catchment into the Clinch River. A long-term strategy is required to monitor the movement of contaminated sediments and to predict the transport of these sediments that could occur during major floods. A monitoring program will provide the information required to (1) evaluate the existing off-site transport of contaminated sediments, (2) evaluate the need for short-term control measures, (3) set priorities for remediation of contaminated areas in White Oak Creek (4) verify the success of completed remedial actions intended to control the movement of contaminated sediments, and (5) develop a computer model to simulate the transport of contaminated sediments in White Oak Creek. A contaminant-transport model will be developed to (1) evaluate the potential for the off-site transport of contaminated sediments during major floods, (2) develop long term control measures and remediation solutions, (3) predict the impact of future land-use changes in White Oak Creek on the transport of contaminated sediment. This report contains a plan for the monitoring and modeling activities required to accomplish these objectives.
Fontaine, T.A.
1991-11-01
Over the past 47 years, operations and waste disposal activities at Oak Ridge National Laboratory have resulted in the contamination of the White Oak Creek drainage system. The containments presenting the highest risk to human health and the environment are particle reactive and are associated with the soils and sediments in White Oak Creek. During floods, the erosion of these sediments results in the transport of contaminants out of the catchment into the Clinch River. A long-term strategy is required to monitor the movement of contaminated sediments and to predict the transport of these sediments that could occur during major floods. A monitoring program will provide the information required to (1) evaluate the existing off-site transport of contaminated sediments, (2) evaluate the need for short-term control measures, (3) set priorities for remediation of contaminated areas in White Oak Creek (4) verify the success of completed remedial actions intended to control the movement of contaminated sediments, and (5) develop a computer model to simulate the transport of contaminated sediments in White Oak Creek. A contaminant-transport model will be developed to (1) evaluate the potential for the off-site transport of contaminated sediments during major floods, (2) develop long term control measures and remediation solutions, (3) predict the impact of future land-use changes in White Oak Creek on the transport of contaminated sediment. This report contains a plan for the monitoring and modeling activities required to accomplish these objectives.
Concentration through large advection
NASA Astrophysics Data System (ADS)
Aleja, D.; López-Gómez, J.
2014-11-01
In this paper we extend the elegant results of Chen, Lam and Lou [6, Section 2], where a concentration phenomenon was established as the advection blows up, to a general class of adventive-diffusive generalized logistic equations of degenerate type. Our improvements are really sharp as we allow the carrying capacity of the species to vanish in some subdomain with non-empty interior. The main technical devices used in the derivation of the concentration phenomenon are Proposition 3.2 of Cano-Casanova and López-Gómez [5], Theorem 2.4 of Amann and López-Gómez [1] and the classical Harnack inequality. By the relevance of these results in spatial ecology, complete technical details seem imperative, because the proof of Theorem 2.2 of [6] contains some gaps originated by an “optimistic” use of Proposition 3.2 of [5]. Some of the general assumptions of [6] are substantially relaxed.
Romero-González, J; Walton, J C; Peralta-Videa, J R; Rodríguez, E; Romero, J; Gardea-Torresdey, J L
2009-01-15
The biosorption of Cr(III) onto packed columns of Agave lechuguilla was analyzed using an advective-dispersive (AD) model and its analytical solution. Characteristic parameters such as axial dispersion coefficients, retardation factors, and distribution coefficients were predicted as functions of inlet ion metal concentration, time, flow rate, bed density, cross-sectional column area, and bed length. The root-mean-square-error (RMSE) values 0.122, 0.232, and 0.285 corresponding to the flow rates of 1, 2, and 3 (10(-3))dm3min(-1), respectively, indicated that the AD model provides an excellent approximation of the simulation of lumped breakthrough curves for the adsorption of Cr(III) by lechuguilla biomass. Therefore, the model can be used for design purposes to predict the effect of varying operational conditions. PMID:18462882
Zhong, Lirong; Szecsody, James E.; Truex, Michael J.; Williams, Mark D.; Liu, Yuanyuan
2015-05-01
Use of gas-phase amendments for in situ remediation of inorganic contaminants in unsaturated sediments of the vadose zone may be advantageous, but there has been limited development and testing of gas remediation technologies. Treatment with ammonia gas has been studied and has a potential for use in treating inorganic contaminants such as uranium because it induces a high pore-water pH causing mineral dissolution and subsequent formation of stable precipitates that decrease the mobility of some contaminants. For field application, knowledge of ammonia transport and the geochemical reactions induced by ammonia is needed. Laboratory studies were conducted to support calculations needed for field treatment design, to quantify advective and diffusive ammonia transport in unsaturated sediments, to evaluate reactions among gas, sediment, and water, and to study reaction-induced pore-water chemistry changes as a function of ammonia delivery conditions. Ammonia gas quickly partitions into sediment pore water and increases pH up to 13.2. Injected ammonia gas front movement can be reasonably predicted by gas flow rate and equilibrium partitioning. The ammonia gas diffusion rate is a function of the water content in the sediment. Measured diffusion front movement was 0.05, 0.03, and 0.02 cm/hr. in sediments with 2.0%, 8.7%, and 13.0% water content, respectively. Sodium, aluminum, and silica pore-water concentrations increase on exposure to ammonia and then decline as aluminosilicates precipitate with declining pH. When uranium is present in the sediment and pore water, up to 85% of the water-leachable uranium was immobilized by ammonia treatment.
Hydrophobic organic contaminant transport property heterogeneity in the Borden Aquifer
NASA Astrophysics Data System (ADS)
Allen-King, Richelle M.; Kalinovich, Indra; Dominic, David F.; Wang, Guohui; Polmanteer, Reid; Divine, Dana
2015-03-01
We determined that the spatial heterogeneity in aquifer properties governing the reactive transport of volatile organic contaminants is defined by the arrangement of lithofacies. We measured permeability (k) and perchloroethene sorption distribution coefficient (Kd) for lithofacies that we delineated for samples from the Canadian Forces Base Borden Aquifer. We compiled existing data and collected 57 new cores to characterize a 30 m section of the aquifer near the test location of Mackay et al. (1986). The k and Kd were measured for samples taken at six elevations from all cores to create a data set consisting of nearly 400 colocated measurements. Through analysis of variance (corrected for multiple comparisons), we determined that the 12 originally mapped lithofacies could be grouped into five relatively distinct chemohydrofacies that capture the variability of both transport properties. The mean of ln k by lithofacies was related to the grain size and the variance was relatively consistent. In contrast, both the mean and variance of ln Kd were greater for more poorly sorted lithofacies, which were also typically more coarse-grained. Half of the aquifer sorption capacity occurred in the three highest-sorbing lithofacies but comprised only 20% of its volume. The model of the aquifer that emerged is that of discontinuous scour-fill deposits of medium sand, generally characterized by greater Kd and k, within laterally extensive fine-grained to very fine-grained sands of lower Kd and k. Our findings demonstrate the importance of considering source rock composition, transport, and deposition processes when constructing conceptual models of chemohydrofacies.
Pore Water PAH Transport in Amended Sediment Caps
NASA Astrophysics Data System (ADS)
Gidley, P. T.; Kwon, S.; Ghosh, U.
2009-05-01
Capping is a common remediation strategy for contaminated sediments that creates a physical barrier between contaminated sediments and the water column. Diffusive flux of contaminants through a sediment cap is small. However, under certain hydrodynamic conditions such as groundwater potential and tidal pumping, groundwater advection can accelerate contaminant transport. Hydrophobic organic contaminants such as polycyclic aromatic hydrocarbons (PAHs) could be transported through the cap under advective conditions. To better understand PAH migration under these conditions, physical models of sediment caps were evaluated in the laboratory through direct measurement of pore water using solid phase micro-extraction with gas chromatography and mass spectrometry. Contaminated sediment and capping material was obtained from an existing Superfund site that was capped at Eagle Harbor, Washington. A PAH dissolution model linked to an advection-dispersion equation with retardation using published organic carbon-water partitioning coefficients (Koc) was compared to measured PAHs in the sediment and cap porewater of the physical model.
EVALUATION OF THE STATE-OF-THE-ART CONTAMINATED SEDIMENT TRANSPORT AND FATE MODELING SYSTEM
Modeling approaches for evaluating the transport and fate of sediment and associated contaminants are briefly reviewed. The main emphasis is on: 1) the application of EFDC (Environmental Fluid Dynamics Code), the state-of-the-art contaminated sediment transport and fate public do...
Stock, Naomi L; Furdui, Vasile I; Muir, Derek C G; Mabury, Scott A
2007-05-15
Perfluorosulfonates (PFSAs) and perfluorocarboxylates (PFCAs) have been hypothesized to reach remote locations such as the Canadian Arctic either indirectly as volatile precursor chemicals that undergo atmospheric transport and subsequent degradation, or directly via oceanic and atmospheric transport of the PFSAs and PFCAs themselves. Water, sediment, and air samples were collected from three Arctic lakes (Amituk, Char, and Resolute) on Cornwallis Island, Nunavut, Canada. Samples were analyzed for PFSAs and PFCAs, precursor chemicals including the fluorotelomer alcohols (FTOHs) and polyfluorinated sulfonamides (FSAs), and precursor degradation products such as the fluorotelomer unsaturated carboxylates (FTUCAs). PFSAs and PFCAs were detected in water and sediment of all three Arctic lakes (concentrations ranged from nondetect to 69 ng/L and nondetect to 85 ng/g dry weight, respectively). FTOHs and FSAs were observed in air samples (mean concentrations ranged from 2.8 to 29 pg/m3), and confirm that volatile precursors are reaching Arctic latitudes. The observation of degradation products, including FTUCAs observed in sediment and atmospheric particles, and N-ethyl perfluorooctanesulfonamide (NEtFOSA) and perfluorooctanesulfonamide (PFOSA) in air samples, indicate that degradation of the FTOHs and FSAs is occurring in the Arctic environment. PFSAs and PFCAs were also observed on atmospheric particles (mean concentrations ranged from < 0.1 to 5.9 pg/m3). In addition, results of this study also indicate that local perfluoroalkyl contamination of Resolute Lake, which is located downstream of an airport wastewater input, has occurred. PMID:17547174
Enhanced subsurface characterization for prediction of contaminant transport using co-kriging
Benson, C.H.; Rashad, S.M.
1996-12-31
Groundwater flow and advective transport were simulated in a heterogeneous synthetic aquifer. These simulations were conducted when the aquifer was fully defined and when it was characterized using a limited amount of hard and soft data (hydraulic conductivity data and soil classifications). Co-kriging was used to combine the data types when estimating the hydraulic conductivity field throughout the aquifer. Results of the flow and transport simulations showed that soil classifications were useful in characterizing the hydraulic conductivity field and reducing errors in statistics describing the plume.
LAYER DEPENDENT ADVECTION IN CMAQ
The advection methods used in CMAQ require that the Courant-Friedrichs-Lewy (CFL) condition be satisfied for numerical stability and accuracy. In CMAQ prior to version 4.3, the ADVSTEP algorithm established CFL-safe synchronization and advection timesteps that were uniform throu...
Contaminant transport model validation: The Oak Ridge Reservation
Lee, R.R.; Ketelle, R.H.
1988-09-01
In the complex geologic setting on the Oak Ridge Reservation, hydraulic conductivity is anisotropic and flow is strongly influenced by an extensive and largely discontinuous fracture network. Difficulties in describing and modeling the aquifer system prompted a study to obtain aquifer property data to be used in a groundwater flow model validation experiment. Characterization studies included the performance of an extensive suite of aquifer test within a 600-square-meter area to obtain aquifer property values to describe the flow field in detail. Following aquifer test, a groundwater tracer test was performed under ambient conditions to verify the aquifer analysis. Tracer migration data in the near-field were used in model calibration to predict tracer arrival time and concentration in the far-field. Despite the extensive aquifer testing, initial modeling inaccurately predicted tracer migration direction. Initial tracer migration rates were consistent with those predicted by the model; however, changing environmental conditions resulted in an unanticipated decay in tracer movement. Evaluation of the predictive accuracy of groundwater flow and contaminant transport models on the Oak Ridge Reservation depends on defining the resolution required, followed by field testing and model grid definition at compatible scales. The use of tracer tests, both as a characterization method and to verify model results, provides the highest level of resolution of groundwater flow characteristics. 3 refs., 4 figs.
Modeling uranium transport in acidic contaminated groundwater with base addition
Zhang, Fan; Luo, Wensui; Parker, Jack C.; Brooks, Scott C; Watson, David B; Jardine, Philip; Gu, Baohua
2011-01-01
This study investigates reactive transport modeling in a column of uranium(VI)-contaminated sediments with base additions in the circulating influent. The groundwater and sediment exhibit oxic conditions with low pH, high concentrations of NO{sub 3}{sup -}, SO{sub 4}{sup 2-}, U and various metal cations. Preliminary batch experiments indicate that additions of strong base induce rapid immobilization of U for this material. In the column experiment that is the focus of the present study, effluent groundwater was titrated with NaOH solution in an inflow reservoir before reinjection to gradually increase the solution pH in the column. An equilibrium hydrolysis, precipitation and ion exchange reaction model developed through simulation of the preliminary batch titration experiments predicted faster reduction of aqueous Al than observed in the column experiment. The model was therefore modified to consider reaction kinetics for the precipitation and dissolution processes which are the major mechanism for Al immobilization. The combined kinetic and equilibrium reaction model adequately described variations in pH, aqueous concentrations of metal cations (Al, Ca, Mg, Sr, Mn, Ni, Co), sulfate and U(VI). The experimental and modeling results indicate that U(VI) can be effectively sequestered with controlled base addition due to sorption by slowly precipitated Al with pH-dependent surface charge. The model may prove useful to predict field-scale U(VI) sequestration and remediation effectiveness.
CONTAMINATION OF U. S. ARCTIC ECOSYSTEMS BY LONG-RANGE TRANSPORT OF ATMOSPHERIC CONTAMINANTS
Various kinds of atmospheric pollutants are found in Arctic environments, including organic contaminants, radionuclides, and pollutants associated with fossil fuel combustion, smelting, and industrial development. hile some of these contaminants originate in the Arctic itself, ot...
Chen, Kuo-Fu
1996-11-01
The health risks for an individual exposed to contaminants released from SRS outfalls from 1989 to 1995 were estimated. The exposure pathways studied are ingestion of drinking water, ingestion of contaminated fish and dermal contact with contaminants in water while swimming. The estimated incremental risks for an individual developing cancer vary from 3.E-06 to 1.0E-05. The estimated total exposure chronic noncancer hazard indices vary from 6.E-02 to 1.E-01. The critical contaminants were ranked based on their cancer risks and chronic noncarcinogenic hazard quotients. For cancer risks, the critical contaminants released from SRS outfalls are arsenic, tetrachloroethylene, and benzene. For chronic noncarcinogenic risks, the critical contaminants released from srs outfalls are cadmium, arsenic, silver, chromium, mercury, selenium, nitrate, manganese, zinc, nickel, uranium, barium, copper, tetrachloroethylene, cyanide, and phenol. The critical pathways in decreasing risk order are ingestion of contaminated fish, ingestion of drinking water and dermal contact with contaminants in water while swimming.
Advection around ventilated U-shaped burrows: A model study
NASA Astrophysics Data System (ADS)
Brand, Andreas; Lewandowski, JöRg; Hamann, Enrico; Nützmann, Gunnar
2013-05-01
Advective transport in the porous matrix of sediments surrounding burrows formed by fauna such as Chironomus plumosus has been generally neglected. A positron emission tomography study recently revealed that the pumping activity of the midge larvae can indeed induce fluid flow in the sediment. We present a numerical model study which explores the conditions at which advective transport in the sediment becomes relevant. A 0.15 m deep U-shaped burrow with a diameter of 0.002 m within the sediment was represented in a 3-D domain. Fluid flow in the burrow was calculated using the Navier-Stokes equation for incompressible laminar flow in the burrow, and flow in the sediment was described by Darcy's law. Nonreactive and reactive transport scenarios were simulated considering diffusion and advection. The pumping activity of the model larva results in considerable advective flow in the sediment at reasonable high permeabilities with flow velocities of up to 7.0 × 10-6 m s-1 close to the larva for a permeability of 3 × 10-12 m2. At permeabilities below 7 × 10-13 m2 advection is negligible compared to diffusion. Reactive transport simulations using first-order kinetics for oxygen revealed that advective flux into the sediment downstream of the pumping larva enhances sedimentary uptake, while the advective flux into the burrow upstream of the larvae inhibits diffusive sedimentary uptake. Despite the fact that both effects cancel each other with respect to total solute uptake, the advection-induced asymmetry in concentration distribution can lead to a heterogeneous solute and redox distribution in the sediment relevant to complex reaction networks.
Advection, diffusion and delivery over a network
Heaton, Luke L.M.; López, Eduardo; Maini, Philip K.; Fricker, Mark D.; Jones, Nick S.
2014-01-01
Many biological, geophysical and technological systems involve the transport of resource over a network. In this paper we present an algorithm for calculating the exact concentration of resource at any point in space or time, given that the resource in the network is lost or delivered out of the network at a given rate, while being subject to advection and diffusion. We consider the implications of advection, diffusion and delivery for simple models of glucose delivery through a vascular network, and conclude that in certain circumstances, increasing the volume of blood and the number of glucose transporters can actually decrease the total rate of glucose delivery. We also consider the case of empirically determined fungal networks, and analyze the distribution of resource that emerges as such networks grow over time. Fungal growth involves the expansion of fluid filled vessels, which necessarily involves the movement of fluid. In three empirically determined fungal networks we found that the minimum currents consistent with the observed growth would effectively transport resource throughout the network over the time-scale of growth. This suggests that in foraging fungi, the active transport mechanisms observed in the growing tips may not be required for long range transport. PMID:23005783
Chemical contaminants in the Wadden Sea: Sources, transport, fate and effects
NASA Astrophysics Data System (ADS)
Laane, R. W. P. M.; Vethaak, A. D.; Gandrass, J.; Vorkamp, K.; Köhler, A.; Larsen, M. M.; Strand, J.
2013-09-01
The Wadden Sea receives contaminants from various sources and via various transport routes. The contaminants described in this overview are various metals (Cd, Cu, Hg, Pb and Zn) and various organic contaminants (polychlorinated biphenyls (PCBs), polycyclic aromatic hydrocarbons (PAHs) and lindane (hexachlorocyclohexane, γ-HCH)). In addition, information is presented about other and emerging contaminants such as antifouling biocides (e.g. TBT and Irgarol), brominated flame retardants (BFRs), poly- and perfluorinated compounds (PFCs) and pharmaceutical and personal care products (PPCPs). Special attention is given to biogeochemical processes that contribute to the mobilization of contaminants in the surface sediments of the Wadden Sea. Finally, the effects on organisms of contaminants are reviewed and discussed. The main source of contaminants in the Wadden Sea are the rivers Rhine (via de Dutch coastal zone), Elbe and Weser. The Wadden Sea is not a sink for contaminants and adsorbed contaminants are transported from east to west. The surface sediments of the Wadden Sea are an important source for contaminants to the water above. The input and concentration of most contaminants have significantly decreased in water, sediments, organisms (e.g., mussel, flounder and bird eggs) in various parts of the Wadden Sea in the last three decades. Remarkably, the Cd concentration in mussels is increasing the last decades. In recent decades, the effects of contaminants on organisms (e.g., flounder, seal) have fallen markedly. Most of the affected populations have recovered, except for TBT induced effects in snails. Little is known about the concentration and effects of most emerging contaminants and the complex environmental mixtures of contaminants. It is recommended to install an international coordinated monitoring programme for contaminants and their effects in the whole Wadden Sea and to identify the chemical contaminants that really cause the effect.
Varank, Gamze; Demir, Ahmet; Yetilmezsoy, Kaan; Bilgili, M Sinan; Top, Selin; Sekman, Elif
2011-11-01
One-dimensional (1D) advection-dispersion transport modeling was conducted as a conceptual approach for the estimation of the transport parameters of fourteen different phenolic compounds (phenol, 2-CP, 2-MP, 3-MP, 4-MP, 2-NP, 4-NP, 2,4-DNP, 2,4-DCP, 2,6-DCP, 2,4,5-TCP, 2,4,6-TCP, 2,3,4,6-TeCP, PCP) and three different inorganic contaminants (Cu, Zn, Fe) migrating downward through the several liner systems. Four identical pilot-scale landfill reactors (0.25 m3) with different composite liners (R1: 0.10+0.10 m of compacted clay liner (CCL), L(e) = 0.20 m, k(e) = 1 × 10(-8) m/s, R2: 0.002-m-thick damaged high-density polyethylene (HDPE) geomembrane overlying 0.10+0.10 m of CCL, L(e) = 0.20 m, k(e) = 1 × 10(-8) m/s, R3: 0.002-m-thick damaged HDPE geomembrane overlying a 0.02-m-thick bentonite layer encapsulated between 0.10+0.10 m CCL, L(e) = 0.22 m, k(e) = 1 × 10(-8) m/s, R4: 0.002-m-thick damaged HDPE geomembrane overlying a 0.02-m-thick zeolite layer encapsulated between 0.10+0.10 m CCL, L(e) = 0.22 m, k(e) = 4.24 × 10(-7) m/s) were simultaneously run for a period of about 540 days to investigate the nature of diffusive and advective transport of the selected organic and inorganic contaminants. The results of 1D transport model showed that the highest molecular diffusion coefficients, ranging from 4.77×10(-10) to 10.67 × 10(-10)m2/s, were estimated for phenol (R4), 2-MP (R1), 2,4-DNP (R2), 2,4-DCP (R1), 2,6-DCP (R2), 2,4,5-TCP (R2) and 2,3,4,6-TeCP (R1). For all reactors, dispersion coefficients of Cu, ranging from 3.47 × 10(-6) m(2)/s to 5.37 × 10(-2) m2/s, was determined to be higher than others obtained for Zn and Fe. Average molecular diffusion coefficients of phenolic compounds were estimated to be about 5.64 × 10(-10) m2/s, 5.37 × 10(-10) m2/s, 2.69 × 10(-10) m2/s and 3.29 × 10(-10) m2/s for R1, R2, R3 and R4 systems, respectively. The findings of this study clearly indicated that about 35-50% of transport of phenolic compounds to the groundwater
The three-dimensional movement of a tracer plume containing bromide and chloride is investigated using the data base from a large-scale natural gradient field experiment on groundwater solute transport. The analysis focuses on the zeroth-, first-, and second-order spatial moments...
Technology Transfer Automated Retrieval System (TEKTRAN)
Improving strategies for monitoring subsurface contaminant transport includes performance comparison of competing models, developed independently or obtained via model abstraction. Model comparison and parameter discrimination involve specific performance indicators selected to better understand s...
Investigations at Pinal Creek, Arizona, evaluated routine sampling procedures for determination of aqueous inorganic geochemistry and assessment of contaminant transport by colloidal mobility. Sampling variables included pump type and flow rate, collection under air or nitrogen,...
Effects of lag and maximum growth in contaminant transport and biodegradation modeling
Wood, B.D.; Dawson, C.N.
1992-06-01
The effects of time lag and maximum microbial growth on biodegradation in contaminant transport are discussed. A mathematical model is formulated that accounts for these effects, and a numerical case study is presented that demonstrates how lag influences biodegradation.
Configuration Management Plan for Long Length Contaminated Equipment Receiver and Transport Trailers
DALE, R.N.
2000-05-01
Long Length Contaminated Equipment Removal System Receiver Trailer and Transport Trailer require a configuration management plan for design, requirements and operations baseline documents. This report serves as the plan for the Trailers.
Configuration Management Plan for Long Length Contaminated Equipment Receiver and Transport Trailers
DALE, R.N.
2000-12-18
Long Length Contaminated Equipment Removal System Receiver Trailers and Transport Trailers require identification and control for the design, requirements and operations baseline documents. This plan serves as those controls for the subject trailers.
Predicting subsurface contaminant transport and transformation requires mathematical models based on a variety of physical, chemical, and biological processes. The mathematical model is an attempt to quantitatively describe observed processes in order to permit systematic forecas...
FIELD-DRIVEN APPROACHES TO SUBSURFACE CONTAMINANT TRANSPORT MODELING.
Observations from field sites provide a means for prioritizing research activities. In the case of petroleum releases, observations may include spiking of concentration distributions that may be related to water table fluctuation, co-location of contaminant plumes with geochemi...
BIOPLUME MODEL FOR CONTAMINANT TRANSPORT AFFECTED BY OXYGEN LIMITED BIODEGRADATION
Many of the organic pollutants entering ground water are potentially biodegradable in the subsurface. This potential has been demonstrated in aquifers contaminated by wood-creosoting process wastes. The persistence of many of these organic compounds in the subsurface indicated ...
NASA Technical Reports Server (NTRS)
Kritz, Mark A.; Le Roulley, Jean-Claude; Danielsen, Edwin F.
1990-01-01
A series of upper tropospheric radon concentration measurements made over the eastern Pacific and west coast of the U.S. during the summers of 1983 and 1984 has revealed the occurrence of unexpectedly high radon concentrations for 9 of the 61 measurements. A frequency distribution plot of the set of 61 observations shows a distinct bimodal distribution, with approximately 2/5 of the observations falling close to 1 pCi/SCM, and 3/5 falling in a high concentration mode centered at about 11 pCi/SCM. Trajectory and synoptic analyses for two of the flights on which such high radon concentrations were observed indicate that this radon-rich air originated in the Asian boundary layer, ascended in cumulus updrafts, and was carried eastward in the fast moving air on the anticyclonic side of the upper tropospheric jet. The results suggest that the combination of rapid vertical transport from the surface boundary layer to the upper troposphere, followed by rapid horizontal transport eastward represents an efficient mode of long-transport for other, chemically reactive atmospheric trace constituents.
Bai, Chunmei; Li, Yusong
2014-08-01
Accurately predicting the transport of contaminants in the field is subject to multiple sources of uncertainty due to the variability of geological settings, the complexity of field measurements, and the scarcity of data. Such uncertainties can be amplified when modeling some emerging contaminants, such as engineered nanomaterials, when a fundamental understanding of their fate and transport is lacking. Typical field work includes collecting concentration at a certain location for an extended period of time, or measuring the movement of plume for an extended period time, which would result in a time series of observation data. This work presents an effort to evaluate the possibility of applying time series analysis, particularly, autoregressive integrated moving average (ARIMA) models, to forecast contaminant transport and distribution in the subsurface environment. ARIMA modeling was first assessed in terms of its capability to forecast tracer transport at two field sites, which had different levels of heterogeneity. After that, this study evaluated the applicability of ARIMA modeling to predict the transport of engineered nanomaterials at field sites, including field measured data of nanoscale zero valent iron and (nZVI) and numerically generated data for the transport of nano-fullerene aggregates (nC60). This proof-of-concept effort demonstrates the possibility of applying ARIMA to predict the contaminant transport in the subsurface environment. Like many other statistical models, ARIMA modeling is only descriptive and not explanatory. The limitation and the challenge associated with applying ARIMA modeling to contaminant transport in the subsurface are also discussed. PMID:24987973
NASA Astrophysics Data System (ADS)
Bai, Chunmei; Li, Yusong
2014-08-01
Accurately predicting the transport of contaminants in the field is subject to multiple sources of uncertainty due to the variability of geological settings, the complexity of field measurements, and the scarcity of data. Such uncertainties can be amplified when modeling some emerging contaminants, such as engineered nanomaterials, when a fundamental understanding of their fate and transport is lacking. Typical field work includes collecting concentration at a certain location for an extended period of time, or measuring the movement of plume for an extended period time, which would result in a time series of observation data. This work presents an effort to evaluate the possibility of applying time series analysis, particularly, autoregressive integrated moving average (ARIMA) models, to forecast contaminant transport and distribution in the subsurface environment. ARIMA modeling was first assessed in terms of its capability to forecast tracer transport at two field sites, which had different levels of heterogeneity. After that, this study evaluated the applicability of ARIMA modeling to predict the transport of engineered nanomaterials at field sites, including field measured data of nanoscale zero valent iron and (nZVI) and numerically generated data for the transport of nano-fullerene aggregates (nC60). This proof-of-concept effort demonstrates the possibility of applying ARIMA to predict the contaminant transport in the subsurface environment. Like many other statistical models, ARIMA modeling is only descriptive and not explanatory. The limitation and the challenge associated with applying ARIMA modeling to contaminant transport in the subsurface are also discussed.
Truex, Michael J.; Oostrom, Martinus; Tartakovsky, Guzel D.
2015-09-01
An approach was developed for evaluating vadose zone transport and attenuation of aqueous wastes containing inorganic (non-volatile) contaminants that were disposed of at the land surface (i.e., directly to the ground in cribs, trenches, tile fields, etc.) and their effect on the underlying groundwater. The approach provides a structured method for estimating transport of contaminants through the vadose zone and the resulting temporal profile of groundwater contaminant concentrations. The intent of the approach is also to provide a means for presenting and explaining the results of the transport analysis in the context of the site-specific waste disposal conditions and site properties, including heterogeneities and other complexities. The document includes considerations related to identifying appropriate monitoring to verify the estimated contaminant transport and associated predictions of groundwater contaminant concentrations. While primarily intended for evaluating contaminant transport under natural attenuation conditions, the approach can also be applied to identify types of, and targets for, mitigation approaches in the vadose zone that would reduce the temporal profile of contaminant concentrations in groundwater, if needed.
Characterization of Contaminant Transport Using Naturally-Occurring U-Series Disequilibria
TEH-LUNG KU
2001-06-01
Study the migration of nuclear waste contaminants in subsurface fractured systems using naturally occurring uranium and thorium-series radionuclides as tracers under in-situ physico-chemical and hydrogeologic conditions. Radioactive disequilibria among members of these decay-series nuclides can provide information on the rates of adsorption-desorption and transport of contaminants as well as on fluid transport and rock dissolution in a natural setting.
NASA Astrophysics Data System (ADS)
Zamani, K.; Bombardelli, F.
2011-12-01
Almost all natural phenomena on Earth are highly nonlinear. Even simplifications to the equations describing nature usually end up being nonlinear partial differential equations. Transport (ADR) equation is a pivotal equation in atmospheric sciences and water quality. This nonlinear equation needs to be solved numerically for practical purposes so academicians and engineers thoroughly rely on the assistance of numerical codes. Thus, numerical codes require verification before they are utilized for multiple applications in science and engineering. Model verification is a mathematical procedure whereby a numerical code is checked to assure the governing equation is properly solved as it is described in the design document. CFD verification is not a straightforward and well-defined course. Only a complete test suite can uncover all the limitations and bugs. Results are needed to be assessed to make a distinction between bug-induced-defect and innate limitation of a numerical scheme. As Roache (2009) said, numerical verification is a state-of-the-art procedure. Sometimes novel tricks work out. This study conveys the synopsis of the experiences we gained during a comprehensive verification process which was done for a transport solver. A test suite was designed including unit tests and algorithmic tests. Tests were layered in complexity in several dimensions from simple to complex. Acceptance criteria defined for the desirable capabilities of the transport code such as order of accuracy, mass conservation, handling stiff source term, spurious oscillation, and initial shape preservation. At the begining, mesh convergence study which is the main craft of the verification is performed. To that end, analytical solution of ADR equation gathered. Also a new solution was derived. In the more general cases, lack of analytical solution could be overcome through Richardson Extrapolation and Manufactured Solution. Then, two bugs which were concealed during the mesh convergence
Advection and diffusion in shoreline change prediction
NASA Astrophysics Data System (ADS)
Anderson, T. R.; Frazer, L. N.
2010-12-01
We added longshore advection and diffusion to the simple cross-shore rate calculation method, as used widely by the USGS and others, to model historic shorelines and to predict future shoreline positions; and applied this to Hawaiian Island beach data. Aerial photographs, sporadically taken throughout the past century, yield usable, albeit limited, historic shoreline data. These photographs provide excellent spatial coverage, but poor temporal resolution, of the shoreline. Due to the sparse historic shoreline data, and the many natural and anthropogenic events influencing coastlines, we constructed a simplistic shoreline change model that can identify long-term behavior of a beach. Our new, two-dimensional model combines the simple rate method to accommodate for cross-shore sediment transport with the classic Pelnard-Considère model for diffusion, as well as a longshore advection speed term. Inverse methods identify cross-shore rate, longshore advection speed, and longshore diffusivity down a sandy coastline. A spatial averaging technique then identifies shoreline segments where one parameter can reasonably account for the cross-shore and longshore transport rates in that area. This produces model results with spatial resolution more appropriate to the temporal spacing of the data. Because changes in historic data can be accounted for by varying degrees of cross-shore and longshore sediment transport - for example, beach erosion can equally be explained by sand moving either off-shore or laterally - we tested several different model scenarios on the data: allowing only cross-shore sediment movement, only longshore movement, and a combination of the two. We used statistical information criteria to determine both the optimal spatial resolution and best-fitting scenario. Finally, we employed a voting method predicting the relaxed shoreline position over time.
Soil and biosolid nano- and macro-colloid properties and contaminant transport behavior
NASA Astrophysics Data System (ADS)
Ghezzi, Jessique L.
Despite indications that they are potential contaminant transport systems and threats to groundwater quality, very little effort has been invested in comparing contaminant transport behavior of natural environmental nanocolloids and their corresponding macrocolloid fractions in the presence of As, Se, Pb, and Cu contaminants. This study involved physico-chemical, mineralogical, stability and contaminant-transport characterizations of nano- (< 100 nm) and macro-colloids (100-2000 nm) fractionated from three Kentucky soils and one biosolid waste. Particle size was investigated with SEM/TEM and dynamic light scattering. Surface reactivity was estimated using CEC and zeta potential. Mineralogical composition was determined by XRD, FTIR, and thermogravimetric analyses. Sorption isotherms assessed affinities for Cu2+, Pb2+, AsO3-, and SeO4 -2 contaminants, while settling kinetics experiments of suspensions at 0, 2 and 10 mg/L contaminants determined stability and transportability potential. Undisturbed 18x30 cm KY Ashton Loam soil monoliths were also used for transport experiments, involving infusion of 50 mg L-1 colloid suspensions spiked with 2 mg L-1 mixed contaminant loads in unsaturated, steady state, unit gradient downward percolation experiments. Overall, nanocolloids exhibited greater stability over corresponding macrocolloids in the presence and absence of contaminants following specific mineralogy trends. Physicochemical characterizations indicated that extensive organic carbon surface coatings and higher Al/Fe:Si ratios may have induced higher stability in the nanocolloid fractions, in spite of some hindrance by nano-aggregation phenomena. In the transport experiments, nanocolloids eluted significantly higher concentrations of colloids, total, and colloid-bound metals than corresponding macrocolloids. Contaminant elutions varied by colloid type, mineralogy and contaminant, with the following sequences: soil-colloids>bio-colloids, smectitic
Prediction of contaminant fate and transport in potable water systems using H2OFate
NASA Astrophysics Data System (ADS)
Devarakonda, Venkat; Manickavasagam, Sivakumar; VanBlaricum, Vicki; Ginsberg, Mark
2009-05-01
BlazeTech has recently developed a software called H2OFate to predict the fate and transport of chemical and biological contaminants in water distribution systems. This software includes models for the reactions of these contaminants with residual disinfectant in bulk water and at the pipe wall, and their adhesion/reactions with the pipe walls. This software can be interfaced with sensors through SCADA systems to monitor water distribution networks for contamination events and activate countermeasures, as needed. This paper presents results from parametric calculations carried out using H2OFate for a simulated contaminant release into a sample water distribution network.
Feed gas contaminant removal in ion transport membrane systems
Carolan, Michael Francis; Miller, Christopher Francis
2008-09-16
Method for gas purification comprising (a) obtaining a feed gas stream containing one or more contaminants selected from the group consisting of volatile metal oxy-hydroxides, volatile metal oxides, and volatile silicon hydroxide; (b) contacting the feed gas stream with a reactive solid material in a guard bed and reacting at least a portion of the contaminants with the reactive solid material to form a solid reaction product in the guard bed; and (c) withdrawing from the guard bed a purified gas stream.
FACT (Version 2.0) - Subsurface Flow and Contaminant Transport Documentation and User's Guide
Aleman, S.E.
2000-05-05
This report documents a finite element code designed to model subsurface flow and contaminant transport, named FACT. FACT is a transient three-dimensional, finite element code designed to simulate isothermal groundwater flow, moisture movement, and solute transport in variably saturated and fully saturated subsurface porous media.
Interest in coupled biodegradation and transport of organic contaminants has expanded greatly in the past several years. In a system in which biodegradation is coupled with solute transport, the magnitude and rate of biodegradation is influenced not only by pr...
Contaminant Transport in the Unsaturated Zone Theory and Modeling
Technology Transfer Automated Retrieval System (TEKTRAN)
Mathematical models are increasingly used to better understand and quantify site-specific subsurface water flow and solute transport processes. This chapter reviews mathematical models for solute transport in predominantly variably-saturated media. After a brief description of the physics and mathe...
WELL IMPACT: A CONTAMINANT FATE AND TRANSPORT MODEL
A chemical reactor analogy is used to derive analytical solutions for contaminant source dissolution in ground water (analytic overlay model). This allows the user to develop a source function based on a conceptual site model. The resulting analytical source function can then b...
Contaminant Transport in Two Central Missouri Karst Recharge Areas
Technology Transfer Automated Retrieval System (TEKTRAN)
Karst watersheds with significant losing streams represent a particularly vulnerable setting for ground water contamination because of the direct connection to surface water. Because of the existing agricultural land-use and future threat of heavy urbanization, two losing stream karst basins were ch...
Contaminant transport in two central Missouri karst recharge areas
Technology Transfer Automated Retrieval System (TEKTRAN)
Karst watersheds with significant losing streams represent a particularly vulnerable setting for ground water contamination because of the direct connection to surface water. Because of the existing agricultural land-use and future threat of heavy urbanization, two losing stream karst basins were ch...
EFFECTS OF NATURAL CYCLIC VARIATIONS ON CONTAMINATED FATE AND TRANSPORT
The studies provide the scientific community with a greater understanding of the physiochemical processes of sediment-contaminant interaction. A primary consideration in sediment clean-up is when to stop, or how clean is acceptable. Present mathematical models assume that ...
Contaminant transport in two central Missouri karst recharge areas
Technology Transfer Automated Retrieval System (TEKTRAN)
Karst watersheds with significant losing streams represent a particularly vulnerable setting for ground water contamination because of the direct connection to surface water. Because of the existing agricultural land-use and future likelihood of urbanization, two losing stream karst basins were chos...
EMERGING CONTAMINANTS IN THE WATER CYCLE: FATE AND TRANSPORT
In the past decade, the scientific community and general public have become increasingly aware of the potential for the presence of unregulated, and generally unmonitored contaminants, found at low concentrations in surface, ground and drinking water. The most common pathway for...
Flow hydrodynamics and contaminant transport in the flow past a lateral square cavity
NASA Astrophysics Data System (ADS)
Escauriaza, Cristian; Polanco, Juan Ignacio; August, Olivia; Bolster, Diogo
2015-11-01
Turbulent flows past lateral cavities play an important role in the transport of contaminants in rivers and streams. Cavities are surface storage zones, where large-scale unsteady coherent structures are the leading mechanisms that produce longer residence times and control the fate of contaminants in the river. In this work we study the recirculating flow and mass transport in a lateral square cavity, by performing numerical simulations with a hybrid URANS/LES turbulence model (DES-LR). We focus on the dynamics of the coherent structures and their impacts on the transport and storage of a passive scalar. In addition, we use the numerical results to develop new 1D models that improve the description of the evolution of the averaged concentration inside the cavity. By transferring the information to larger spatial scales, we provide new insights on the mechanisms of contaminant transport and analyze the overall effects of surface storage zones in open channel flows. supported by Fondecyt grant 1130940.
Influencing factors on particle-bound contaminant transport in the Elbe estuary
NASA Astrophysics Data System (ADS)
Kleisinger, Carmen; Haase, Holger; Schubert, Birgit
2016-04-01
Particulate matter, i.e. suspended particulate matter and sediments in rivers and estuaries, often are contaminated with trace metals and selected organic contaminants and are mainly associated with fine-grained fractions. Transport processes and fate of particles in estuaries are influenced by several factors, e.g. freshwater discharge, tide, flow velocity and dredging activities (Kappenberg et al., 2007). Understanding the transport processes in estuaries may help to achieve the objectives of the Water Framework Directive and the Marine Strategy Framework Directive. The German Federal Institute of Hydrology (BfG) operates for more than 20 years five monitoring sites in the Elbe estuary in order to monitor the development of particle-bound contaminant concentrations over time and to understand their transport mechanisms. Results of the monitoring revealed freshwater discharge as an important influencing factor on the transport of contaminated particulate matter (Ackermann et al., 2007). The bidirectional transport of marine and fluvial water and particulate matter in estuaries results in a turbidity zone where large amounts of particulate matter are temporarily retained and thus in a delayed transport of particulate matter towards the sea. The extent and the location of the turbidity zone as well as the ratio of highly contaminated fluvial and less contaminated marine sediments at a given location are mainly influenced by the freshwater discharge (Kowalewska et al., 2011). Furthermore, at high freshwater discharge conditions the highly contaminated particulate matter from fluvial origin are transported downstream the estuary, whereas at low freshwater discharges, upstream transport of less contaminated marine sediments prevails. Hence, residence times of particulate matter in the estuary are difficult to estimate. Furthermore, sedimentation areas with flow reduced conditions, e.g. wadden areas or branches of the Elbe estuary, may act as sinks for particle bound
Rasschaert, G; Houf, K; De Zutter, L
2007-01-01
The possible colonization of the intestines and contamination of broilers after transport to the slaughterhouse with Campylobacter strains present in cleaned and disinfected transport containers was investigated. Seven broiler flocks with a Campylobacter-free status were sampled once just before loading at the farm and once just before slaughter. On both occasions, samples were also taken from the exterior of the birds and from the intestinal content. Transport containers used to transport the flock were sampled on the farm just before loading the birds. Campylobacters were enumerated and genotyped by flagellin gene A PCR-restriction fragment length polymorphism and pulsed-field gel electrophoresis. In total, 25 of the 35 sampled containers were Campylobacter contaminated, and 30 genotypes were found. Three broiler flocks became colonized on the farm between initial status determination and transport to the slaughterhouse, and three Campylobacter-free flocks were externally contaminated after transport. In none of the seven flocks was evidence found of intestinal colonization or cocolonization due to transport in Campylobacter-contaminated containers. PMID:17265858
NASA Astrophysics Data System (ADS)
Padilla, I. Y.; Steele, K.
2008-05-01
The karst aquifer of the North Coast of Puerto Rico represents a significant source of water for drinking purposes, as well as ecosystem sustainability. The same characteristics making this aquifer the most productive in the island, fast infiltration and rapid flow in karst conduits, make the aquifers vulnerable highly vulnerable to contamination. Once in the ground water, organic contaminants move through the karst aquifers by complex pathways dictated by system characteristics and flow regimes. Ground water flow in karst aquifers is subscribed to two types of flow systems: conduit flow and diffuse flow. Transport in conduit-flow dominated systems tends to convey solutes rapidly through the system to a discharge or point without much attenuation. Transport in diffuse- flow systems, on the other hand, causes significant solute retardation and serves as a long-term source of contamination. Although it is common to attribute one type of predominant flow regime, most carbonate aquifers are characterized by a mixture of both flow systems. The north coast aquifer of Puerto Rico has been impacted by a large number of contaminates sites. During the last 25 years, 10 Superfund sites have been declared in the zone and others are being evaluated for inclusion in the National Priority List. The work presented herein addresses the potential impact of these sites on the extent of contamination and discusses the transport mechanisms affecting the transport and persistence of organic contaminants in the north coast aquifer of Puerto Rico. Preliminary evaluation indicates that fate and transport of these contaminants is controlled by a combinations of conduit- and diffuse-flow mechanisms, where conduits tend to concentrate water and contaminants and convey it rapidly or to "trapping" diffusive-flow zones of smaller pore-size zones.
Concentration polarization, surface currents, and bulk advection in a microchannel
NASA Astrophysics Data System (ADS)
Nielsen, Christoffer P.; Bruus, Henrik
2014-10-01
We present a comprehensive analysis of salt transport and overlimiting currents in a microchannel during concentration polarization. We have carried out full numerical simulations of the coupled Poisson-Nernst-Planck-Stokes problem governing the transport and rationalized the behavior of the system. A remarkable outcome of the investigations is the discovery of strong couplings between bulk advection and the surface current; without a surface current, bulk advection is strongly suppressed. The numerical simulations are supplemented by analytical models valid in the long channel limit as well as in the limit of negligible surface charge. By including the effects of diffusion and advection in the diffuse part of the electric double layers, we extend a recently published analytical model of overlimiting current due to surface conduction.
NASA Astrophysics Data System (ADS)
Yin, Y.; Sykes, J. F.
2006-12-01
Transport parameter estimation and contaminant source identification are critical steps in the development of a physically based groundwater contaminant transport model. For most transient field scale problems, the high computational burden required by parameter identification algorithms combined with sparse data sets often limits calibration. However, when data are available, a high performance computing system and parallel computing may make the calibration process feasible. The selection of the optimization algorithm is also critical. In this paper, the contaminant transport and source parameters were estimated and compared using optimization with two heuristic search algorithms (a dynamically dimensioned search and a parallelized micro genetic algorithm) and a gradient based multi-start PEST algorithm which were implemented on the Shared Hierarchical Academic Research Computing Network (Sharcnet). The case study is located in New Jersey where improper waste disposal resulted in the contamination of down gradient public water supply wells. Using FRAC3DVS, a physically based transient three-dimensional groundwater flow model with spatially and temporally varying recharge was developed and calibrated using both approximately 9 years of head data from continuous well records and data over a period of approximately 30 years from traditional monitoring wells. For the contaminant system, the parameters that were estimated include source leaching rate, source concentration, dispersivities, and retardation coefficient. The groundwater domain was discretized using 214,520 elements. With highly changing pump rates at the 7 municipal wells, time increments over the approximately 30 year simulation period varied dynamically between several days and 3 months. On Sharcnet, one forward simulation on a single processor of both transient flow and contaminant transport takes approximately 3 to 4 hours. The contaminant transport model calibration results indicate that overall
Wastewater contaminant transport and treatment in a nutrient limited ribbed fen
NASA Astrophysics Data System (ADS)
McCarter, C. P. R.; Price, J. S.; Branfireun, B. A.
2015-12-01
To minimize the discharge of wastewater contaminants from remote northern communities and mining operations, fen peatlands in sub-arctic regions are used for tertiary wastewater treatment to detain, transform, and remove these contaminants. However, there is a limited understanding of contaminant transport and treatment in fen peatlands, particularly in sub-arctic Canada. To better characterize wastewater contaminant transport and treatment in these systems, approximately 44 m3 day-1 of simulated wastewater, concentrated custom-blend fertilizer (NO3-, PO33-, and SO42-) and Cl- diluted with water, was pumped into a small 0.5 ha sub-arctic ribbed fen continuously for 47 days (July 15th -August 31st 2014). Contaminant concentration of 3 similar ribbed fens varied between 0.0-3.0 mg L-1 over the study period (May - September 2014). An exponential increase in transmissivity (2.4 to 16.8 m2 day-1) as the water table rose (~0.16 m) increased the average linear groundwater velocity (0.5 to 3.4 m day-1) and resulted in rapid SO42- (0.8 m day-1) and Cl- (1.9 m day-1) transport. Notwithstanding the rapid transport of Cl-, diffusion into inactive pores still retarded Cl- transport by a factor of 1.8. Contrary to the rapid transport of SO42- and Cl-, the other contaminants were rapidly removed from the pore water (likely through biological uptake or adsorption) and minimal transport was observed (0.29 and 0.04 m day-1 for PO33- and NO3-, respectively). Northern ribbed fens have a large capacity to detain certain wastewater contaminants (e.g., NO3- and PO33-), yet allow rapid transport of others (e.g., SO42- and Cl-). Thus, these peatlands have the potential to significantly decrease wastewater contamination in northern aquatic environment by both biogeochemical and physical processes but careful management of the hydrology is required to prevent the release of mobile contaminants.
POWELL, KIMBERLYR.
2004-05-25
Implementation of monitored natural attenuation (MNA) as a remediation method requires a mechanistic understanding of the natural attenuation processes occurring at a given site. For inorganic contaminants, natural attenuation typically involves a decrease in metal toxicity and/or mobility. These natural processes include dilution, dispersion, sorption (including adsorption, absorption, and precipitation), and redox processes. In order to better quantify these processes in terms of metal availability, sequential extraction experiments were carried out on subsurface soil samples impacted by a low pH, high sulfate, metals (Be, Ni, U, As) plume associated with the long-term operation of a coal plant at the Savannah River Site. These laboratory scale studies provide mechanistic information regarding the solid phases in the soils associated with natural attenuation of the contaminant metals. This data provides input to be evaluated in the definition of the contaminant source term as well as transport of contaminants for site transport models.
BUOYANT ADVECTION OF GASES IN UNSATURATED SOIL
Seely, Gregory E.; Falta, Ronald W.; Hunt, James R.
2010-01-01
In unsaturated soil, methane and volatile organic compounds can significantly alter the density of soil gas and induce buoyant gas flow. A series of laboratory experiments was conducted in a two-dimensional, homogeneous sand pack with gas permeabilities ranging from 110 to 3,000 darcy. Pure methane gas was injected horizontally into the sand and steady-state methane profiles were measured. Experimental results are in close agreement with a numerical model that represents the advective and diffusive components of methane transport. Comparison of simulations with and without gravitational acceleration permits identification of conditions where buoyancy dominates methane transport. Significant buoyant flow requires a Rayleigh number greater than 10 and an injected gas velocity sufficient to overcome dilution by molecular diffusion near the source. These criteria allow the extension of laboratory results to idealized field conditions for methane as well as denser-than-air vapors produced by volatilizing nonaqueous phase liquids trapped in unsaturated soil. PMID:20396624
Contaminant transport through a coal washery discard reactive wall
Gray, S.C.; Indraratna, B.; Yassini, I.
1999-07-01
This study examines the utilization of coal washery discard (CWD) as an inexpensive and readily available reactive wall material for the treatment of contaminated groundwater within the Illawarra region of New South Wales, Australia. The technology has been applied to a blast furnace slag (BFS) emplacement to attenuate an alkaline plume migrating towards a sensitive marine water body. Preliminary field performance data indicates that the CWD wall is reducing the pH of the plume to acceptable levels for marine ecosystems. The primary removal mechanisms within the CWD have been identified, however, have not be quantified at this stage.
Carrigan, C.R.; Hudson, G.B.
1998-06-01
'The intent of this research program is to obtain an improved understanding of vadose zone transport processes and to develop field and modeling techniques required to characterize contaminant transport in the unsaturated zone at DOE sites. For surface spills and near-surface leaks of chemicals, the vadose zone may well become a long-term source of contamination for the underlying water table. Transport of contaminants can occur in both the liquid and gas phases of the unsaturated zone. This transport occurs naturally as a result of diffusion, buoyancy forces (gravity), capillarity and barometric pressure variations. In some cases transport can be enhanced by anisotropies present in hydrologic regimes. This is particularly true for gas-phase transport which may be subject to vertical pumping resulting from atmospheric pressure changes. For liquid-phase flows, heterogeneity may enhance the downward transport of contaminants to the water table depending on soil properties and the scale of the surface spill or near-surface leak. Characterization techniques based upon the dynamics of transport processes are likely to yield a better understanding of the potential for contaminant transport at a specific site than methods depending solely on hydrologic properties derived from a borehole. Such dynamic-characterization techniques can be useful for evaluating sites where contamination presently exists as well as for providing an objective basis to evaluate the efficacy of proposed as well as implemented clean-up technologies. The real-time monitoring of processes that may occur during clean-up of tank waste and the mobility of contaminants beneath the Hanford storage tanks during sluicing operations is one example of how techniques developed in this effort can be applied to current remediation problems. In the future, such dynamic-characterization methods might also be used as part of the site-characterization process for determining suitable locations of new DOE facilities
Gorelick, S.M.; Voss, C.I.; Gill, P.E.; Murray, W.; Saunders, M.A.; Wright, M.H.
1984-01-01
A simulation-management methodology is demonstrated for the rehabilitation of aquifers that have been subjected to chemical contamination. Finite element groundwater flow and contaminant transport simulation are combined with nonlinear optimization. The model is capable of determining well locations plus pumping and injection rates for groundwater quality control. Examples demonstrate linear or nonlinear objective functions subject to linear and nonlinear simulation and water management constraints. -from Authors
TRANSPORT OF CHEMICAL CONTAMINANTS IN KARST TERRANES: OUTLINE AND SUMMARY
Chemical spills that reach an aquifer in karst terranes do not behave like those in granular or highly fractured aquifers. pills reaching diffuse-flow aquifers display relatively slow transport, are radially dispersive, and can be tracked through the use of monitoring wells. pill...
Analysis of heat transfer and contaminant transport in fume hoods
Pathanjali, C.; Rahman, M.M.
1996-12-31
The paper presents the analysis of three-dimensional flow patterns and the associated heat and mass transfer mechanisms in a fume hood enclosure. The flow enters the hood through the front window opening (positive x-direction) and leaves the cupboard through an opening on the top of the hood (positive z-direction). The flow was assumed to be fully turbulent. The flow pattern for different sash openings were studied. The flow pattern around an object located at the bottom of the hood was studied for different locations of the object. It was found that air entering the hood proceeds directly to the back wall, impinges it and turns upward toward the top wall and exits through the outlet. The flow finds its way around any object forming a recirculating region at its training surface. With an increase in the sash opening, the velocity becomes higher and the fluid traces the path to the outlet more quickly. The volume occupied by recirculating flow decreases with increase in sash opening. Both temperature and concentration were found to be maximum near the source and gradually decreased as the heated air or gaseous contaminant entrained with incoming air. The local concentration decreased with increase in sash opening area. The results will be very useful to design experiments with optimum sash opening providing adequate disposal of contaminants with minimum use of conditioned air inside the room.
Advances in Dynamic Transport of Organic Contaminants in Karst Groundwater Systems
NASA Astrophysics Data System (ADS)
Padilla, I. Y.; Vesper, D.; Alshawabkeh, A.; Hellweger, F.
2011-12-01
Karst groundwater systems develop in soluble rocks such as limestone, and are characterized by high permeability and well-developed conduit porosity. These systems provide important freshwater resources for human consumption and ecological integrity of streams, wetlands, and coastal zones. The same characteristics that make karst aquifers highly productive make them highly vulnerable to contamination. As a result, karst aquifers serve as an important route for contaminants exposure to humans and wildlife. Transport of organic contaminants in karst ground-water occurs in complex pathways influenced by the flow mechanism predominating in the aquifer: conduit-flow dominated systems tend to convey solutes rapidly through the system to a discharge point without much attenuation; diffuse-flow systems, on the other hand, can cause significant solute retardation and slow movement. These two mechanisms represent end members of a wide spectrum of conditions found in karst areas, and often a combination of conduit- and diffuse-flow mechanisms is encountered, where both flow mechanisms can control the fate and transport of contaminants. This is the case in the carbonate aquifers of northern Puerto Rico. This work addresses advances made on the characterization of fate and transport processes in karst ground-water systems characterized by variable conduit and/or diffusion dominated flow under high- and low-flow conditions. It involves laboratory-scale physical modeling and field-scale sampling and historical analysis of contaminant distribution. Statistical analysis of solute transport in Geo-Hydrobed physical models shows the heterogeneous character of transport dynamics in karstic units, and its variability under different flow regimes. Field-work analysis of chlorinated volatile organic compounds and phthalates indicates a large capacity of the karst systems to store and transmit contaminants. This work is part of the program "Puerto Rico Testsite for Exploring Contamination
Advective Mechanisms in Tree Island Formation
NASA Astrophysics Data System (ADS)
Stothoff, S.
2002-05-01
Tree islands are important landscape features in the Florida Everglades. Tres islands are formed of peat deposited on the shallow limestone bedrock, and have been stressed as the system has changed in response to anthropogenic activities due to the sensitivity of organic soils to hydrologic cycles. The plume shape aligned with flow direction for typical tree islands is characteristic of advective transport, despite the rather low flow velocities in the system. Hypothesized mechanisms for the plume shape include sediment transport downstream from the head of the island (often anchored by a bedrock rise), or nutrient transport downstream allowing plants to produce more sediments in situ. Understanding mechanisms controlling tree island shape will aid in understanding the response of tree islands to hydrologic management. An integrated system of field, laboratory, and modeling studies is underway, with the first effort aimed at bounding the importance of the simpler sediment transport processes before tackling more-complex nutrient transport processes. The numerical model integrating the field and laboratory efforts is a 3D finite volume model considering water flow in the shallow groundwater/surface-water system together with sediment transport. The model can account for variable vegetative resistance through the flow column, including the important case where a dense mat forms at the surface. Model components specific for this system and associated data requirements are presented.
Color image analysis of contaminants and bacteria transport in porous media
Rashidi, M.; Dehmeshid, J.; Dacini, F.; Cole, L.; Dickenson, E.
1997-07-01
Transport of contaminants and bacteria in aqueous heterogeneous saturated porous systems have been studies experimentally using a novel fluorescent microscopic imaging (FMI) technique. The approach involves color visualization and quantification of bacterium and contaminant distributions within a transparent porous column. By introducing stained bacteria and an organic dye as a contaminant into the column and illuminating the porous regions with a planar sheet of laser beam, contaminant and bacterial transport processes through the porous medium can be observed and measured microscopically. A computer controlled CCD camera is used to record the fluorescent images as a function of time. These images are recorded by a frame accurate high resolution VCR and are then analyzed using a color image analysis code written in our laboratories. The color images are digitalized this way and simultaneously concentration and velocity distributions of both contaminant and bacterium are evaluated as a function of time and pore characteristics. The approach provides a unique dynamic probe to observe these transport processes microscopically. These results are extremely valuable in in-situ bioremediation problems since microscopic particle-contaminant-bacterium interactions are the key to understanding and optimization of these processes.
Bourigault, C; Daniel, L; Jourdain, S; Hardy, E; Heriaud, K; Virmaux, M; Eniafe-Eveillard, B; Lejeune, B
2007-03-01
From September to December 2004, contaminations were found in fifteen conservation and transport mediums of corneal grafts at the tissue bank of Brest, including seven by Sphingomonas paucimobilis. The pulsed-field gel electrophoresis made it possible to establish the genotypic profiles of each strain and to compare them. Similarities were found between certain strains of the contaminated mediums and those of the thermostated double boiler of the tissue bank. The link between the contamination and the defrosting of the mediums in the double boiler was thus established. Measures of prevention are currently proposed to defrost the bottles like the use of a dry bath to replace the current one. PMID:16631318
NASA Astrophysics Data System (ADS)
Ostrowski, J.; Shlomi, S.; Michalak, A.
2007-12-01
The process of estimating the release history of a contaminant in groundwater relies on coupling a limited number of concentration measurements with a groundwater flow and transport model in an inverse modeling framework. The information provided by available measurements is generally not sufficient to fully characterize the unknown release history; therefore, an accurate assessment of the estimation uncertainty is required. The modeler's level of confidence in the transport parameters, expressed as pdfs, can be incorporated into the inverse model to improve the accuracy of the release estimates. In this work, geostatistical inverse modeling is used in conjunction with Monte Carlo sampling of transport parameters to estimate groundwater contaminant release histories. Concentration non-negativity is enforced using a Gibbs sampling algorithm based on a truncated normal distribution. The method is applied to two one-dimensional test cases: a hypothetical dataset commonly used in validating contaminant source identification methods, and data collected from a tetrachloroethylene and trichloroethylene plume at the Dover Air Force Base in Delaware. The estimated release histories and associated uncertainties are compared to results from a geostatistical inverse model where uncertainty in transport parameters is ignored. Results show that the a posteriori uncertainty associated with the model that accounts for parameter uncertainty is higher, but that this model provides a more realistic representation of the release history based on available data. This modified inverse modeling technique has many applications, including assignment of liability in groundwater contamination cases, characterization of groundwater contamination, and model calibration.
INFLUENCE OF COUPLED PROCESSES ON CONTAMINANT FATE AND TRANSPORT IN SUBSURFACE ENVIRONMENTS
Jardine, Philip M
2008-01-01
The following chapter emphasizes subsurface environmental research investigations over the past 10 to 15 years that couple hydrological, geochemical, and biological processes as related to contaminant fate and transport. An attempt is made to focus on field-scale studies with possible reference to laboratory-scale endeavors. Much of the research discussed reflects investigations of the influence of coupled processes on the fate and transport of inorganic, radionuclide, and organic contaminants in subsurface environments as a result of natural processes or energy and weapons production endeavors that required waste disposal. The chapter provides on overview of the interaction between hydro-bio-geochemical processes in structured, heterogeneous subsurface environments and how these interactions control contaminant fate and transport, followed by experimental and numerical subsurface science research and case studies involving specific classes of inorganic and organic contaminants. Lastly, thought provoking insights are highlighted on why the study of subsurface coupled processes is paramount to understanding potential future contaminant fate and transport issues of global concern.
Modeling the emission, transport and deposition of contaminated dust from a mine tailing site.
Stovern, Michael; Betterton, Eric A; Sáez, A Eduardo; Villar, Omar Ignacio Felix; Rine, Kyle P; Russell, Mackenzie R; King, Matt
2014-01-01
Mining operations are potential sources of airborne particulate metal and metalloid contaminants through both direct smelter emissions and wind erosion of mine tailings. The warmer, drier conditions predicted for the Southwestern US by climate models may make contaminated atmospheric dust and aerosols increasingly important, due to potential deleterious effects on human health and ecology. Dust emissions and dispersion of contaminants from the Iron King Mine tailings in Dewey-Humboldt, Arizona, a Superfund site, are currently being investigated through in situ field measurements and computational fluid dynamics modeling. These tailings are significantly contaminated with lead and arsenic with an average soil concentration of 1616 and 1420 ppm, respectively. Similar levels of these contaminants have also been measured in soil samples taken from the area surrounding the mine tailings. Using a computational fluid dynamics model, we have been able to model dust transport from the mine tailings to the surrounding region. The model includes a distributed Eulerian model to simulate fine aerosol transport and a Lagrangian approach to model fate and transport of larger particles. In order to improve the accuracy of the dust transport simulations both regional topographical features and local weather patterns have been incorporated into the model simulations. PMID:24552963
Impact of scaling and body movement on contaminant transport in airliner cabins
NASA Astrophysics Data System (ADS)
Mazumdar, Sagnik; Poussou, Stephane B.; Lin, Chao-Hsin; Isukapalli, Sastry S.; Plesniak, Michael W.; Chen, Qingyan
2011-10-01
Studies of contaminant transport have been conducted using small-scale models. This investigation used validated Computational Fluid Dynamics (CFD) to examine if a small-scale water model could reveal the same contaminant transport characteristics as a full-scale airliner cabin. But due to similarity problems and the difficulty of scaling the geometry, a perfect scale up from a small water model to an actual air model was found to be impossible. The study also found that the seats and passengers tended to obstruct the lateral transport of the contaminants and confine their spread to the aisle of the cabin. The movement of a crew member or a passenger could carry a contaminant in its wake to as many rows as the crew member or passenger passed. This could be the reason why a SARS infected passenger could infect fellow passengers who were seated seven rows away. To accurately simulate the contaminant transport, the shape of the moving body should be a human-like model.
Organic contaminant transport and fate in the subsurface: evolution of knowledge and understanding
Essaid, Hedeff I.; Bekins, Barbara A.; Cozzarelli, Isabelle M.
2015-01-01
Toxic organic contaminants may enter the subsurface as slightly soluble and volatile nonaqueous phase liquids (NAPLs) or as dissolved solutes resulting in contaminant plumes emanating from the source zone. A large body of research published in Water Resources Research has been devoted to characterizing and understanding processes controlling the transport and fate of these organic contaminants and the effectiveness of natural attenuation, bioremediation, and other remedial technologies. These contributions include studies of NAPL flow, entrapment, and interphase mass transfer that have advanced from the analysis of simple systems with uniform properties and equilibrium contaminant phase partitioning to complex systems with pore-scale and macroscale heterogeneity and rate-limited interphase mass transfer. Understanding of the fate of dissolved organic plumes has advanced from when biodegradation was thought to require oxygen to recognition of the importance of anaerobic biodegradation, multiple redox zones, microbial enzyme kinetics, and mixing of organic contaminants and electron acceptors at plume fringes. Challenges remain in understanding the impacts of physical, chemical, biological, and hydrogeological heterogeneity, pore-scale interactions, and mixing on the fate of organic contaminants. Further effort is needed to successfully incorporate these processes into field-scale predictions of transport and fate. Regulations have greatly reduced the frequency of new point-source contamination problems; however, remediation at many legacy plumes remains challenging. A number of fields of current relevance are benefiting from research advances from point-source contaminant research. These include geologic carbon sequestration, nonpoint-source contamination, aquifer storage and recovery, the fate of contaminants from oil and gas development, and enhanced bioremediation.
Organic contaminant transport and fate in the subsurface: Evolution of knowledge and understanding
NASA Astrophysics Data System (ADS)
Essaid, Hedeff I.; Bekins, Barbara A.; Cozzarelli, Isabelle M.
2015-07-01
Toxic organic contaminants may enter the subsurface as slightly soluble and volatile nonaqueous phase liquids (NAPLs) or as dissolved solutes resulting in contaminant plumes emanating from the source zone. A large body of research published in Water Resources Research has been devoted to characterizing and understanding processes controlling the transport and fate of these organic contaminants and the effectiveness of natural attenuation, bioremediation, and other remedial technologies. These contributions include studies of NAPL flow, entrapment, and interphase mass transfer that have advanced from the analysis of simple systems with uniform properties and equilibrium contaminant phase partitioning to complex systems with pore-scale and macroscale heterogeneity and rate-limited interphase mass transfer. Understanding of the fate of dissolved organic plumes has advanced from when biodegradation was thought to require oxygen to recognition of the importance of anaerobic biodegradation, multiple redox zones, microbial enzyme kinetics, and mixing of organic contaminants and electron acceptors at plume fringes. Challenges remain in understanding the impacts of physical, chemical, biological, and hydrogeological heterogeneity, pore-scale interactions, and mixing on the fate of organic contaminants. Further effort is needed to successfully incorporate these processes into field-scale predictions of transport and fate. Regulations have greatly reduced the frequency of new point-source contamination problems; however, remediation at many legacy plumes remains challenging. A number of fields of current relevance are benefiting from research advances from point-source contaminant research. These include geologic carbon sequestration, nonpoint-source contamination, aquifer storage and recovery, the fate of contaminants from oil and gas development, and enhanced bioremediation.
Feed gas contaminant control in ion transport membrane systems
Carolan, Michael Francis; Minford, Eric; Waldron, William Emil
2009-07-07
Ion transport membrane oxidation system comprising an enclosure having an interior and an interior surface, inlet piping having an internal surface and adapted to introduce a heated feed gas into the interior of the enclosure, and outlet piping adapted to withdraw a product gas from the interior of the enclosure; one or more planar ion transport membrane modules disposed in the interior of the enclosure, each membrane module comprising mixed metal oxide material; and a preheater adapted to heat a feed gas to provide the heated feed gas to the inlet piping, wherein the preheater comprises an interior surface. Any of the interior surfaces of the enclosure, the inlet piping, and the preheater may be lined with a copper-containing metal lining. Alternatively, any of the interior surfaces of the inlet piping and the preheater may be lined with a copper-containing metal lining and the enclosure may comprise copper.
Consistency problem with tracer advection in the Atmospheric Model GAMIL
NASA Astrophysics Data System (ADS)
Zhang, Kai; Wan, Hui; Wang, Bin; Zhang, Meigen
2008-03-01
The radon transport test, which is a widely used test case for atmospheric transport models, is carried out to evaluate the tracer advection schemes in the Grid-Point Atmospheric Model of IAP-LASG (GAMIL). Two of the three available schemes in the model are found to be associated with significant biases in the polar regions and in the upper part of the atmosphere, which implies potentially large errors in the simulation of ozone-like tracers. Theoretical analyses show that inconsistency exists between the advection schemes and the discrete continuity equation in the dynamical core of GAMIL and consequently leads to spurious sources and sinks in the tracer transport equation. The impact of this type of inconsistency is demonstrated by idealized tests and identified as the cause of the aforementioned biases. Other potential effects of this inconsistency are also discussed. Results of this study provide some hints for choosing suitable advection schemes in the GAMIL model. At least for the polar-region-concentrated atmospheric components and the closely correlated chemical species, the Flux-Form Semi-Lagrangian advection scheme produces more reasonable simulations of the large-scale transport processes without significantly increasing the computational expense.
ADVECTION INFLUENCES ON EVAPOTRANSPIRATION OF ALFALFA IN A SEMIARID ENVIRONMENT
Technology Transfer Automated Retrieval System (TEKTRAN)
Advective enhancement of crop evapotranspiration (ET) occurs when drier, hotter air is transported into the crop by wind and can be an important factor in the water balance of irrigated crops in a semiarid climate. Thirteen days of moderate to extremely high ET rates of irrigated alfalfa (Medicago ...
NASA Astrophysics Data System (ADS)
Thorbjarnarson, Kathryn W.; Mackay, Douglas M.
1997-01-01
A tracer-elution experiment was conducted in a 9-m-thick alluvial sand aquifer at the Rocky Mountain Arsenal, Denver, Colorado, within an extensive 1,1,1-trichloroethene and trichloroethene plume. The forced-gradient flow field was controlled by an injection well and an extraction well separated by 8.4 m and aligned in the direction of the natural-gradient flow. Upon extraction, the contaminant-laden water was treated by air stripping and reinjected into the injection well. Iodide tracer was added to the injection flow during the initial 27.5 h of the experiment. Tracer transport and organic contaminant elution were monitored by four 0.15-m-screened drive points and a fully penetrating monitoring well. Relative permeabilities, dispersivities and retardation factors were estimated from tracer breakthrough and contaminant elution curves by the moment method and by curve-fitting with an advection-dispersion model. Tracer transport through the four strata sampled by the drive points indicated a permeability variation of three orders of magnitude. Contaminant elution was not observed in the lowest-permeability stratum monitored during the experiment. In all monitored strata, contaminant elution was controlled primarily by permeability effects on water flow and exhibited minimal retardation or desorption effects. The fully penetrating monitoring well exhibited a tracer response primarily from the more permeable strata with the addition of tracer from the less permeable strata producing an increased breakthrough spreading. This increased spreading or dispersion was reflected in a higher longitudinal dispersivity estimate (1.2 m assuming a homogeneous aquifer) than dispersivity estimates from the drive-point sampler tracer curves (ranging from 5 to 21 cm). Contaminant elution curves from the fully penetrating monitoring well exhibited an initial response primarily from the more permeable strata (rapid elution of contaminants) and provided no insight into the elution
NASA Astrophysics Data System (ADS)
Mosthaf, K.; Rosenberg, L.; Balbarini, N.; Broholm, M. M.; Bjerg, P. L.; Binning, P. J.
2014-12-01
It is important to understand the fate and transport of contaminants in limestone aquifers because they are a major drinking water resource. This is challenging because they are highly heterogeneous; with micro-porous grains, flint inclusions, and being heavily fractured. Several modeling approaches have been developed to describe contaminant transport in fractured media, such as the discrete fracture (with various fracture geometries), equivalent porous media (with and without anisotropy), and dual porosity models. However, these modeling concepts are not well tested for limestone geologies. Given available field data and model purpose, this paper therefore aims to develop, examine and compare modeling approaches for transport of contaminants in fractured limestone aquifers. The model comparison was conducted for a contaminated site in Denmark, where a plume of a dissolved contaminant (PCE) has migrated through a fractured limestone aquifer. Multilevel monitoring wells have been installed at the site and available data includes information on spill history, extent of contamination, geology and hydrogeology. To describe the geology and fracture network, data from borehole logs was combined with an analysis of heterogeneities and fractures from a nearby excavation (analog site). Methods for translating the geological information and fracture mapping into each of the model concepts were examined. Each model was compared with available field data, considering both model fit and measures of model suitability. An analysis of model parameter identifiability and sensitivity is presented. Results show that there is considerable difference between modeling approaches, and that it is important to identify the right one for the actual scale and model purpose. A challenge in the use of field data is the determination of relevant hydraulic properties and interpretation of aqueous and solid phase contaminant concentration sampling data. Traditional water sampling has a bias
Transport of organic contaminants in geomembranes under stress
Xiao, S. ); Moresoli, C. . Dept. of Chemical Engineering); Burczyk, A. ); Pintauro, P.; De Kee, D. . Dept. of Chemical Engineering)
1999-07-01
The transport properties of aqueous solutions of benzene, dichloromethane, and trichloroethylene through extended polyvinyl chloride and high-density polyethylene geomembranes are investigated. It is found that extension enhances the permeation rates of the penetrants through polyvinyl chloride geomembranes, and the opposite effect is found in the case of high-density polyethylene. This difference in response is attributed to the type of structural change, which occurred as a result of the extension. The diffusivities of a mixture of the three organic solvents through the geomembranes are also determined.
Through the processes of atmospheric transport and deposition, many anthropogenic contaminants such as industrial organics, pesticides, and trace metals have become widely distributed around the globe. ue to the phenomenon of long-range atmospheric transport, even the most remote...
Simonds, F. William
2002-01-01
A three-dimensional finite-difference steady-state ground-water flow model was developed to simulate hydraulic conditions at the Area 6 Landfill, Naval Air Station Whidbey Island, near Oak Harbor, Washington. Remediation efforts were started in 1995 in an attempt to contain trichloroethene and other contaminants in the ground water. The model was developed as a tool to test the effectiveness of the pump-and-treat remediation efforts as well as alternative remediation strategies. The model utilized stratigraphic data from approximately 76 Navy and 19 private wells to define the geometry of the shallow, intermediate, and deep aquifers and the intervening confining layers. Initial aquifer parameters and recharge estimates from aquifer tests and published remedial investigation reports were used in the model and then adjusted until simulated water levels closely matched observed water-level data collected prior to the onset of remediation in 1995. The calibrated model was then modified to depict the remedial pump-and-treat system, in which contaminated ground water is extracted, treated, and returned to the ground surface for infiltration. The water levels simulated by the modified model were compared with observed water levels for the 1998 calendar year, during which time the pump-and-treat system was in nearly continuous operation and the ground-water system had equilibrated to steady-state conditions. Although artificial boundaries were used in the model, the choice of model boundary conditions was simulation in the area of primary concern surrounding the western contaminant plume and extraction wells. Particle tracking results indicate that the model can effectively simulate the advective transport of contaminants from the source area to the pumping wells and thus be used to test alternative remedial pumping strategies.