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...
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...
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...
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...
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...
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...
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
Solution of the advection-dispersion equation: Continuous load of finite duration
Runkel, R.L.
1996-01-01
Field studies of solute fate and transport in streams and rivers often involve an. experimental release of solutes at an upstream boundary for a finite period of time. A review of several standard references on surface-water-quality modeling indicates that the analytical solution to the constant-parameter advection-dispersion equation for this type of boundary condition has been generally overlooked. Here an exact analytical solution that considers a continuous load of unite duration is compared to an approximate analytical solution presented elsewhere. Results indicate that the exact analytical solution should be used for verification of numerical solutions and other solute-transport problems wherein a high level of accuracy is required. ?? ASCE.
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...
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....
İbiş, Birol
2014-01-01
This paper aims to obtain the approximate solution of time-fractional advection-dispersion equation (FADE) involving Jumarie's modification of Riemann-Liouville derivative by the fractional variational iteration method (FVIM). FVIM provides an analytical approximate solution in the form of a convergent series. Some examples are given and the results indicate that the FVIM is of high accuracy, more efficient, and more convenient for solving time FADEs. PMID:24578662
Ibiş, Birol; Bayram, Mustafa
2014-01-01
This paper aims to obtain the approximate solution of time-fractional advection-dispersion equation (FADE) involving Jumarie's modification of Riemann-Liouville derivative by the fractional variational iteration method (FVIM). FVIM provides an analytical approximate solution in the form of a convergent series. Some examples are given and the results indicate that the FVIM is of high accuracy, more efficient, and more convenient for solving time FADEs. PMID:24578662
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)
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...
Healy, R.W.; Russell, T.F.
1992-01-01
A finite-volume Eulerian-Lagrangian local adjoint method for solution of the advection-dispersion equation is developed and discussed. The method is mass conservative and can solve advection-dominated ground-water solute-transport problems accurately and efficiently. An integrated finite-difference approach is used in the method. A key component of the method is that the integral representing the mass-storage term is evaluated numerically at the current time level. Integration points, and the mass associated with these points, are then forward tracked up to the next time level. The number of integration points required to reach a specified level of accuracy is problem dependent and increases as the sharpness of the simulated solute front increases. Integration points are generally equally spaced within each grid cell. For problems involving variable coefficients it has been found to be advantageous to include additional integration points at strategic locations in each well. These locations are determined by backtracking. Forward tracking of boundary fluxes by the method alleviates problems that are encountered in the backtracking approaches of most characteristic methods. A test problem is used to illustrate that the new method offers substantial advantages over other numerical methods for a wide range of problems.
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
Healy, R.W.; Russell, T.F.
1998-01-01
We extend the finite-volume Eulerian-Lagrangian localized adjoint method (FVELLAM) for solution of the advection-dispersion equation to two dimensions. The method can conserve mass globally and is not limited by restrictions on the size of the grid Peclet or Courant number. Therefore, it is well suited for solution of advection-dominated ground-water solute transport problems. In test problem comparisons with standard finite differences, FVELLAM is able to attain accurate solutions on much coarser space and time grids. On fine grids, the accuracy of the two methods is comparable. A critical aspect of FVELLAM (and all other ELLAMs) is evaluation of the mass storage integral from the preceding time level. In FVELLAM this may be accomplished with either a forward or backtracking approach. The forward tracking approach conserves mass globally and is the preferred approach. The backtracking approach is less computationally intensive, but not globally mass conservative. Boundary terms are systematically represented as integrals in space and time which are evaluated by a common integration scheme in conjunction with forward tracking through time. Unlike the one-dimensional case, local mass conservation cannot be guaranteed, so slight oscillations in concentration can develop, particularly in the vicinity of inflow or outflow boundaries. Published by Elsevier Science Ltd.
Healy, R.W.; Russell, T.F.
1993-01-01
Test results demonstrate that the finite-volume Eulerian-Lagrangian localized adjoint method (FVELLAM) outperforms standard finite-difference methods for solute transport problems that are dominated by advection. FVELLAM systematically conserves mass globally with all types of boundary conditions. Integrated finite differences, instead of finite elements, are used to approximate the governing equation. This approach, in conjunction with a forward tracking scheme, greatly facilitates mass conservation. The mass storage integral is numerically evaluated at the current time level, and quadrature points are then tracked forward in time to the next level. Forward tracking permits straightforward treatment of inflow boundaries, thus avoiding the inherent problem in backtracking of characteristic lines intersecting inflow boundaries. FVELLAM extends previous results by obtaining mass conservation locally on Lagrangian space-time elements. -from Authors
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.
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
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...
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.
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...
FRACTIONAL SOLUTE TRANSPORT EQUATION EVALUATED WITH THE MISCIBLE DISPLACEMENT EXPERIMENTAL DATA
Technology Transfer Automated Retrieval System (TEKTRAN)
A new solute transport model has been recently developed assuming that the movements of solute particles in hierarchically-structured porous media belongs to the family of Lévy motions rather than to the Brownian motion. The one-dimensional fractional advective-dispersive transport equation, or FADE...
Technology Transfer Automated Retrieval System (TEKTRAN)
Solute transport in soils and sediments is commonly simulated with the parabolic advective-dispersive equation, or ADE. In the last decades, it has been reported that this model cannot take in account several important features of solute movement through soil. Recently, a new model base on the assu...
A Lagrangian-Eulerian finite element method with adaptive gridding for advection-dispersion problems
Ijiri, Y.; Karasaki, K.
1994-02-01
In the present paper, a Lagrangian-Eulerian finite element method with adaptive gridding for solving advection-dispersion equations is described. The code creates new grid points in the vicinity of sharp fronts at every time step in order to reduce numerical dispersion. The code yields quite accurate solutions for a wide range of mesh Peclet numbers and for mesh Courant numbers well in excess of 1.
NASA Astrophysics Data System (ADS)
Dean, A. M.; Benson, D. A.; Major, E.
2010-12-01
By adding a fractional-in-time term to the traditional advection dispersion equation, a model is able to simulate a late-time heavy-tailed contaminant breakthrough curve. This heavy-tailed breakthrough curve is observed in data collected during a conservative tracer “push-pull” test at the Macrodispersion Experiment (MADE) site. A time fractional advection dispersion equation (fADE) is able to predict power law tailing of conservative solutes by accounting for solutes transferring between the mobile and relatively immobile phases. Solutes can become trapped in a low permeability zone where the transport is controlled by diffusion instead of advection. It has been observed that the late-time heavy-tailed breakthrough curve may follow a power law due to the movement into these low flow zones. By solving the time fADE in a particle tracking program (SLIM-FAST) the model accounts for mass transfer between various phases and produces the same power law tail as observed in field data. For the implementation of the time fADE, in SLIM-FAST, the particles move based on a random-walk motion but have the ability to transition into a relatively immobile phase after (exponentially) random mobile times. Following a period in the immobile phase, the particle re-enters the mobile phase to be moved by advection and Fickian dispersion. To test the fADE approach, a recent single-well push-pull tracer test at the MADE site is reproduced using a groundwater flow code (ParFlow) and a particle tracking code (SLIM-FAST) using various immobile residence-time distributions.
NASA Astrophysics Data System (ADS)
Benson, D. A.; Zhang, Y.
2006-12-01
Conservative solute transport through natural media is typically "anomalous" or non-Fickian. The anomalous transport may be characterized by faster than linear growth of the centered second moment, or non-Gaussian leading or trailing edges of a plume emanating from a point source. These characteristics develop because of non-local dependence on either past (time) or far upstream (space) concentrations. Non-local equations developed to describe anomalous dispersion usually focus on constant transport parameters and/or independence of the transport on space dimension. These simplifications have been useful for fitting simple transport processes, such as laboratory column tests or 1-D projections of field data. However, they may be insufficient for real field settings, where direction-dependent depositional processes and nonstationary heterogeneity can occur. We develop a generalized, multi-dimensional, spatiotemporal fractional advection- dispersion equation (fADE) with variable parameters to characterize regional-scale anomalous dispersion processes including trapping in immobile zones and/or super-Fickian rapid transport. A Lagrangian numerical model of the space-time fractional transport equation is developed in which solute particles can disperse in both space and time, depending on the medium heterogeneity properties, such as the connectivity and statistical distributions of high versus low-permeability deposits. In the generalized fADE, the range of the order of fractional time derivative is (0 2], representing a wide range of possible trapping behavior. The extension of the order to the range (1 2] is novel to transport theory. We apply the numerical model in 1-D and 2-D to the MADE site tritium plumes, and results indicate that this method can capture the main behaviors of realistic plumes, including local variations of spreading, direction-dependent scaling rates, and arbitrary rapid transport along preferential flow paths. Since the governing equation
NASA Astrophysics Data System (ADS)
Tertre, E.; Hubert, F.; Bruzac, S.; Pacreau, M.; Ferrage, E.; Prêt, D.
2013-07-01
The present study aims at testing the validity of using an Na+/Ca2+ ion-exchange model, derived from batch data to interpret experimental Ca2+-for-Na+ exchange breakthrough curves obtained on vermiculite (a common swelling clay mineral in surface environments). The ion-exchange model was constructed considering the multi-site nature of the vermiculite surface as well as the exchange of all aqueous species (Mg2+ derived from the dissolution of the solid and H+). The proposed ion-exchange model was then coupled with a transport model, and the predicted breakthrough curves were compared with the experimental ones obtained using a well stirred flow-through reactor. For a given solute residence time in the reactor (typically 50 min), our thermodynamic model based on instantaneous equilibrium was found to accurately reproduce several of the experimental breakthrough curves, depending on the Na+ and Ca2+ concentrations of the influents pumped through the reactor. However the model failed to reproduce experimental breakthrough curves obtained at high flow rates and low chemical gradient between the exchanger phase and the solution. An alternative model based on a hybrid equilibrium/kinetic approach was thus used and allowed predicting experimental data. Based on these results, we show that a simple parameter can be used to differentiate between thermodynamic and kinetic control of the exchange reaction with water flow. The results of this study are relevant for natural systems where two aquatic environments having contrasted chemistries interact. Indeed, the question regarding the attainment of a full equilibrium in such a system during the contact time of the aqueous phase with the particle/colloid remains most often open. In this context, we show that when a river (a flow of fresh water) encounters marine colloids, a systematic full equilibrium can be assumed (i.e., the absence of kinetic effects) when the residence time of the solute in 1 m3 of the system is ⩾6200 h.
Solute transport in solution conduits exhibiting multi-peaked breakthrough curves
NASA Astrophysics Data System (ADS)
Field, Malcolm S.; Leij, Feike J.
2012-05-01
SummarySolute transport in karst aquifers is primarily constrained to solution conduits where transport is rapid, turbulent, and relatively unrestrictive. Breakthrough curves generated from tracer tests are typically positively-skewed and may exhibit multiple peaks. In order to understand the circumstances under which multi-peaked positively skewed breakthrough curves occur, physical experiments utilizing single- and multiple-flow channels were conducted. Experiments also included waterfalls, short-term solute detention in pools, and flow obstructions. Results demonstrated that breakthrough curve skewness nearly always occurs to some degree but is magnified as immobile-flow regions are encountered. Multi-peaked breakthrough curves occurred when flow in the main channel became partially occluded from blockage in the main channel that forced divergence of solute into auxiliary channels and when waterfalls and detention in pools occurred. Currently, multi-peaked breakthrough curves are fitted by a multi-dispersion model in which a series of curves generated by the advection-dispersion equation are fitted to each measured peak by superimposing the measured breakthrough curve to obtain a combined model fit with a consequent set of estimated velocities and dispersions. In this paper, a dual-advection dispersion equation with first-order mass transfer between conduits was derived. The dual-advection dispersion equation was then applied to the multi-peaked breakthrough curves obtained from the physical experiments in order to obtain some insight into the operative solute-transport processes through the acquisition of a consequent set of velocities, dispersions, and related parameters. Successful application of the dual-advection, dispersion equation to a tracer test that exhibited dual peaks for a karst aquifer known to consist of two connected but mostly separate conduits confirmed the appropriateness of using a multi-dispersion type model when conditions warrant.
Knopman, Debra S.; Voss, Clifford I.
1987-01-01
The spatial and temporal variability of sensitivities has a significant impact on parameter estimation and sampling design for studies of solute transport in porous media. Physical insight into the behavior of sensitivities is offered through an analysis of analytically derived sensitivities for the one-dimensional form of the advection-dispersion equation. When parameters are estimated in regression models of one-dimensional transport, the spatial and temporal variability in sensitivities influences variance and covariance of parameter estimates. Several principles account for the observed influence of sensitivities on parameter uncertainty. (1) Information about a physical parameter may be most accurately gained at points in space and time. (2) As the distance of observation points from the upstream boundary increases, maximum sensitivity to velocity during passage of the solute front increases. (3) The frequency of sampling must be 'in phase' with the S shape of the dispersion sensitivity curve to yield the most information on dispersion. (4) The sensitivity to the dispersion coefficient is usually at least an order of magnitude less than the sensitivity to velocity. (5) The assumed probability distribution of random error in observations of solute concentration determines the form of the sensitivities. (6) If variance in random error in observations is large, trends in sensitivities of observation points may be obscured by noise. (7) Designs that minimize the variance of one parameter may not necessarily minimize the variance of other parameters.
Technology Transfer Automated Retrieval System (TEKTRAN)
Field tests were conducted to obtain irrigation evaluation and solute transport data that were used to calibrate and validate an advection-dispersion model for furrow irrigation. Empirical infiltration equation and roughness parameters were estimated from the field data. These estimates were used a...
Huang, Y.H.; Saiers, J.E.; Harvey, J.W.; Noe, G.B.; Mylon, S.
2008-01-01
The movement of particulate matter within wetland surface waters affects nutrient cycling, contaminant mobility, and the evolution of the wetland landscape. Despite the importance of particle transport in influencing wetland form and function, there are few data sets that illuminate, in a quantitative way, the transport behavior of particulate matter within surface waters containing emergent vegetation. We report observations from experiments on the transport of 1 ??m latex microspheres at a wetland field site located in Water Conservation Area 3A of the Florida Everglades. The experiments involved line source injections of particles inside two 4.8-m-long surface water flumes constructed within a transition zone between an Eleocharis slough and Cladium jamaicense ridge and within a Cladium jamaicense ridge. We compared the measurements of particle transport to calculations of two-dimensional advection-dispersion model that accounted for a linear increase in water velocities with elevation above the ground surface. The results of this analysis revealed that particle spreading by longitudinal and vertical dispersion was substantially greater in the ridge than within the transition zone and that particle capture by aquatic vegetation lowered surface water particle concentrations and, at least for the timescale of our experiments, could be represented as an irreversible, first-order kinetics process. We found generally good agreement between our field-based estimates of particle dispersion and water velocity and estimates determined from published theory, suggesting that the advective-dispersive transport of particulate matter within complex wetland environments can be approximated on the basis of measurable properties of the flow and aquatic vegetation. Copyright 2008 by the American Geophysical Union.
NASA Astrophysics Data System (ADS)
Bodin, Jacques
2015-03-01
In this study, new multi-dimensional time-domain random walk (TDRW) algorithms are derived from approximate one-dimensional (1-D), two-dimensional (2-D), and three-dimensional (3-D) analytical solutions of the advection-dispersion equation and from exact 1-D, 2-D, and 3-D analytical solutions of the pure-diffusion equation. These algorithms enable the calculation of both the time required for a particle to travel a specified distance in a homogeneous medium and the mass recovery at the observation point, which may be incomplete due to 2-D or 3-D transverse dispersion or diffusion. The method is extended to heterogeneous media, represented as a piecewise collection of homogeneous media. The particle motion is then decomposed along a series of intermediate checkpoints located on the medium interface boundaries. The accuracy of the multi-dimensional TDRW method is verified against (i) exact analytical solutions of solute transport in homogeneous media and (ii) finite-difference simulations in a synthetic 2-D heterogeneous medium of simple geometry. The results demonstrate that the method is ideally suited to purely diffusive transport and to advection-dispersion transport problems dominated by advection. Conversely, the method is not recommended for highly dispersive transport problems because the accuracy of the advection-dispersion TDRW algorithms degrades rapidly for a low Péclet number, consistent with the accuracy limit of the approximate analytical solutions. The proposed approach provides a unified methodology for deriving multi-dimensional time-domain particle equations and may be applicable to other mathematical transport models, provided that appropriate analytical solutions are available.
Purely Lagrangian Simulation of Advection, Dispersion, Precipitation, and Dissolution
NASA Astrophysics Data System (ADS)
Benson, D.; Zhang, Y.; Reeves, D. M.
2008-05-01
We extend the advantages of Lagrangian random walk particle tracking (RWPT) methods that have long been used to simulate advection and dispersion in highly heterogeneous media. By formulating dissolution as a random, independent decay process, the classical continuum rate law is recovered. Formulating the random precipitation process requires a consideration of the probability that two nearby particles will coincide in a given time period. This depends on local mixing (as by diffusion) and the total domain particle number density, which are fixed and therefore easy to calculate. The result is that the classical law of mass action for equilibrium reactions can be reproduced in an ensemble sense. The same number of parameters for A+B ⇌ C are needed in a probabilistic versus continuum reaction simulation-- —one each for forward and backward probabilities that correspond to rates. The random nature of the simulations allows for significant disequilibrium in any given region at any time that is independent of the numerical details such as time stepping or particle density. This is exemplified by nearby or intermingled groups of reactants and little or no product--—a result that is often noted in the field that is difficult to reconcile with continuum methods or coarse-grained Eulerian models. Our results support recent results of perturbed advection-dispersion-reaction continuum models (Luo et al., WRR 44, 2008), and suggest that many different kinds of reactions can be easily added to existing RWPT codes.
Reactive solute transport in acidic streams
Broshears, R.E.
1996-01-01
Spatial and temporal profiles of Ph and concentrations of toxic metals in streams affected by acid mine drainage are the result of the interplay of physical and biogeochemical processes. This paper describes a reactive solute transport model that provides a physically and thermodynamically quantitative interpretation of these profiles. The model combines a transport module that includes advection-dispersion and transient storage with a geochemical speciation module based on MINTEQA2. Input to the model includes stream hydrologic properties derived from tracer-dilution experiments, headwater and lateral inflow concentrations analyzed in field samples, and a thermodynamic database. Simulations reproduced the general features of steady-state patterns of observed pH and concentrations of aluminum and sulfate in St. Kevin Gulch, an acid mine drainage stream near Leadville, Colorado. These patterns were altered temporarily by injection of sodium carbonate into the stream. A transient simulation reproduced the observed effects of the base injection.
Simulation models for conservative and nonconservative solute transport in streams
Runkel, R.L.
1995-01-01
Solute transport in streams is governed by a suite of hydrologic and chemical processes. Interactions between hydrologic processes and chemical reactions may be quantified through a combination of field-scale experimentation and simulation modeling. Two mathematical models that simulate conservative and nonconservative solute transport in streams are presented. A model for conservative solutes that considers One Dimensional Transport with Inflow and Storage (OTIS) may be used in conjunction with tracer-dilution methods to quantify hydrologic transport processes (advection, dispersion, lateral inflow and transient storage). For nonconservative solutes, a model known as OTEQ may be used to quantify chemical processes within the context of hydrologic transport. OTEQ combines the transport mechanisms in OTIS with a chemical equilibrium sub-model that considers complexation, precipitation/dissolution and sorption. OTEQ has been used to quantify processes affecting trace metals in two streams in the Rocky Mountains of Colorado, USA.
Solute transport in dual-permeability porous media
NASA Astrophysics Data System (ADS)
Leij, Feike J.; Toride, Nobuo; Field, Malcolm S.; Sciortino, Antonella
2012-04-01
A dual-advection dispersion equation (DADE) is presented and solved to describe solute transport in structured or layered porous media with different nonzero flow rates in two distinct pore domains with linear solute transfer between them. This dual-permeability model constitutes a generalized version of the advection-dispersion equation (ADE) for transport in uniform porous media and the mobile-immobile model (MIM) for transport in media with a mobile and an immobile pore domain. Analytical tools for the DADE have mostly been lacking. An analytical solution has therefore been derived using Laplace transformation with time and modal decomposition based on matrix diagonalization, assuming the same dispersivity for both domains. Temporal moments are derived for the DADE and contrasted with those for the ADE and the MIM. The effective dispersion coefficient for the DADE approaches that of the ADE for a similar velocity in both pore domains and large values for the first-order transfer parameter, and approaches that of the MIM for the opposite conditions. The solution of the DADE is used to illustrate how differences in pore water velocity between the domains and low transfer rates will lead to double peaks in the volume- or flux-averaged concentration profiles versus time or position. The DADE is applied to optimize experimental breakthrough curves for an Andisol with a distinct intra- and interaggregate porosity. The DADE improved the description of the breakthrough data compared to the ADE and the MIM.
Probabilistic exposure risk assessment with advective-dispersive well vulnerability criteria
NASA Astrophysics Data System (ADS)
Enzenhoefer, Rainer; Nowak, Wolfgang; Helmig, Rainer
2012-02-01
Time-related advection-based well-head protection zones are commonly used to manage the contamination risk of drinking water wells. According to current water safety plans advanced risk management schemes are needed to better control and monitor all possible hazards within catchments. The goal of this work is to cast the four advective-dispersive intrinsic well vulnerability criteria by Frind et al. [1] into a framework of probabilistic risk assessment framework. These criteria are: (i) arrival time, (ii) level of peak concentration, (iii) time until first arrival of critical concentrations and (iv) exposure time. Our probabilistic framework yields catchment-wide maps of probabilities to not comply with these criteria. This provides indispensable information for catchment managers to perform probabilistic exposure risk assessment and thus improves the basis for risk-informed well-head management. We resolve heterogeneity with high-resolution Monte Carlo simulations and use a new reverse formulation of temporal moment transport equations to keep computational costs low. Our method is independent of dimensionality and boundary conditions, and can account for arbitrary sources of uncertainty. It can be coupled with any method for conditioning on available data. For simplicity, we demonstrate the concept on a 2D example that includes conditioning on synthetic data.
NASA Astrophysics Data System (ADS)
Bernabé, Y.; Wang, Y.; Qi, T.; Li, M.
2016-02-01
The main purpose of this work is to investigate the relationship between passive advection-dispersion and permeability in porous materials presumed to be statistically homogeneous at scales larger than the pore scale but smaller than the reservoir scale. We simulated fluid flow through pipe network realizations with different pipe radius distributions and different levels of connectivity. The flow simulations used periodic boundary conditions, allowing monitoring of the advective motion of solute particles in a large periodic array of identical network realizations. In order to simulate dispersion, we assumed that the solute particles obeyed Taylor dispersion in individual pipes. When a particle entered a pipe, a residence time consistent with local Taylor dispersion was randomly assigned to it. When exiting the pipe, the particle randomly proceeded into one of the pipes connected to the original one according to probabilities proportional to the outgoing volumetric flow in each pipe. For each simulation we tracked the motion of at least 6000 solute particles. The mean fluid velocity was 10-3 ms-1, and the distance traveled was on the order of 10 m. Macroscopic dispersion was quantified using the method of moments. Despite differences arising from using different types of lattices (simple cubic, body-centered cubic, and face-centered cubic), a number of general observations were made. Longitudinal dispersion was at least 1 order of magnitude greater than transverse dispersion, and both strongly increased with decreasing pore connectivity and/or pore size variability. In conditions of variable hydraulic radius and fixed pore connectivity and pore size variability, the simulated dispersivities increased as power laws of the hydraulic radius and, consequently, of permeability, in agreement with previously published experimental results. Based on these observations, we were able to resolve some of the complexity of the relationship between dispersivity and permeability.
Simulating water, solute, and heat transport in the subsurface with the VS2DI software package
Healy, R.W.
2008-01-01
The software package VS2DI was developed by the U.S. Geological Survey for simulating water, solute, and heat transport in variably saturated porous media. The package consists of a graphical preprocessor to facilitate construction of a simulation, a postprocessor for visualizing simulation results, and two numerical models that solve for flow and solute transport (VS2DT) and flow and heat transport (VS2DH). The finite-difference method is used to solve the Richards equation for flow and the advection-dispersion equation for solute or heat transport. This study presents a brief description of the VS2DI package, an overview of the various types of problems that have been addressed with the package, and an analysis of the advantages and limitations of the package. A review of other models and modeling approaches for studying water, solute, and heat transport also is provided. ?? Soil Science Society of America. All rights reserved.
Wexler, Eliezer J.
1989-01-01
Analytical solutions to the advective-dispersive solute transport equation are useful in predicting the fate of solutes in groundwater. Analytical solutions compiled from available literature or derived by the author are presented in this report for a variety of boundary condition types and solute-source configuration in one-, two-, and three-dimensional systems with uniform groundwater flow. A set of user-oriented computer programs was created to evaluate these solutions and to display the results in tabular and computer-graphics format. These programs incorporate many features that enhance their accuracy, ease of use, and versatility. Documentation for the programs describes their operation and required input data, and presents the results of sample problems. Derivations of select solutions, source codes for the computer programs, and samples of program input and output also are described. (USGS)
Wexler, Eliezer J.
1992-01-01
Analytical solutions to the advective-dispersive solute-transport equation are useful in predicting the fate of solutes in ground water. Analytical solutions compiled from available literature or derived by the author are presented for a variety of boundary condition types and solute-source configurations in one-, two-, and three-dimensional systems having uniform ground-water flow. A set of user-oriented computer programs was created to evaluate these solutions and to display the results in tabular and computer-graphics format. These programs incorporate many features that enhance their accuracy, ease of use, and versatility. Documentation for the programs describes their operation and required input data, and presents the results of sample problems. Derivations of selected solutions, source codes for the computer programs, and samples of program input and output also are included.
NASA Astrophysics Data System (ADS)
Cartwright, Ian
Advection-dispersion fluid flow models implicitly assume that the infiltrating fluid flows through an already fluid-saturated medium. However, whether rocks contain a fluid depends on their reaction history, and whether any initial fluid escapes. The behaviour of different rocks may be illustrated using hypothetical marble compositions. Marbles with diverse chemistries (e.g. calcite + dolomite + quartz) are relatively reactive, and will generally produce a fluid during heating. By contrast, marbles with more restricted chemistries (e.g. calcite + quartz or calcite-only) may not. If the rock is not fluid bearing when fluid infiltration commences, mineralogical reactions may produce a reaction-enhanced permeability in calcite + dolomite + quartz or calcite + quartz, but not in calcite-only marbles. The permeability production controls the pattern of mineralogical, isotopic, and geochemical resetting during fluid flow. Tracers retarded behind the mineralogical fronts will probably be reset as predicted by the advection-dispersion models; however, tracers that are expected to be reset ahead of the mineralogical fronts cannot progress beyond the permeability generating reaction. In the case of very unreactive lithologies (e.g. pure calcite marbles, cherts, and quartzites), the first reaction to affect the rocks may be a metasomatic one ahead of which there is little pervasive resetting of any tracer. Centimetre-scale layering may lead to the formation of self-perpetuating fluid channels in rocks that are not fluid saturated due to the juxtaposition of reactants. Such layered rocks may show patterns of mineralogical resetting that are not predicted by advection-dispersion models. Patterns of mineralogical and isotopic resetting in marbles from a number of terrains, for example: Chillagoe, Marulan South, Reynolds Range (Australia); Adirondack Mountains, Old Woman Mountains, Notch Peak (USA); and Stephen Cross Quarry (Canada) vary as predicted by these models.
Anwar, S.; Cortis, A.; Sukop, M.
2008-10-20
Lattice Boltzmann models simulate solute transport in porous media traversed by conduits. Resulting solute breakthrough curves are fitted with Continuous Time Random Walk models. Porous media are simulated by damping flow inertia and, when the damping is large enough, a Darcy's Law solution instead of the Navier-Stokes solution normally provided by the lattice Boltzmann model is obtained. Anisotropic dispersion is incorporated using a direction-dependent relaxation time. Our particular interest is to simulate transport processes outside the applicability of the standard Advection-Dispersion Equation (ADE) including eddy mixing in conduits. The ADE fails to adequately fit any of these breakthrough curves.
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.
Analytically-derived sensitivities in one-dimensional models of solute transport in porous media
Knopman, D.S.
1987-01-01
Analytically-derived sensitivities are presented for parameters in one-dimensional models of solute transport in porous media. Sensitivities were derived by direct differentiation of closed form solutions for each of the odel, and by a time integral method for two of the models. Models are based on the advection-dispersion equation and include adsorption and first-order chemical decay. Boundary conditions considered are: a constant step input of solute, constant flux input of solute, and exponentially decaying input of solute at the upstream boundary. A zero flux is assumed at the downstream boundary. Initial conditions include a constant and spatially varying distribution of solute. One model simulates the mixing of solute in an observation well from individual layers in a multilayer aquifer system. Computer programs produce output files compatible with graphics software in which sensitivities are plotted as a function of either time or space. (USGS)
NASA Astrophysics Data System (ADS)
Chen, Jui-Sheng; Jang, Cheng-Shin; Cheng, Chung-Ting; Liu, Chen-Wuing
2010-09-01
SummaryThis study presents a novel mathematical model for describing the transport of the remedial reagent in a vertical circulation flow field in an anisotropic aquifer. To develop the mathematical model, the radial and vertical components of the pore water velocity are calculated first by using an analytical solution for steady-state drawdown distribution near a vertical circulation well. Next, the obtained radial and vertical components of the pore water velocity are then incorporated into a three-dimensional axisymmetrical advection-dispersion equation in cylindrical coordinates from which to build the reagent transport equation. The Laplace transform finite difference technique is applied to solve the three-dimensional axisymmetrical advection-dispersion equation with spatial variable-dependent coefficients. The developed mathematical model is used to investigate the effects of various parameters such as hydraulic conductivity anisotropy, longitudinal and transverse dispersivities, the placement of the extraction and injection screened intervals of the vertical circulation well and the injection modes on the transport regime of the remedial reagent. Results show that those parameters have different degrees of impacts on the distribution of the remedial reagent. The mathematical model provides an effective tool for designing and operating an enhanced groundwater remediation in an anisotropic aquifer using the vertical circulation well technology.
Lancaster, Jill; Downes, Barbara J
2014-12-01
Many communities comprise species that select resources that are patchily distributed in an environment that is otherwise unsuitable or suboptimal. Effects of this patchiness can depend on the characteristics of patch arrays and animal movements, and produce non-intuitive outcomes in which population densities are unrelated to resource abundance. Resource mosaics are predicted to have only weak effects, however, where patches are ephemeral or organisms are transported advectively. The running waters of streams and benthic invertebrates epitomize such systems, but empirical tests of resource mosaics are scarce. We sampled 15 common macroinvertebrates inhabiting distinct detritus patches at four sites within a sand-bed stream, where detritus formed a major resource of food and living space. At each site, environmental variables were measured for 100 leaf packs; invertebrates were counted in 50 leaf packs. Sites differed in total abundance of detritus, leaf pack sizes and invertebrate densities. Multivariate analysis indicated that patch size was the dominant environmental variable, but invertebrate densities differed significantly between sites even after accounting for patch size. Leaf specialists showed positive and strong density-area relationships, except where the patch size range was small and patches were aggregated. In contrast, generalist species had weaker and variable responses to patch sizes. Population densities were not associated with total resource abundance, with the highest densities of leaf specialists in sites with the least detritus. Our results demonstrate that patchy resources can affect species even in communities where species are mobile, have advective dispersal, and patches are relatively ephemeral. PMID:25190216
NASA Astrophysics Data System (ADS)
Pang, Liping; Goltz, Mark; Close, Murray
2003-01-01
In this note, we applied the temporal moment solutions of [Das and Kluitenberg, 1996. Soil Sci. Am. J. 60, 1724] for one-dimensional advective-dispersive solute transport with linear equilibrium sorption and first-order degradation for time pulse sources to analyse soil column experimental data. Unlike most other moment solutions, these solutions consider the interplay of degradation and sorption. This permits estimation of a first-order degradation rate constant using the zeroth moment of column breakthrough data, as well as estimation of the retardation factor or sorption distribution coefficient of a degrading solute using the first moment. The method of temporal moment (MOM) formulae was applied to analyse breakthrough data from a laboratory column study of atrazine, hexazinone and rhodamine WT transport in volcanic pumice sand, as well as experimental data from the literature. Transport and degradation parameters obtained using the MOM were compared to parameters obtained by fitting breakthrough data from an advective-dispersive transport model with equilibrium sorption and first-order degradation, using the nonlinear least-square curve-fitting program CXTFIT. The results derived from using the literature data were also compared with estimates reported in the literature using different equilibrium models. The good agreement suggests that the MOM could provide an additional useful means of parameter estimation for transport involving equilibrium sorption and first-order degradation. We found that the MOM fitted breakthrough curves with tailing better than curve fitting. However, the MOM analysis requires complete breakthrough curves and relatively frequent data collection to ensure the accuracy of the moments obtained from the breakthrough data.
Pang, Liping; Goltz, Mark; Close, Murray
2003-01-01
In this note, we applied the temporal moment solutions of [Das and Kluitenberg, 1996. Soil Sci. Am. J. 60, 1724] for one-dimensional advective-dispersive solute transport with linear equilibrium sorption and first-order degradation for time pulse sources to analyse soil column experimental data. Unlike most other moment solutions, these solutions consider the interplay of degradation and sorption. This permits estimation of a first-order degradation rate constant using the zeroth moment of column breakthrough data, as well as estimation of the retardation factor or sorption distribution coefficient of a degrading solute using the first moment. The method of temporal moment (MOM) formulae was applied to analyse breakthrough data from a laboratory column study of atrazine, hexazinone and rhodamine WT transport in volcanic pumice sand, as well as experimental data from the literature. Transport and degradation parameters obtained using the MOM were compared to parameters obtained by fitting breakthrough data from an advective-dispersive transport model with equilibrium sorption and first-order degradation, using the nonlinear least-square curve-fitting program CXTFIT. The results derived from using the literature data were also compared with estimates reported in the literature using different equilibrium models. The good agreement suggests that the MOM could provide an additional useful means of parameter estimation for transport involving equilibrium sorption and first-order degradation. We found that the MOM fitted breakthrough curves with tailing better than curve fitting. However, the MOM analysis requires complete breakthrough curves and relatively frequent data collection to ensure the accuracy of the moments obtained from the breakthrough data. PMID:12498577
Bad behavior of Godunov mixed methods for strongly anisotropic advection-dispersion equations
NASA Astrophysics Data System (ADS)
Mazzia, Annamaria; Manzini, Gianmarco; Putti, Mario
2011-09-01
We study the performance of Godunov mixed methods, which combine a mixed-hybrid finite element solver and a Godunov-like shock-capturing solver, for the numerical treatment of the advection-dispersion equation with strong anisotropic tensor coefficients. It turns out that a mesh locking phenomenon may cause ill-conditioning and reduce the accuracy of the numerical approximation especially on coarse meshes. This problem may be partially alleviated by substituting the mixed-hybrid finite element solver used in the discretization of the dispersive (diffusive) term with a linear Galerkin finite element solver, which does not display such a strong ill conditioning. To illustrate the different mechanisms that come into play, we investigate the spectral properties of such numerical discretizations when applied to a strongly anisotropic diffusive term on a small regular mesh. A thorough comparison of the stiffness matrix eigenvalues reveals that the accuracy loss of the Godunov mixed method is a structural feature of the mixed-hybrid method. In fact, the varied response of the two methods is due to the different way the smallest and largest eigenvalues of the dispersion (diffusion) tensor influence the diagonal and off-diagonal terms of the final stiffness matrix. One and two dimensional test cases support our findings.
Gureghian, A.B.
1990-08-01
Analytical solutions based on the Laplace transforms are presented for the one-dimensional, transient, advective-dispersive transport of a reacting radionuclide through a discrete planar fracture with constant aperture subject to diffusion in the surrounding rock matrix where both regions of solute migration display residual concentrations. The dispersion-free solutions, which are of closed form, are also reported. The solution assumes that the ground-water flow regime is under steady-state and isothermal conditions and that the rock matrix is homogeneous, isotropic, and saturated with stagnant water. The verification of the solution was performed by means of related analytical solutions dealing with particular aspects of the transport problem under investigation on the one hand, and a numerical solution capable of handling the complete problem on the other. The integrals encountered in the general solution are evaluated by means of a composite Gauss-Legendre quadrature scheme. 9 refs., 8 figs., 32 tabs.
Scott, D.T.; Gooseff, M.N.; Bencala, K.E.; Runkel, R.L.
2003-01-01
The hydrologic processes of advection, dispersion, and transient storage are the primary physical mechanisms affecting solute transport in streams. The estimation of parameters for a conservative solute transport model is an essential step to characterize transient storage and other physical features that cannot be directly measured, and often is a preliminary step in the study of reactive solutes. Our study used inverse modeling to estimate parameters of the transient storage model OTIS (One dimensional Transport with Inflow and Storage). Observations from a tracer injection experiment performed on Uvas Creek, California, USA, are used to illustrate the application of automated solute transport model calibration to conservative and nonconservative stream solute transport. A computer code for universal inverse modeling (UCODE) is used for the calibrations. Results of this procedure are compared with a previous study that used a trial-and-error parameter estimation approach. The results demonstrated 1) importance of the proper estimation of discharge and lateral inflow within the stream system; 2) that although the fit of the observations is not much better when transient storage is invoked, a more randomly distributed set of residuals resulted (suggesting non-systematic error), indicating that transient storage is occurring; 3) that inclusion of transient storage for a reactive solute (Sr2+) provided a better fit to the observations, highlighting the importance of robust model parameterization; and 4) that applying an automated calibration inverse modeling estimation approach resulted in a comprehensive understanding of the model results and the limitation of input data.
NONUNIFORM AND UNSTEADY SOLUTE TRANSPORT IN FURROW IRRRIGATION: I. MODEL DEVELOPMENT
Technology Transfer Automated Retrieval System (TEKTRAN)
A model for solving a cross-section-averaged Advection-Dispersion Equation (ADE) was developed to simulate the transport of fertilizer in furrow irrigation. The advection and dispersion processes were solved separately at each time step by implementing a method of characteristics with cubic spline i...
Reactive solute transport in streams. 1. Development of an equilibrium- based model
Runkel, R.L.; Bencala, K.E.; Broshears, R.E.; Chapra, S.C.
1996-01-01
An equilibrium-based solute transport model is developed for the simulation of trace metal fate and transport in streams. The model is formed by coupling a solute transport model with a chemical equilibrium submodel based on MINTEQ. The solute transport model considers the physical processes of advection, dispersion, lateral inflow, and transient storage, while the equilibrium submodel considers the speciation and complexation of aqueous species, precipitation/dissolution and sorption. Within the model, reactions in the water column may result in the formation of solid phases (precipitates and sorbed species) that are subject to downstream transport and settling processes. Solid phases on the streambed may also interact with the water column through dissolution and sorption/desorption reactions. Consideration of both mobile (water-borne) and immobile (streambed) solid phases requires a unique set of governing differential equations and solution techniques that are developed herein. The partial differential equations describing physical transport and the algebraic equations describing chemical equilibria are coupled using the sequential iteration approach.
NASA Astrophysics Data System (ADS)
Kelkar, Sharad; Ding, Mei; Chu, Shaoping; Robinson, Bruce A.; Arnold, Bill; Meijer, Arend; Eddebbarh, Al-Aziz
2010-09-01
This paper presents a study of solute transport through ground water in the saturated zone and the resulting breakthrough curves (BTCs), using a field-scale numerical model that incorporates the processes of advection, dispersion, matrix diffusion in fractured volcanic formations, sorption, and colloid-facilitated transport. Such BTCs at compliance boundaries are often used as performance measures for a site. The example considered here is that of the saturated zone study prepared for the Yucca Mountain license application. The saturated zone at this site occurs partly in volcanic, fractured rock formations and partly in alluvial formations. This paper presents a description of the site and the ground water flow model, the development of the conceptual model of transport, model uncertainties, model validation, and the influence of uncertainty in input parameters on the downstream BTCs at the Yucca Mountain site.
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.
Dynamic typology of hydrothermal systems: competing effects of advection, dispersion and reactivity
NASA Astrophysics Data System (ADS)
Dolejs, David
2016-04-01
Genetic interpretation hydrothermal systems relies on recognition of (i) hydrothermal fluid source, (ii) fluid migration pathways, and (iii) deposition site identified by hydrothermal alteration and/or mineralization. Frequently, only the last object is of interest or accessible to direct observation, but constraints on the fluid source (volume) and pathways can be obtained from evaluation of the time-integrated fluid flux during hydrothermal event. Successful interpretation of the petrological record, that is, progress of alteration reactions, relies on identification of individual contributions arising from solute advection (to the deposition site), its lateral dispersion, and reaction efficiency. Although these terms are all applicable in a mass-conservation relationship within the framework of the transport theory, they are rarely considered simultaneously and their relative magnitudes evaluated. These phenomena operate on variable length and time scales, and may in turn provide insight into the system dynamics such as flow, diffusion and reaction rates, or continuous vs. episodic behavior of hydrothermal events. In addition, here we demonstrate that they also affect estimate of the net fluid flux, frequently by several orders of magnitude. The extent of alteration and mineralization reactions between the hydrothermal fluid and the host environment is determined by: (i) temperature, pressure or any other gradients across the mineralization site, (ii) magnitude of disequilibrium at inflow to the mineralization site, which is related to physico-chemical gradient between the fluid source and the mineralization site, and (iii) chemical redistribution (dispersion) within the mineralization site. We introduce quantitative mass-transport descriptors - Péclet and Damköhler II numbers - to introduce division into dispersion-dominated, advection-dominated and reaction-constrained systems. Dispersive systems are characterized by lateral solute redistribution, driven by
NASA Astrophysics Data System (ADS)
Bouchelaghem, F.; Vulliet, L.
2001-10-01
The development of a predictive model of behaviour of porous media during injection of miscible grout, taking into account convection, dilution and filtration of grout solution with interstitial water, as well as consolidation aspects, is presented. Model assumptions are reviewed and discussed first. During the establishment of the model, we insist on surface terms and their physical relevance in expressing adsorption effects. Constitutive laws such as Fick's law for diffusive mass transport, hydrodynamic dispersion tensor dealing with miscibility, are modified by taking into account filtration effects. A new surface term appears in mass balance equations as a consequence of filtration. According to the filtration laws used, an initial filtration rate is estimated on the basis of a one-dimensional experimental campaign. The field equations are discretized by using Galerkin finite element and -scheme standard method. For transport equation, Streamline Upwind Petrov Galerkin method is employed to prevent numerical oscillations. Lastly, confrontation of numerical results with laboratory experiments constitutes a first step to validate the model on a realistic basis.
NASA Astrophysics Data System (ADS)
Zaramella, M.; Marion, A.; Lewandowski, J.; Nützmann, G.
2016-07-01
Solute transport in rivers is controlled by surface flow hydrodynamics and by transient storage in dead zones, pockets of vegetation and hyporheic sediments where mass exchange and retention are governed by complex mechanisms. The physics of these processes are generally investigated by optimization of transient storage models (TSMs) to experimental data often yielding inconsistent and equifinal parameter sets. Uncertainty on parameters estimation is found to depend not only on the rates of exchange between the stream and storage zones, the stream-water velocity and the stream reach length according to the experimental Damkohler number (DaI), but also on the relative significance between transient storage and longitudinal dispersion on breakthrough curves (BTCs). An optimization strategy was developed and applied to an experimental dataset obtained from tracer tests in a small lowland river, analyzing BTCs generated through tracer injections under different conditions. The method supplies a tool to estimate model parameters from observed data through the analysis of the relative parameter significance. To analyze model performance a double compartment TSM was optimized by a regular fit procedure based on simple root mean square error minimization and by a fit based on a relative significance analysis of mechanism signatures. As a result consistent longitudinal dispersion and transient storage parameters were obtained when the signature targeted optimization was used.
NASA Astrophysics Data System (ADS)
Swanson, R. D.; Binley, A. M.; Keating, K.; France, S.; Osterman, G. K.; Day-Lewis, F. D.; Singha, K.
2013-12-01
The advection-dispersion equation fails to describe non-Fickian solute transport in saturated porous media, necessitating the use of other models. The dual-domain mass transfer (DDMT) model partitions the total porosity into mobile and less-mobile domains with solute exchange between the domains; consequently, the DDMT model can produce a better fit to breakthrough curves (BTCs) in systems defined by more- and less-mobile components. However, direct experimental estimation of DDMT model parameters such as rate of exchange and the mobile and less-mobile porosities remains elusive. Consequently, model parameters are often calculated purely as a model fitting exercise. There is a clear need for material characterization techniques that can offer some insight into the pore space geometrical arrangement, particularly if such techniques can be extended to the field scale. Here, we interpret static direct-current (DC) resistivity, complex resistivity (CR) and nuclear magnetic resonance (NMR) geophysical measurements in the characterization of mass transfer parameters. We use two different samples of the zeolite clinoptilolite, a material shown to demonstrate solute mass transfer due to a significant intragranular porosity, along with glass beads as a control. We explore the relation between geophysical and DDMT parameters in conjunction with supporting material characterization methods. Our results reveal how these geophysical measurements can offer some insight into the pore structures controlling the observed anomalous transport behavior.
The Mechanism of Field-Scale Solute Transport: An insight from Numerical Simulations
NASA Astrophysics Data System (ADS)
Russo, David
2014-05-01
Field-scale transport of conservative (chloride) and reactive (nitrate) solutes was analyzed by means of two different model processes for the local description of the transport. The first is the classical, one-region advection dispersion equation (ADE) model, while the second is the two-region, mobile-immobile (MIM) model. The analyses were performed by means of detailed three-dimensional (3-D), numerical simulations of the flow and the transport considering realistic features of the soil-water-plant-atmosphere system, pertinent to a turf field located in the Glil Yam site, Israel, irrigated with treated waste water (TWW). Simulated water content and concentration profiles were compared with available measurements of their counterparts. Results of the analyses suggest that the behavior of both the conservative and the reactive solutes in the Glil Yam site is quantified better when the transport on the local scale is modeled as a two-region, MIM model, than when a single-region, ADE model is used. Reconstruction of the shape of the measured solute concentration profiles using the MIM transport model, required relatively large immobile water content fraction and relatively small mass transfer coefficient. These results suggest that in the case of initially non-zero solute concentration profile (e.g., chloride and nitrate), the 3-D ADE transport model may significantly overestimate the groundwater contamination hazard posed by the solutes moving through the vadose zone, as compared with the 3-D MIM transport model, while the opposite is true in the case of initially zero solute concentration profile (e.g., carbamazepine). These findings stem from the combination of relatively large immobile water content fraction and relatively small mass transfer coefficient taken into account in the MIM transport model. In the first case, this combination forces a considerable portion of the solute mass to remain in the immobile region of the water-filled pores, while the opposite
NASA Astrophysics Data System (ADS)
Geiger, S.; Cortis, A.; Birkholzer, J. T.
2010-12-01
Solute transport in fractured porous media is typically "non-Fickian"; that is, it is characterized by early breakthrough and long tailing and by nonlinear growth of the Green function-centered second moment. This behavior is due to the effects of (1) multirate diffusion occurring between the highly permeable fracture network and the low-permeability rock matrix, (2) a wide range of advection rates in the fractures and, possibly, the matrix as well, and (3) a range of path lengths. As a consequence, prediction of solute transport processes at the macroscale represents a formidable challenge. Classical dual-porosity (or mobile-immobile) approaches in conjunction with an advection-dispersion equation and macroscopic dispersivity commonly fail to predict breakthrough of fractured porous media accurately. It was recently demonstrated that the continuous time random walk (CTRW) method can be used as a generalized upscaling approach. Here we extend this work and use results from high-resolution finite element-finite volume-based simulations of solute transport in an outcrop analogue of a naturally fractured reservoir to calibrate the CTRW method by extracting a distribution of retention times. This procedure allows us to predict breakthrough at other model locations accurately and to gain significant insight into the nature of the fracture-matrix interaction in naturally fractured porous reservoirs with geologically realistic fracture geometries.
Geiger, S.; Cortis, A.; Birkholzer, J.T.
2010-04-01
Solute transport in fractured porous media is typically 'non-Fickian'; that is, it is characterized by early breakthrough and long tailing and by nonlinear growth of the Green function-centered second moment. This behavior is due to the effects of (1) multirate diffusion occurring between the highly permeable fracture network and the low-permeability rock matrix, (2) a wide range of advection rates in the fractures and, possibly, the matrix as well, and (3) a range of path lengths. As a consequence, prediction of solute transport processes at the macroscale represents a formidable challenge. Classical dual-porosity (or mobile-immobile) approaches in conjunction with an advection-dispersion equation and macroscopic dispersivity commonly fail to predict breakthrough of fractured porous media accurately. It was recently demonstrated that the continuous time random walk (CTRW) method can be used as a generalized upscaling approach. Here we extend this work and use results from high-resolution finite element-finite volume-based simulations of solute transport in an outcrop analogue of a naturally fractured reservoir to calibrate the CTRW method by extracting a distribution of retention times. This procedure allows us to predict breakthrough at other model locations accurately and to gain significant insight into the nature of the fracture-matrix interaction in naturally fractured porous reservoirs with geologically realistic fracture geometries.
Weissmann, Gary S
2013-12-06
The objective of this project was to characterize the influence that naturally complex geologic media has on anomalous dispersion and to determine if the nature of dispersion can be estimated from the underlying heterogeneous media. The UNM portion of this project was to provide detailed representations of aquifer heterogeneity through producing highly-resolved models of outcrop analogs to aquifer materials. This project combined outcrop-scale heterogeneity characterization (conducted at the University of New Mexico), laboratory experiments (conducted at Sandia National Laboratory), and numerical simulations (conducted at Sandia National Laboratory and Colorado School of Mines). The study was designed to test whether established dispersion theory accurately predicts the behavior of solute transport through heterogeneous media and to investigate the relationship between heterogeneity and the parameters that populate these models. The dispersion theory tested by this work was based upon the fractional advection-dispersion equation (fADE) model. Unlike most dispersion studies that develop a solute transport model by fitting the solute transport breakthrough curve, this project explored the nature of the heterogeneous media to better understand the connection between the model parameters and the aquifer heterogeneity. We also evaluated methods for simulating the heterogeneity to see whether these approaches (e.g., geostatistical) could reasonably replicate realistic heterogeneity. The UNM portion of this study focused on capturing realistic geologic heterogeneity of aquifer analogs using advanced outcrop mapping methods.
CFEST Coupled Flow, Energy & Solute Transport Version CFEST005 Theory Guide
Freedman, Vicky L.; Chen, Yousu; Gupta, Sumant K.
2005-11-01
This document presents the mathematical theory implemented in the CFEST (Coupled Flow, Energy, and Solute Transport) simulator. The simulator is a three-dimensional finite element model that can be used for evaluating flow and solute mass transport. Although the theory for thermal transport is presented in this guide, it has not yet been fully implemented in the simulator. The flow module is capable of simulating both confined and unconfined aquifer systems, as well as constant and variable density fluid flows. For unconfined aquifers, the model uses a moving boundary for the water table, deforming the numerical mesh so that the uppermost nodes are always at the water table. For solute transport, changes in concentration of a single dissolved chemical constituent are computed for advective and hydrodynamic transport, linear sorption represented by a retardation factor, and radioactive decay. Once fully implemented, transport of thermal energy in the groundwater and solid matrix of the aquifer can also be used to model aquifer thermal regimes. Mesh construction employs “collapsible”, hexahedral finite elements in a three-dimensional coordinate system. CFEST uses the Galerkin finite element method to convert the partial differential equations to algebraic form. To solve the coupled equations for momentum, solute and heat transport, either Picard or Newton-Raphson iterative schemes are used to treat nonlinearities. An upstream weighted residual finite-element method is used to solve the advective-dispersive transport and energy transfer equations, which circumvents problems of numerical oscillation problems. Matrix solutions of the flow and transport problems are performed using efficient iterative solvers available in ITPACK and PETSc, solvers that are available in the public domain. These solvers are based on the preconditioned conjugate gradient and ORTHOMIN methods for symmetric and a nonsymmetric matrices, respectively.
NASA Astrophysics Data System (ADS)
Cornaton, F.; Perrochet, P.
2006-09-01
Groundwater age and life expectancy probability density functions (pdf) have been defined, and solved in a general three-dimensional context by means of forward and backward advection-dispersion equations [Cornaton F, Perrochet P. Groundwater age, life expectancy and transit time distributions in advective-dispersive systems; 1. Generalized reservoir theory. Adv Water Res (xxxx)]. The discharge and recharge zones transit time pdfs were then derived by applying the reservoir theory (RT) to the global system, thus considering as ensemble the union of all inlet boundaries on one hand, and the union of all outlet boundaries on the other hand. The main advantages in using the RT to calculate the transit time pdf is that the outlet boundary geometry does not represent a computational limiting factor (e.g. outlets of small sizes), since the methodology is based on the integration over the entire domain of each age, or life expectancy, occurrence. In the present paper, we extend the applicability of the RT to sub-drainage basins of groundwater reservoirs by treating the reservoir flow systems as compartments which transfer the water fluxes to a particular discharge zone, and inside which mixing and dispersion processes can take place. Drainage basins are defined by the field of probability of exit at outlet. In this way, we make the RT applicable to each sub-drainage system of an aquifer of arbitrary complexity and configuration. The case of the well-head protection problem is taken as illustrative example, and sensitivity analysis of the effect of pore velocity variations on the simulated ages is carried out.
Mehl, S.; Hill, M.C.
2001-01-01
Five common numerical techniques for solving the advection-dispersion equation (finite difference, predictor corrector, total variation diminishing, method of characteristics, and modified method of characteristics) were tested using simulations of a controlled conservative tracer-test experiment through a heterogeneous, two-dimensional sand tank. The experimental facility was constructed using discrete, randomly distributed, homogeneous blocks of five sand types. This experimental model provides an opportunity to compare the solution techniques: the heterogeneous hydraulic-conductivity distribution of known structure can be accurately represented by a numerical model, and detailed measurements can be compared with simulated concentrations and total flow through the tank. The present work uses this opportunity to investigate how three common types of results - simulated breakthrough curves, sensitivity analysis, and calibrated parameter values - change in this heterogeneous situation given the different methods of simulating solute transport. The breakthrough curves show that simulated peak concentrations, even at very fine grid spacings, varied between the techniques because of different amounts of numerical dispersion. Sensitivity-analysis results revealed: (1) a high correlation between hydraulic conductivity and porosity given the concentration and flow observations used, so that both could not be estimated; and (2) that the breakthrough curve data did not provide enough information to estimate individual values of dispersivity for the five sands. This study demonstrates that the choice of assigned dispersivity and the amount of numerical dispersion present in the solution technique influence estimated hydraulic conductivity values to a surprising degree.
Impact of space-time mesh adaptation on solute transport modeling in porous media
NASA Astrophysics Data System (ADS)
Esfandiar, Bahman; Porta, Giovanni; Perotto, Simona; Guadagnini, Alberto
2015-02-01
We implement a space-time grid adaptation procedure to efficiently improve the accuracy of numerical simulations of solute transport in porous media in the context of model parameter estimation. We focus on the Advection Dispersion Equation (ADE) for the interpretation of nonreactive transport experiments in laboratory-scale heterogeneous porous media. When compared to a numerical approximation based on a fixed space-time discretization, our approach is grounded on a joint automatic selection of the spatial grid and the time step to capture the main (space-time) system dynamics. Spatial mesh adaptation is driven by an anisotropic recovery-based error estimator which enables us to properly select the size, shape, and orientation of the mesh elements. Adaptation of the time step is performed through an ad hoc local reconstruction of the temporal derivative of the solution via a recovery-based approach. The impact of the proposed adaptation strategy on the ability to provide reliable estimates of the key parameters of an ADE model is assessed on the basis of experimental solute breakthrough data measured following tracer injection in a nonuniform porous system. Model calibration is performed in a Maximum Likelihood (ML) framework upon relying on the representation of the ADE solution through a generalized Polynomial Chaos Expansion (gPCE). Our results show that the proposed anisotropic space-time grid adaptation leads to ML parameter estimates and to model results of markedly improved quality when compared to classical inversion approaches based on a uniform space-time discretization.
NASA Astrophysics Data System (ADS)
Botter, Gianluca; Benettin, Paolo; Rinaldo, Andrea
2015-04-01
Advection-dispersion equations have been extensively used to model flow and transport processes through heterogeneous media like hillslopes and groundwater systems. Therein, the spreading of solute plumes and the shape of the breakthrough curve is known to be controlled by the macrodispersion coefficient, which embeds the underlying heterogeneity of velocities and flowpaths. On a nearly parallel track, the use of travel time distributions (TTDs) has become increasingly widespread in catchment hydrology, to establish a formal linkage between input and output chemographs through suitable transfer functions. Recent theoretical advances and real-world applications have shown that the structure of travel time distributions in time variable flow systems like watersheds is strongly related to the time variability of the water storage and input/output fluxes. The dynamical structure of TTDs has been proved to be effectively parametrized through suitable StorAge Selection (SAS) functions, that express in a spatially integrated fashion how the set of ages available within a control volume are selected and removed by the output fluxes. In this contribution, we analyze the relationship between Advection-Dispersion Models and StorAge Selection Functions, with examples for one-dimensional transport in a finite domain with constant convection and dispersion coefficient. Our results show that when the dispersion is high (say, Pe < 10), the distribution of ages leaving the system through the control plane is similar to the distribution of ages available within the storage, thereby leading to uniform SAS functions (random sampling). Implications for the interpretation and the prediction of the chemical response of rivers are discussed through the application of the SAS functions to model solute circulation in highly monitored watersheds belonging to diverse regions of the world. We suggest that the use of Storage Selection functions in different fields of hydrology may bring
NASA Astrophysics Data System (ADS)
Furbish, David Jon; Childs, Elise M.; Haff, Peter K.; Schmeeckle, Mark W.
2009-09-01
We formulate soil grain transport by rain splash as a stochastic advection-dispersion process. By taking into account the intermittency of grain motions activated by raindrop impacts, the formulation indicates that gradients in raindrop intensity, and thus grain activity (the volume of grains in motion per unit area) can be as important as gradients in grain concentration and surface slope in effecting transport. This idea is confirmed by rain splash experiments and manifest in topographic roughening via mound growth beneath desert shrubs. The formulation provides a framework for describing transport and dispersal of any soil material moveable by rain splash, including soil grains, soil-borne pathogens and nutrients, seeds, or debitage. As such it shows how classic models of topographic "diffusion" reflect effects of slope-dependent grain drift, not diffusion, and it highlights the role of rain splash in the ecological behavior of desert shrubs as "resource islands." Specifically, the growth of mounds beneath shrub canopies, where differential rain splash initially causes more grains to be splashed inward beneath the protective canopy than outward, involves the "harvesting" of nearby soil material, including nutrients. Mounds thus represent temporary storage of soil derived from areas surrounding the shrubs. As the inward grain flux associated with differential rain splash is sustained over the shrub lifetime, mound material is effectively sequestered from erosional processes that might otherwise move this material downslope. With shrub death and loss of the protective canopy, differential rain splash vanishes and the mound material is dispersed to the surrounding area, again subject to downslope movement.
Webb, S.W.
1996-05-01
Two models for gas-phase diffusion and advection in porous media, the Advective-Dispersive Model (ADM) and the Dusty-Gas Model (DGM), are reviewed. The ADM, which is more widely used, is based on a linear addition of advection calculated by Darcy`s Law and ordinary diffusion using Fick`s Law. Knudsen diffusion is often included through the use of a Klinkenberg factor for advection, while the effect of a porous medium on the diffusion process is through a porosity-tortuosity-gas saturation multiplier. Another, more comprehensive approach for gas-phase transport in porous media has been formulated by Evans and Mason, and is referred to as the Dusty- Gas Model (DGM). This model applies the kinetic theory of gases to the gaseous components and the porous media (or ``dust``) to develop an approach for combined transport due to ordinary and Knudsen diffusion and advection including porous medium effects. While these two models both consider advection and diffusion, the formulations are considerably different, especially for ordinary diffusion. The various components of flow (advection and diffusion) are compared for both models. Results from these two models are compared to isothermal experimental data for He-Ar gas diffusion in a low-permeability graphite. Air-water vapor comparisons have also been performed, although data are not available, for the low-permeability graphite system used for the helium-argon data. Radial and linear air-water heat pipes involving heat, advection, capillary transport, and diffusion under nonisothermal conditions have also been considered.
NASA Astrophysics Data System (ADS)
Cardenas, M. Bayani
2009-12-01
The transition from non-Fickian to Fickian macroscale transport is explicitly demonstrated for an increasing array of three-dimensional pores with vortices in between a lattice of cubic packed spheres by microscale finite element Navier-Stokes flow and transport simulations. Solute residence time distribution begins with a power law for one pore but gradually and eventually transforms to an exponential distribution typical of classic dispersive transport after about ten pores. Parameter fitting of an analytical solution to the 1-D advection-dispersion equation using the simulated breakthrough curves leads to fitted pore velocities within 1% of actual values and an asymptotic fitted dispersion coefficient after a few pores. Therefore, after dozens of pores, bulk transport can be described by the advection-dispersion equation. Persistent vortices in similarly structured porous media subjected to similar grain-scale Reynolds and Peclet numbers may have minimal contribution to anomalous transport observed at larger scales.
NASA Astrophysics Data System (ADS)
Gao, Guangyao; Zhan, Hongbin; Feng, Shaoyuan; Huang, Guanhua; Mao, Xiaomin
2009-10-01
SummaryThis study compared five different models for evaluating solute transport in a 1250-cm long, saturated and highly heterogeneous soil column. The five models were: the convection-dispersion equation (CDE), the mobile-immobile model (MIM), the convective lognormal transfer function model (CLT), the spatial fractional advection-dispersion equation (FADE) and the continuous time random walk model (CTRW). Each of these models was used to fit the breakthrough curve (BTC) at each distance individually and was also used to fit the BTCs at different distances simultaneously. Dependence of estimated parameters on distance was investigated. The estimated parameters at 200 cm were used to make predictions at subsequent distances. Highly anomalous transport behavior was observed in the column as the BTCs demonstrated significantly irregular shape and long tailing. This study indicated that CDE, CLT and FADE were unable to describe the anomalous BTCs adequately and their parameters changed with transport distance significantly. Compared to CDE, CLT and FADE, MIM better captured the evolution of anomalous BTCs. However, MIM did not explain the distinct BTC tailing satisfactorily. In contrast to MIM, CTRW better simulated the long tails of BTCs. The spreading parameter ( β) of CTRW was close to one and remained approximately constant at different travel distances. To make the comparison of these five models more general beyond the specific transport condition in the soil column, a generic evaluation of the advantages and disadvantages of these five models was presented in terms of their theory framework and a priori knowledge of the model behaviors.
Embry, Irucka; Roland, Victor; Agbaje, Oluropo; Watson, Valetta; Martin, Marquan; Painter, Roger; Byl, Tom; Sharpe, Lonnie
2013-01-01
A new residence-time distribution (RTD) function has been developed and applied to quantitative dye studies as an alternative to the traditional advection-dispersion equation (AdDE). The new method is based on a jointly combined four-parameter gamma probability density function (PDF). The gamma residence-time distribution (RTD) function and its first and second moments are derived from the individual two-parameter gamma distributions of randomly distributed variables, tracer travel distance, and linear velocity, which are based on their relationship with time. The gamma RTD function was used on a steady-state, nonideal system modeled as a plug-flow reactor (PFR) in the laboratory to validate themore » effectiveness of the model. The normalized forms of the gamma RTD and the advection-dispersion equation RTD were compared with the normalized tracer RTD. The normalized gamma RTD had a lower mean-absolute deviation (MAD) (0.16) than the normalized form of the advection-dispersion equation (0.26) when compared to the normalized tracer RTD. The gamma RTD function is tied back to the actual physical site due to its randomly distributed variables. The results validate using the gamma RTD as a suitable alternative to the advection-dispersion equation for quantitative tracer studies of non-ideal flow systems.« less
Modeling variably saturated subsurface solute transport with MODFLOW-UZF and MT3DMS
Morway, Eric D.; Niswonger, Richard G.; Langevin, Christian D.; Bailey, Ryan T.; Healy, Richard W.
2013-01-01
The MT3DMS groundwater solute transport model was modified to simulate solute transport in the unsaturated zone by incorporating the unsaturated-zone flow (UZF1) package developed for MODFLOW. The modified MT3DMS code uses a volume-averaged approach in which Lagrangian-based UZF1 fluid fluxes and storage changes are mapped onto a fixed grid. Referred to as UZF-MT3DMS, the linked model was tested against published benchmarks solved analytically as well as against other published codes, most frequently the U.S. Geological Survey's Variably-Saturated Two-Dimensional Flow and Transport Model. Results from a suite of test cases demonstrate that the modified code accurately simulates solute advection, dispersion, and reaction in the unsaturated zone. Two- and three-dimensional simulations also were investigated to ensure unsaturated-saturated zone interaction was simulated correctly. Because the UZF1 solution is analytical, large-scale flow and transport investigations can be performed free from the computational and data burdens required by numerical solutions to Richards' equation. Results demonstrate that significant simulation runtime savings can be achieved with UZF-MT3DMS, an important development when hundreds or thousands of model runs are required during parameter estimation and uncertainty analysis. Three-dimensional variably saturated flow and transport simulations revealed UZF-MT3DMS to have runtimes that are less than one tenth of the time required by models that rely on Richards' equation. Given its accuracy and efficiency, and the wide-spread use of both MODFLOW and MT3DMS, the added capability of unsaturated-zone transport in this familiar modeling framework stands to benefit a broad user-ship.
Modeling variably saturated subsurface solute transport with MODFLOW-UZF and MT3DMS.
Morway, Eric D; Niswonger, Richard G; Langevin, Christian D; Bailey, Ryan T; Healy, Richard W
2013-03-01
The MT3DMS groundwater solute transport model was modified to simulate solute transport in the unsaturated zone by incorporating the unsaturated-zone flow (UZF1) package developed for MODFLOW. The modified MT3DMS code uses a volume-averaged approach in which Lagrangian-based UZF1 fluid fluxes and storage changes are mapped onto a fixed grid. Referred to as UZF-MT3DMS, the linked model was tested against published benchmarks solved analytically as well as against other published codes, most frequently the U.S. Geological Survey's Variably-Saturated Two-Dimensional Flow and Transport Model. Results from a suite of test cases demonstrate that the modified code accurately simulates solute advection, dispersion, and reaction in the unsaturated zone. Two- and three-dimensional simulations also were investigated to ensure unsaturated-saturated zone interaction was simulated correctly. Because the UZF1 solution is analytical, large-scale flow and transport investigations can be performed free from the computational and data burdens required by numerical solutions to Richards' equation. Results demonstrate that significant simulation runtime savings can be achieved with UZF-MT3DMS, an important development when hundreds or thousands of model runs are required during parameter estimation and uncertainty analysis. Three-dimensional variably saturated flow and transport simulations revealed UZF-MT3DMS to have runtimes that are less than one tenth of the time required by models that rely on Richards' equation. Given its accuracy and efficiency, and the wide-spread use of both MODFLOW and MT3DMS, the added capability of unsaturated-zone transport in this familiar modeling framework stands to benefit a broad user-ship. PMID:22834908
Influence of karst evolution on solute transport evaluated by process-based numerical modelling
NASA Astrophysics Data System (ADS)
Hubinger, Bernhard; Birk, Steffen
2010-05-01
Karst waters are of major interest in water resources management. Because of their inherent properties karst systems show great vulnerability with regard to contaminants. Karst systems include highly permeable solution conduit networks formed by chemical aggressive water embedded in a fissured matrix. Small initial voids are widened and thus act as preferential passages, where flow is rapid and often turbulent. Water discharging at karst spring originates from different pathways with different residence times. Contaminant transport through conduit pathways is very rapid, whereas flow through the fissured porous matrix is much slower. Thus, on the one hand, pollutants may be rapidly transported and reach high concentrations at the karst spring shortly after their release; on the other hand, the existence of slow flow components may cause the pollution to last for long times. In this work, solute transport properties of karst aquifers are investigated using generic conduit networks of hydraulically connected proto-conduits with initially log-normally distributed apertures in the millimetre range and below. Conduit evolution is modelled by coupling flow, transport, and dissolution processes, whereby single conduits are widened up to the metre range. Thus, different stages of karst evolution can be distinguished. The resulting flow systems provide the basis for modelling advective-dispersive transport of non-reactive solutes through the network of more or less widened (proto-)conduits. The general transport characteristics in karst systems as well as the influence of heterogeneities and structures on solute transport are illustrated for cases of direct injection into the conduit systems at different evolutionary stages. The resulting breakthrough curves typically show several distinct, chronologically shifted peaks with long tailings, which appears to be similar to data from field tracer experiments.
NASA Astrophysics Data System (ADS)
Swanson, Ryan D.; Binley, Andrew; Keating, Kristina; France, Samantha; Osterman, Gordon; Day-Lewis, Frederick D.; Singha, Kamini
2015-02-01
The advection-dispersion equation (ADE) fails to describe commonly observed non-Fickian solute transport in saturated porous media, necessitating the use of other models such as the dual-domain mass-transfer (DDMT) model. DDMT model parameters are commonly calibrated via curve fitting, providing little insight into the relation between effective parameters and physical properties of the medium. There is a clear need for material characterization techniques that can provide insight into the geometry and connectedness of pore spaces related to transport model parameters. Here, we consider proton nuclear magnetic resonance (NMR), direct-current (DC) resistivity, and complex conductivity (CC) measurements for this purpose, and assess these methods using glass beads as a control and two different samples of the zeolite clinoptilolite, a material that demonstrates non-Fickian transport due to intragranular porosity. We estimate DDMT parameters via calibration of a transport model to column-scale solute tracer tests, and compare NMR, DC resistivity, CC results, which reveal that grain size alone does not control transport properties and measured geophysical parameters; rather, volume and arrangement of the pore space play important roles. NMR cannot provide estimates of more-mobile and less-mobile pore volumes in the absence of tracer tests because these estimates depend critically on the selection of a material-dependent and flow-dependent cutoff time. Increased electrical connectedness from DC resistivity measurements are associated with greater mobile pore space determined from transport model calibration. CC was hypothesized to be related to length scales of mass transfer, but the CC response is unrelated to DDMT.
Simulation and analysis of solute transport in 2D fracture/pipe networks: The SOLFRAC program
NASA Astrophysics Data System (ADS)
Bodin, Jacques; Porel, Gilles; Delay, Fred; Ubertosi, Fabrice; Bernard, Stéphane; de Dreuzy, Jean-Raynald
2007-01-01
The Time Domain Random Walk (TDRW) method has been recently developed by Delay and Bodin [Delay, F. and Bodin, J., 2001. Time domain random walk method to simulate transport by advection-dispersion and matrix diffusion in fracture networks. Geophys. Res. Lett., 28(21): 4051-4054.] and Bodin et al. [Bodin, J., Porel, G. and Delay, F., 2003c. Simulation of solute transport in discrete fracture networks using the time domain random walk method. Earth Planet. Sci. Lett., 6566: 1-8.] for simulating solute transport in discrete fracture networks. It is assumed that the fracture network can reasonably be represented by a network of interconnected one-dimensional pipes (i.e. flow channels). Processes accounted for are: (1) advection and hydrodynamic dispersion in the channels, (2) matrix diffusion, (3) diffusion into stagnant zones within the fracture planes, (4) sorption reactions onto the fracture walls and in the matrix, (5) linear decay, and (6) mass sharing at fracture intersections. The TDRW method is handy and very efficient in terms of computation costs since it allows for the one-step calculation of the particle residence time in each bond of the network. This method has been programmed in C++, and efforts have been made to develop an efficient and user-friendly software, called SOLFRAC. This program is freely downloadable at the URL http://labo.univ-poitiers.fr/hydrasa/intranet/telechargement.htm. It calculates solute transport into 2D pipe networks, while considering different types of injections and different concepts of local dispersion within each flow channel. Post-simulation analyses are also available, such as the mean velocity or the macroscopic dispersion at the scale of the entire network. The program may be used to evaluate how a given transport mechanism influences the macroscopic transport behaviour of fracture networks. It may also be used, as is the case, e.g., with analytical solutions, to interpret laboratory or field tracer test experiments
NASA Astrophysics Data System (ADS)
Gallo, Claudio; Paniconi, Claudio; Gambolati, Giuseppe
The two-domain concept is widely used in modelling transport in heterogeneous porous media and transport of rate-limited sorbing contaminants. When a first-order kinetic relationship is used to represent the transfer of mass between domains, the model can be expressed as a modified advection-dispersion equation describing general transport coupled to a first-order ordinary differential equation accounting for mass transfer. Different approaches can be used to solve the resulting system, including: simultaneously solving the coupled transport and kinetic equations; discretising and algebraically solving the mass transfer equation and substituting it into the transport equation; solving the mass transfer equation analytically and substituting the integral solution into the transport equation to obtain a single integro-differential equation; and solving the system in Laplace space and back-transforming the solution into the time domain. These four approaches — coupled, algebraic substitution, integro-differential, and finite element Laplace transform (FELT) — are evaluated on the basis of their general features and on their performance in two test cases. The results indicate that the algebraic substitution approach is robust and, on scalar computers, verr efficient. The FELT approach is easily parallelised and achieves good speed-up on supercomputers, but the method is restricted to time-invariant velocity and saturation fields, and is only useful for obtaining the solution at or not too far from the maximum simulation time. The integro-differential method is as efficient as but less robust than the algebraic substitution approach, requiring a small time step size when the mass transfer coefficient is very large. Finally, the coupled approach is robust and flexible, but requires the solution of a system of equations twice as large as the other methods. On balance, the algebraic substitution and, to a lesser extent, the integro-differential methods appear to be the
NASA Astrophysics Data System (ADS)
Bouchelaghem, F.; Vulliet, L.; Leroy, D.; Laloui, L.; Descoeudres, F.
2001-10-01
A model was developed, to describe miscible grout propagation in a saturated deformable porous medium, based on Bear's statistical model with spatial volume averaging. In a previous paper, the model was first successfully confronted to one-dimensional laboratory experiments.In the present paper, the numerical model is used to simulate practical grouting operation in a cylindrical injection model. The cylindrical injection model lends itself to study main flow and propagation character istics for a dispersed suspension-type grout, under axisymmetric conditions close to real scale conditions.Comparison between numerical solutions and experimental results is essential to confirm the validity and accuracy of the proposed model from a phenomenological standpoint. The numerical model performances show that the underlying mathematical model constitutes a realistic predictive model reproducing most prominent features during injection of a suspension-type grout into a deformable porous medium. The basic mechanism by which injected miscible grout permeates a soil mass is discussed in detail. Such a tool leads to quality control criteria for grouting on a theoretical basis, which complements existing criteria acquired through engineering practice.
Cox, T.J.; Runkel, R.L.
2008-01-01
Past applications of one-dimensional advection, dispersion, and transient storage zone models have almost exclusively relied on a central differencing, Eulerian numerical approximation to the nonconservative form of the fundamental equation. However, there are scenarios where this approach generates unacceptable error. A new numerical scheme for this type of modeling is presented here that is based on tracking Lagrangian control volumes across a fixed (Eulerian) grid. Numerical tests are used to provide a direct comparison of the new scheme versus nonconservative Eulerian numerical methods, in terms of both accuracy and mass conservation. Key characteristics of systems for which the Lagrangian scheme performs better than the Eulerian scheme include: nonuniform flow fields, steep gradient plume fronts, and pulse and steady point source loadings in advection-dominated systems. A new analytical derivation is presented that provides insight into the loss of mass conservation in the nonconservative Eulerian scheme. This derivation shows that loss of mass conservation in the vicinity of spatial flow changes is directly proportional to the lateral inflow rate and the change in stream concentration due to the inflow. While the nonconservative Eulerian scheme has clearly worked well for past published applications, it is important for users to be aware of the scheme's limitations. ?? 2008 ASCE.
Boundary effects on solute transport in finite soil columns
NASA Astrophysics Data System (ADS)
Schwartz, R. C.; McInnes, K. J.; Juo, A. S. R.; Wilding, L. P.; Reddell, D. L.
1999-03-01
This study investigates the influence of inlet and outlet disturbances and formulated boundary conditions on the estimation of the dispersion coefficient and retardation factor for short soil columns. Unsaturated miscible displacement experiments utilizing a Br- tracer were carried out on undisturbed columns of a fine-textured Ultisol. Solutions were applied using either a fritted plate or an array of dispensing tips that produced droplets at a prescribed flow rate. One- and two-layer analytical solutions of the advective-dispersive equation were fitted to effluent concentrations using nonlinear least squares parameter optimization. Comparison of two-layer simulations with experimental data indicated that the analytical solution with a semi-infinite interface boundary best approximated effluent concentrations under the conditions of this study. This solution corresponds to a continuous flux concentration and a macroscopically discontinuous resident concentration at the interface between the soil and porous plates. Parameter estimates were not significantly different with respect to the application method used at the inlet. This may be attributed to a less uniform distribution of solution onto the soil surface by the drip apparatus and/or by the presence of stagnant regions within the inlet reservoir and hence increased dispersion within the inlet platen apparatus. Two-layer simulations indicated that the dispersion coefficient was underestimated by 14-27% when the influence of the inlet and outlet apparatus were not included in the fitted solution of the advective-dispersive equation. In addition, use of one-layer analytical solutions caused the retardation factor to be overestimated by no more than the fractional increase in pore volume imparted by the platen apparatus.
Continuous time random walks for non-local radial solute transport
NASA Astrophysics Data System (ADS)
Dentz, Marco; Kang, Peter K.; Le Borgne, Tanguy
2015-08-01
This study formulates and analyzes continuous time random walk (CTRW) models in radial flow geometries for the quantification of non-local solute transport induced by heterogeneous flow distributions and by mobile-immobile mass transfer processes. To this end we derive a general CTRW framework in radial coordinates starting from the random walk equations for radial particle positions and times. The particle density, or solute concentration is governed by a non-local radial advection-dispersion equation (ADE). Unlike in CTRWs for uniform flow scenarios, particle transition times here depend on the radial particle position, which renders the CTRW non-stationary. As a consequence, the memory kernel characterizing the non-local ADE, is radially dependent. Based on this general formulation, we derive radial CTRW implementations that (i) emulate non-local radial transport due to heterogeneous advection, (ii) model multirate mass transfer (MRMT) between mobile and immobile continua, and (iii) quantify both heterogeneous advection in a mobile region and mass transfer between mobile and immobile regions. The expected solute breakthrough behavior is studied using numerical random walk particle tracking simulations. This behavior is analyzed by explicit analytical expressions for the asymptotic solute breakthrough curves. We observe clear power-law tails of the solute breakthrough for broad (power-law) distributions of particle transit times (heterogeneous advection) and particle trapping times (MRMT model). The combined model displays two distinct time regimes. An intermediate regime, in which the solute breakthrough is dominated by the particle transit times in the mobile zones, and a late time regime that is governed by the distribution of particle trapping times in immobile zones. These radial CTRW formulations allow for the identification of heterogeneous advection and mobile-immobile processes as drivers of anomalous transport, under conditions relevant for field tracer
NASA Astrophysics Data System (ADS)
Gao, G.; Zhan, H.; Feng, S.; Huang, G.; Mao, X.
2008-12-01
This study compared five different models for evaluating solute transport in a 1,250-cm long, saturated and highly heterogeneous soil column. The five models were: the convection-dispersion equation (CDE), the mobile-immobile model (MIM), the convective lognormal transfer function model (CLT), the spatial fractional advection-dispersion equation (FADE) and the continuous time random walk model (CTRW). These models were used to fit each breakthrough curve (BTC) and also fitted to the measured BTCs at different distances simultaneously. In addition, the estimated parameters at 200 cm were used to predict the BTCs at subsequent distances. Non-Fickian transport behavior was found to dominate as the measured BTCs demonstrated nonsigmoidal shape and distinct tailing, and solute transport underwent a transition from notable non-Fickian to Fickian behavior with the increase of transport scale in terms of the CTRW theory. The CDE, CLT and FADE were all unable to describe the measured non-Fickian BTCs adequately although the FADE provided better simulation results at the tailing parts of BTCs than CDE and CLT. Compared to the CDE, CLT and FADE, both the MIM and CTRW better captured the full evolution of the measured BTCs. However, the modeling results of MIM at the tails of BTCs were somewhat smaller than the measured results, while the modeling results of CTRW were over the measured results at the BTCs tails. A generic analysis of the advantages and disadvantages of these five models under various circumstances was also provided.
Younes, A; Delay, F; Fajraoui, N; Fahs, M; Mara, T A
2016-08-01
The concept of dual flowing continuum is a promising approach for modeling solute transport in porous media that includes biofilm phases. The highly dispersed transit time distributions often generated by these media are taken into consideration by simply stipulating that advection-dispersion transport occurs through both the porous and the biofilm phases. Both phases are coupled but assigned with contrasting hydrodynamic properties. However, the dual flowing continuum suffers from intrinsic equifinality in the sense that the outlet solute concentration can be the result of several parameter sets of the two flowing phases. To assess the applicability of the dual flowing continuum, we investigate how the model behaves with respect to its parameters. For the purpose of this study, a Global Sensitivity Analysis (GSA) and a Statistical Calibration (SC) of model parameters are performed for two transport scenarios that differ by the strength of interaction between the flowing phases. The GSA is shown to be a valuable tool to understand how the complex system behaves. The results indicate that the rate of mass transfer between the two phases is a key parameter of the model behavior and influences the identifiability of the other parameters. For weak mass exchanges, the output concentration is mainly controlled by the velocity in the porous medium and by the porosity of both flowing phases. In the case of large mass exchanges, the kinetics of this exchange also controls the output concentration. The SC results show that transport with large mass exchange between the flowing phases is more likely affected by equifinality than transport with weak exchange. The SC also indicates that weakly sensitive parameters, such as the dispersion in each phase, can be accurately identified. Removing them from calibration procedures is not recommended because it might result in biased estimations of the highly sensitive parameters. PMID:27182791
NASA Astrophysics Data System (ADS)
Younes, A.; Delay, F.; Fajraoui, N.; Fahs, M.; Mara, T. A.
2016-08-01
The concept of dual flowing continuum is a promising approach for modeling solute transport in porous media that includes biofilm phases. The highly dispersed transit time distributions often generated by these media are taken into consideration by simply stipulating that advection-dispersion transport occurs through both the porous and the biofilm phases. Both phases are coupled but assigned with contrasting hydrodynamic properties. However, the dual flowing continuum suffers from intrinsic equifinality in the sense that the outlet solute concentration can be the result of several parameter sets of the two flowing phases. To assess the applicability of the dual flowing continuum, we investigate how the model behaves with respect to its parameters. For the purpose of this study, a Global Sensitivity Analysis (GSA) and a Statistical Calibration (SC) of model parameters are performed for two transport scenarios that differ by the strength of interaction between the flowing phases. The GSA is shown to be a valuable tool to understand how the complex system behaves. The results indicate that the rate of mass transfer between the two phases is a key parameter of the model behavior and influences the identifiability of the other parameters. For weak mass exchanges, the output concentration is mainly controlled by the velocity in the porous medium and by the porosity of both flowing phases. In the case of large mass exchanges, the kinetics of this exchange also controls the output concentration. The SC results show that transport with large mass exchange between the flowing phases is more likely affected by equifinality than transport with weak exchange. The SC also indicates that weakly sensitive parameters, such as the dispersion in each phase, can be accurately identified. Removing them from calibration procedures is not recommended because it might result in biased estimations of the highly sensitive parameters.
NASA Astrophysics Data System (ADS)
Swanson, R. D.; Binley, A.; Keating, K.; Haggerty, R.; Day-Lewis, F. D.; Singha, K.
2012-12-01
The advection-dispersion equation cannot describe non-Fickian solute transport in saturated porous media and often fails to match breakthrough curve (BTC) history. The multirate mass transfer (MRMT) model partitions the total porosity into immobile and mobile domains with a distribution of exchange rates between the two domains; consequently, the MRMT model produces a better fit to BTCs. However, direct experimental support for the MRMT model parameters remains elusive and model parameters are often estimated a posteriori by an optimization procedure. Complex and direct-current electrical resistivity methods have been used to monitor non-Fickian solute transport in groundwater, but the electrical response has yet to be interpreted within a multirate framework. Here, we investigate electrical geophysical methods to improve our characterization of MRMT parameters. We explore the electrical response in two separate steps: (1) we simulate the direct current electrical response within a multirate framework in order to estimate, from temporal moments, an effective, single rate of mass transfer, and; (2) we develop an empirical link between length scales of multirate mass transfer and length scales of relaxation time distributions measured from complex resistivity at the laboratory scale for the zeolite clinoptilolite which has previously demonstrated MRMT behavior. We use nuclear magnetic resonance measurements of the zeolite to estimate the mobile and immobile porosity of the sample. This study demonstrates our approach at the laboratory scale and offers future perspectives for field investigations.
Convergent radial tracing of viral and solute transport in gneiss saprolite.
Taylor, Richard; Tindimugaya, Callist; Barker, John; Macdonald, David; Kulabako, Robinah
2010-01-01
Deeply weathered crystalline rock aquifer systems comprising unconsolidated saprolite and underlying fractured bedrock (saprock) underlie 40% of sub-Saharan Africa. The vulnerability of this aquifer system to contamination, particularly in rapidly urbanizing areas, remains poorly understood. In order to assess solute and viral transport in saprolite derived from Precambrian gneiss, forced-gradient tracer experiments using chloride and Escherichia coli phage PhiX174 were conducted in southeastern Uganda. The bacteriophage tracer was largely unrecovered; adsorption to the weathered crystalline rock matrix is inferred and enabled by the low pH (5.7) of site ground water and the bacteriophage's relatively high isoelectric point (pI = 6.6). Detection of the applied PhiX174 phage in the pumping well discharge at early times during the experiment traces showed, however, that average ground water flow velocities exceed that of the inert solute tracer, chloride. This latter finding is consistent with observations in other hydrogeological environments where statistically extreme sets of microscopic flow velocities are considered to transport low numbers of fecal pathogens and their proxies along a selected range of linked ground water pathways. Application of a radial advection-dispersion model with an exponentially decaying source term to the recovered chloride tracer estimates a dispersivity (alpha) of 0.8 +/- 0.1 m over a distance of 4.15 m. Specific yield (S(y)) is estimated to be 0.02 from volume balance calculations based on tracer experiments. As single-site observations, our estimates of saprolite S(y) and alpha are tentative but provide a starting point for assessing the vulnerability of saprolite aquifers in sub-Saharan Africa to contamination and estimating quantitatively the impact of climate and abstraction on ground water storage. PMID:19245374
NASA Astrophysics Data System (ADS)
Liang, Ching-Ping; Hsu, Shao-Yiu; Chen, Jui-Sheng
2016-09-01
It is recommended that an in-situ infiltration tracer test is considered for simultaneously determining the longitudinal and transverse dispersion coefficients in soil. Analytical solutions have been derived for two-dimensional advective-dispersive transport in a radial geometry in the literature which can be used for interpreting the result of such a tracer test. However, these solutions were developed for a transport domain with an unbounded-radial extent and an infinite thickness of vadose zone which might not be realistically manifested in the actual solute transport during a field infiltration tracer test. Especially, the assumption of infinite thickness of vadose zone should be invalid for infiltration tracer tests conducted in soil with a shallow groundwater table. This paper describes an analytical model for interpreting the results of an infiltration tracer test based on improving the transport domain with a bounded-radial extent and a finite thickness of vadose zone. The analytical model is obtained with the successive application of appropriate integral transforms and their corresponding inverse transforms. A comparison of the newly derived analytical solution against the previous analytical solutions in which two distinct sets of radial extent and thickness of vadose zone are considered is conducted to determine the influence of the radial and exit boundary conditions on the solute transport. The results shows that both the radial and exit boundary conditions substantially affect the trailing segment of the breakthrough curves for a soil medium with large dispersion coefficients. Previous solutions derived for a transport domain with an unbounded-radial and an infinite thickness of vadose zone boundary conditions give lower concentration predictions compared with the proposed solution at late times. Moreover, the differences between two solutions are amplified when the observation positions are near the groundwater table. In addition, we compare our
Transport solutions for cleaner air.
Kelly, Frank J; Zhu, Tong
2016-05-20
In cities across the globe, road transport remains an important source of air pollutants that are linked with acute and chronic health effects. Decreasing vehicle emissions--while maintaining or increasing commuter journeys--remains a major challenge for city administrators. In London, congestion-charging and a citywide low-emission zone failed to bring nitrogen dioxide concentrations under control. In Beijing, controls on the purchase and use of cars have not decreased transport emissions to a sufficient extent. As cities continue to grow, not even zero-emission vehicles are the solution. Moving increasingly large numbers of people efficiently around a city can only be achieved by expanding mass transit systems. PMID:27199415
Solute transport at fracture intersections
NASA Astrophysics Data System (ADS)
Mourzenko, V. V.; Yousefian, F.; Kolbah, B.; Thovert, J.-F.; Adler, P. M.
2002-01-01
A numerical study of three-dimensional solute transportat fracture intersections by using a particle tracking technique is presented.Two models of orthogonal fracture intersection are considered, namely, twoparallel-walled channels and two rough-walled Gaussian fractures. The fluidvelocity is calculated by solving the three-dimensional Stokes equation withno-slip boundary condition at the solid wall. Examples of individual trajectoriesof particles are first given in order to illustrate the main features of thephenomenon. Solute mass partitioning between outgoing fracture branches isconsidered for various transport regimes, characterized by the local Pécletnumber, and for various ratios of the flow rates in the intersecting channels.Generally speaking, it can be said that at dominant diffusion the influenceof the flow rates ratio is weak, while it is important in the opposite situation.Validity of the classical models of solute mixing, stream tube routing, andperfect mixing is analyzed by comparing their predictions with the numericaldata. Preliminary recommendations are made for the use of these results inlarge-scale modeling.
NASA Astrophysics Data System (ADS)
Wissmeier, L. C.; Barry, D. A.
2009-12-01
Computer simulations of water availability and quality play an important role in state-of-the-art water resources management. However, many of the most utilized software programs focus either on physical flow and transport phenomena (e.g., MODFLOW, MT3DMS, FEFLOW, HYDRUS) or on geochemical reactions (e.g., MINTEQ, PHREEQC, CHESS, ORCHESTRA). In recent years, several couplings between both genres of programs evolved in order to consider interactions between flow and biogeochemical reactivity (e.g., HP1, PHWAT). Software coupling procedures can be categorized as ‘close couplings’, where programs pass information via the memory stack at runtime, and ‘remote couplings’, where the information is exchanged at each time step via input/output files. The former generally involves modifications of software codes and therefore expert programming skills are required. We present a generic recipe for remotely coupling the PHREEQC geochemical modeling framework and flow and solute transport (FST) simulators. The iterative scheme relies on operator splitting with continuous re-initialization of PHREEQC and the FST of choice at each time step. Since PHREEQC calculates the geochemistry of aqueous solutions in contact with soil minerals, the procedure is primarily designed for couplings to FST’s for liquid phase flow in natural environments. It requires the accessibility of initial conditions and numerical parameters such as time and space discretization in the input text file for the FST and control of the FST via commands to the operating system (batch on Windows; bash/shell on Unix/Linux). The coupling procedure is based on PHREEQC’s capability to save the state of a simulation with all solid, liquid and gaseous species as a PHREEQC input file by making use of the dump file option in the TRANSPORT keyword. The output from one reaction calculation step is therefore reused as input for the following reaction step where changes in element amounts due to advection/dispersion
Coupled Fluid Energy Solute Transport
Energy Science and Technology Software Center (ESTSC)
1992-02-13
CFEST is a Coupled Fluid, Energy, and Solute Transport code for the study of a multilayered, nonisothermal ground-water system. It can model discontinuous as well as continuous layers, time-dependent and constant source/sinks, and transient as well as steady-state flow. The finite element method is used for analyzing isothermal and nonisothermal events in a confined aquifer system. Only single-phase Darcian flow is considered. In the Cartesian coordinate system, flow in a horizontal plane, in a verticalmore » plane, or in a fully three-dimensional region can be simulated. An option also exists for the axisymmetric analysis of a vertical cross section. The code employs bilinear quadrilateral elements in all two dimensional analyses and trilinear quadrilateral solid elements in three dimensional simulations. The CFEST finite element formulation can approximate discontinuities, major breaks in slope or thickness, and fault zones in individual hydrogeologic units. The code accounts for heterogeneity in aquifer permeability and porosity and accommodates anisotropy (collinear with the Cartesian coordinates). The variation in the hydraulic properties is described on a layer-by-layer basis for the different hydrogeologic units. Initial conditions can be prescribed hydraulic head or pressure, temperature, or concentration. CFEST can be used to support site, repository, and waste package subsystem assessments. Some specific applications are regional hydrologic characterization; simulation of coupled transport of fluid, heat, and salinity in the repository region; consequence assessment due to natural disruption or human intrusion scenarios in the repository region; flow paths and travel-time estimates for transport of radionuclides; and interpretation of well and tracer tests.« less
NASA Astrophysics Data System (ADS)
Cremer, Clemens; Neuweiler, Insa; Bechtold, Michel; Vanderborght, Jan
2016-04-01
behavior depends on the magnitude of the flow rates and hydraulic conductivity curves of the materials. Based on the unsaturated hydraulic conductivity at the intersection point of conductivity curves, we are able to define an estimate of flow rates at which the dynamic of the upper boundary condition significantly alters preferential flow paths through the system. If flow rates are low, with regard to the materials hydraulic conductivity at the intersection point, the influence of dynamic boundary conditions is small. If flow rates are in the range of the unsaturated hydraulic conductivity at intersection, solute is trapped in the fine material during upwards transport, which results in a more pronounced tailing. For flow rates exceeding the intersection conductivity, a redistribution at the soil surface can occur. References: Bechtold, M., S. Haber-Pohlmeier, J. Vanderborght, A. Pohlmeier, T.P.A. Ferré and H. Veerecken. 2011a. Near-surface solute redistribution during evaporation. Geophys. Res. Lett., 38, L17404, doi:10.1029/2011GL048147. Bechtold, M., J. Vanderborght, O. Ippisch and H. Vereecken. 2011b. Efficient random walk particle tracking algorithm for advective dispersive transport in media with discontinuous dispersion coefficients and water contents. Water Resour. Res., 47, W10526, doi: 10.1029/2010WR010267. Ippisch O., H.-J. Vogel and P. Bastian. 2006. Validity limits fort he van Genuchten-Mualem model and implications for parameter estimation and numerical simulation. Adv. Water Resour., 29, 1780-1789, doi: 10.1016/j.advwateres.2005.12.011. Lehmann, P. and D. Or. 2009. Evaporation and capillary coupling across vertical textural contrasts in porous media. Phys. Rev. E, 80, 046318, doi:10.1103/PhysRevE.80.046318.
Water and solute transport in agricultural soils predicted by volumetric clay and silt contents.
Karup, Dan; Moldrup, Per; Paradelo, Marcos; Katuwal, Sheela; Norgaard, Trine; Greve, Mogens H; de Jonge, Lis W
2016-09-01
Solute transport through the soil matrix is non-uniform and greatly affected by soil texture, soil structure, and macropore networks. Attempts have been made in previous studies to use infiltration experiments to identify the degree of preferential flow, but these attempts have often been based on small datasets or data collected from literature with differing initial and boundary conditions. This study examined the relationship between tracer breakthrough characteristics, soil hydraulic properties, and basic soil properties. From six agricultural fields in Denmark, 193 intact surface soil columns 20cm in height and 20cm in diameter were collected. The soils exhibited a wide range in texture, with clay and organic carbon (OC) contents ranging from 0.03 to 0.41 and 0.01 to 0.08kgkg(-1), respectively. All experiments were carried out under the same initial and boundary conditions using tritium as a conservative tracer. The breakthrough characteristics ranged from being near normally distributed to gradually skewed to the right along with an increase in the content of the mineral fines (particles ≤50μm). The results showed that the mineral fines content was strongly correlated to functional soil structure and the derived tracer breakthrough curves (BTCs), whereas the OC content appeared less important for the shape of the BTC. Organic carbon was believed to support the stability of the soil structure rather than the actual formation of macropores causing preferential flow. The arrival times of 5% and up to 50% of the tracer mass were found to be strongly correlated with volumetric fines content. Predicted tracer concentration breakthrough points as a function of time up to 50% of applied tracer mass could be well fitted to an analytical solution to the classical advection-dispersion equation. Both cumulative tracer mass and concentration as a function of time were well predicted from the simple inputs of bulk density, clay and silt contents, and applied tracer mass
Water and solute transport in agricultural soils predicted by volumetric clay and silt contents
NASA Astrophysics Data System (ADS)
Karup, Dan; Moldrup, Per; Paradelo, Marcos; Katuwal, Sheela; Norgaard, Trine; Greve, Mogens H.; de Jonge, Lis W.
2016-09-01
Solute transport through the soil matrix is non-uniform and greatly affected by soil texture, soil structure, and macropore networks. Attempts have been made in previous studies to use infiltration experiments to identify the degree of preferential flow, but these attempts have often been based on small datasets or data collected from literature with differing initial and boundary conditions. This study examined the relationship between tracer breakthrough characteristics, soil hydraulic properties, and basic soil properties. From six agricultural fields in Denmark, 193 intact surface soil columns 20 cm in height and 20 cm in diameter were collected. The soils exhibited a wide range in texture, with clay and organic carbon (OC) contents ranging from 0.03 to 0.41 and 0.01 to 0.08 kg kg- 1, respectively. All experiments were carried out under the same initial and boundary conditions using tritium as a conservative tracer. The breakthrough characteristics ranged from being near normally distributed to gradually skewed to the right along with an increase in the content of the mineral fines (particles ≤ 50 μm). The results showed that the mineral fines content was strongly correlated to functional soil structure and the derived tracer breakthrough curves (BTCs), whereas the OC content appeared less important for the shape of the BTC. Organic carbon was believed to support the stability of the soil structure rather than the actual formation of macropores causing preferential flow. The arrival times of 5% and up to 50% of the tracer mass were found to be strongly correlated with volumetric fines content. Predicted tracer concentration breakthrough points as a function of time up to 50% of applied tracer mass could be well fitted to an analytical solution to the classical advection-dispersion equation. Both cumulative tracer mass and concentration as a function of time were well predicted from the simple inputs of bulk density, clay and silt contents, and applied tracer
Lin, Guang; Tartakovsky, Alexandre M.
2009-05-01
In this study, a probabilistic collocation method (PCM) on sparse grids was used to solve stochastic equations describing flow and transport in three-dimensional in saturated, randomly heterogeneous porous media. Karhunen-Lo\\`{e}ve (KL) decomposition was used to represent the three-dimensional log hydraulic conductivity $Y=\\ln K_s$. The hydraulic head $h$ and average pore-velocity $\\bf v$ were obtained by solving the three-dimensional continuity equation coupled with Darcy's law with random hydraulic conductivity field. The concentration was computed by solving a three-dimensional stochastic advection-dispersion equation with stochastic average pore-velocity $\\bf v$ computed from Darcy's law. PCM is an extension of the generalized polynomial chaos (gPC) that couples gPC with probabilistic collocation. By using the sparse grid points, PCM can handle a random process with large number of random dimensions, with relatively lower computational cost, compared to full tensor products. Monte Carlo (MC) simulations have also been conducted to verify accuracy of the PCM. By comparing the MC and PCM results for mean and standard deviation of concentration, it is evident that the PCM approach is computational more efficient than Monte Carlo simulations. Unlike the conventional moment-equation approach, there is no limitation on the amplitude of random perturbation in PCM. Furthermore, PCM on sparse grids can efficiently simulate solute transport in randomly heterogeneous porous media with large variances.
Solute transport by a volatile solvent
NASA Astrophysics Data System (ADS)
Brown, Glenn O.; McWhorter, David B.
1990-05-01
In relatively dry porous media, water is transported as both liquid and vapor. Exact knowledge of this two-phase transport, and the phase transfer of water associated with it, is required for the prediction of solute transport. Combined liquid and vapor transport is examined starting from basic principles. An analytic solution is presented for the case of isothermal, transient, one-dimensional sorption of water with constant liquid content boundaries. A relation is also obtained for the evaporation and condensation within the flow field. A numerical solution for the solute transport is obtained which takes maximum advantage of the analytical flow solution. Using the properties of Lurgi retorted oil shale, several special cases are examined which show the relative importance of the separate phases in the total transport of water, the effects on the phase transfer, and the solute transport. It is expected that these methods and results can be applied to other problems in multiple phase transport, such as hazardous waste disposal and pesticide transport.
NASA Astrophysics Data System (ADS)
Gotovac, Hrvoje; Srzic, Veljko
2014-05-01
Contaminant transport in natural aquifers is a complex, multiscale process that is frequently studied using different Eulerian, Lagrangian and hybrid numerical methods. Conservative solute transport is typically modeled using the advection-dispersion equation (ADE). Despite the large number of available numerical methods that have been developed to solve it, the accurate numerical solution of the ADE still presents formidable challenges. In particular, current numerical solutions of multidimensional advection-dominated transport in non-uniform velocity fields are affected by one or all of the following problems: numerical dispersion that introduces artificial mixing and dilution, grid orientation effects, unresolved spatial and temporal scales and unphysical numerical oscillations (e.g., Herrera et al, 2009; Bosso et al., 2012). In this work we will present Eulerian Lagrangian Adaptive Fup Collocation Method (ELAFCM) based on Fup basis functions and collocation approach for spatial approximation and explicit stabilized Runge-Kutta-Chebyshev temporal integration (public domain routine SERK2) which is especially well suited for stiff parabolic problems. Spatial adaptive strategy is based on Fup basis functions which are closely related to the wavelets and splines so that they are also compactly supported basis functions; they exactly describe algebraic polynomials and enable a multiresolution adaptive analysis (MRA). MRA is here performed via Fup Collocation Transform (FCT) so that at each time step concentration solution is decomposed using only a few significant Fup basis functions on adaptive collocation grid with appropriate scales (frequencies) and locations, a desired level of accuracy and a near minimum computational cost. FCT adds more collocations points and higher resolution levels only in sensitive zones with sharp concentration gradients, fronts and/or narrow transition zones. According to the our recent achievements there is no need for solving the large
Mammalian ion-coupled solute transporters.
Hediger, M A; Kanai, Y; You, G; Nussberger, S
1995-01-01
Active transport of solutes into and out of cells proceeds via specialized transporters that utilize diverse energy-coupling mechanisms. Ion-coupled transporters link uphill solute transport to downhill electrochemical ion gradients. In mammals, these transporters are coupled to the co-transport of H+, Na+, Cl- and/or to the countertransport of K+ or OH-. By contrast, ATP-dependent transporters are directly energized by the hydrolysis of ATP. The development of expression cloning approaches to select cDNA clones solely based on their capacity to induce transport function in Xenopus oocytes has led to the cloning of several ion-coupled transporter cDNAs and revealed new insights into structural designs, energy-coupling mechanisms and physiological relevance of the transporter proteins. Different types of mammalian ion-coupled transporters are illustrated by discussing transporters isolated in our own laboratory such as the Na+/glucose co-transporters SGLT1 and SGLT2, the H(+)-coupled oligopeptide transporters PepT1 and PepT2, and the Na(+)- and K(+)-dependent neuronal and epithelial high affinity glutamate transporter EAAC1. Most mammalian ion-coupled organic solute transporters studied so far can be grouped into the following transporter families: (1) the predominantly Na(+)-coupled transporter family which includes the Na+/glucose co-transporters SGLT1, SGLT2, SGLT3 (SAAT-pSGLT2) and the inositol transporter SMIT, (2) the Na(+)- and Cl(-)-coupled transporter family which includes the neurotransmitter transporters of gamma-amino-butyric acid (GABA), serotonin, dopamine, norepinephrine, glycine and proline as well as transporters of beta-amino acids, (3) the Na(+)- and K(+)-dependent glutamate/neurotransmitter family which includes the high affinity glutamate transporters EAAC1, GLT-1, GLAST, EAAT4 and the neutral amino acid transporters ASCT1 and SATT1 reminiscent of system ASC and (4) the H(+)-coupled oligopeptide transporter family which includes the intestinal H
NASA Astrophysics Data System (ADS)
Schuetz, Tobias; Lange, Jens; Weiler, Markus
2010-05-01
Constructed wetlands in agricultural headwater catchments may serve as simple treatment systems to improve retention and mitigation of agricultural non-point-source pollution. To calculate and predict retention capacities of 6 different constructed wetland systems concerning micro-pollutants, we used a one-dimensional solute transport model to compare the results of a series of multi-tracer experiments. The investigated wetland systems consisted of two surface flow wetlands with permanent through flow, two vegetated ditches, a forest buffer zone and a flood detention pond. Transport behaviour was investigated using different tracers: salt and two differently sorptive fluorescent dyes (Sulphorhodamine B and fluoresceine). The hypothesis that shallow and vegetated systems offer the highest sorption capacity for sorptive but mobile pollutants was tested applying a solute transport model to the observed tracer breakthrough. The transport model OTIS (Runkel, 1998) which includes advection, dispersion and lateral exchange to a transient storage was optimized to observed breakthrough of applied tracers at defined cross-sections along the wetlands. Optimized model parameters include dispersivity, cross-sectional areas of both stream and transient storage, as well as an exchange coefficient. Sorption was included based on the KD value, mass of accessible sediment and a sorption coefficient. We assumed that each measurable cross-section is a combination of dead zones and flowing parts. For three of the wetland systems we could exclude lateral in- and outflows. For the other systems, a quantification of lateral flows was possible. We used the set of conservative tracer data to calculate conservative transport characteristics and cross-sections. Then we applied the calibrated model on the sorptive tracer data only using sorption capacity in the storage zone as a calibration parameter and observed KD values and mass of accessible sediment. The results for the different tracer
Bencala, K.E.; McKnight, Diane M.; Zellweger, G.W.
1990-01-01
Physical parameters characterizing solute transport in the Snake River were variable along a 5.2-km study reach. Stream cross-sectional area and volumetric inflow each varied by a factor of 3. Because of transient storage, the residence time of injected tracers in the Snake River was longer than would be calculated by consideration of convective travel time alone. Distributed inflows along the stream were a significant source of in-stream chemical variations. These transport characteristics of the Snake River were established on the basis of the assumption of lithium as an ideally conservative tracer and use of simulations of advection, dispersion, and transient storage. -from Authors
Electrofuels: Versatile Transportation Energy Solutions
2010-07-01
Electrofuels Project: ARPA-E’s Electrofuels Project is using microorganisms to create liquid transportation fuels in a new and different way that could be up to 10 times more energy efficient than current biofuel production methods. ARPA-E is the only U.S. government agency currently funding research on Electrofuels.
Performance Analysis of Solution Transportation Absorption Chiller
NASA Astrophysics Data System (ADS)
Kiani, Behdad; Hamamoto, Yoshinori; Akisawa, Atsushi; Kashiwagi, Takao
Thermally activated advanced absorption cycles are considered promising candidates to replace CFCs, HCFCs and HFCs for residential and commercial applications. In such absorption systems, it is desirable to utilize the waste heat from industries for heating and cooling applications in commercial and residential sectors. For this purpose, it is necessary to transport energy over some distance because the waste heat source and demand are generally located apart from each other. Transportation of steam, hot water or chilled water requires high construction costs for insulation. There is an efficient method of energy transportation using absorption system called “ Solution Transportation Absorption System (STA)”. The solution is transported at an ambient temperature so that tube-insulations not required. This paper shows the simulation of the abovementioned system and the optimal result, using mathematical optimization. The optimum system with industry‧s waste heat utilization is obtained. At the end, the effect on the pollution emission and energy conservation is obtained.
Guyonvarch, Estelle; Ramin, Elham; Kulahci, Murat; Plósz, Benedek Gy
2015-10-15
The present study aims at using statistically designed computational fluid dynamics (CFD) simulations as numerical experiments for the identification of one-dimensional (1-D) advection-dispersion models - computationally light tools, used e.g., as sub-models in systems analysis. The objective is to develop a new 1-D framework, referred to as interpreted CFD (iCFD) models, in which statistical meta-models are used to calculate the pseudo-dispersion coefficient (D) as a function of design and flow boundary conditions. The method - presented in a straightforward and transparent way - is illustrated using the example of a circular secondary settling tank (SST). First, the significant design and flow factors are screened out by applying the statistical method of two-level fractional factorial design of experiments. Second, based on the number of significant factors identified through the factor screening study and system understanding, 50 different sets of design and flow conditions are selected using Latin Hypercube Sampling (LHS). The boundary condition sets are imposed on a 2-D axi-symmetrical CFD simulation model of the SST. In the framework, to degenerate the 2-D model structure, CFD model outputs are approximated by the 1-D model through the calibration of three different model structures for D. Correlation equations for the D parameter then are identified as a function of the selected design and flow boundary conditions (meta-models), and their accuracy is evaluated against D values estimated in each numerical experiment. The evaluation and validation of the iCFD model structure is carried out using scenario simulation results obtained with parameters sampled from the corners of the LHS experimental region. For the studied SST, additional iCFD model development was carried out in terms of (i) assessing different density current sub-models; (ii) implementation of a combined flocculation, hindered, transient and compression settling velocity function; and (iii
NASA Astrophysics Data System (ADS)
Swanson, Ryan David
The advection-dispersion equation (ADE) fails to describe non-Fickian solute transport breakthrough curves (BTCs) in saturated porous media in both laboratory and field experiments, necessitating the use of other models. The dual-domain mass transfer (DDMT) model partitions the total porosity into mobile and less-mobile domains with an exchange of mass between the two domains, and this model can reproduce better fits to BTCs in many systems than ADE-based models. However, direct experimental estimation of DDMT model parameters remains elusive and model parameters are often calculated a posteriori by an optimization procedure. Here, we investigate the use of geophysical tools (direct-current resistivity, nuclear magnetic resonance, and complex conductivity) to estimate these model parameters directly. We use two different samples of the zeolite clinoptilolite, a material shown to demonstrate solute mass transfer due to a significant internal porosity, and provide the first evidence that direct-current electrical methods can track solute movement into and out of a less-mobile pore space in controlled laboratory experiments. We quantify the effects of assuming single-rate DDMT for multirate mass transfer systems. We analyze pore structures using material characterization methods (mercury porosimetry, scanning electron microscopy, and X-ray computer tomography), and compare these observations to geophysical measurements. Nuclear magnetic resonance in conjunction with direct-current resistivity measurements can constrain mobile and less-mobile porosities, but complex conductivity may have little value in relation to mass transfer despite the hypothesis that mass transfer and complex conductivity lengths scales are related. Finally, we conduct a geoelectrical monitored tracer test at the Macrodispersion Experiment (MADE) site in Columbus, MS. We relate hydraulic and electrical conductivity measurements to generate a 3D hydraulic conductivity field, and compare to
Liu, Gaisheng; Lu, Zhiming; Zhang, Dongxiao
2007-01-01
A new approach has been developed for solving solute transport problems in randomly heterogeneous media using the Karhunen-Loève-based moment equation (KLME) technique proposed by Zhang and Lu (2004). The KLME approach combines the Karhunen-Loève decomposition of the underlying random conductivity field and the perturbative and polynomial expansions of dependent variables including the hydraulic head, flow velocity, dispersion coefficient, and solute concentration. The equations obtained in this approach are sequential, and their structure is formulated in the same form as the original governing equations such that any existing simulator, such as Modular Three-Dimensional Multispecies Transport Model for Simulation of Advection, Dispersion, and Chemical Reactions of Contaminants in Groundwater Systems (MT3DMS), can be directly applied as the solver. Through a series of two-dimensional examples, the validity of the KLME approach is evaluated against the classical Monte Carlo simulations. Results indicate that under the flow and transport conditions examined in this work, the KLME approach provides an accurate representation of the mean concentration. For the concentration variance, the accuracy of the KLME approach is good when the conductivity variance is 0.5. As the conductivity variance increases up to 1.0, the mismatch on the concentration variance becomes large, although the mean concentration can still be accurately reproduced by the KLME approach. Our results also indicate that when the conductivity variance is relatively large, neglecting the effects of the cross terms between velocity fluctuations and local dispersivities, as done in some previous studies, can produce noticeable errors, and a rigorous treatment of the dispersion terms becomes more appropriate.
Multilevel transport solution of LWR reactor cores
Jose Ignacio Marquez Damian; Cassiano R.E. de Oliveira; HyeonKae Park
2008-09-01
This work presents a multilevel approach for the solution of the transport equation in typical LWR assemblies and core configurations. It is based on the second-order, even-parity formulation of the transport equation, which is solved within the framework provided by the finite element-spherical harmonics code EVENT. The performance of the new solver has been compared with that of the standard conjugate gradient solver for diffusion and transport problems on structured and unstruc-tured grids. Numerical results demonstrate the potential of the multilevel scheme for realistic reactor calculations.
Solute transport through a deforming porous medium
NASA Astrophysics Data System (ADS)
Peters, Glen P.; Smith, David W.
2002-06-01
Solute transport through a porous medium is typically modelled assuming the porous medium is rigid. However, many applications exist where the porous medium is deforming, including, municipal landfill liners, mine tailings dams, and land subsidence. In this paper, mass balance laws are used to derive the flow and transport equations for a deforming porous medium. The equations are derived in both spatial and material co-ordinate systems. Solute transport through an engineered landfill liner is used as an illustrative example to show the differences between the theory for a rigid porous medium, and small and large deformation analysis of a deforming porous medium. It is found that the large deformation model produces shorter solute breakthrough times, followed by the small deformation model, and then the rigid porous medium model. It is also found that it is important to include spatial and temporal void ratio variations in the large deformation analysis. It is shown that a non-linear large deformation model may greatly reduce the solute breakthrough time, compared to a standard transport analysis typically employed by environmental engineers.
Linear transport models for adsorbing solutes
NASA Astrophysics Data System (ADS)
Roth, K.; Jury, W. A.
1993-04-01
A unified linear theory for the transport of adsorbing solutes through soils is presented and applied to analyze movement of napropamide through undisturbed soil columns. The transport characteristics of the soil are expressed in terms of the travel time distribution of the mobile phase which is then used to incorporate local interaction processes. This approach permits the analysis of all linear transport processes, not only the small subset for which a differential description is known. From a practical point of view, it allows the direct use of measured concentrations or fluxes of conservative solutes to characterize the mobile phase without first subjecting them to any model. For complicated flow regimes, this may vastly improve the identification of models and estimation of their parameters for the local adsorption processes.
Quasi 3D modeling of water flow and solute transport in vadose zone and groundwater
NASA Astrophysics Data System (ADS)
Yakirevich, A.; Kuznetsov, M.; Weisbrod, N.; Pachepsky, Y. A.
2013-12-01
The complexity of subsurface flow systems calls for a variety of concepts leading to the multiplicity of simplified flow models. One commonly used simplification is based on the assumption that lateral flow and transport in unsaturated zone is insignificant unless the capillary fringe is involved. In such cases the flow and transport in the unsaturated zone above groundwater level can be simulated as a 1D phenomenon, whereas through groundwater they are viewed as 2D or 3D phenomena. A new approach for a numerical scheme for 3D variably saturated flow and transport is presented. A Quasi-3D approach allows representing flow in the 'vadose zone - aquifer' system by a series of 1D Richards' equations solved in variably-saturated zone and by 3D-saturated flow equation in groundwater (modified MODFLOW code). The 1D and 3D equations are coupled at the phreatic surface in a way that aquifer replenishment is calculated using the Richards' equation, and solving for the moving water table does not require definition of the specific yield parameter. The 3D advection-dispersion equation is solved in the entire domain by the MT3D code. Using implicit finite differences approximation to couple processes in the vadose zone and groundwater provides mass conservation and increase of computational efficiency. The above model was applied to simulate the impact of irrigation on groundwater salinity in the Alto Piura aquifer (Northern Peru). Studies on changing groundwater quality in arid and semi-arid lands show that irrigation return flow is one of the major factors contributing to aquifer salinization. Existing mathematical models do not account explicitly for the solute recycling during irrigation on a daily scale. Recycling occurs throughout the unsaturated and saturated zones, as function of the solute mass extracted from pumping wells. Salt concentration in irrigation water is calculated at each time step as a function of concentration of both surface water and groundwater
Solute transport in heterogeneous porous formations
NASA Astrophysics Data System (ADS)
Demmy, George Gary, Jr.
1999-10-01
This work quantifies relationships between the spatial, or Eulerian, distribution of the properties of a chemically and physically heterogeneous porous medium and those as observed along the natural, or Lagrangian, trajectories that a fluid particle traces in a steady and irrotational flow. From these relationships, expressions that relate the transport of solutes through the porous medium along the natural trajectories to the aforementioned Eulerian distributions are developed. The effects of injection mode upon global measures of transport as reflected by the temporal moments of breakthrough curves and spatial moments of a solute plume are developed. The coupled effects of correlation of a linear equilibrium sorption to the underlying log hydraulic conductivity field and injection mode on the evolving temporal moments of mass breakthrough curve and the coupled effects of correlation of a first-order decay coefficient and injection mode upon the spatial moments of a solute plume are examined.
Probing Nanoscale Thermal Transport in Surfactant Solutions
Cao, Fangyu; Liu, Ying; Xu, Jiajun; He, Yadong; Hammouda, B.; Qiao, Rui; Yang, Bao
2015-01-01
Surfactant solutions typically feature tunable nanoscale, internal structures. Although rarely utilized, they can be a powerful platform for probing thermal transport in nanoscale domains and across interfaces with nanometer-size radius. Here, we examine the structure and thermal transport in solution of AOT (Dioctyl sodium sulfosuccinate) in n-octane liquids using small-angle neutron scattering, thermal conductivity measurements, and molecular dynamics simulations. We report the first experimental observation of a minimum thermal conductivity occurring at the critical micelle concentration (CMC): the thermal conductivity of the surfactant solution decreases as AOT is added till the onset of micellization but increases as more AOT is added. The decrease of thermal conductivity with AOT loading in solutions in which AOT molecules are dispersed as monomers suggests that even the interfaces between individual oleophobic headgroup of AOT molecules and their surrounding non-polar octane molecules can hinder heat transfer. The increase of thermal conductivity with AOT loading after the onset of micellization indicates that the thermal transport in the core of AOT micelles and across the surfactant-oil interfaces, both of which span only a few nanometers, are efficient. PMID:26534840
Probing Nanoscale Thermal Transport in Surfactant Solutions
NASA Astrophysics Data System (ADS)
Cao, Fangyu; Liu, Ying; Xu, Jiajun; He, Yadong; Hammouda, B.; Qiao, Rui; Yang, Bao
2015-11-01
Surfactant solutions typically feature tunable nanoscale, internal structures. Although rarely utilized, they can be a powerful platform for probing thermal transport in nanoscale domains and across interfaces with nanometer-size radius. Here, we examine the structure and thermal transport in solution of AOT (Dioctyl sodium sulfosuccinate) in n-octane liquids using small-angle neutron scattering, thermal conductivity measurements, and molecular dynamics simulations. We report the first experimental observation of a minimum thermal conductivity occurring at the critical micelle concentration (CMC): the thermal conductivity of the surfactant solution decreases as AOT is added till the onset of micellization but increases as more AOT is added. The decrease of thermal conductivity with AOT loading in solutions in which AOT molecules are dispersed as monomers suggests that even the interfaces between individual oleophobic headgroup of AOT molecules and their surrounding non-polar octane molecules can hinder heat transfer. The increase of thermal conductivity with AOT loading after the onset of micellization indicates that the thermal transport in the core of AOT micelles and across the surfactant-oil interfaces, both of which span only a few nanometers, are efficient.
Probing Nanoscale Thermal Transport in Surfactant Solutions.
Cao, Fangyu; Liu, Ying; Xu, Jiajun; He, Yadong; Hammouda, B; Qiao, Rui; Yang, Bao
2015-01-01
Surfactant solutions typically feature tunable nanoscale, internal structures. Although rarely utilized, they can be a powerful platform for probing thermal transport in nanoscale domains and across interfaces with nanometer-size radius. Here, we examine the structure and thermal transport in solution of AOT (Dioctyl sodium sulfosuccinate) in n-octane liquids using small-angle neutron scattering, thermal conductivity measurements, and molecular dynamics simulations. We report the first experimental observation of a minimum thermal conductivity occurring at the critical micelle concentration (CMC): the thermal conductivity of the surfactant solution decreases as AOT is added till the onset of micellization but increases as more AOT is added. The decrease of thermal conductivity with AOT loading in solutions in which AOT molecules are dispersed as monomers suggests that even the interfaces between individual oleophobic headgroup of AOT molecules and their surrounding non-polar octane molecules can hinder heat transfer. The increase of thermal conductivity with AOT loading after the onset of micellization indicates that the thermal transport in the core of AOT micelles and across the surfactant-oil interfaces, both of which span only a few nanometers, are efficient. PMID:26534840
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.
Lagrangian simulation of multidimensional anomalous transport at the MADE site
NASA Astrophysics Data System (ADS)
Zhang, Yong; Benson, David A.
2008-04-01
Contaminant transport through regional-scale natural geological formations typically exhibits several ``anomalous'' features, including direction-dependent spreading rates, channeling along preferential flow paths, trapping of solute in relatively immobile domains, and/or the local variation of transport speed. Simulating these plume characteristics can be computationally intensive using a traditional advection-dispersion equation (ADE) because anomalous features of transport generally depend on local-scale subsurface properties. Here we develop an alternative simulation approach that solves the full nonlocal, multidimensional, spatiotemporal fractional-order ADE with variable coefficients in a Lagrangian framework using a novel non-Markovian random walk method. This model allows us to simulate anomalous plumes without the need to explicitly define local-scale heterogeneity. The simple model accurately simulates the tritium plume measured at the extensively characterized MADE test site.
Compositional transport in solidifying aqueous binary solution
NASA Astrophysics Data System (ADS)
Zhong, Jin-Qiang; Yin, Zuo-Chao; Xue, Qiwei; Wettlaufer, John
2013-11-01
We observe the formation of double-diffusive layers adjacent to mushy layers that form during the directional solidification of aqueous ammonium chloride. The plumes emerging from chimney's in the mushy layers continuously supply a buoyancy flux in the (finite) liquid region above, driving downward motion of double-diffusive layers. The downward velocity of the layers is found to be in good agreement with a filling box model that captures the crucial hydrodynamics of the entraining buoyant plumes and compositional transport. We demonstrate that the buoyancy flux through the system decays according to a similarity solution. We note that the experimental findings provide some insight into the brine transport in growing sea ice.
NASA Astrophysics Data System (ADS)
Karwan, D. L.; Saiers, J. E.; Gravelle, J. A.
2006-12-01
Transport of suspended particulate material (SPM) plays a fundamental role in biogeochemical cycling within forested watershed and is a primary water quality concern in managed forests. When present in excess, SPM can degrade aquatic habitat for fish communities and disrupt the connection between surface and ground water. Although monthly and annual suspended loads are often monitored, little is known about the hydrologic transport of these particles from suspension to monitoring location. The downstream transport of suspended particles can be influenced by transient storage mechanisms, such as settling and resuspension, stagnation in side pools, exchange with the streambed or hyporheic zone, and entrapment on stream vegetation and coarse woody debris. A tracer injection experiment was performed in order to compare transport and transient storage of suspended clay-sized particles, comprised of titanium dioxide (1 - 2 μm diameter), with that of a conservative solute, bromide, under baseflow conditions. The solute and particle tracers were applied to a second-order North Idaho stream for four hours and water samples were collected at four locations downstream before, during, and after the tracer injections for analysis of bromide and TiO2 concentrations. A one- dimensional numerical model was applied in inverse mode to the measured breakthrough curve data to quantify the processes that governed solute and particle transport. The results of this analysis indicate that transient- storage processes exerted only a minor influence on the advective-dispersive transport of bromide, while TiO2 transport was influenced by advection, dispersion, sedimentation (on the stream bed and on aquatic vegetation), and slow resuspension. Results of our analysis illustrate the mechanisms and timescale of SPM transport in this watershed and provide insight into the potential response of SPM concentrations to elevated sediment inputs.
NASA Astrophysics Data System (ADS)
Drummond, J. D.; Covino, T. P.; Aubeneau, A. F.; Patil, S.; Leong, D. N.; Ran, L.; Packman, A. I.; Schumer, R.
2010-12-01
Solute transport and hydrologic retention strongly affect biogeochemical processes that are critical to stream ecosystems. Tracer injections are used to characterize solute transport and storage in stream reaches, but the range of processes accurately resolved using this approach is not clear. The solute residence time distribution (RTD) depends on both in-stream mixing and exchange with the hyporheic zone. For shorter residence times, in-stream breakthrough curves (BTCs) can be modeled well with the classical advection-dispersion equation, whereas longer RTDs produce highly skewed in-stream BTCs for which traditional solute models are inappropriate. Observed BTCs have most commonly been modeled with in-stream advection-dispersion plus an exponential RTD, but process-based models suggest that hyporheic retention extends to much longer times and a power-law RTD is more appropriate. We synthesized results from a variety of tracer-injection studies to investigate how experimental design and tracer sensitivity influence the interpretation of tailing behavior and RTDs. We found that BTC tails are often not well observed in stream tracer experiments. The two main reasons for this are: 1) experimental truncation, which occurs when sampling ends before all tracer mass reaches the sampling location, and 2) sensitivity truncation, when tracer concentrations in the tail are too low to be detected reliably above background levels. Continuous Time Random Walk (CTRW) theory was used to determine the effects of tail truncation on tracer mass recovery and tailing behavior. Tail truncation due to both experimental and sensitivity truncation decreased mass recovery and obscured assessment of BTC tailing. Failure to consider tail truncation leads to underestimation of the retention of solutes in the streambed and subsurface (i.e., transient storage). Based on these findings, we propose criteria for stream tracer experiments to minimize tail truncation and improve inverse modeling of
Analytical solutions for anomalous dispersion transport
NASA Astrophysics Data System (ADS)
O'Malley, D.; Vesselinov, V. V.
2014-06-01
Groundwater flow and transport often occur in a highly heterogeneous environment (potentially heterogeneous at multiple spatial scales) and is impacted by geochemical reactions, advection, diffusion, and other pore scale processes. All these factors can give rise to large-scale anomalous dispersive behavior that can make complex model representation and prediction of plume concentrations challenging due to difficulties unraveling all the complexities associated with the governing processes, flow medium, and their parameters. An alternative is to use upscaled stochastic models of anomalous dispersion, and this is the approach used here. Within a probabilistic framework, we derive a number of analytical solutions for several anomalous dispersion models. The anomalous dispersion models are allowed to be either non-Gaussian (α-stable Lévy), correlated, or nonstationary from the Lagrangian perspective. A global sensitivity analysis is performed to gain a greater understanding of the extent to which uncertainty in the parameters associated with the anomalous behavior can be narrowed by examining concentration measurements from a network of monitoring wells and to demonstrate the computational speed of the solutions. The developed analytical solutions are encoded and available for use in the open source computational framework MADS (http://mads.lanl.gov).
Percolation and transport in a sandy soil under a natural hydraulic gradient
Green, C.T.; Stonestrom, D.A.; Bekins, B.A.; Akstin, K.C.; Schulz, M.S.
2005-01-01
[1] Unsaturated flow and transport under a natural hydraulic gradient in a Mediterranean climate were investigated with a field tracer experiment combined with laboratory analyses and numerical modeling. Bromide was applied to the surface of a sandy soil during the dry season. During the subsequent rainy season, repeated sediment sampling tracked the movement of bromide through the profile. Analysis of data on moisture content, matric pressure, unsaturated hydraulic conductivity, bulk density, and soil texture and structure provides insights into parameterization and use of the advective-dispersive modeling approach. Capturing the gross features of tracer and moisture movement with model simulations required an order-of-magnitude increase in laboratory-measured hydraulic conductivity. Wetting curve characteristics better represented field results, calling into question the routine estimation of hydraulic characteristics based only on drying conditions. Measured increases in profile moisture exceeded cumulative precipitation in early winter, indicating that gains from dew drip can exceed losses from evapotranspiration during periods of heavy ("Tule") fog. A single-continuum advective-dispersive modeling approach could not reproduce a peak of bromide that was retained near the soil surface for over 3 years. Modeling of this feature required slow exchange of solute at a transfer rate of 0.5-1 ?? 10-4 d-1 with an immobile volume approaching the residual moisture content.
Multiscale solute transport upscaling for a three-dimensional hierarchical porous medium
NASA Astrophysics Data System (ADS)
Zhang, Mingkan; Zhang, Ye
2015-03-01
A laboratory-generated hierarchical, fully heterogeneous aquifer model (FHM) provides a reference for developing and testing an upscaling approach that integrates large-scale connectivity mapping with flow and transport modeling. Based on the FHM, three hydrostratigraphic models (HSMs) that capture lithological (static) connectivity at different resolutions are created, each corresponding to a sedimentary hierarchy. Under increasing system lnK variances (0.1, 1.0, 4.5), flow upscaling is first conducted to calculate equivalent hydraulic conductivity for individual connectivity (or unit) of the HSMs. Given the computed flow fields, an instantaneous, conservative tracer test is simulated by all models. For the HSMs, two upscaling formulations are tested based on the advection-dispersion equation (ADE), implementing space versus time-dependent macrodispersivity. Comparing flow and transport predictions of the HSMs against those of the reference model, HSMs capturing connectivity at increasing resolutions are more accurate, although upscaling errors increase with system variance. Results suggest: (1) by explicitly modeling connectivity, an enhanced degree of freedom in representing dispersion can improve the ADE-based upscaled models by capturing non-Fickian transport of the FHM; (2) when connectivity is sufficiently resolved, the type of data conditioning used to model transport becomes less critical. Data conditioning, however, is influenced by the prediction goal; (3) when aquifer is weakly-to-moderately heterogeneous, the upscaled models adequately capture the transport simulation of the FHM, despite the existence of hierarchical heterogeneity at smaller scales. When aquifer is strongly heterogeneous, the upscaled models become less accurate because lithological connectivity cannot adequately capture preferential flows; (4) three-dimensional transport connectivities of the hierarchical aquifer differ quantitatively from those analyzed for two-dimensional systems
An analytical model for predicting transport in a coupled vadose/phreatic system
Tomasko, D.
1997-05-01
A simple analytical model is presented for predicting the transport of a contaminant in both the unsaturated (vadose) and saturated (phreatic) zones following a surficial spill. The model incorporates advection, dispersion, adsorption, and first-order decay in both zones and couples the transport processes at the water table. The governing equation is solved by using the method of Laplace transforms, with numerical inversion of the Laplace space equation for concentration. Because of the complexity of the functional form for the Laplace space solution, a numerical methodology using the real and imaginary parts of a Fourier series was implemented. To reduce conservatism in the model, dilution at the water table was also included. Verification of the model is demonstrated by its ability to reproduce the source history at the surface and to replicate appropriate one-dimensional transport through either the vadose or phreatic zone. Because of its simplicity and lack of detailed input data requirements, the model is recommended for scoping calculations.
Modeling flow and solute transport in irrigation furrows
Technology Transfer Automated Retrieval System (TEKTRAN)
This paper presents an internally coupled flow and solute transport model for free-draining irrigation furrows. Furrow hydraulics is simulated with a numerical zero-inertia model and solute transport is computed with a model based on a numerical solution of the cross-section averaged advection-dispe...
Analytical Solutions for Sequentially Reactive Transport with Different Retardation Factors
Sun, Y; Buscheck, T A; Mansoor, K; Lu, X
2001-08-01
Integral transforms have been widely used for deriving analytical solutions for solute transport systems. Often, analytical solutions can only be written in closed form in frequency domains and numerical inverse-transforms have to be involved to obtain semi-analytical solutions in the time domain. For this reason, previously published closed form solutions are restricted either to a small number of species or to the same retardation assumption. In this paper, we applied the solution scheme proposed by Bauer et al. in the time domain. Using available analytical solutions of a single species transport with first-order decay without coupling with its parent species concentration as fundamental solutions, a daughter species concentration can be expressed as a linear function of those fundamental solutions. The implementation of the solution scheme is straight forward and exact analytical solutions are derived for one- and three-dimensional transport systems.
NASA Astrophysics Data System (ADS)
Laine-Kaulio, H.; Karvonen, T.; Koivusalo, H.; Lauren, A.; Saastamoinen, S.
2009-04-01
) with Guelph permeameter in the field, iii) and by means of inverse modelling. The inverse model application was based on calibration of a one-dimensional groundwater model against data on groundwater levels in the study slope. Conductivities of the different soil horizons were adjusted to reproduce the measured groundwater levels of a recession period after artificial irrigation. Conductivity results, together with soil physical and water retention data were applied to parameterise a three-dimensional flow and advection-dispersion model. The model was used to simulate the transport of a chloride tracer plume in the study slope during artificial irrigation. A line-type irrigation source was installed upslope from the study section of the slope. Changes in groundwater levels and chloride concentrations within the study section were observed through well screens. Chloride as a conservative tracer provided an indicator for subsurface flow in the study slope. Intensive irrigation rates were applied to initiate fast lateral preferential flow. Saturated hydraulic conductivities obtained with the three methods were remarkably different. Conductivities obtained with the Guelph permeameter and the groundwater model reduced clearly with soil depth. Higher conductivities near soil surface were due to loose soil structure and preferential flowpaths. Soil core samples yielded the lowest estimates for the saturated hydraulic conductivity, as they represented the small-scale conductivity of the soil texture and soil matrix. The hillslope-scale groundwater model produced the highest estimates that characterised the large-scale structural properties and their impact on lateral preferential flow. Average saturated hydraulic conductivities in the soil core samples were 6E-6 m/s in the eluvial horizon, transition zone and subsoil, and 1E-5 m/s in the illuvial horizon. The average conductivities based on the Guelph measurements varied from 2E-5 m/s in the subsoil to 5E-5 m/s in the eluvial
DCPT v1.0 - New particle tracker for modeling transport in dual-continuum - User's Manual
Pan, Lehua; Liu, Hui Hai; Cushey, Mark; Bodvarsson, Gudmundur
2001-04-01
DCPT (Dual-Continuum Particle Tracker) V1.0 is a new software for simulating solute transport in the subsurface. It is based on the random-walk method for modeling transport processes such as advection, dispersion/diffusion, linear sorption, radioactive decay, and fracture-matrix mass exchange (in fractured porous media). The user shall provide flow-field and other parameters in the form of input files. In Comparison to several analytical and numerical solutions for a number of test cases, DCPT shows excellent performance in both accuracy and efficiency. This report serves as a user's manual of DCPT V1.0. It includes theoretical basis, numerical methods, software structure, input/output description, and examples.
Runkel, Robert L.
2010-01-01
OTEQ is a mathematical simulation model used to characterize the fate and transport of waterborne solutes in streams and rivers. The model is formed by coupling a solute transport model with a chemical equilibrium submodel. The solute transport model is based on OTIS, a model that considers the physical processes of advection, dispersion, lateral inflow, and transient storage. The equilibrium submodel is based on MINTEQ, a model that considers the speciation and complexation of aqueous species, acid-base reactions, precipitation/dissolution, and sorption. Within OTEQ, reactions in the water column may result in the formation of solid phases (precipitates and sorbed species) that are subject to downstream transport and settling processes. Solid phases on the streambed may also interact with the water column through dissolution and sorption/desorption reactions. Consideration of both mobile (waterborne) and immobile (streambed) solid phases requires a unique set of governing differential equations and solution techniques that are developed herein. The partial differential equations describing physical transport and the algebraic equations describing chemical equilibria are coupled using the sequential iteration approach. The model's ability to simulate pH, precipitation/dissolution, and pH-dependent sorption provides a means of evaluating the complex interactions between instream chemistry and hydrologic transport at the field scale. This report details the development and application of OTEQ. Sections of the report describe model theory, input/output specifications, model applications, and installation instructions. OTEQ may be obtained over the Internet at http://water.usgs.gov/software/OTEQ.
NASA Astrophysics Data System (ADS)
Zhu, Chen
2004-08-01
Coprecipitation of barite with trace constituents was simulated with consideration of aqueous speciation and complexation, mixing properties for the binary solid solutions (Zhu, this issue), precipitation and dissolution kinetics, and advective-dispersive transport. Speciation-solubility modeling was used to reproduce BaSO 4-RaSO 4 coprecipitation experimental results, and to calculate CrO 42- aqueous concentrations in equilibrium with a Ba(SO 4,CrO 4) solid solution. Kinetic reaction path modeling was used to simulate the coprecipitation of barite with RaSO 4 to form an onion-like chemically zoned solid upon the cooling of oil field brine. A one-dimensional coupled reactive mass transport model shows a strikingly different transport pattern for the tracer Ra 2+, when the dominant attenuation reaction is with solid solution (Ba, Ra) SO 4 as compared to the case when it is controlled by pure RaSO 4 and barite solids under local equilibrium conditions. A self-enrichment of Ra 2+ in the groundwater and aquifer solid matrix—higher concentrations of Ra 2+ downstream from the reaction front—results from the coprecipitation reaction and advective-dispersive transport. This self-enrichment process generates a secondary tracer source, which has tracer concentrations higher than that of the original source. On the other hand, coprecipitation reactions can reduce Ra 2+ concentrations in groundwater to a much lower level (below ppb) than that of pure RaSO 4(c) solubility (near ppm), which has been used to establish the Ra 2+ concentration limits in groundwater, soil, and nuclear waste repositories.
Modeling two-dimensional reactive transport using a Godunov-mixed finite element method
NASA Astrophysics Data System (ADS)
James, Andrew I.; Jawitz, James W.
2007-05-01
SummaryThe development of a model to simulate transport of materials in variable-depth flows is discussed. The model numerically approximates solutions to the advection-dispersion-reaction equation using a time-splitting technique where the advective, dispersive, and reactive parts of the equation are solved separately. An explicit finite-volume Godunov method is used to approximate the advective part while a hybridized mixed finite element method is used to solve for the dispersive step. A backward Euler method is used to solve the reactive component. Rather than solving each component once at each time step, the advective and reactive steps are fractionally and symmetrically split around the dispersive step, so that half of a reactive and advective step are solved before and after each dispersive step. Since the dispersive step is implicit, but computationally expensive, while the advective step is explicit but has time step constraints, this allows stable and more efficient schemes to be implemented in contrast to non-split or simple time-split algorithms. This technique allows problems with high grid Peclet numbers, such as transport problems with sharp solute fronts, to be solved without oscillations in the solution and with virtually no artificial diffusion. By applying the technique to variable depth flows, a variety of applications to transport and reaction problems in surface water and unconfined aquifers can be undertaken. Numerical results for several non-reactive and reactive transport problems in one- and two-dimensions are presented. Observed convergence rates are up to second-order for these simulations.
A quasilinear model for solute transport under unsaturated flow
Houseworth, J.E.; Leem, J.
2009-05-15
We developed an analytical solution for solute transport under steady-state, two-dimensional, unsaturated flow and transport conditions for the investigation of high-level radioactive waste disposal. The two-dimensional, unsaturated flow problem is treated using the quasilinear flow method for a system with homogeneous material properties. Dispersion is modeled as isotropic and is proportional to the effective hydraulic conductivity. This leads to a quasilinear form for the transport problem in terms of a scalar potential that is analogous to the Kirchhoff potential for quasilinear flow. The solutions for both flow and transport scalar potentials take the form of Fourier series. The particular solution given here is for two sources of flow, with one source containing a dissolved solute. The solution method may easily be extended, however, for any combination of flow and solute sources under steady-state conditions. The analytical results for multidimensional solute transport problems, which previously could only be solved numerically, also offer an additional way to benchmark numerical solutions. An analytical solution for two-dimensional, steady-state solute transport under unsaturated flow conditions is presented. A specific case with two sources is solved but may be generalized to any combination of sources. The analytical results complement numerical solutions, which were previously required to solve this class of problems.
Innovative Solutions to Challenges in Pupil Transportation.
ERIC Educational Resources Information Center
Ross, Jonathan; Burkybile, Sharon
2000-01-01
States have had to budget increasing amounts for mandated pupil-transportation services as their state transportation aid has been slashed dramatically. Among school districts, cooperation and coordination through shared services (consortia) have resulted in safer, more reliable, and more efficient transportation. Implementation advice is…
Solute transporters in plant thylakoid membranes
Schoefs, Benoît
2010-01-01
Plants utilize sunlight to drive photosynthetic energy conversion in the chloroplast thylakoid membrane. Here are located four major photosynthetic complexes, about which we have great knowledge in terms of structure and function. However, much less we know about auxiliary proteins, such as transporters, ensuring an optimum function and turnover of these complexes. The most prominent thylakoid transporter is the proton-translocating ATP-synthase. Recently, four additional transporters have been identified in the thylakoid membrane of Arabidopsis thaliana, namely one copper-transporting P-ATPase, one chloride channel, one phosphate transporter, and one ATP/ADP carrier. Here, we review the current knowledge on the function and physiological role of these transporters during photosynthesis and light stress in plants. Subsequently, we make a survey on the outlook of thylakoid activities awaiting identification of responsible proteins. Such knowledge is necessary to understand the thylakoid network of transporters, and to design strategies for bioengineering crop plants in the future. PMID:20585503
CHROMIUM TRANSPORT, OXIDATION, AND ADSORPTION IN MANGANESE-COATED SAND
We examine how the processes of advection, dispersion, oxidation-reduction, and adsorption combine to affect the transport of chromium through columns packed with pyrolusite (P-MnO$-coated sand. We find that P-Mn02 effectively oxidizes Cr@I) to Cr(VI) and that the extent of oxida...
Biological solutions to transport network design.
Bebber, Daniel P; Hynes, Juliet; Darrah, Peter R; Boddy, Lynne; Fricker, Mark D
2007-09-22
Transport networks are vital components of multicellular organisms, distributing nutrients and removing waste products. Animal and plant transport systems are branching trees whose architecture is linked to universal scaling laws in these organisms. In contrast, many fungi form reticulated mycelia via the branching and fusion of thread-like hyphae that continuously adapt to the environment. Fungal networks have evolved to explore and exploit a patchy environment, rather than ramify through a three-dimensional organism. However, there has been no explicit analysis of the network structures formed, their dynamic behaviour nor how either impact on their ecological function. Using the woodland saprotroph Phanerochaete velutina, we show that fungal networks can display both high transport capacity and robustness to damage. These properties are enhanced as the network grows, while the relative cost of building the network decreases. Thus, mycelia achieve the seemingly competing goals of efficient transport and robustness, with decreasing relative investment, by selective reinforcement and recycling of transport pathways. Fungal networks demonstrate that indeterminate, decentralized systems can yield highly adaptive networks. Understanding how these relatively simple organisms have found effective transport networks through a process of natural selection may inform the design of man-made networks. PMID:17623638
Solute Transport Across a Contact Interface in Deformable Porous Media
Ateshian, Gerard A.; Maas, Steve; Weiss, Jeffrey A.
2012-01-01
A finite element formulation of neutral solute transport across a contact interface between deformable porous media is implemented and validated against analytical solutions. By reducing the integral statements of external virtual work on the two contacting surfaces into a single contact integral, the algorithm automatically enforces continuity of solute molar flux across the contact interface, whereas continuity of the effective solute concentration (a measure of the solute mechano-chemical potential) is achieved using a penalty method. This novel formulation facilitates the analysis of problems in biomechanics where the transport of metabolites across contact interfaces of deformable tissues may be of interest. This contact algorithm is the first to address solute transport across deformable interfaces, and is made available in the public domain, open-source finite element code FEBio (http://mrl.sci.utah.edu/software). PMID:22281406
Solute transport across a contact interface in deformable porous media.
Ateshian, Gerard A; Maas, Steve; Weiss, Jeffrey A
2012-04-01
A finite element formulation of neutral solute transport across a contact interface between deformable porous media is implemented and validated against analytical solutions. By reducing the integral statements of external virtual work on the two contacting surfaces into a single contact integral, the algorithm automatically enforces continuity of solute molar flux across the contact interface, whereas continuity of the effective solute concentration (a measure of the solute mechano-chemical potential) is achieved using a penalty method. This novel formulation facilitates the analysis of problems in biomechanics where the transport of metabolites across contact interfaces of deformable tissues may be of interest. This contact algorithm is the first to address solute transport across deformable interfaces, and is made available in the public domain, open-source finite element code FEBio (http://www.febio.org). PMID:22281406
Electrokinetic induced solute dispersion in porous media; pore network modeling
NASA Astrophysics Data System (ADS)
Li, Shuai; Schotting, Ruud; Raoof, Amir
2013-04-01
Electrokinetic flow plays an important role in remediation process, separation technique, and chromatography. The solute dispersion is a key parameter to determine transport efficiency. In this study, we present the electrokinetic effects on solute dispersion in porous media at the pore scale, using a pore network model. The analytical solution of the electrokinetic coupling coefficient was obtained to quantity the fluid flow velocity in a cylinder capillary. The effect of electrical double layer on the electrokinetic coupling coefficient was investigated by applying different ionic concentration. By averaging the velocity over cross section within a single pore, the average flux was obtained. Applying such single pore relationships, in the thin electrical double layer limit, to each and every pore within the pore network, potential distribution and the induced fluid flow was calculated for the whole domain. The resulting pore velocities were used to simulate solute transport within the pore network. By averaging the results, we obtained the breakthrough curve (BTC) of the average concentration at the outlet of the pore network. Optimizing the solution of continuum scale advection-dispersion equation to such a BTC, solute dispersion coefficient was estimated. We have compared the dispersion caused by electrokinetic flow and pure pressure driven flow under different Peclet number values. In addition, the effect of microstructure and topological properties of porous media on fluid flow and solute dispersion is presented, mainly based on different pore coordination numbers.
Transport of Organic Solutes in Clay Formations
The research is a pilot investigation for the SERDP (Strategic Environmental Research and Development Program, DoD) founded project, Impact of Clay-DNAPL Interactions on Transport and Storage of Chlorinated Solvents in Low Permeability Zones, from 2010-2012. The report tries to s...
Towards a transport approach that acknowledges mixing and dispersion.
NASA Astrophysics Data System (ADS)
Carrera, J.; soler Sagarra, J.; de Dreuzy, J. R.; Dentz, M.
2014-12-01
It is generally accepted that the Advection-Dispersion Equation (ADE) is a poor representation of transport for problems beyond assessing the extent of a solute plume. Specifically, mixing must be honored for proper assessment of chemical reactions. Therefore, it is necessary to develop a transport approach that acknowledges dispersion (for adequate representation of solute spreading) and mixing (for adequate representation of chemical reactions). Non-local in time solute transport formulations have been considered a hopeful alternative to the ADE because they overcome many of its limitations. We have computed the deviation from gaussian mixing obtained in transport through highly heterogeneous media and compared it with that of non-local in time formulations. We find that these underestimate such deviation. Therefore, they are not sufficient; more sophisticated approaches are needed. An appealing option is to extend non-locality also to space, but this opens a broad range of possibilities. We explore some non-local in space and time formulations, so as to define the constraints that these must meet in order to be valid representations (valid in the sense of reproducing the actual spreading and mixing rates) of solute transport through heterogeneous media.
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.
Future variability of solute transport in a macrotidal estuary
NASA Astrophysics Data System (ADS)
Robins, Peter E.; Lewis, Matt J.; Simpson, John H.; Howlett, Eleanor R.; Malham, Shelagh K.
2014-12-01
The physical controls on salt distribution and river-sourced conservative solutes, including the potential implications of climate change, are investigated referring to model simulations of a macrotidal estuary. In the UK, such estuaries typically react rapidly to rainfall events and, as such, are often in a state of non-equilibrium in terms of solute transport; hence are particularly sensitive to climate extremes. Sea levels are projected to rise over the 21st century, extending the salinity maximum upstream in estuaries, which will also affect downstream solute transport, promoting estuarine trapping and reducing offshore dispersal of material. Predicted 'drier summers' and 'wetter winters' in the UK will influence solute transport further still; we found that projected river flow climate changes were more influential than sea-level rise, especially for low flow conditions. Our simulations show that projected climate change for the UK is likely to increase variability in estuarine solute transport and, specifically, increase the likelihood of estuarine trapping during summer, mainly due to drier weather conditions. Future changes in solute transport were less certain during winter, since increased river flow will to some extent counter-act the effects of sea-level rise. Our results have important implications for non-conservative nutrient transport, water quality, coastal management and ecosystem resilience.
The mass transport (advection-dispersion) equations allowing coupled second-order reaction (i.e. Omega sub 1, C sub 1) + (omega sub 2, C sub 2) (R sub 12) -> Re) between two constituents are derived and result in a set of coupled nonlinear partial differential equations. Neglecti...
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
Effects of Soil Behavior on Solute Transport in Groundwater
NASA Astrophysics Data System (ADS)
Jeng, Dong-Sheng; Zhang, Huijie
2010-05-01
The evaluation of solute transport in groundwater is particularly important for environmental engineers involved in the design of urban environments. In general, the simulation of solute transport in porous medium has been linked with fluid flow, which has commonly based on Darcy law. Unlike previous work, we use a more generalized fluid flow model with poro-elastic theory, in which Darcy model is one of its special cases. The new feature of the new model is the inclusion of soil characteristics and behavior in the prediction of solute transport in aquifers. Based on the new model, numerical example demonstrates significant influence of poro-elastic soil behavior on the movement of zone of peak concentration of solute in groundwater.
Moridis, George J.
2001-10-10
In this paper, semianalytical solutions are developed for the problem of transport of radioactive or reactive solute tracers through a layered system of heterogeneous fractured media with misaligned fractures. The tracer transport equations in the non-flowing matrix account for (a) diffusion, (b) surface diffusion, (c) mass transfer between the mobile and immobile water fractions, (d) linear kinetic or equilibrium physical, chemical, or combined solute sorption or colloid filtration, and (e) radioactive decay or first-order chemical reactions. The tracer-transport equations in the fractures account for the same processes, in addition to advection and hydrodynamic dispersion. Any number of radioactive decay daughter products (or products of a linear, first-order reaction chain) can be tracked. The solutions, which are analytical in the Laplace space, are numerically inverted to provide the solution in time and can accommodate any number of fractured and/or porous layers. The solutions are verified using analytical solutions for limiting cases of solute and colloid transport through fractured and porous media. The effect of important parameters on the transport of {sup 3}H, {sup 237}Np and {sup 239}Pu (and its daughters) is investigated in several test problems involving layered geological systems of varying complexity.
George J. Moridis
2001-10-01
In this paper, semianalytical solutions are developed for the problem of transport of radioactive or reactive solute tracers through a layered system of heterogeneous fractured media with misaligned fractures. The tracer transport equations in the non-flowing matrix account for (a) diffusion, (b) surface diffusion, (c) mass transfer between the mobile and immobile water fractions, (d) linear kinetic or equilibrium physical, chemical, or combined solute sorption or colloid filtration, and (e) radioactive decay or first-order chemical reactions. The tracer-transport equations in the fractures account for the same processes, in addition to advection and hydrodynamic dispersion. Any number of radioactive decay daughter products (or products of a linear, first-order reaction chain) can be tracked. The solutions, which are analytical in the Laplace space, are numerically inverted to provide the solution in time and can accommodate any number of fractured and/or porous layers. The solutions are verified using analytical solutions for limiting cases of solute and colloid transport through fractured and porous media. The effect of important parameters on the transport of {sup 3}H, {sup 237}Np and {sup 239}Pu (and its daughters) is investigated in several test problems involving layered geological systems of varying complexity.
End-Member Formulation of Solid Solutions and Reactive Transport
Lichtner, Peter C.
2015-09-01
A model for incorporating solid solutions into reactive transport equations is presented based on an end-member representation. Reactive transport equations are solved directly for the composition and bulk concentration of the solid solution. Reactions of a solid solution with an aqueous solution are formulated in terms of an overall stoichiometric reaction corresponding to a time-varying composition and exchange reactions, equivalent to reaction end-members. Reaction rates are treated kinetically using a transition state rate law for the overall reaction and a pseudo-kinetic rate law for exchange reactions. The composition of the solid solution at the onset of precipitation is assumed to correspond to the least soluble composition, equivalent to the composition at equilibrium. The stoichiometric saturation determines if the solid solution is super-saturated with respect to the aqueous solution. The method is implemented for a simple prototype batch reactor using Mathematica for a binary solid solution. Finally, the sensitivity of the results on the kinetic rate constant for a binary solid solution is investigated for reaction of an initially stoichiometric solid phase with an undersaturated aqueous solution.
Insights into non-Fickian solute transport in carbonates
Bijeljic, Branko; Mostaghimi, Peyman; Blunt, Martin J
2013-01-01
[1] We study and explain the origin of early breakthrough and long tailing plume behavior by simulating solute transport through 3-D X-ray images of six different carbonate rock samples, representing geological media with a high degree of pore-scale complexity. A Stokes solver is employed to compute the flow field, and the particles are then transported along streamlines to represent advection, while the random walk method is used to model diffusion. We compute the propagators (concentration versus displacement) for a range of Peclet numbers (Pe) and relate it to the velocity distribution obtained directly on the images. There is a very wide distribution of velocity that quantifies the impact of pore structure on transport. In samples with a relatively narrow spread of velocities, transport is characterized by a small immobile concentration peak, representing essentially stagnant portions of the pore space, and a dominant secondary peak of mobile solute moving at approximately the average flow speed. On the other hand, in carbonates with a wider velocity distribution, there is a significant immobile peak concentration and an elongated tail of moving fluid. An increase in Pe, decreasing the relative impact of diffusion, leads to the faster formation of secondary mobile peak(s). This behavior indicates highly anomalous transport. The implications for modeling field-scale transport are discussed. Citation: Bijeljic, B., P. Mostaghimi, and M. J. Blunt (2013), Insights into non-Fickian solute transport in carbonates, Water Resour. Res., 49, 2714–2728, doi:10.1002/wrcr.20238. PMID:24223444
A KINETIC MODEL FOR CELL DENSITY DEPENDENT BACTERIAL TRANSPORT IN POROUS MEDIA
A kinetic transport model with the ability to account for variations in cell density of the aqueous and solid phases was developed for bacteria in porous media. Sorption kinetics in the advective-dispersive-sorptive equation was described by assuming that adsorption was proportio...
Larsen, Erik Hviid; Sørensen, Jakob Balslev; Sørensen, Jens Nørkær
2000-01-01
A mathematical model of an absorbing leaky epithelium is developed for analysis of solute coupled water transport. The non-charged driving solute diffuses into cells and is pumped from cells into the lateral intercellular space (lis). All membranes contain water channels with the solute passing those of tight junction and interspace basement membrane by convection-diffusion. With solute permeability of paracellular pathway large relative to paracellular water flow, the paracellular flux ratio of the solute (influx/outflux) is small (2–4) in agreement with experiments. The virtual solute concentration of fluid emerging from lis is then significantly larger than the concentration in lis. Thus, in absence of external driving forces the model generates isotonic transport provided a component of the solute flux emerging downstream lis is taken up by cells through the serosal membrane and pumped back into lis, i.e., the solute would have to be recirculated. With input variables from toad intestine (Nedergaard, S., E.H. Larsen, and H.H. Ussing, J. Membr. Biol. 168:241–251), computations predict that 60–80% of the pumped flux stems from serosal bath in agreement with the experimental estimate of the recirculation flux. Robust solutions are obtained with realistic concentrations and pressures of lis, and with the following features. Rate of fluid absorption is governed by the solute permeability of mucosal membrane. Maximum fluid flow is governed by density of pumps on lis-membranes. Energetic efficiency increases with hydraulic conductance of the pathway carrying water from mucosal solution into lis. Uphill water transport is accomplished, but with high hydraulic conductance of cell membranes strength of transport is obscured by water flow through cells. Anomalous solvent drag occurs when back flux of water through cells exceeds inward water flux between cells. Molecules moving along the paracellular pathway are driven by a translateral flow of water, i.e., the model
Molecular level water and solute transport in reverse osmosis membranes
NASA Astrophysics Data System (ADS)
Lueptow, Richard M.; Shen, Meng; Keten, Sinan
2015-11-01
The water permeability and rejection characteristics of six solutes, methanol, ethanol, 2-propanol, urea, Na+, and Cl-, were studied for a polymeric reverse osmosis (RO) membrane using non-equilibrium molecular dynamics simulations. Results indicate that water flux increases with an increasing fraction of percolated free volume in the membrane polymer structure. Solute molecules display Brownian motion and hop from pore to pore as they pass through the membrane. The solute rejection depends on both the size of the solute molecule and the chemical interaction of the solute with water and the membrane. When the open spaces in the polymeric structure are such that solutes have to shed at least one water molecule from their solvation shell to pass through the membrane molecular structure, the water-solute pair interaction energy governs solute rejection. Organic solutes more easily shed water molecules than ions to more readily pass through the membrane. Hydrogen-bonding sites for molecules like urea also lead to a higher rejection. These findings underline the importance of the solute's solvation shell and solute-water-membrane chemistry in solute transport and rejection in RO membranes. Funded by the Institute for Sustainability and Energy at Northwestern with computing resources from XSEDE (NSF grant ACI-1053575).
Analytical model of solute transport by unsteady unsaturated gravitational infiltration.
Lessoff, S C; Indelman, P
2004-08-01
Penetration of reactive solute into a soil during a cycle of water infiltration and redistribution is investigated by deriving analytical closed form solutions for fluid flux, moisture content and contaminant concentration. The solution is developed for gravitational flow and advective transport and is applied to two scenarios of solute applications encountered in the applications: a finite pulse of solute dissolved in irrigation water and an instantaneous pulse broadcasted onto the soil surface. Through comparison to simulations of Richards' flow, capillary suction is shown to have contrasting effects on the upper and lower boundaries of the fluid pulse, speeding penetration of the wetting front and reducing the rate of drying. This leads to agreement between the analytical and numerical solutions for typical field and experimental conditions. The analytical solution is further incorporated into a stochastic column model of flow and transport to compute mean solute concentration in a heterogeneous field. An unusual phenomenon of plume contraction is observed at long times of solute propagation during the drying stage. The mean concentration profiles match those of the Monte-Carlo simulations for capillary length scales typical of sandy soils. PMID:15240168
Influence of surfactants on unsaturated water flow and solute transport
NASA Astrophysics Data System (ADS)
Karagunduz, Ahmet; Young, Michael H.; Pennell, Kurt D.
2015-04-01
Surfactants can reduce soil water retention by changing the surface tension of water and the contact angle between the liquid and solid phases. As a result, water flow and solute transport in unsaturated soil may be altered in the presence of surfactants. In this study, the effects of a representative nonionic surfactant, Triton X-100, on coupled water flow and nonreactive solute transport during unsaturated flow conditions were evaluated. Batch reactor experiments were conducted to measure the surfactant sorption characteristics, while unsaturated transport experiments were performed in columns packed with 40-270 mesh Ottawa sand at five initial water contents. Following the introduction of surfactant solution, the rate of water percolation through the sand increased; however, this period of rapid water drainage was followed by decreased water percolation due to the reduction in soil water content and the corresponding decrease in unsaturated hydraulic conductivity behind the surfactant front. The observed changes in water percolation occurred sequentially, and resulted in faster nonreactive solute transport than was observed in the absence of surfactant. A one-dimensional mathematical model accurately described coupled water flow, surfactant, and solute transport under most experimental conditions. Differences between model predictions and experimental data were observed in the column study performed at the lowest water content (0.115 cm3/cm3), which was attributed to surfactant adsorption at the air-water interface. These findings demonstrate the potential influence of surfactants additives on unsaturated water flow and solute transport in soils, and demonstrate a methodology to couple these processes in a predictive modeling tool.
The secret to successful solute-transport modeling.
Konikow, Leonard F
2011-01-01
Modeling subsurface solute transport is difficult-more so than modeling heads and flows. The classical governing equation does not always adequately represent what we see at the field scale. In such cases, commonly used numerical models are solving the wrong equation. Also, the transport equation is hyperbolic where advection is dominant, and parabolic where hydrodynamic dispersion is dominant. No single numerical method works well for all conditions, and for any given complex field problem, where seepage velocity is highly variable, no one method will be optimal everywhere. Although we normally expect a numerically accurate solution to the governing groundwater-flow equation, errors in concentrations from numerical dispersion and/or oscillations may be large in some cases. The accuracy and efficiency of the numerical solution to the solute-transport equation are more sensitive to the numerical method chosen than for typical groundwater-flow problems. However, numerical errors can be kept within acceptable limits if sufficient computational effort is expended. But impractically long simulation times may promote a tendency to ignore or accept numerical errors. One approach to effective solute-transport modeling is to keep the model relatively simple and use it to test and improve conceptual understanding of the system and the problem at hand. It should not be expected that all concentrations observed in the field can be reproduced. Given a knowledgeable analyst, a reasonable description of a hydrogeologic framework, and the availability of solute-concentration data, the secret to successful solute-transport modeling may simply be to lower expectations. PMID:21039449
Transport in disordered media with spatially nonuniform fields.
Scher, Harvey; Willbrand, Karen; Berkowitz, Brian
2010-03-01
The theoretical treatment of transport in a disordered system in the presence of a system-wide force field F(x) or spatially varying macroscopic velocity field v(x) is developed in the framework of continuous time random walk (CTRW). The physical basis of CTRW and related fractional derivative equations relies on a mapping of the aggregate of transition rates w(s,s'), between sites s and s', in the Master equation describing the system kinetics, onto a joint probability distribution function psi(s,t). This distribution is calculated from the ensemble average of a position-dependent functional of w(s,s'); the procedure is effective when the scale of heterogeneities is much smaller than the system size. However, statistical homogeneity does not hold in the presence of large heterogeneities, which control the macroscopic v(x), or in the case of an interaction of F(x) with the transition rates. The transport equation, incorporating large-scale heterogeneity, involves the use of a local ensemble average to obtain a position-dependent psi(s,t;x); this determines a memory function, M(t;x), which is convoluted with the advection-dispersion operator. A prototype transport equation for a system with statistical inhomogeneity is developed as an integrodifferential equation. It is solved numerically for particles migrating with a steady-state Darcy velocity v(x , determined for different permeability fields and boundary conditions. The nature of the solutions as a function of key transport parameters (e.g., a characteristic time tc) is explored, and solutions are also compared to those of the advection-dispersion equation for v(x) and to a laboratory experiment. This transport equation is in contrast to the fractional Fokker-Planck equation, which is based on a decoupling of F(x) or v(x) with the transition rates w(s,s'). Further, an analytic expression for the effect of a variance of the ensemble average on the solution of the CTRW transport equation is derived. PMID:20365692
On diagonalization of coupled hydrologic transport and geochemical reaction equations
Yeh, Gour-Tsyh; Cheng, Hwai-Ping
1996-12-31
Two basic ingredients present in modeling the transport of reactive multi-components: the transport is described by a set of advection-dispersion-reactive partial differential equations (PDEs) based on the principle of mass balance; the chemical reactions, under the assumptions of local equilibrium, are described by a set of highly nonlinear algebraic equations (AEs) base on the principles of mole balance and mass action. For a typical application, the complete set of nonlinear PDEs and AEs consist of more than one hundred simultaneous equations. Thus, it is impractical to solve this set of equations simultaneously. General practice is to divide this set of equations into two subsets: one is the primary governing equations (PGEs) consisting of mainly the transport equations and the other one is the secondary governing equations consisting of mainly the geochemical reaction equations. The PGEs are solved for the chosen primary dependent variables (PDVs) and the SGEs are used to compute for the secondary dependent variables (SDVs). The major difficulties in simulating the reactive transport is the numerical solution of PGEs. From the computational point of view, the solution of the set of highly nonlinear PDEs are solved either with the direct substitution approach (DSA) or with the sequential iteration approach (SIA). For DSA, geochemical equilibrium reaction equations are substituted into the hydrologic transport equations to results in a set of nonlinear partial differential equations.
A fracture network model for water flow and solute transport
Robinson, B.A.
1989-01-01
This paper summarizes code development work and sample calculations for FRACNET, a two-dimensional steady state simulator of fluid flow and solute transport in fractured porous media. The model analyzes flow and transport by generating a fracture network based on statistical characteristics of fractures obtained from well logs and other data. After a network is generated, flow and tracer transport are computed for appropriate boundary conditions and wellbore source/sink terms. In addition, for a given realization, the code can be used to indicate whether the medium can be treated as an equivalent porous medium. 18 refs., 7 figs.
Solute transport along preferential flow paths in unsaturated fractures
Su, G.W.; Geller, J.T.; Pruess, K.; Hunt, J.R.
2001-01-01
Laboratory experiments were conducted to study solute transport along preferential flow paths in unsaturated, inclined fractures. Qualitative aspects of solute transport were identified in a miscible dye tracer experiment conducted in a transparent replica of a natural granite fracture. Additional experiments were conducted to measure the breakthrough curves of a conservative tracer introduced into an established preferential flow path in two different fracture replicas and a rock-replica combination. The influence of gravity was investigated by varying fracture inclination. The relationship between the travel times of the solute and the relative influence of gravity was substantially affected by two modes of intermittent flow that occurred: the snapping rivulet and the pulsating blob modes. The measured travel times of the solute were evaluated with three transfer function models: the axial dispersion, the reactors-in-series, and the lognormal models. The three models described the solute travel times nearly equally well. A mechanistic model was also formulated to describe transport when the pulsating blob mode occurred which assumed blobs of water containing solute mixed with residual pools of water along the flow path.
Dispersion of solutes in porous media
NASA Astrophysics Data System (ADS)
Hunt, A. G.; Skinner, T. E.; Ewing, R. P.; Ghanbarian-Alavijeh, B.
2011-04-01
A recently introduced theory of solute transport in porous media is tested by comparison with experiment. The solute transport is predicted using an adaptation of the cluster statistics of percolation theory to critical path analysis together with knowledge of how the structure of such percolation clusters affects the time of transport across them. Only the effects of a single scale of medium heterogeneity are incorporated, and a minimal amount of information regarding the structure of the medium is required. This framework is used to find effectively the distributions of solute velocities and travel distances and thus generate arrival time distributions. The comparison with experiment focuses on the dispersivity (the ratio of the second to the first moment of the spatial solute distribution). The predictions of the theory in the absence of diffusion are verified by comparing with over 2200 experiments over length scales from a few microns to 100 km. At larger length scales (centimeters on up) about 95% of the data lie within our predicted bounds. At smaller length scales approximately 99.8% of the data lie where we predict. These comparisons are not trivial as the typical values of the dispersivity increase by ten orders of magnitude over ten orders of magnitude of length scale. Noteworthy is that the classical advection-dispersion (ADE) equation predicts that the dispersivity should be independent of length scale! This agreement with experiment requires rethinking of the relevance of diffusion and multi-scale heterogeneity and would also appear to signal the complete inappropriateness of using the classical ADE or any of its derivatives to model solute transport.
Insights into non-Fickian solute transport in carbonates
NASA Astrophysics Data System (ADS)
Bijeljic, Branko; Mostaghimi, Peyman; Blunt, Martin J.
2013-05-01
We study and explain the origin of early breakthrough and long tailing plume behavior by simulating solute transport through 3-D X-ray images of six different carbonate rock samples, representing geological media with a high degree of pore-scale complexity. A Stokes solver is employed to compute the flow field, and the particles are then transported along streamlines to represent advection, while the random walk method is used to model diffusion. We compute the propagators (concentration versus displacement) for a range of Peclet numbers (Pe) and relate it to the velocity distribution obtained directly on the images. There is a very wide distribution of velocity that quantifies the impact of pore structure on transport. In samples with a relatively narrow spread of velocities, transport is characterized by a small immobile concentration peak, representing essentially stagnant portions of the pore space, and a dominant secondary peak of mobile solute moving at approximately the average flow speed. On the other hand, in carbonates with a wider velocity distribution, there is a significant immobile peak concentration and an elongated tail of moving fluid. An increase in Pe, decreasing the relative impact of diffusion, leads to the faster formation of secondary mobile peak(s). This behavior indicates highly anomalous transport. The implications for modeling field-scale transport are discussed.
Effect of proteolytic enzymes on transepithelial solute transport
Niewoehner, D.E.; Sinha, A.A.; Rice, K.; Cadman, S.; Wangensteen, D.
1986-10-01
The effects of proteases on air-space clearance (AC) of small ((/sup 14/C)sucrose, 342 daltons) and large (/sup 125/I-neutral dextran, 70,000 daltons) solutes were studied in isolated, fluid-filled hamster lungs that were perfused in a nonrecirculating system. When instilled into the air spaces, porcine pancreatic elastase (0.1-0.4 mg/ml) and bovine pancreatic trypsin (BPT) (0.5-2.0 mg/ml), but neither Clostridium histolyticum collagenase (5.0 mg/ml) nor phenylmethylsulfonyl fluoride-inactivated BPT caused large increases in the AC of both tracer molecules. BPT-induced solute clearance was further characterized functionally and morphologically. The functional characteristics of solute AC under steady-state conditions did not indicate that transepithelial transport was diffusion-limited. Inhibition by millimolar concentrations of Zn/sup 2 +/ and by lung cooling, along with electron microscopic studies employing horseradish peroxidase as a macromolecule tracer, were consistent with epithelial solute transport by a vesicular mechanism (transcytosis). Solute transport from the interstitial compartment to the lung exterior was shown to occur via two pathways. By unknown mechanisms BPT caused small amounts of water to flow through an incompletely identified, extravascular pathway. In BPT-exposed lungs efflux of /sup 125/I-dextran 70 occurred almost exclusively through this pathway, whereas (/sup 14/C)sucrose was transported to the lung exterior partly through this same pathway and partly through the vasculature. The large differences in the diffusion coefficients of the two tracers may have accounted for these observed patterns of solute efflux from the lung. The possible significance of our findings to the pathogenesis of experimental emphysema are discussed.
Simulation of transportation of low enriched uranium solutions
Hope, E.P.; Ades, M.J.
1996-08-01
A simulation of the transportation by truck of low enriched uranium solutions has been completed for NEPA purposes at the Savannah River Site. The analysis involves three distinct source terms, and establishes the radiological risks of shipment to three possible destinations. Additionally, loading accidents were analyzed to determine the radiological consequences of mishaps during handling and delivery. Source terms were developed from laboratory measurements of chemical samples from low enriched uranium feed materials being stored at SRS facilities, and from manufacturer data on transport containers. The transportation simulations were accomplished over the INTERNET using the DOE TRANSNET system at Sandia National Laboratory. The HIGHWAY 3.3 code was used to analyze routing scenarios, and the RADTRAN 4 code was used to analyze incident free and accident risks of transporting radiological materials. Loading accidents were assessed using the Savannah River Site AXAIR89Q and RELEASE 2 codes.
Modeling reactive geochemical transport of concentrated aqueous solutions
NASA Astrophysics Data System (ADS)
Zhang, Guoxiang; Zheng, Zuoping; Wan, Jiamin
2005-02-01
Aqueous solutions with ionic strength larger than 1 M are usually considered concentrated aqueous solutions. These solutions can be found in some natural systems and are also industrially produced and released into accessible natural environments, and as such, they pose a big environmental problem. Concentrated aqueous solutions have unique thermodynamic and physical properties. They are usually strongly acidic or strongly alkaline, with the ionic strength possibly reaching 30 M or higher. Chemical components in such solutions are incompletely dissociated. The thermodynamic activities of both ionic and molecular species in these solutions are determined by the ionic interactions. In geological media the problem is further complicated by the interactions between the solutions and sediments and rocks. The chemical composition of concentrated aqueous solutions when migrating through the geological media may be drastically altered by these strong fluid-rock interactions. To effectively model reactive transport of concentrated aqueous solutions, we must take into account the ionic interactions. For this purpose we substantially extended an existing reactive transport code, BIO-CORE2D©, by incorporating a Pitzer ion interaction model to calculate the ionic activity. In the present paper, the model and two test cases of the model are briefly introduced. We also simulate a laboratory column experiment in which the leakage of highly alkaline waste fluid stored at Hanford (a U.S. Department of Energy site, located in Washington State) was studied. Our simulation captures the measured pH evolution and indicates that all the reactions controlling the pH evolution, including cation exchanges and mineral dissolution/precipitation, are coupled.
Heterogeneous solute transport in a tile-drained field
NASA Astrophysics Data System (ADS)
Basile, A.; Comegna, A.; Coppola, A.; Hassan, S.; Haikal, M. A.; Kassab, M.; Lamaddalena, N.
2009-04-01
Preferential flow and its diverse attributes: i) macropore flow; ii) fingered flow; iii) funnel flow, cannot be described by a single process hypothesis and are unpredictable from a priori analysis of field characteristics due to the inability of sampling methods to capture minute features triggering such flows. Most solute transport techniques are expensive and require extensive soil disturbance. Moreover, solute transport in heterogeneous porous media cannot always be conceptualized as being either a convective-dispersive or a stochastic-convective process. One approach to predict subsurface leaching could be the coupling of near surface measurements with a generalized transport model. A steady state field tracer experiment was conducted on a tile-drained "Terra Rossa" plot located in Valenzano (Bari - Italy), to test whether TDR BTCs measured 1 m a part along a transect of 40 m can be used in such a way for accurate prediction of tile's BTC. A Generalized Transfer Function (GTF) (Zhang, 2000) was fitted to the observed concentration a three depths for each site along the transect to identify the transfer function parameters. To account for vertical transport in the unsaturated zone and lateral divergence near the tile, these parameters were used in a 2D model (Utermann, 1990) to predict earlier breakthrough of tile flux concentration. The 2D model predictions of the flux concentrations were similar to the observed values, nearly reproducing the channel-like nature of solute flow.
Conservative and reactive solute transport in constructed wetlands
Keefe, S.H.; Barber, L.B.; Runkel, R.L.; Ryan, J.N.; McKnight, Diane M.; Wass, R.D.
2004-01-01
The transport of bromide, a conservative tracer, and rhodamine WT (RWT), a photodegrading tracer, was evaluated in three wastewater-dependent wetlands near Phoenix, Arizona, using a solute transport model with transient storage. Coupled sodium bromide and RWT tracer tests were performed to establish conservative transport and reactive parameters in constructed wetlands with water losses ranging from (1) relatively impermeable (15%), (2) moderately leaky (45%), and (3) significantly leaky (76%). RWT first-order photolysis rates and sorption coefficients were determined from independent field and laboratory experiments. Individual wetland hydraulic profiles influenced the extent of transient storage interaction in stagnant water areas and consequently RWT removal. Solute mixing and transient storage interaction occurred in the impermeable wetland, resulting in 21% RWT mass loss from main channel and storage zone photolysis (10%) and sorption (11%) reactions. Advection and dispersion governed solute transport in the leaky wetland, limiting RWT photolysis removal (1.2%) and favoring main channel sorption (3.6%). The moderately leaky wetland contained islands parallel to flow, producing channel flow and minimizing RWT losses (1.6%).
Long-Term Transport of Cryptosporidium Parvum
NASA Astrophysics Data System (ADS)
Andrea, C.; Harter, T.; Hou, L.; Atwill, E. R.; Packman, A.; Woodrow-Mumford, K.; Maldonado, S.
2005-12-01
The protozoan pathogen Cryptosporidium parvum is a leading cause of waterborne disease. Subsurface transport and filtration in natural and artificial porous media are important components of the environmental pathway of this pathogen. It has been shown that the oocysts of C. parvum show distinct colloidal properties. We conducted a series of laboratory studies on sand columns (column length: 10 cm - 60 cm, flow rates: 0.7 m/d - 30 m/d, ionic strength: 0.01 - 100 mM, filter grain size: 0.2 - 2 mm, various solution chemistry). Breakthrough curves were measured over relatively long time-periods (hundreds to thousands of pore volumes). We show that classic colloid filtration theory is a reasonable tool for predicting the initial breakthrough, but it is inadequate to explain the significant tailing observed in the breakthrough of C. parvum oocyst through sand columns. We discuss the application of the Continuous Time Random Walk approach to account for the strong tailing that was observed in our experiments. The CTRW is generalized transport modeling framework, which includes the classic advection-dispersion equation (ADE), the fractional ADE, and the multi-rate mass transfer model as special cases. Within this conceptual framework, it is possible to distinguish between the contributions of pore-scale geometrical (physical) disorder and of pore-scale physico-chemical heterogeneities (e.g., of the filtration, sorption, desorption processes) to the transport of C. parvum oocysts.
Stochastic analysis of transport of conservative solutes in caisson experiments
Dagan, G.
1995-02-01
The Los Alamos National Laboratory has conducted in the past a series of experiments of transport of conservative and reactive solutes. The experimental setup and the experimental results are presented in a series of reports. The main aim of the experiments was to validate models of transport of solutes in unsaturated flow at the caisson intermediate scale, which is much larger than the one pertaining to laboratory columns. First attempts to analyze the experimental results were by one-dimensional convective-dispersion models. These models could not explain the observed solute breakthrough curves and particularly the large solute dispersion in the caisson effluent Since there were some question marks about the uniformity of water distribution at the caisson top, the transport experiments were repeated under conditions of saturated flow. In these experiments constant heads were applied at the top and the bottom of the caisson and the number of concentration monitoring stations was quadrupled. The analysis of the measurements by the same one-dimensional model indicated clearly that the fitted dispersivity is much larger than the pore-sole dispersivity and that it grows with the distance in an approximately linear fashion. This led to the conclusion, raised before, that transport in the caisson is dominated by heterogeneity effects, i.e. by spatial variability of the material Such effects cannot be captured by traditional one-dimensional models. In order to account for the effect of heterogeneity, the saturated flow experiments have been analyzed by using stochastic transport modeling. The apparent linear growth of dispersivity with distance suggested that the system behaves like a stratified one. Consequently, the model of Dagan and Bresier has been adopted in order to interpret concentration measurements. In this simple model the caisson is viewed as a bundle of columns of different permeabilities, which are characterized by a p.d.f. (probability denasity function).
Generic transport coefficients of a confined electrolyte solution.
Yoshida, Hiroaki; Mizuno, Hideyuki; Kinjo, Tomoyuki; Washizu, Hitoshi; Barrat, Jean-Louis
2014-11-01
Physical parameters characterizing electrokinetic transport in a confined electrolyte solution are reconstructed from the generic transport coefficients obtained within the classical nonequilibrium statistical thermodynamic framework. The electro-osmotic flow, the diffusio-osmotic flow, the osmotic current, as well as the pressure-driven Poiseuille-type flow, the electric conduction, and the ion diffusion are described by this set of transport coefficients. The reconstruction is demonstrated for an aqueous NaCl solution between two parallel charged surfaces with a nanoscale gap, by using the molecular dynamic (MD) simulations. A Green-Kubo approach is employed to evaluate the transport coefficients in the linear-response regime, and the fluxes induced by the pressure, electric, and chemical potential fields are compared with the results of nonequilibrium MD simulations. Using this numerical scheme, the influence of the salt concentration on the transport coefficients is investigated. Anomalous reversal of diffusio-osmotic current, as well as that of electro-osmotic flow, is observed at high surface charge densities and high added-salt concentrations. PMID:25493746
Generic transport coefficients of a confined electrolyte solution
NASA Astrophysics Data System (ADS)
Yoshida, Hiroaki; Mizuno, Hideyuki; Kinjo, Tomoyuki; Washizu, Hitoshi; Barrat, Jean-Louis
2014-11-01
Physical parameters characterizing electrokinetic transport in a confined electrolyte solution are reconstructed from the generic transport coefficients obtained within the classical nonequilibrium statistical thermodynamic framework. The electro-osmotic flow, the diffusio-osmotic flow, the osmotic current, as well as the pressure-driven Poiseuille-type flow, the electric conduction, and the ion diffusion are described by this set of transport coefficients. The reconstruction is demonstrated for an aqueous NaCl solution between two parallel charged surfaces with a nanoscale gap, by using the molecular dynamic (MD) simulations. A Green-Kubo approach is employed to evaluate the transport coefficients in the linear-response regime, and the fluxes induced by the pressure, electric, and chemical potential fields are compared with the results of nonequilibrium MD simulations. Using this numerical scheme, the influence of the salt concentration on the transport coefficients is investigated. Anomalous reversal of diffusio-osmotic current, as well as that of electro-osmotic flow, is observed at high surface charge densities and high added-salt concentrations.
Technology Transfer Automated Retrieval System (TEKTRAN)
Traditionally, a two-parameter partial differential equation has been used to describe the one-dimensional convective-dispersive transport of chemicals in field soils. The parameters in this equation include the dispersion coefficient and a distribution coefficient, the latter accounting for interac...
Scaling and predicting solute transport processes in riverine ecosystems
NASA Astrophysics Data System (ADS)
Gonzalez-Pinzon, R.; Haggerty, R.; Camacho Botero, L. A.
2012-12-01
In the last three decades, research on solute transport and nutrient processing has revealed complex interactions between landscapes and stream ecosystems, and numerous attempts to scale and predict these processes have been primarily limited by the difficulty of measuring and extrapolating hydrodynamic and geomorphic characteristics. We hypothesize that there should be predictable patterns in the way that streams interact with their landscapes, because those interactions are in the form of energy, mass and momentum, which are conservative and interrelated properties. Therefore, despite local hydrogeomorphic characteristics define the actual extent of solute transport processes in a given riverine ecosystem, the physical imprints marked-up in breakthrough curves (BTCs) should have scaling properties. To evaluate our hypothesis we created an extensive database that includes 133 BTCs from conservative tracer experiments conducted under different hydrologic conditions (1 lt/s to 1197 m3/s), different experimental conditions (10s of meters to 10s of kilometers), different geographic positions (South and North America, Europe, Australia, Antarctica), and different types of lotic environments, i.e., urban manmade channels, forested headwater streams, desert-like streams, hyporheic wells, and major rivers. We investigated the existence of patterns in conservative solute transport using a model-independent approach, i.e., temporal moments of the histories of tracer experiments. Our results show that the normalized first absolute moment is correlated with the second and third moments with R2>0.99 for all riverine ecosystems. Most importantly, the first central temporal moment of the distributions (mean travel time) is correlated with the second (variance) with an R2>0.93, and the correlation between the second central moment and the third central moment (skewness) takes the form of the coefficient of skewness (CSK) with an R2>0.98, defining a statistically averaged CSK= 1
Radial reactive solute transport in an aquifer-aquitard system
NASA Astrophysics Data System (ADS)
Wang, Quanrong; Zhan, Hongbin
2013-11-01
Radial reactive transport is investigated in an aquifer-aquitard system considering the important processes such as advection, radial and vertical dispersions for the aquifer, vertical advection and dispersion for the aquitards, and first-order biodegradation or radioactive decay. We solved the coupled governing equations of transport in the aquifer and the aquitards by honoring the continuity of concentration and mass flux across the aquifer-aquitard interfaces and recognizing the concentration variation along the aquifer thickness. This effort improved the averaged-approximation (AA) model, which dealt with radial dispersion in an aquifer-aquitard system by excluding the aquitard advection. To compare with our new solution, we expanded the AA model by including the aquitard advection. The expanded AA model considerably overestimated the mass in the upper aquitard when an upward advection existed there. The rates of mass change in the upper aquitard from the new solution and the AA model solution increased with time following sub-linear fashions. The times corresponding to the peak values of the residence time distributions for the AA model, the expanded AA model, and the new model were almost the same. The residence time distributions seemed to follow the Maxwell-Boltzmann distribution closely when plotting the time in logarithmic scale. In addition, we developed a finite-element COMSOL Multiphysics simulation of the problem, and found that the COMSOL solution agreed with the new solution well.
EFFECTIVE POROSITY IMPLIES EFFECTIVE BULK DENSITY IN SORBING SOLUTE TRANSPORT
Flach, G.
2012-02-27
The concept of an effective porosity is widely used in solute transport modeling to account for the presence of a fraction of the medium that effectively does not influence solute migration, apart from taking up space. This non-participating volume or ineffective porosity plays the same role as the gas phase in single-phase liquid unsaturated transport: it increases pore velocity, which is useful towards reproducing observed solute travel times. The prevalent use of the effective porosity concept is reflected by its prominent inclusion in popular texts, e.g., de Marsily (1986), Fetter (1988, 1993) and Zheng and Bennett (2002). The purpose of this commentary is to point out that proper application of the concept for sorbing solutes requires more than simply reducing porosity while leaving other material properties unchanged. More specifically, effective porosity implies the corresponding need for an effective bulk density in a conventional single-porosity model. The reason is that the designated non-participating volume is composed of both solid and fluid phases, both of which must be neglected for consistency. Said another way, if solute does not enter the ineffective porosity then it also cannot contact the adjoining solid. Conceptually neglecting the fluid portion of the non-participating volume leads to a lower (effective) porosity. Likewise, discarding the solid portion of the non-participating volume inherently leads to a lower or effective bulk density. In the author's experience, practitioners virtually never adjust bulk density when adopting the effective porosity approach.
Some Exact Solutions in Energy Dependent Transport Theory
NASA Astrophysics Data System (ADS)
Williams, M. M. R.
1980-01-01
Some exact solutions are obtained for energy dependent slowing down problems with energy dependent cross sections. The transport equation is solved using the backward-forward model of Fermi. Also studied is the energy dependent diffusion equation. Using these models, and a novel technique involving difference equations, it has been possible to find explicit, and numerically useful, solutions for slowing down from a plane, monoenergetic source in an infinite medium. The slowing down density and the energy deposition function are obtained which are of value in reactor physics and radiation damage calculations.
Upscaling transport with mass transfer models: Mean behavior and propagation of uncertainty
NASA Astrophysics Data System (ADS)
Fernã Ndez-Garcia, D.; Llerar-Meza, G.; Gómez-HernáNdez, J. Jaime
2009-10-01
The choice of an adequate large-scale conceptual transport model constitutes a major challenge associated with the upscaling of solute transport. Among the different alternatives to the classical advection-dispersion model, the (multirate) mass transfer model has been proposed as a valuable and convenient alternative to model the large-scale behavior of solute transport. This paper evaluates the use of mass transfer models as a constitutive equation for upscaling solute transport. To achieve this, we compare Monte Carlo simulations of solute transport at two different support scales. Transport simulations performed at the smallest scale represent a set of reference transport solutions described at a high resolution, which are contrasted against transport simulations obtained using an upscaled model (low resolution). Several formulations of the multirate mass transfer model, which differ in the type of memory function (single rate, double rate, and truncated power law), are used as a constitutive transport equation. The large-scale scenario represents a simplified model obtained by partially homogenizing the reference solution. Results show that the double-rate and the truncated power law mass transfer models are capable of properly describing the ensemble average behavior of the main features associated with the integrated breakthrough curves. However, the uncertainty associated with the upscaled mass transfer models was substantially smaller than that attributed to the reference solution. Importantly, the cumulative distribution function of concentrations associated with the upscaled model follows a distribution similar to the reference solution but with smaller statistical dispersion. The reason is that while appropriate memory functions can be used to preserve the residence time distribution of mass particles during upscaling, the lack of memory in space prevents the model from reproducing mass fluxes in all directions. Specifically, the reproduction of mass
Water and chloride transport in a fine-textured soil in a feedlot pen
NASA Astrophysics Data System (ADS)
Veizaga, E. A.; Rodríguez, L.; Ocampo, C. J.
2015-11-01
Cattle feeding in feedlot pens produces large amounts of manure and animal urine. Manure solutions resulting from surface runoff are composed of numerous chemical constituents whose leaching causes salinization of the soil profile. There is a relatively large number of studies on preferential flow characterization and modeling in clayed soils. However, research on water flow and solute transport derived from cattle feeding operations in fine-textured soils under naturally occurring precipitation events is less frequent. A field monitoring and modeling investigation was conducted at two plots on a fine-textured soil near a feedlot pen in Argentina to assess the potential of solute leaching into the soil profile. Soil pressure head and chloride concentration of the soil solution were used in combination with HYDRUS-1D numerical model to simulate water flow and chloride transport resorting to the concept of mobile/immobile-MIM water for solute transport. Pressure head sensors located at different depths registered a rapid response to precipitation suggesting the occurrence of preferential flow-paths for infiltrating water. Cracks and small fissures were documented at the field site where the % silt and % clay combined is around 94%. Chloride content increased with depth for various soil pressure head conditions, although a dilution process was observed as precipitation increased. The MIM approach improved numerical results at one of the tested sites where the development of cracks and macropores is likely, obtaining a more dynamic response in comparison with the advection-dispersion equation.
Water and chloride transport in a fine-textured soil in a feedlot pen.
Veizaga, E A; Rodríguez, L; Ocampo, C J
2015-11-01
Cattle feeding in feedlot pens produces large amounts of manure and animal urine. Manure solutions resulting from surface runoff are composed of numerous chemical constituents whose leaching causes salinization of the soil profile. There is a relatively large number of studies on preferential flow characterization and modeling in clayed soils. However, research on water flow and solute transport derived from cattle feeding operations in fine-textured soils under naturally occurring precipitation events is less frequent. A field monitoring and modeling investigation was conducted at two plots on a fine-textured soil near a feedlot pen in Argentina to assess the potential of solute leaching into the soil profile. Soil pressure head and chloride concentration of the soil solution were used in combination with HYDRUS-1D numerical model to simulate water flow and chloride transport resorting to the concept of mobile/immobile-MIM water for solute transport. Pressure head sensors located at different depths registered a rapid response to precipitation suggesting the occurrence of preferential flow-paths for infiltrating water. Cracks and small fissures were documented at the field site where the % silt and % clay combined is around 94%. Chloride content increased with depth for various soil pressure head conditions, although a dilution process was observed as precipitation increased. The MIM approach improved numerical results at one of the tested sites where the development of cracks and macropores is likely, obtaining a more dynamic response in comparison with the advection-dispersion equation. PMID:26348833
Numerical error in groundwater flow and solute transport simulation
NASA Astrophysics Data System (ADS)
Woods, Juliette A.; Teubner, Michael D.; Simmons, Craig T.; Narayan, Kumar A.
2003-06-01
Models of groundwater flow and solute transport may be affected by numerical error, leading to quantitative and qualitative changes in behavior. In this paper we compare and combine three methods of assessing the extent of numerical error: grid refinement, mathematical analysis, and benchmark test problems. In particular, we assess the popular solute transport code SUTRA [Voss, 1984] as being a typical finite element code. Our numerical analysis suggests that SUTRA incorporates a numerical dispersion error and that its mass-lumped numerical scheme increases the numerical error. This is confirmed using a Gaussian test problem. A modified SUTRA code, in which the numerical dispersion is calculated and subtracted, produces better results. The much more challenging Elder problem [Elder, 1967; Voss and Souza, 1987] is then considered. Calculation of its numerical dispersion coefficients and numerical stability show that the Elder problem is prone to error. We confirm that Elder problem results are extremely sensitive to the simulation method used.
Long-term tritium transport through field-scale compacted soil liner
Toupiol, C.; Willingham, T.W.; Valocchi, A.J.; Werth, C.J.; Krapac, I.G.; Stark, T.D.; Daniel, D.E.
2002-01-01
A 13-year study of tritium transport through a field-scale earthen liner was conducted by the Illinois State Geological Survey to determine the long-term performance of compacted soil liners in limiting chemical transport. Two field-sampling procedures (pressure-vacuum lysimeter and core sampling) were used to determine the vertical tritium concentration profiles at different times and locations within the liner. Profiles determined by the two methods were similar and consistent. Analyses of the concentration profiles showed that the tritium concentration was relatively uniformly distributed horizontally at each sampling depth within the liner and thus there was no apparent preferential transport. A simple one-dimensional analytical solution to the advective-dispersive solute transport equation was used to model tritium transport through the liner. Modeling results showed that diffusion was the dominant contaminant transport mechanism. The measured tritium concentration profiles were accurately modeled with an effective diffusion coefficient of 6 ?? 10-4 mm2/s, which is in the middle of the range of values reported in the literature.
Model prediction uncertainty of bromide and pesticides transport in laboratory column
NASA Astrophysics Data System (ADS)
Dusek, Jaromir; Dohnal, Michal; Snehota, Michal; Sobotkova, Martina; Ray, Chittaranjan; Vogel, Tomas
2016-04-01
Knowledge of transport parameters of reactive solutes such as pesticides is a prerequisite for reliable predictions of their fate and transport in soil porous systems. Water flow and transport of bromide tracer and five pesticides (atrazine, imazaquin, sulfometuron methyl, S-metolachlor, and imidacloprid) through an undisturbed soil column of tropical Oxisol were analyzed using a one-dimensional numerical model. Laboratory column leaching experiment with three flow interruptions was conducted. The applied numerical model is based on Richards' equation for solving water flow and the advection-dispersion equation for solving solute transport. A global optimization method was used to evaluate the model's sensitivity to transport parameters and the uncertainty of model predictions. Within the Monte Carlo modeling framework, multiple forward simulations searching through the parametric space, were executed to describe the observed breakthrough curves. All pesticides were found to be relatively mobile. Experimental data indicated significant non-conservative behavior of bromide tracer. All pesticides, with the exception of imidacloprid, were found less persistent. Three of the five pesticides (atrazine, sulfometuron methyl, and S-metolachlor) were better described by the linear kinetic sorption model, while the breakthrough curves of imazaquin and imidacloprid were more appropriately approximated using nonlinear instantaneous sorption. Sensitivity analysis suggested that the model is most sensitive to sorption distribution coefficient. The prediction limits contained most of the measured points of the experimental breakthrough curves, indicating adequate model concept and model structure for the description of transport processes in the soil column under study.
Nick, H M; Paluszny, A; Blunt, M J; Matthai, S K
2011-11-01
A second order in space accurate implicit scheme for time-dependent advection-dispersion equations and a discrete fracture propagation model are employed to model solute transport in porous media. We study the impact of the fractures on mass transport and dispersion. To model flow and transport, pressure and transport equations are integrated using a finite-element, node-centered finite-volume approach. Fracture geometries are incrementally developed from a random distributions of material flaws using an adoptive geomechanical finite-element model that also produces fracture aperture distributions. This quasistatic propagation assumes a linear elastic rock matrix, and crack propagation is governed by a subcritical crack growth failure criterion. Fracture propagation, intersection, and closure are handled geometrically. The flow and transport simulations are separately conducted for a range of fracture densities that are generated by the geomechanical finite-element model. These computations show that the most influential parameters for solute transport in fractured porous media are as follows: fracture density and fracture-matrix flux ratio that is influenced by matrix permeability. Using an equivalent fracture aperture size, computed on the basis of equivalent permeability of the system, we also obtain an acceptable prediction of the macrodispersion of poorly interconnected fracture networks. The results hold for fractures at relatively low density. PMID:22181492
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.
In dealing with problems related to land-based nuclear waste management, a number of analytical and approximate solutions were developed to quantify radionuclide transport through fractures contained in the porous formation. t has been reported that by treating the radioactive de...
SOLUTIONS APPROXIMATING SOLUTE TRANSPORT IN A LEAKY AQUIFER RECEIVING WASTEWATER INJECTION
A mathematical model amenable to analytical solution techniques is developed for the investigation of contaminant transport from an injection well into a leaky aquifer system, which comprises a pumped and an unpumped aquifer connected to each other by an aquitard. A steady state ...
SOLUTIONS APPROXIMATING SOLUTE TRANSPORT IN A LEAKY AQUIFER RECEIVING WASTEWATER INJECTION
A mathematical model amenable to analytical solution techniques is developed for the investigation of contaminant transport from an injection well into a leaky aquifer system, which comprises a pumped and an unpumped aquifer connected to each other by an aquitard. teady state gro...
Modeling Multi-process Transport of Pathogens in Porous Media
NASA Astrophysics Data System (ADS)
Cheng, L.; Brusseau, M. L.
2004-12-01
The transport behavior of microorganisms in porous media is of interest with regard to the fate of pathogens associated with wastewater recharge, riverbank filtration, and land application of biosolids. This interest has fomented research on the transport of pathogens in the subsurface environment. The factors influencing pathogen transport within the subsurface environment include advection, dispersion, filtration, and inactivation. The filtration process, which mediates the magnitude and rate of pathogen retention, comprises several mechanisms such as attachment to porous-medium surfaces, straining, and sedimentation. We present a mathematical model wherein individual filtration mechanisms are explicitly incorporated along with advection, dispersion, and inactivation. The performance of the model is evaluated by applying it to several data sets obtained from miscible-displacement experiments conducted using various pathogens. Input parameters are obtained to the extent possible from independent means.
Hydrophilic solute transport across the rat blood-brain barrier
Lucchesi, K.J.
1987-01-01
Brain capillary permeability-surface area products (PS) of hydrophilic solutes ranging in size from 180 to 5,500 Daltons were measured in rats according to the method of Ohno, Pettigrew and Rapoport. The distribution volume of 70 KD dextran at 10 minutes after i.v. injection was also measured to determine the residual volume of blood in brain tissue at the time of sacrifice. Small test solutes were injected in pairs in order to elucidate whether their transfer into the brain proceeds by diffusion through water- or lipid-filled channels or by vesicular transport. This issue was examined in rats whose blood-brain barrier (BBB) was presumed to be intact (untreated) and in rats that received intracarotid infusions to open the BBB (isosmotic salt (ISS) and hyperosmolar arabinose). Ohno PS values of {sup 3}H-inulin and {sup 14}C-L-glucose in untreated rats were found to decrease as the labelling time was lengthened. This was evidence that a rapidly equilibrating compartment exists between blood and brain that renders the Ohno two-compartment model inadequate for computing true transfer rate constants. When the data were reanalyzed using a multi-compartment graphical analysis, solutes with different molecular radii were found to enter the brain at approximately equal rates. Furthermore, unidirectional transport is likely to be initiated by solute adsorption to a glycocalyx coat on the luminal surface of brain capillary endothelium. Apparently, more inulin than L-glucose was adsorbed, which may account for its slightly faster transfer across the BBB. After rats were treated with intracarotid infusions of ISS or hyperosmolar arabinose, solute PS values were significantly increased, but the ratio of PS for each of the solute pairs approached that of their free-diffusion coefficients.
Ground-water solute transport with hydrogeochemical reactions
Garcia-Delgado, R.A.; Koussis, A.D.
1997-03-01
Chemical contamination of ground water is typically associated with multicomponent solutions of reactive substances, the mobility of which is affected by their reactivity. In predicting geochemical transport, it is therefore important that the liquid and solid phase reactions be modeled, along with the flow-controlled processes. This demanding task is typically carried out on powerful computers. Frequently, however, field data are available for a limited number of species, or, a small number of species suffices to characterize ground-water quality. In such cases it is desirable to be able to model the transport on a widely available class of inexpensive computers. The authors report on the development of a 2-D model for the transport of reactive species that runs efficiently on PCs. The model follows a modified one-step procedure that adopts total (aqueous and adsorbed) concentrations and aqueous concentrations of components, and accounts for aqueous complexation and for competitive sorption via isotherms or selectivity coefficients. The use of principal directions of transport coordinates, dimensional splitting, and a specialized algorithm for handling advection-dominated transport render it compact and efficient. Mass conservation is satisfied with high accuracy.
Biotic controls on solute distribution and transport in headwater catchments
NASA Astrophysics Data System (ADS)
Herndon, E. M.; Dere, A. L.; Sullivan, P. L.; Norris, D.; Reynolds, B.; Brantley, S. L.
2015-01-01
Solute concentrations in stream water vary with discharge in patterns that record complex feedbacks between hydrologic and biogeochemical processes. In a comparison of headwater catchments underlain by shale in Pennsylvania, USA (Shale Hills) and Wales, UK (Plynlimon), dissimilar concentration-discharge behaviors are best explained by contrasting landscape distributions of soil solution chemistry - especially dissolved organic carbon (DOC) - that have been established by patterns of vegetation. Specifically, elements that are concentrated in organic-rich soils due to biotic cycling (Mn, Ca, K) or that form strong complexes with DOC (Fe, Al) are spatially heterogeneous in pore waters because organic matter is heterogeneously distributed across the catchments. These solutes exhibit non-chemostatic "bioactive" behavior in the streams, and solute concentrations either decrease (Shale Hills) or increase (Plynlimon) with increasing discharge. In contrast, solutes that are concentrated in soil minerals and form only weak complexes with DOC (Na, Mg, Si) are spatially homogeneous in pore waters across each catchment. These solutes are chemostatic in that their stream concentrations vary little with stream discharge, likely because these solutes are released quickly from exchange sites in the soils during rainfall events. Differences in the hydrologic connectivity of organic-rich soils to the stream drive differences in concentration behavior between catchments. As such, in catchments where soil organic matter (SOM) is dominantly in lowlands (e.g., Shale Hills), bioactive elements are released to the stream early during rainfall events, whereas in catchments where SOM is dominantly in uplands (e.g., Plynlimon), bioactive elements are released later during rainfall events. The distribution of vegetation and SOM across the landscape is thus a key component for predictive models of solute transport in headwater catchments.
Analytical solution of two-dimensional solute transport in an aquifer-aquitard system.
Zhan, Hongbin; Wen, Zhang; Huang, Guanhua; Sun, Dongmin
2009-07-21
This study deals with two-dimensional solute transport in an aquifer-aquitard system by maintaining rigorous mass conservation at the aquifer-aquitard interface. Advection, longitudinal dispersion, and transverse vertical dispersion are considered in the aquifer. Vertical advection and diffusion are considered in the aquitards. The first-type and the third-type boundary conditions are considered in the aquifer. This study differs from the commonly used averaged approximation (AA) method that treats the mass flux between the aquifer and aquitard as an averaged volumetric source/sink term in the governing equation of transport in the aquifer. Analytical solutions of concentrations in the aquitards and aquifer and mass transported between the aquifer and upper or lower aquitard are obtained in the Laplace domain, and are subsequently inverted numerically to yield results in the real time domain (the Zhan method). The breakthrough curves (BTCs) and distribution profiles in the aquifer obtained in this study are drastically different from those obtained using the AA method. Comparison of the numerical simulation using the model MT3DMS and the Zhan method indicates that the numerical result differs from that of the Zhan method for an asymmetric case when aquitard advections are at the same direction. The AA method overestimates the mass transported into the upper aquitard when an upward advection exists in the upper aquitard. The mass transported between the aquifer and the aquitard is sensitive to the aquitard Peclet number, but less sensitive to the aquitard diffusion coefficient. PMID:19477033
JOVIAN STRATOSPHERE AS A CHEMICAL TRANSPORT SYSTEM: BENCHMARK ANALYTICAL SOLUTIONS
Zhang Xi; Shia Runlie; Yung, Yuk L.
2013-04-20
We systematically investigated the solvable analytical benchmark cases in both one- and two-dimensional (1D and 2D) chemical-advective-diffusive systems. We use the stratosphere of Jupiter as an example but the results can be applied to other planetary atmospheres and exoplanetary atmospheres. In the 1D system, we show that CH{sub 4} and C{sub 2}H{sub 6} are mainly in diffusive equilibrium, and the C{sub 2}H{sub 2} profile can be approximated by modified Bessel functions. In the 2D system in the meridional plane, analytical solutions for two typical circulation patterns are derived. Simple tracer transport modeling demonstrates that the distribution of a short-lived species (such as C{sub 2}H{sub 2}) is dominated by the local chemical sources and sinks, while that of a long-lived species (such as C{sub 2}H{sub 6}) is significantly influenced by the circulation pattern. We find that an equator-to-pole circulation could qualitatively explain the Cassini observations, but a pure diffusive transport process could not. For slowly rotating planets like the close-in extrasolar planets, the interaction between the advection by the zonal wind and chemistry might cause a phase lag between the final tracer distribution and the original source distribution. The numerical simulation results from the 2D Caltech/JPL chemistry-transport model agree well with the analytical solutions for various cases.
Konikow, L.F.; Hornberger, G.Z.
2003-01-01
This report describes modifications to a U.S. Geological Survey (USGS) threedimensional solute-transport model (MODFLOWGWT), which is incorporated into the USGS MODFLOW ground-water model as the Ground- Water Transport (GWT) Process. The modifications improve the capability of MODFLOW-GWT to accurately simulate solute transport in simulations that represent a nonzero flux across an aquifer boundary. In such situations, the new Boundary Flux Package (BFLX) will allow the user flexibility to assign the flux to specific cell faces, although that flexibility is limited for certain types of fluxes (such as recharge and evapotranspiration, which can only be assigned to the top face if either is to be represented as a boundary flux). The approach is consistent with that used in the MODPATH model. The application of the BFLX Package was illustrated using a test case in which the Lake Package was active. The results using the BFLX Package showed noticeably higher magnitudes of velocity in the cells adjacent to the lake than previous results without the BFLX Package. Consequently, solute was transported slightly faster through the lake-aquifer system when the BFLX Package is active. However, the overall solute distributions did not differ greatly from simulations made without using the BFLX Package.
James, Andrew I.; Jawitz, James W.; Munoz-Carpena, Rafael
2009-01-01
A model to simulate transport of materials in surface water and ground water has been developed to numerically approximate solutions to the advection-dispersion equation. This model, known as the Transport and Reaction Simulation Engine (TaRSE), uses an algorithm that incorporates a time-splitting technique where the advective part of the equation is solved separately from the dispersive part. An explicit finite-volume Godunov method is used to approximate the advective part, while a mixed-finite element technique is used to approximate the dispersive part. The dispersive part uses an implicit discretization, which allows it to run stably with a larger time step than the explicit advective step. The potential exists to develop algorithms that run several advective steps, and then one dispersive step that encompasses the time interval of the advective steps. Because the dispersive step is computationally most expensive, schemes can be implemented that are more computationally efficient than non-time-split algorithms. This technique enables scientists to solve problems with high grid Peclet numbers, such as transport problems with sharp solute fronts, without spurious oscillations in the numerical approximation to the solution and with virtually no artificial diffusion.
Reactive transport in porous media: a comparison of model prediction with laboratory visualization.
Gramling, Carolyn M; Harvey, Charles F; Meigs, Lucy C
2002-06-01
Groundwater transport models that accurately describe spreading of nonreactive solutes in an aquifer can poorly predict concentrations of reactive solutes. The dispersive term in the advection-dispersion equation can overpredict pore-scale mixing, and thereby overpredict homogeneous chemical reaction. We quantified this experimentally by imaging instantaneous colorimetric reactions between solutions of aqueous CuSO4 and EDTA4- within a 30-cm long translucent chamber packed with cryolite sand that closely matched the optical index of refraction of water. A charge-coupled device camera was used to quantify concentrations of blue CuEDTA2- within the chamber as it was produced by mixing of the two reactants at different flow rates. We compared these experimental results with a new analytic solution for instantaneous bimolecular reaction coupled with advection and dispersion of the product and reactants. For all flow rates, the concentrations of CuEDTA2- recorded in the experiments were about 20% less than predicted by the analytic solution, thereby demonstrating that models assuming complete mixing at the pore scale can overpredict reaction during transport. PMID:12075812
NASA Astrophysics Data System (ADS)
Rogiers, Bart
2015-04-01
Since a few years, an increasing number of contributed R packages is becoming available, in the field of hydrology. Hydrological time series analysis packages, lumped conceptual rainfall-runoff models, distributed hydrological models, weather generators, and different calibration and uncertainty estimation methods are all available. Also a few packages are available for solving partial differential equations. Subsurface hydrological modelling is however still seldomly performed in R, or with codes interfaced with R, despite the fact that excellent geostatistical packages, model calibration/inversion options and state-of-the-art visualization libraries are available. Moreover, other popular scientific programming languages like matlab and python have packages for pre- and post-processing files of MODFLOW (Harbaugh 2005) and MT3DMS (Zheng 2010) models. To fill this gap, we present here the development versions of the RMODFLOW and RMT3DMS packages, which allow pre- and post-processing MODFLOW and MT3DMS input and output files from within R. File reading and writing functions are currently available for different packages, and plotting functions are foreseen making use of the ggplot2 package (plotting system based on the grammar of graphics; Wickham 2009). The S3 generic-function object oriented programming style is used for this. An example is provided, making modifications to an existing model, and visualization of the model output. References Harbaugh, A. (2005). MODFLOW-2005: The US Geological Survey Modular Ground-water Model--the Ground-water Flow Process, U.S. Geological Survey Techniques and Methods 6-A16 (p. 253). Wickham, H. (2009). ggplot2: elegant graphics for data analysis. Springer New York, 2009. Zheng, C. (2010). MT3DMS v5.3, a modular three-dimensional multispecies transport model for simulation of advection, dispersion and chemical reactions of contaminants in groundwater systems. Supplemental User's Guide. (p. 56).
Willis, Catherine; Rubin, Jacob
1987-01-01
In this paper we consider examples of chemistry-affected transport processes in porous media. A moving boundary problem which arises during transport with precipitation-dissolution reactions is solved by three different numerical methods. Two of these methods (one explicit and one implicit) are based on an integral formulation of mass balance and lead to an approximation of a weak solution. These methods are compared to a front-tracking scheme. Although the two approaches are conceptually different, the numerical solutions showed good agreement. As the ratio of dispersion to convection decreases, the methods based on the integral formulation become computationally more efficient. Specific reactions were modeled to examine the dependence of the system on the physical and chemical parameters.
Phononic heat transport in the transient regime: An analytic solution
NASA Astrophysics Data System (ADS)
Tuovinen, Riku; Säkkinen, Niko; Karlsson, Daniel; Stefanucci, Gianluca; van Leeuwen, Robert
2016-06-01
We investigate the time-resolved quantum transport properties of phonons in arbitrary harmonic systems connected to phonon baths at different temperatures. We obtain a closed analytic expression of the time-dependent one-particle reduced density matrix by explicitly solving the equations of motion for the nonequilibrium Green's function. This is achieved through a well-controlled approximation of the frequency-dependent bath self-energy. Our result allows for exploring transient oscillations and relaxation times of local heat currents, and correctly reduces to an earlier known result in the steady-state limit. We apply the formalism to atomic chains, and benchmark the validity of the approximation against full numerical solutions of the bosonic Kadanoff-Baym equations for the Green's function. We find good agreement between the analytic and numerical solutions for weak contacts and baths with a wide energy dispersion. We further analyze relaxation times from low to high temperature gradients.
Coarse grained modeling of transport properties in monoclonal antibody solution
NASA Astrophysics Data System (ADS)
Swan, James; Wang, Gang
Monoclonal antibodies and their derivatives represent the fastest growing segment of the bio pharmaceutical industry. For many applications such as novel cancer therapies, high concentration, sub-cutaneous injections of these protein solutions are desired. However, depending on the peptide sequence within the antibody, such high concentration formulations can be too viscous to inject via human derived force alone. Understanding how heterogenous charge distribution and hydrophobicity within the antibodies leads to high viscosities is crucial to their future application. In this talk, we explore a coarse grained computational model of therapeutically relevant monoclonal antibodies that accounts for electrostatic, dispersion and hydrodynamic interactions between suspended antibodies to predict assembly and transport properties in concentrated antibody solutions. We explain the high viscosities observed in many experimental studies of the same biologics.
Modeling solute transport by DLA in soils of northeastern Egypt.
Hamed, Yasser Ahmed; Yasuda, Hiroshi; Persson, Magnus; Berndtsson, Ronny; Wang, Xin-ping
2015-01-01
Arid soils in Egypt display large variability in solute transport properties, causing problems in soil management. To characterize this variability, dye infiltration experiments were conducted on four plots representing three main soil types in northeastern Egypt. The plots represented both cultivated and uncultivated land use. The observed dye patterns displayed a large variability and especially the clay soils indicated a high degree of preferential flow. The loamy sand and sandy soils displayed a more uniform dye distribution indicating more homogeneous soil properties. The observed dye patterns were modeled using a diffusion limited aggregation (DLA) model. The DLA is a random walk model where model parameters can be optimized using genetic algorithms (GA). The DLA model reproduced the observed dye patterns for all soils in an excellent way. The best fit was obtained with a specific combination of directional random walk probabilities Pu, Pd, Pr, and Pl for each plot (correlation 0.97-0.99). To account for soil layers with different hydraulic properties a two layer DLA model was developed. For all plots the Pu (upward random walk probability) was higher for the upper more homogeneous soil layer. The overall results showed that spatial variability resulting from solute transport for the investigated soils can be modeled using a DLA approach. PMID:25790463
Modeling Solute Transport by DLA in Soils of Northeastern Egypt
Hamed, Yasser Ahmed; Yasuda, Hiroshi; Persson, Magnus; Berndtsson, Ronny; Wang, Xin-ping
2015-01-01
Arid soils in Egypt display large variability in solute transport properties, causing problems in soil management. To characterize this variability, dye infiltration experiments were conducted on four plots representing three main soil types in northeastern Egypt. The plots represented both cultivated and uncultivated land use. The observed dye patterns displayed a large variability and especially the clay soils indicated a high degree of preferential flow. The loamy sand and sandy soils displayed a more uniform dye distribution indicating more homogeneous soil properties. The observed dye patterns were modeled using a diffusion limited aggregation (DLA) model. The DLA is a random walk model where model parameters can be optimized using genetic algorithms (GA). The DLA model reproduced the observed dye patterns for all soils in an excellent way. The best fit was obtained with a specific combination of directional random walk probabilities Pu, Pd, Pr, and Pl for each plot (correlation 0.97–0.99). To account for soil layers with different hydraulic properties a two layer DLA model was developed. For all plots the Pu (upward random walk probability) was higher for the upper more homogeneous soil layer. The overall results showed that spatial variability resulting from solute transport for the investigated soils can be modeled using a DLA approach. PMID:25790463
Record setting during dispersive transport in porous media
NASA Astrophysics Data System (ADS)
Edery, Yaniv; Kostinski, Alex; Berkowitz, Brian
2011-08-01
How often does a contaminant ‘particle’ migrating in a porous medium set a distance record, i.e., advance farther from the origin than at all previous time steps? This question is of fundamental importance in characterizing the nature of the leading edge of a contaminant plume as it is transported through an aquifer. It was proven theoretically by Majumdar and Ziff (2008) that, in the 1d case for pure diffusion, record setting of a random walker scales with n1/2, where n is the number of steps, regardless of the length and time distribution of steps. Here, we use numerical simulations, benchmarked against the 1d analytical solution, to extend this result also for pure diffusion in 2d and 3d domains. We then consider transport in the presence of a drift (i.e., advective-dispersive transport), and show that the record-setting pace of random walkers changes abruptly from $\\propto$ n1/2 to $\\propto$ n1. We explore the dependence of the prefactor on the distribution of step length and number of spatial dimensions. The key implication is that when, after a brief transitional period, the scaling regime commences, the maximum distance reached by the leading edge of a migrating contaminant plume scales linearly with n, regardless of the drift magnitude.
Soil properties and preferential solute transport at the field scale
NASA Astrophysics Data System (ADS)
Koestel, J. K.; Luong, N. M.; Nørgaard, T.; Vendelboe, A. L.; Moldrup, P.; Jarvis, N. J.; Lamandé, M.; Iversen, B. V.; Wollesen de Jonge, L.
2012-04-01
An important fraction of water flow and solute transport through soil takes place through preferential flow paths. Although this had been already observed in the nineteenth century, it had been forgotten by the scientific community until it was rediscovered during the 1970s. The awareness of the relevance of preferential flow was broadly re-established in the community by the early 1990s. However, since then, the notion remains widespread among soil scientists that the occurrence and strength of preferential flow cannot be predicted from measurable proxy variables such as soil properties or land management practices (e.g. Beven, K., 1991, Modeling preferential flow - an uncertain future, Preferential Flow, 1-11). In our study, we present evidence that disproves this notion. We evaluated breakthrough curve experiments under a constant irrigation rate of 1 cm/h conducted on 65 soil columns (20 cm diameter and 20 height) which had been sampled from an approximately 1 ha large loamy field-site in Silstrup, Denmark. We show that the holdback factor, which is an indicator for the strength of preferential transport, is strongly correlated to the bulk density, which in turn is correlated to the organic matter content. By applying multiple linear regression in a bootstrapping framework, we could estimate the holdback factor from the bulk density and the very fine sand fraction with a coefficient of determination of 0.65. Our results raise hopes that it is indeed possible to establish pedotransfer functions for soil susceptibility to preferential flow and transport.
Polymer Dynamics Effects on Solute Transport in Hairy Nanoparticle Membranes
NASA Astrophysics Data System (ADS)
Buenning, Eileen; Bilchak, Connor; Durning, Christopher; Benicewicz, Brian; Sokolov, Alexei; Kumar, Sanat
Molecular transport measurements in matrix-free grafted nanoparticle (MFGNP) films have shown remarkable enhancement of permeability and ideal selectivity of small condensable molecules and simple gases over the neat polymer melts and conventional, dispersed nanoparticle composites. Films comprised of covalently-attached poly(methyl acrylate) PMA chains to the surface of 14nm silica particles self-assemble into ordered arrays, and we postulate this structure plays an important role in regulating solute transport. This self-assembly creates interstitial spaces between the nanoparticle cores, which the polymer chains can only fill by stretching. Here we use small-angle neutron scattering (SANS), broadband dielectric spectroscopy (BDS), rheology and temperature-modulated differential scanning calorimetry (TMDSC) to probe polymer chain and segmental dynamics and investigate this hypothesis of chain stretching in MFGNP materials. We found that grafting slows both chain and segmental relaxation, and increases fragility, indicating that the chains are more ``frustrated'' in the grafted systems. We propose that the effects of the chain/surface interactions on chain dynamics leads to an increase in available free volume and thus enhances transport properties in MFGNP systems. Special thanks to the NSF GRFP and the DOE SCGSR programs.
NASA Astrophysics Data System (ADS)
Gimmi, T.; Waber, H. N.; Gautschi, A.; Rübel, A.
2007-04-01
In order to characterize the large-scale transport properties of the Opalinus Clay formation, the pore water isotope composition (δ18O and δ2H) was determined on samples from the deep borehole Benken (northeastern Switzerland) across Jurassic argillaceous rocks. The sequence of claystones and marls, delimited by two aquifers, is located at depth from about 400 to 700 m and exhibits very low hydraulic conductivities (below 10-13 m s-1). The isotope data of the pore water were obtained from core samples by diffusive vapor equilibration, vacuum distillation, and squeezing. Compared with the other methods, vacuum distillation led to too low values. To evaluate the large-scale transport properties of the formation, we performed a series of advective-dispersive model calculations and compared them with the experimental data. In accordance with the hydrogeological history, we varied initial and boundary conditions as well as model parameters. The main results can be summarized as follows: (1) Molecular diffusion to the underlying aquifer can explain the general features of the isotope profiles, (2) no signatures of advective flow could be detected, (3) the evolution time is of the order of 0.5-1 Ma (relying on laboratory diffusion coefficients) with a possible range of about 0.2-2 Ma, which is geologically plausible, and (4) parameters measured on small scales (centimeters or meters and months) are also plausible at the formation scale (tens of meters and millions of years) for the sediments investigated.
Stochastic analysis of a field-scale unsaturated transport experiment
NASA Astrophysics Data System (ADS)
Severino, G.; Comegna, A.; Coppola, A.; Sommella, A.; Santini, A.
2010-10-01
Modelling of field-scale transport of chemicals is of deep interest to public as well as private sectors, and it represents an area of active theoretical research in many environmentally-based disciplines. However, the experimental data needed to validate field-scale transport models are very limited due to the numerous logistic difficulties that one faces out. In the present paper, the migration of a tracer (Cl -) was monitored during its movement in the unsaturated zone beneath the surface of 8 m × 50 m sandy soil. Under flux-controlled, steady-state water flow ( Jw = 10 mm/day) was achieved by bidaily sprinkler irrigation. A pulse of 105 g/m 2 KCl was applied uniformly to the surface, and subsequently leached downward by the same (chloride-free) flux Jw over the successive two months. Chloride concentration monitoring was carried out in seven measurement campaigns (each one corresponding to a given time) along seven (parallel) transects. The mass recovery was near 100%, therefore underlining the very good-quality of the concentration data-set. The chloride concentrations are used to test two field-scale models of unsaturated transport: (i) the Advection-Dispersion Equation (ADE), which models transport far from the zone of solute entry, and (ii) the Stochastic- Convective Log- normal (CLT) transfer function model, which instead accounts for transport near the release zone. Both the models provided an excellent representation of the solute spreading at z > 0.45 m (being z = 0.45 m the calibration depth). As a consequence, by the depth z ≈ 50 cm one can regard transport as Fickian. The ADE model dramatically underestimates solute spreading at shallow depths. This is due to the boundary effects which are not captured by the ADE. The CLT model appears to be a more robust tool to mimic transport at every depth.
Bencala, K.E.
1984-01-01
Solute transport in streams is determined by the interaction of physical and chemical processes. Data from an injection experiment for chloride and several cations indicate significant influence of solute-streambed processes on transport in a mountain stream. These data are interpreted in terms of transient storage processes for all tracers and sorption processes for the cations. Process parameter values are estimated with simulations based on coupled quasi-two-dimensional transport and first-order mass transfer sorption. Comparative simulations demonstrate the relative roles of the physical and chemical processes in determining solute transport. -from Author
NASA Astrophysics Data System (ADS)
Sittig, S.; Vrugt, J. A.; Kasteel, R.; Groeneweg, J.; Vereecken, H.
2011-12-01
Persistent antibiotics in the soil potentially contaminate the groundwater and affect the quality of drinking water. To improve our understanding of antibiotic transport in soils, we performed laboratory transport experiments in soil columns under constant irrigation conditions with repeated applications of chloride and radio-labeled SDZ. The tracers were incorporated in the first centimeter, either with pig manure or with solution. Breakthrough curves and concentration profiles of the parent compound and the main transformation products were measured. The goal is to describe the observed nonlinear and kinetic transport behavior of SDZ. Our analysis starts with synthetic transport data for the given laboratory flow conditions for tracers which exhibit increasingly complex interactions with the solid phase. This first step is necessary to benchmark our inverse modeling approach for ideal situations. Then we analyze the transport behavior using the column experiments in the laboratory. Our analysis uses a Markov chain Monte Carlo sampler (Differential Evolution Adaptive Metropolis algorithm, DREAM) to efficiently search the parameter space of an advective-dispersion model. Sorption of the antibiotics to the soil was described using a model regarding reversible as well as irreversible sorption. This presentation will discuss our initial findings. We will present the data of our laboratory experiments along with an analysis of parameter uncertainty.
A comparison of solute-transport solution techniques based on inverse modelling results
Mehl, S.; Hill, M.C.
2000-01-01
Five common numerical techniques (finite difference, predictor-corrector, total-variation-diminishing, method-of-characteristics, and modified-method-of-characteristics) were tested using simulations of a controlled conservative tracer-test experiment through a heterogeneous, two-dimensional sand tank. The experimental facility was constructed using randomly distributed homogeneous blocks of five sand types. This experimental model provides an outstanding opportunity to compare the solution techniques because of the heterogeneous hydraulic conductivity distribution of known structure, and the availability of detailed measurements with which to compare simulated concentrations. The present work uses this opportunity to investigate how three common types of results-simulated breakthrough curves, sensitivity analysis, and calibrated parameter values-change in this heterogeneous situation, given the different methods of simulating solute transport. The results show that simulated peak concentrations, even at very fine grid spacings, varied because of different amounts of numerical dispersion. Sensitivity analysis results were robust in that they were independent of the solution technique. They revealed extreme correlation between hydraulic conductivity and porosity, and that the breakthrough curve data did not provide enough information about the dispersivities to estimate individual values for the five sands. However, estimated hydraulic conductivity values are significantly influenced by both the large possible variations in model dispersion and the amount of numerical dispersion present in the solution technique.Five common numerical techniques (finite difference, predictor-corrector, total-variation-diminishing, method-of-characteristics, and modified-method-of-characteristics) were tested using simulations of a controlled conservative tracer-test experiment through a heterogeneous, two-dimensional sand tank. The experimental facility was constructed using randomly
Coupling of solute transport and cell expansion in pea stems
NASA Technical Reports Server (NTRS)
Schmalstig, J. G.; Cosgrove, D. J.
1990-01-01
As cells expand and are displaced through the elongation zone of the epicotyl of etiolated pea (Pisum sativum L. var Alaska) seedlings, there is little net dilution of the cell sap, implying a coordination between cell expansion and solute uptake from the phloem. Using [14C] sucrose as a phloem tracer (applied to the hypogeous cotyledons), the pattern of label accumulation along the stem closely matched the growth rate pattern: high accumulation in the growing zone, little accumulation in nongrowing regions. Several results suggest that a major portion of phloem contents enters elongating cells through the symplast. We propose that the coordination between phloem transport and cell expansion is accomplished via regulatory pathways affecting both plasmodesmata conductivity and cell expansion.
Effects of isotope selection on solution convergence in HZE transport
NASA Technical Reports Server (NTRS)
Kim, Myung-Hee Y.; Wilson, John W.; Kiefer, Richard L.; Thibeault, Sheila A.
1994-01-01
A fragmenting iron ion produces hundreds of isotopes during nuclear reactions. These isotopes are represented in the solution of the transport problem. A reduced set of isotopes is selected to minimize the computational burden but introduces error in the final result. A minimum list of 122 isotopes is required for adequate representation of the mass and charge distributions of the secondary radiation fields. A reduced set of 80 isotopes is sufficient to represent the charge distribution alone and represents reasonably well the linear energy transfer properties of the iron beam. Because iron fragmentation produces nearly every isotope lighter than iron, the resulting 122-isotope list should be adequate for ion beams with charges equal to or less than 26.
Uncertainty and Sensitivity Analyses of Model Predictions of Solute Transport
NASA Astrophysics Data System (ADS)
Skaggs, T. H.; Suarez, D. L.; Goldberg, S. R.
2012-12-01
Soil salinity reduces crop production on about 50% of irrigated lands worldwide. One roadblock to increased use of advanced computer simulation tools for better managing irrigation water and soil salinity is that the models usually do not provide an estimate of the uncertainty in model predictions, which can be substantial. In this work, we investigate methods for putting confidence bounds on HYDRUS-1D simulations of solute leaching in soils. Uncertainties in model parameters estimated with pedotransfer functions are propagated through simulation model predictions using Monte Carlo simulation. Generalized sensitivity analyses indicate which parameters are most significant for quantifying uncertainty. The simulation results are compared with experimentally observed transport variability in a number of large, replicated lysimeters.
An exact solution of solute transport by one-dimensional random velocity fields
Cvetkovic, V.D.; Dagan, G.; Shapiro, A.M.
1991-01-01
The problem of one-dimensional transport of passive solute by a random steady velocity field is investigated. This problem is representative of solute movement in porous media, for example, in vertical flow through a horizontally stratified formation of variable porosity with a constant flux at the soil surface. Relating moments of particle travel time and displacement, exact expressions for the advection and dispersion coefficients in the Focker-Planck equation are compared with the perturbation results for large distances. The first- and second-order approximations for the dispersion coefficient are robust for a lognormal velocity field. The mean Lagrangian velocity is the harmonic mean of the Eulerian velocity for large distances. This is an artifact of one-dimensional flow where the continuity equation provides for a divergence free fluid flux, rather than a divergence free fluid velocity. ?? 1991 Springer-Verlag.
Hierarchical Adaptive Solution of Radiation Transport Problems on Unstructured Grids
Dr. Cassiano R. E de Oliveira
2008-06-30
Computational radiation transport has steadily gained acceptance in the last decade as a viable modeling tool due to the rapid advancements in computer software and hardware technologies. It can be applied for the analysis of a wide range of problems which arise in nuclear reactor physics, medical physics, atmospheric physics, astrophysics and other areas of engineering physics. However, radiation transport is an extremely chanllenging computational problem since the governing equation is seven-deimensional (3 in space, 2 in direction, 1 in energy, and 1 in time) with a high degree of coupleing betwen these variables. If not careful, this relatively large number of independent variables when discretized can potentially lead to sets of linear equations of intractable size. Though parallel computing has allowed the solution of very large problems, avaliable computational resources will always be finite due to the fact that every more sophisticated multiphysics models are being demanded by industry. There is thus the pressing requirement to optimize the discretizations so as to minimize the effort and maximize the accuracy.
NASA Astrophysics Data System (ADS)
Chrysikopoulos, Constantinos V.; Manariotis, Ioannis D.; Syngouna, Vasiliki I.
2014-05-01
Accurate prediction of colloid and biocolloid transport in porous media relies heavily on usage of suitable dispersion coefficients. The widespread procedure for dispersion coefficient determination consists of conducting conservative tracer experiments and subsequently fitting the collected breakthrough data with a selected advection-dispersion transport model. The fitted dispersion coefficient is assumed to characterize the porous medium and is often used thereafter to analyze experimental results obtained from the same porous medium with other solutes, colloids, and biocolloids. The classical advection-dispersion equation implies that Fick's first law of diffusion adequately describes the dispersion process, or that the dispersive flux is proportional to the concentration gradient. Therefore, the above-described procedure inherently assumes that the dispersive flux of all solutes, colloids and biocolloids under the same flow field conditions is exactly the same. Furthermore, the available mathematical models for colloid and biocoloid transport in porous media do not adequately account for gravity effects. Here an extensive laboratory study was undertaken in order to assess whether the dispersivity, which traditionally has been considered to be a property of the porous medium, is dependent on colloid particle size, interstitial velocity and length scale. The breakthrough curves were successfully simulated with a mathematical model describing colloid and biocolloid transport in homogeneous, water saturated porous media. The results demonstrated that the dispersivity increases very slowly with increasing interstitial velocity, and increases with column length. Furthermore, contrary to earlier results, which were based either on just a few experimental observations or experimental conditions leading to low mass recoveries, dispersivity was positively correlated with colloid particle size. Also, transport experiments were performed with biocolloids (bacteriophages:
Improved parallel solution techniques for the integral transport matrix method
Zerr, Robert J; Azmy, Yousry Y
2010-11-23
Alternative solution strategies to the parallel block Jacobi (PBJ) method for the solution of the global problem with the integral transport matrix method operators have been designed and tested. The most straightforward improvement to the Jacobi iterative method is the Gauss-Seidel alternative. The parallel red-black Gauss-Seidel (PGS) algorithm can improve on the number of iterations and reduce work per iteration by applying an alternating red-black color-set to the subdomains and assigning multiple sub-domains per processor. A parallel GMRES(m) method was implemented as an alternative to stationary iterations. Computational results show that the PGS method can improve on the PBJ method execution by up to {approx}50% when eight sub-domains per processor are used. However, compared to traditional source iterations with diffusion synthetic acceleration, it is still approximately an order of magnitude slower. The best-performing case are opticaUy thick because sub-domains decouple, yielding faster convergence. Further tests revealed that 64 sub-domains per processor was the best performing level of sub-domain division. An acceleration technique that improves the convergence rate would greatly improve the ITMM. The GMRES(m) method with a diagonal block preconditioner consumes approximately the same time as the PBJ solver but could be improved by an as yet undeveloped, more efficient preconditioner.
Incorporating Super-Diffusion due to Sub-Grid Heterogeneity to Capture Non-Fickian Transport.
Baeumer, Boris; Zhang, Yong; Schumer, Rina
2015-01-01
Numerical transport models based on the advection-dispersion equation (ADE) are built on the assumption that sub-grid cell transport is Fickian such that dispersive spreading around the average velocity is symmetric and without significant tailing on the front edge of a solute plume. However, anomalous diffusion in the form of super-diffusion due to preferential pathways in an aquifer has been observed in field data, challenging the assumption of Fickian dispersion at the local scale. This study develops a fully Lagrangian method to simulate sub-grid super-diffusion in a multidimensional regional-scale transport model by using a recent mathematical model allowing super-diffusion along the flow direction given by the regional model. Here, the time randomizing procedure known as subordination is applied to flow field output from MODFLOW simulations. Numerical tests check the applicability of the novel method in mapping regional-scale super-diffusive transport conditioned on local properties of multidimensional heterogeneous media. PMID:25214174
NASA Astrophysics Data System (ADS)
Morales-Casique, E.; Lezama-Campos, J. L.; Guadagnini, A.; Neuman, S. P.
2013-05-01
Modeling tracer transport in geologic porous media suffers from the corrupt characterization of the spatial distribution of hydrogeologic properties of the system and the incomplete knowledge of processes governing transport at multiple scales. Representations of transport dynamics based on a Fickian model of the kind considered in the advection-dispersion equation (ADE) fail to capture (a) the temporal variation associated with the rate of spreading of a tracer, and (b) the distribution of early and late arrival times which are often observed in field and/or laboratory scenarios and are considered as the signature of anomalous transport. Elsewhere we have presented exact stochastic moment equations to model tracer transport in randomly heterogeneous aquifers. We have also developed a closure scheme which enables one to provide numerical solutions of such moment equations at different orders of approximations. The resulting (ensemble) average and variance of concentration fields were found to display a good agreement against Monte Carlo - based simulation results for mildly heterogeneous (or well-conditioned strongly heterogeneous) media. Here we explore the ability of the moment equations approach to describe the distribution of early arrival times and late time tailing effects which can be observed in Monte-Carlo based breakthrough curves (BTCs) of the (ensemble) mean concentration. We show that BTCs of mean resident concentration calculated at a fixed space location through higher-order approximations of moment equations display long tailing features of the kind which is typically associated with anomalous transport behavior and are not represented by an ADE model with constant dispersive parameter, such as the zero-order approximation.
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.
Modeling colloid and microorganism transport and release with transients in solution ionic strength
Technology Transfer Automated Retrieval System (TEKTRAN)
The transport and fate of colloids, microorganisms, and nanoparticles in subsurface environments is strongly influenced by transients in solution ionic strength (IS). A sophisticated dual-permeability transport model that is capable of simulating exponential, hyperexponential, uniform, and nonmonot...
Product Lifecycle Management and the Quest for Sustainable Space Transportation Solutions
NASA Technical Reports Server (NTRS)
Caruso, Pamela W.
2009-01-01
This viewgraph presentation reviews NASA Marshall's effort to sustain space transportation solutions through product lines that include: 1) Propulsion and Transportation Systems; 2) Life Support Systems; and 3) and Earth and Space Science Spacecraft Systems, and Operations.
NASA Astrophysics Data System (ADS)
Drummond, J. D.; Covino, T. P.; Aubeneau, A. F.; Leong, D.; Patil, S.; Schumer, R.; Packman, A. I.
2012-09-01
Hydrologic transport and retention strongly affect biogeochemical processes that are critical to stream ecosystems. Tracer injection studies are often used to characterize solute transport and retention in stream reaches, but the range of processes accurately resolved with this approach is not clear. Solute residence time distributions depend on both in-stream mixing and exchange with the hyporheic zone and the larger groundwater system. Observed in-stream breakthrough curves have most commonly been modeled with in-stream advection-dispersion plus an exponential residence time distribution, but process-based models suggest that hyporheic exchange is a fractal process, and that hyporheic residence time distributions are more appropriately characterized by power law tailing. We synthesized results from a variety of tracer-injection studies to investigate the information content of tracer breakthrough curves. We found that breakthrough curve tails are often not well characterized in stream tracer experiments. The two main reasons for this are: 1) experimental truncation of breakthrough curves, which occurs when sampling ends before all tracer mass reaches the sampling location, and 2) sensitivity truncation of breakthrough curves, when tracer concentrations in the tail are too low to be detected reliably above background levels. Tail truncation reduces observed mass recovery and obscures assessment of breakthrough curve tailing and solute residence time. Failure to consider tail truncation leads to underestimation of hyporheic exchange and solute retention and to corresponding overestimation of hyporheic biogeochemical transformation rates. Based on these findings, we propose criteria for improved design of in-stream tracer injection experiments to improve assessment of solute tailing behavior.
Long range transport of colloids in aqueous solutions
NASA Astrophysics Data System (ADS)
Florea, Daniel; Musa, Sami; Huyghe, Jacques M. R. J.; Wyss, Hans M.
2013-03-01
Colloids in aqueous suspensions can experience strong, extremely long range repulsive forces near interfaces such as biological tissues, gels, ion exchange resins or metals. As a result exclusion zones extending over several millimeters can be formed. While this phenomenon has been previously described, a physical understanding of this process is still lacking. This exclusion zone formation is puzzling because the typical forces acting on colloidal particles are limited to much shorter distances and external fields that could drive the particles are absent. Here we study the exclusion zone formation in detail by following the time and distance-dependent forces acting on the particles. We present a simple model that accounts for our experimental data and directly links the exclusion zone formation to an already known physical transport phenomenon. We show that the effect can be tuned by changing the zeta potential of the particles or by varying the species present in the aqueous solution. We thus provide a direct physical explanation for the intriguing exclusion zone formation and we illustrate how this effect can be exploited in a range of industrial applications.
Semianalytical Solutions of Radioactive or Reactive Tracer Transport in Layered Fractured Media
G.J. Moridis; G. S. Bodvarsson
2001-10-01
In this paper, semianalytical solutions are developed for the problem of transport of radioactive or reactive tracers (solutes or colloids) through a layered system of heterogeneous fractured media with misaligned fractures. The tracer transport equations in the matrix account for (a) diffusion, (b) surface diffusion (for solutes only), (c) mass transfer between the mobile and immobile water fractions, (d) linear kinetic or equilibrium physical, chemical, or combined solute sorption or colloid filtration, and (e) radioactive decay or first order chemical reactions. Any number of radioactive decay daughter products (or products of a linear, first-order reaction chain) can be tracked. The tracer-transport equations in the fractures account for the same processes, in addition to advection and hydrodynamic dispersion. Additionally, the colloid transport equations account for straining and velocity adjustments related to the colloidal size. The solutions, which are analytical in the Laplace space, are numerically inverted to provide the solution in time and can accommodate any number of fractured and/or porous layers. The solutions are verified using analytical solutions for limiting cases of solute and colloid transport through fractured and porous media. The effect of important parameters on the transport of {sup 3}H, {sup 237}Np and {sup 239}Pu (and its daughters) is investigated in several test problems involving layered geological systems of varying complexity. {sup 239}Pu colloid transport problems in multilayered systems indicate significant colloid accumulations at straining interfaces but much faster transport of the colloid than the corresponding strongly sorbing solute species.
Semianalytical solutions of radioactive or reactive tracer transport in layered fractured media
Moridis, G.J.; Bodvarsson, G.S.
2001-10-10
In this paper, semianalytical solutions are developed for the problem of transport of radioactive or reactive tracers (solutes or colloids) through a layered system of heterogeneous fractured media with misaligned fractures. The tracer transport equations in the matrix account for (a) diffusion, (b) surface diffusion (for solutes only), (c) mass transfer between the mobile and immobile water fractions, (d) linear kinetic or equilibrium physical, chemical, or combined solute sorption or colloid filtration, and (e) radioactive decay or first order chemical reactions. Any number of radioactive decay daughter products (or products of a linear, first-order reaction chain) can be tracked. The tracer-transport equations in the fractures account for the same processes, in addition to advection and hydrodynamic dispersion. Additionally, the colloid transport equations account for straining and velocity adjustments related to the colloidal size. The solutions, which are analytical in the Laplace space, are numerically inverted to provide the solution in time and can accommodate any number of fractured and/or porous layers. The solutions are verified using analytical solutions for limiting cases of solute and colloid transport through fractured and porous media. The effect of important parameters on the transport of {sup 3}H, {sup 237}Np and {sup 239}Pu (and its daughters) is investigated in several test problems involving layered geological systems of varying complexity. {sup 239}Pu colloid transport problems in multilayered systems indicate significant colloid accumulations at straining interfaces but much faster transport of the colloid than the corresponding strongly sorbing solute species.
A Functional-Phylogenetic Classification System for Transmembrane Solute Transporters
Saier, Milton H.
2000-01-01
A comprehensive classification system for transmembrane molecular transporters has been developed and recently approved by the transport panel of the nomenclature committee of the International Union of Biochemistry and Molecular Biology. This system is based on (i) transporter class and subclass (mode of transport and energy coupling mechanism), (ii) protein phylogenetic family and subfamily, and (iii) substrate specificity. Almost all of the more than 250 identified families of transporters include members that function exclusively in transport. Channels (115 families), secondary active transporters (uniporters, symporters, and antiporters) (78 families), primary active transporters (23 families), group translocators (6 families), and transport proteins of ill-defined function or of unknown mechanism (51 families) constitute distinct categories. Transport mode and energy coupling prove to be relatively immutable characteristics and therefore provide primary bases for classification. Phylogenetic grouping reflects structure, function, mechanism, and often substrate specificity and therefore provides a reliable secondary basis for classification. Substrate specificity and polarity of transport prove to be more readily altered during evolutionary history and therefore provide a tertiary basis for classification. With very few exceptions, a phylogenetic family of transporters includes members that function by a single transport mode and energy coupling mechanism, although a variety of substrates may be transported, sometimes with either inwardly or outwardly directed polarity. In this review, I provide cross-referencing of well-characterized constituent transporters according to (i) transport mode, (ii) energy coupling mechanism, (iii) phylogenetic grouping, and (iv) substrates transported. The structural features and distribution of recognized family members throughout the living world are also evaluated. The tabulations should facilitate familial and functional
Solute transport in eroded and rehabilitated prairie landforms. 1. Nonreactive solute.
Papiernik, Sharon K; Koskinen, William C; Yates, Scott R
2009-08-26
Information regarding solute and water transport as affected by soil properties, topography, and climatic conditions is required to improve and validate transport models. This study evaluated the dissipation of bromide applied to the soil surface in the fall and spring to undisturbed (eroded) and rehabilitated landforms, in which topsoil was moved from depositional areas to the eroded upper slope. Despite large changes in soil properties, the amount and center of mass of bromide remaining in the top 1 m of soil was the same in undisturbed and rehabilitated plots. Approximately 60% of the fall-applied bromide was lost during the winter and early spring, presumably due to leaching and runoff. The center of mass of spring-applied bromide remained at depths of <30 cm. At the end of the experiment, 33% of the spring-applied bromide was detected in soil and 56% in corn plants. These results suggest that little bromide was leached out of the root zone in the spring and that plant uptake was a major route of bromide dissipation during the growing season. PMID:19653694
Waniewski, Jacek; Antosiewicz, Stefan; Baczynski, Daniel; Poleszczuk, Jan; Pietribiasi, Mauro; Lindholm, Bengt; Wankowicz, Zofia
2016-01-01
During peritoneal dialysis (PD), the peritoneal membrane undergoes ageing processes that affect its function. Here we analyzed associations of patient age and dialysis vintage with parameters of peritoneal transport of fluid and solutes, directly measured and estimated based on the pore model, for individual patients. Thirty-three patients (15 females; age 60 (21-87) years; median time on PD 19 (3-100) months) underwent sequential peritoneal equilibration test. Dialysis vintage and patient age did not correlate. Estimation of parameters of the two-pore model of peritoneal transport was performed. The estimated fluid transport parameters, including hydraulic permeability (LpS), fraction of ultrasmall pores (α u), osmotic conductance for glucose (OCG), and peritoneal absorption, were generally independent of solute transport parameters (diffusive mass transport parameters). Fluid transport parameters correlated whereas transport parameters for small solutes and proteins did not correlate with dialysis vintage and patient age. Although LpS and OCG were lower for older patients and those with long dialysis vintage, αu was higher. Thus, fluid transport parameters--rather than solute transport parameters--are linked to dialysis vintage and patient age and should therefore be included when monitoring processes linked to ageing of the peritoneal membrane. PMID:26989432
Waniewski, Jacek; Antosiewicz, Stefan; Baczynski, Daniel; Poleszczuk, Jan; Pietribiasi, Mauro; Lindholm, Bengt; Wankowicz, Zofia
2016-01-01
During peritoneal dialysis (PD), the peritoneal membrane undergoes ageing processes that affect its function. Here we analyzed associations of patient age and dialysis vintage with parameters of peritoneal transport of fluid and solutes, directly measured and estimated based on the pore model, for individual patients. Thirty-three patients (15 females; age 60 (21–87) years; median time on PD 19 (3–100) months) underwent sequential peritoneal equilibration test. Dialysis vintage and patient age did not correlate. Estimation of parameters of the two-pore model of peritoneal transport was performed. The estimated fluid transport parameters, including hydraulic permeability (LpS), fraction of ultrasmall pores (αu), osmotic conductance for glucose (OCG), and peritoneal absorption, were generally independent of solute transport parameters (diffusive mass transport parameters). Fluid transport parameters correlated whereas transport parameters for small solutes and proteins did not correlate with dialysis vintage and patient age. Although LpS and OCG were lower for older patients and those with long dialysis vintage, αu was higher. Thus, fluid transport parameters—rather than solute transport parameters—are linked to dialysis vintage and patient age and should therefore be included when monitoring processes linked to ageing of the peritoneal membrane. PMID:26989432
LI, MING-HSU; SIEGEL, MALCOLM D.; YEH, GOUR-TSYH
1999-09-20
The couplings among chemical reaction rates, advective and diffusive transport in fractured media or soils, and changes in hydraulic properties due to precipitation and dissolution within fractures and in rock matrix are important for both nuclear waste disposal and remediation of contaminated sites. This paper describes the development and application of LEHGC2.0, a mechanistically-based numerical model for simulation of coupled fluid flow and reactive chemical transport including both fast and slow reactions invariably saturated media. Theoretical bases and numerical implementations are summarized, and two example problems are demonstrated. The first example deals with the effect of precipitation-dissolution on fluid flow and matrix diffusion in a two-dimensional fractured media. Because of the precipitation and decreased diffusion of solute from the fracture into the matrix, retardation in the fractured medium is not as large as the case wherein interactions between chemical reactions and transport are not considered. The second example focuses on a complicated but realistic advective-dispersive-reactive transport problem. This example exemplifies the need for innovative numerical algorithms to solve problems involving stiff geochemical reactions.
Hydrogeochemistry and simulated solute transport, Piceance Basin, northwestern Colorado
Robson, S.G.; Saulnier, G.J., Jr.
1981-01-01
Oil-shale mining activities in Piceance basin in northwestern Colorado could adversely affect the ground- and surface-water quality in the basin. This study of the hydrology and geochemistry of the area used ground-water solute-transport-modeling techniques to investigate the possible impact of the mines on water quality. Maps of the extent and structure of the aquifer were prepared and show that a saturated thickness of 2,000 feet occurs in the northeast part of the basin. Ground-water recharge in the upland areas in the east, south, and west parts of the basin moves down into deeper zones in the aquifer and laterally to the discharge areas along Piceance and Yellow Creeks. The saline zone and the unsaturated zone provide the majority of the dissolved solids found in the ground water. Precipitation, ion-exchange, and oxidation-reduction reactions are also occuring in the aquifer. Model simulations of ground-water pumpage in tracts C-a and C-b indicate that the altered direction of ground-water movement near the pumped mines will cause an improvement in ground-water quality near the mines and a degradation of water quality downgradient from the tracts. Model simulations of mine leaching in tract C-a and C-b indicate that equal rates of mine leaching in the tracts will produce much different effects on the water quality in the basin. Tract C-a, by virtue of its remote location from perennial streams, will primarily degrade the ground-water quality over a large area to the northeast of the tract. Tract C-b, by contrast, will primarily degrade the surface-water quality in Piceance Creek, with only localized effects on the ground-water quality. (USGS)
Hydrogeochemistry and simulated solute transport, Piceance Basin, northwestern Colorado
Robson, Stanley G.; Saulnier, George J.
1980-01-01
Oil-shale mining activities in Piceance basin in northwestern Colorado could adversely affect the ground- and surface-water quality in the basin. This study of the hydrology and geochemistry of the area used groundwater solute-transport-modeling techniques to investigate the possible impact of the mines on water quality. Maps of the extent and structure of the aquifer were prepared and show that a saturated thickness of 2,000 feet occurs in the northeast part of the basin. Ground-water recharge in the upland areas in the east, south, and west parts of the basin moves down into deeper zones in the aquifer and laterally to the discharge areas along Piceance and Yellow Creeks. The saline zone and the unsaturated zone provide the majority of the dissolved solids found in the ground water. Precipitation, ion-exchange, and oxidation-reduction reactions are also occurring in the aquifer. Model simulations of groundwater pumpage in tracts C-a and C-b indicate that the altered direction of groundwater movement near the pumped mines will cause an improvement in groundwater quality near the mines and a degradation of water quality downgradient from the tracts. Model simulations of mine leaching in tract C-a and C-b indicate that equal rates of mine leaching in the tracts will produce much different effects on the water quality in the basin. Tract C-a, by virtue of its remote location from perennial streams, will primarily degrade the groundwater quality over a large area to the northeast of the tract. Tract C-b, by contrast, will primarily degrade the surface-water quality in Piceance Creek, with only localized effects on the groundwater quality. (USGS)
NASA Astrophysics Data System (ADS)
Charbonneau, A. M.; Novakowski, K. S.; Ross, N.
2004-05-01
In fractured rock environments, the process of matrix diffusion strongly influences the persistence of groundwater contamination. Biofilm growth, which occurs naturally in fractured rock, can have a significant effect on the solute transport properties of the matrix, and thus the mobility of contaminants in the system. A semi-analytical model that accounts for the presence of a biofilm has been developed from a solution of the advection-dispersion equation for solute transport in radial coordinates. The model is intended to aid in simulating radial diffusion experiments. It describes solute transport by diffusion from a cylindrical reservoir into a layered medium of finite diameter. The model accounts for linear adsorption and decay in each layer, as well as the periodic addition and/or removal of fluid of known solute concentration from the reservoir. The boundary value problem is solved using the Laplace transform method and Crout's method of LU Decomposition, and numerically inverted using the De Hoog algorithm. According to a sensitivity analysis, compared to the no-biofilm case, diffusion of conservative tracers is most sensitive to the depth of penetration of the biofilm into the matrix, as well as the effective porosity and diffusion coefficient of that penetrated zone. The model developed in this study has been used to interpret the results of diffusion experiments undertaken to investigate mass transport into intact rock samples in the presence of a biofilm and to assign average mass transport parameters to the system. Laboratory data demonstrates that the presence of a biofilm acts to limit diffusive transfer between a discrete fracture and the matrix.
Demonstrations in Solute Transport Using Dyes: Part II. Modeling.
ERIC Educational Resources Information Center
Butters, Greg; Bandaranayake, Wije
1993-01-01
A solution of the convection-dispersion equation is used to describe the solute breakthrough curves generated in the demonstrations in the companion paper. Estimation of the best fit model parameters (solute velocity, dispersion, and retardation) is illustrated using the method of moments for an example data set. (Author/MDH)
An exploration of coupled surface-subsurface solute transport in a fully integrated catchment model
NASA Astrophysics Data System (ADS)
Liggett, Jessica E.; Partington, Daniel; Frei, Sven; Werner, Adrian D.; Simmons, Craig T.; Fleckenstein, Jan H.
2015-10-01
Coupling surface and subsurface water flow in fully integrated hydrological codes is becoming common in hydrological research; however, the coupling of surface-subsurface solute transport has received much less attention. Previous studies on fully integrated solute transport focus on small scales, simple geometric domains, and have not utilised many different field data sources. The objective of this study is to demonstrate the inclusion of both flow and solute transport in a 3D, fully integrated catchment model, utilising high resolution observations of dissolved organic carbon (DOC) export from a wetland complex during a rainfall event. A sensitivity analysis is performed to span a range of transport conditions for the surface-subsurface boundary (e.g. advective exchange only, advection plus diffusion, advection plus full mechanical dispersion) and subsurface dispersivities. The catchment model captures some aspects of observed catchment behaviour (e.g. solute discharge at the catchment outlet, increasing discharge from wetlands with increased stream discharge, and counter-clockwise concentration-discharge relationships), although other known behaviours are not well represented in the model (e.g. slope of concentration-discharge plots). Including surface-subsurface solute transport aids in evaluating internal model processes, however there are challenges related to the influence of dispersion across the surface-subsurface interface, and non-uniqueness of the solute transport solution. This highlights that obtaining solute field data is especially important for constraining integrated models of solute transport.
Water flow and solute transport in floating fen root mats
NASA Astrophysics Data System (ADS)
Stofberg, Sija F.; EATM van der Zee, Sjoerd
2015-04-01
be very similar and likely functionally related. Our experimental field data were used for modelling water flow and solute transport in floating fens, using HYDRUS 2D. Fluctuations of surface water and root mat, as well as geometry and unsaturated zone parameters can have a major influence on groundwater fluctuations and the exchange between rain and surface water and the water in the root mats. In combination with the duration of salt pulses in surface water, and sensitivity of fen plants to salinity (Stofberg et al. 2014, submitted), risks for rare plants can be anticipated.
Houseworth, J.E.
2004-09-16
Exact analytical solutions are presented for solute transport in an unsaturated fracture and porous rock matrix. The problem includes advective transport in the fracture and rock matrix as well as advective and diffusive fracture-matrix exchange. Linear sorption in the fracture and matrix and radioactive decay are also treated. The solution is for steady, uniform transport velocities within the fracture and matrix, but allows for independent specification of each of the velocities. The problem is first solved in terms of the solute concentrations that result from an instantaneous point source. Superposition integrals are then used to derive the solute mass flux at a fixed downstream position from an instantaneous point source and for the solute concentrations that result from a continuous point source. Solutions are derived for cases with the solute source in the fracture and the solute source in the matrix. The analytical solutions are closed-form and are expressed in terms of algebraic functions, exponentials, and error functions. Comparisons between the analytical solutions and numerical simulations, as well as sensitivity studies, are presented. Increased sensitivity to cross-flow and solute source location is found for increasing Peclet number. The numerical solutions are found to compare well with the analytical solutions at lower Peclet numbers ,but show greater deviation at higher Peclet numbers.
VERIFICATION OF TRANSPORT CODES BY THE METHOD OF MANUFACTURED SOLUTIONS: THE ATTILA EXPERIENCE
S. D. PAUTZ
2001-03-19
We extend the Method of Manufactured Solutions (MMS) to the verification of transport codes. We derive analytic fixed sources required by the MMS procedure for several types of transport problems and apply the method to the Attila transport code. By means of this method we discover and correct several coding mistakes in Attila and ultimately verify its correct implementation for the problems studied. Our studies reveal that the MMS procedure is a useful tool for transport code development.
Fate and Transport of Graphene Oxide in Granular Porous Media: Experimental Results and Modeling
NASA Astrophysics Data System (ADS)
Gao, Bin
2014-05-01
Although graphene oxide (GO) has been used in many applications to improve human life quality, its environmental fate and behavior are still largely unknown. In this work, a range of laboratory experiments were conducted to explore the aggregation, deposition, and transport mechanisms of GO nano-sheets in porous media under various conditions. Stability experimental data showed that both cation valence and pH showed significant effect on the aggregation of GO sheets. The measured critical coagulation concentrations were in good agreement with the predictions of the extended Schulze-Hardy rule. Sand column experimental results indicated that deposition and transport of GO in porous media were strongly dependent on solution ionic strength. Particularly, GO showed high mobility under low ionic strength conditions in both saturated and unsaturated columns. Increasing ionic strength dramatically increased the retention of GO in porous media, mainly through secondary-minimum deposition. Recovery rates of GO in unsaturated sand columns were lower than that in saturated columns under the same ionic strength conditions, suggesting moisture content also played an important role in the retention of GO in porous media. Findings from the bubble column experiments showed that the GO did not attach to the air-water interface, which is consistent with the XDLVO predictions. Additional retention mechanisms, such as film straining, thus could be responsible to the reduced mobility of GO in unsaturated porous media. The breakthrough curves of GO in saturated and unsaturated columns could be accurately simulated by an advection-dispersion-reaction model.
Smoothed Particle Hydrodynamics Stochastic Model for Flow and Transport in Porous Media
Tartakovsky, Alexandre M.; Tartakovsky, Daniel M.; Meakin, Paul
2008-11-03
A meso-scale stochastic Lagrangian particle model was developed and used to simulate conservative and reactive transport in porous media. In the stochastic model, the fluid flow in a porous continuum is governed by a combination of a Langevin equation and continuity equation. Pore-scale velocity fluctuations, the source of hydrodynamic dispersion, are represented by the white noise. A smoothed particle hydrodynamics method was used to solve the governing equations. Changes in the properties of the fluid particles (e.g., the solute concentration) are governed by the advection-diffusion equation. The separate treatment of advective and diffusive mixing in the stochastic transport model is more realistic than the classical advection-dispersion theory, which uses a single effective diffusion coefficient (the dispersion coefficient) to describe both types of mixing leading to over-prediction of mixing induced effective reaction rates. The stochastic model predicts much lower reaction product concentrations in mixing induced reactions. In addition, the dispersion theory predicts more stable fronts (with a higher effective fractal dimension) than the stochastic model during the growth of Rayleigh-Taylor instabilities.
Numerical solution of the radiation transport equation in disk geometry
NASA Technical Reports Server (NTRS)
Spagna, George F., Jr.; Leung, Chun Ming
1987-01-01
An efficient numerical method for solving the problem of radiation transport in a dusty medium with two dimensional (2-D) disk geometry is described. It is a generalization of the one-dimensional quasi-diffusion method in which the transport equation is cast in diffusion form and then solved as a boundary value problem. The method should be applicable to a variety of astronomical sources, the dynamics of which are angular-momentum dominated and hence not accurately treated by spherical geometry, e.g., protoplanetary nebulae, circumstellar disks, interstellar molecular clouds, accretion disks, and disk galaxies. The computational procedure and practical considerations for implementing the method are described in detail. To illustrate the effects of 2-D radiation transport, some model results (dust temperature distributions and IR flux spectra) for externally heated, interstellar dust clouds with spherically symmetric and disk geometry are compared.
Light-driven solute transport in Halobacterium halobium
NASA Technical Reports Server (NTRS)
Lanyi, J. K.
1979-01-01
The cell membrane of Halobacterium halobium exhibits differential regions which contain crystalline arrays of a single kind of protein, termed bacteriorhodopsin. This bacterial retinal-protein complex resembles the visual pigment and, after the absorption of protons, translocates H(+) across the cell membrane, leading to an electrochemical gradient for protons between the inside and the outside of the cell. Thus, light is an alternate source of energy in these bacteria, in addition to terminal oxidation. The paper deals with work on light-driven transport in H. halobium with cell envelope vesicles. The discussion covers light-driven movements of H(+), Na(+), and K(+); light-driven amino acid transport; and apparent allosteric control of amino acid transport. The scheme of energy coupling in H. halobium vesicles appears simple, its quantitative details are quite complex and reveal regulatory phenomena. More knowledge is required of the way the coupling components are regulated by the ion gradients present.
A dual-porosity model for simulating solute transport in oil shale
Glover, K.C.
1987-01-01
A model is described for simulating three-dimensional groundwater flow and solute transport in oil shale and associated geohydrologic units. The model treats oil shale as a dual-porosity medium by simulating flow and transport within fractures using the finite-element method. Diffusion of solute between fractures and the essentially static water of the shale matrix is simulated by including an analytical solution that acts as a source-sink term to the differential equation of solute transport. While knowledge of fracture orientation and spacing is needed to effectively use the model, it is not necessary to map the locations of individual fractures. The computer program listed in the report incorporates many of the features of previous dual-porosity models while retaining a practical approach to solving field problems. As a result the theory of solute transport is not extended in any appreciable way. The emphasis is on bringing together various aspects of solute transport theory in a manner that is particularly suited to the unusual groundwater flow and solute transport characteristics of oil shale systems. (Author 's abstract)
Influence of spatial and temporal flow variability on solute transport in catchments
NASA Astrophysics Data System (ADS)
Selroos, Jan-Olof; Destouni, Georgia
2015-04-01
The present study quantifies the separate and combined effects of spatial and temporal variability of waterborne solute transport through catchments. The questions addressed are whether, when and why different types of variability may dominate catchment-scale transport. We utilize a versatile numerical solute transport code with a particle-based Monte Carlo time domain random walk method to simulate waterborne transport through a generic catchment. The methodology is exemplified by performing simulations using data on spatiotemporal flow and transport variability from direct stream discharge observations and independently calculated advective solute travel time distributions for catchments within the water management district Northern Baltic Proper (NBP) in Mid-Eastern Sweden. A main conclusion of the study is that projections of catchment mass loading based on spatial variability alone are robust estimates of long-term average solute transport development. This is especially true when annually aggregated mass load rather than finer temporal resolution of mass flux is considered. Temporal variability yields short-term fluctuations around the long-term average solute breakthrough development, and earlier or later arrival than the latter, depending on the timing and duration of solute input relative to the temporal flow variability. The exact temporal characteristics of future solute breakthroughs are thus fundamentally uncertain but their statistical expectation may be well quantified by only spatial variability account.