Measurement and analysis of non-Fickian dispersion in heterogeneous porous media

Contaminant breakthrough behavior in a variety of heterogeneous porous media was measured in laboratory experiments, and evaluated in terms of both the classical advection-dispersion equation (ADE) and the continuous time random walk (CTRW) framework. Heterogeneity can give rise to non-Fickian transport patterns, which are distinguished by "anomalous" early arrival and late time tails in breakthrough curves. Experiments were conducted in two mid-scale laboratory flow cells packed with clean, sieved sand of specified grain sizes. Three sets of experiments were performed, using a "homogeneous" packing, a randomly heterogeneous packing using sand of two grain sizes, and an exponentially correlated structure using sand of three grain sizes. Concentrations of sodium chloride tracer were monitored at the inflow reservoir and measured at the outflow reservoir. Breakthrough curves were then analyzed by comparison to fitted solutions from the ADE and CTRW formulations. In all three systems, including the "homogeneous" one, subtle yet measurable differences between Fickian and non-Fickian transport were observed. Quantitative analysis demonstrated that the CTRW theory characterized the full shape of the breakthrough curves far more effectively than the ADE.     

Modeling solute transport in one-dimensional homogeneous and heterogeneous soil columns with continuous time random walk

In this paper, we used the continuous time random walk (CTRW) framework to characterize the transport process in 1250-cm long one-dimensional homogenous and heterogeneous soil columns at the experiments conducted by Huang et al. [Huang, K., Toride, N., van Genuchten, M.Th., 1995. Experimental investigation of solute transport in large, homogeneous and heterogeneous, saturated soil columns. Trans. Porous Media. 18, 283-302]. The transport process was also simulated by using the advection-dispersion equation (ADE) and the spatial fractional advection-dispersion equation (FADE) for comparison. In the homogeneous soil column, the non-Fickian behavior is found at the distances less than 1000cm with beta values larger than 1.60, but less than 2, and Fickian form transport is obtained at distances larger than 1000cm with beta values larger than 2. In the heterogeneous soil column, we found the most anomalous behavior at distances from 200cm to 700cm with beta values ranging from 0.894 to 0.958, and non-Fickian transport process is observed at distances larger than 800cm with beta values in the range between 1 and 1.3. More significant non-Fickian behavior is found for transport in the heterogeneous soil column than that in the homogeneous soil column. The CTRW fits to the breakthrough curves (BTCs) have lower values of root mean square error (RMSE) and higher values of determination coefficient (r(2)), with respect to the fits of ADE and FADE. The CTRW model also is better captures the full evolution of BTCs, and especially their tails.     

Continuous time random walk model better describes the tailing of atrazine transport in soil

Contaminant transport in soils is complicated and involves some physical and chemical nonequilibrium processes. In this research, the soil column displacement experiments of Cl⁻ and atrazine under different flow velocities were carried out. The data sets of Cl⁻ transport in sandy loam fitted to the convection dispersion equation (CDE) and the two-region model (TRM) indicated that the effects of physical nonequilibrium process produced by immobile water on the breakthrough curves (BTCs) of Cl⁻ and atrazine transport through the repacking soil columns were negligible. The two-site model (TSM) and the continuous time random walk (CTRW) were also used to fit atrazine transport behavior at the flow rate of 19.86 cm h⁻¹. The CTRW convincingly captured the full evolution of atrazine BTC in the soil column, especially for the part of long tailing. However, the TSM failed to characterize the tailing of atrazine BTC in the soil column. The calculated fraction of equilibrium sorption sites, F, ranging from 0.78 to 0.80 for all flow rates suggested the contribution of nonequilibrium sorption sites to the asymmetry of atrazine BTCs. Furthermore, the data sets for the flow rates of 6.68 cm h⁻¹ and 32.81 cm h⁻¹ were predicted by the TSM and the CTRW. As to the flow rate of 6.68 cm h⁻¹, the CTRW predicted the entire BTC of atrazine transport better than the TSM did. For the flow rate of 32.81 cm h⁻¹, the CTRW characterized the late part of the tail better, while the TSM failed to predict the tailings of atrazine BTC.     

Eulerian derivation of the fractional advection-dispersion equation

A fractional advection-dispersion equation (ADE) is a generalization of the classical ADE in which the second-order derivative is replaced with a fractional-order derivative. In contrast to the classical ADE, the fractional ADE has solutions that resemble the highly skewed and heavy-tailed breakthrough curves observed in field and laboratory studies. These solutions, known as alpha-stable distributions, are the result of a generalized central limit theorem which describes the behavior of sums of finite or infinite-variance random variables. We use this limit theorem in a model which sums the length of particle jumps during their random walk through a heterogeneous porous medium. If the length of solute particle jumps is not constrained to a representative elementary volume (REV), dispersive flux is proportional to a fractional derivative. The nature of fractional derivatives is readily visualized and their parameters are based on physical properties that are measurable. When a fractional Fick's law replaces the classical Fick's law in an Eulerian evaluation of solute transport in a porous medium, the result is a fractional ADE. Fractional ADEs are ergodic equations since they occur when a generalized central limit theorem is employed.     

Simplified contaminant source depletion models as analogs of multiphase simulators

Four simplified dense non-aqueous phase liquid (DNAPL) source depletion models recently introduced in the literature are evaluated for the prediction of long-term effects of source depletion under natural gradient flow. These models are simple in form (a power function equation is an example) but are shown here to serve as mathematical analogs to complex multiphase flow and transport simulators. The spill and subsequent dissolution of DNAPLs was simulated in domains having different hydrologic characteristics (variance of the log conductivity field=0.2, 1 and 3) using the multiphase flow and transport simulator UTCHEM. The dissolution profiles were fitted using four analytical models: the equilibrium streamtube model (ESM), the advection dispersion model (ADM), the power law model (PLM) and the Damkohler number model (DaM). All four models, though very different in their conceptualization, include two basic parameters that describe the mean DNAPL mass and the joint variability in the velocity and DNAPL distributions. The variability parameter was observed to be strongly correlated with the variance of the log conductivity field in the ESM and ADM but weakly correlated in the PLM and DaM. The DaM also includes a third parameter that describes the effect of rate-limited dissolution, but here this parameter was held constant as the numerical simulations were found to be insensitive to local-scale mass transfer. All four models were able to emulate the characteristics of the dissolution profiles generated from the complex numerical simulator, but the one-parameter PLM fits were the poorest, especially for the low heterogeneity case.     

Adherence of Sulfur Dioxide Concentrations in the Vicinity of a Steam Plant to Plume Dispersion Models

Air monitoring data for a calendar year at one of the TVA power plants has been used to evaluate the appropriateness of the Sutton, the Bosanquet and Pearson, and the USPHS-TVA atmospheric dispersion models to predict ground level concentrations of sulfur dioxide from emission and meterological data. Aerometric data included one half hourly average sulfur dioxide concentrations, recorded by four Thomas autometers, and the necessary meterological parameters for the solving of atmospheric dispersion models. Based on these meterological parameters and observed plume rise data, over 4000 one half hourly average maximum and minimum expected ground line sulfur dioxide concentrations were predicted for each of the above dispersion models by the use of computer techniques. The plant is a line source; however, an empirical correction was applied to emission data to reduce them to emissions for an equivalent point source. The predicted sulfur dioxide levels for each of the dispersion models were compared to the measured levels throughout the year. Three different sets of diffusion coefficients were applied to the Sutton model and successful predictions, according to a criterion utilizing an acceptable range of concentration, varied from 66 to 93%. The Bosanquet and Pearson model produced successful predictions 90% of the time, while the USPHS-TVA model was successful 94% of the time.Unsuccessful predictions were primarily overestimates.     

One-dimensional simulation of solute transfer in saturated–unsaturated porous media using the discontinuous finite elements method

A one-dimensional transport model for simulating water flow and solute transport in homogeneous–heterogeneous, saturated–unsaturated porous media is presented. The model is composed of a combination of accurate numerical algorithms for solving the nonlinear Richard's and advection–dispersion equations (ADE). The mixed form of Richard's equation is solved using a standard finite element method (FEM) with primary variable switching. The transport equation is solved using operator splitting, with the discontinuous finite element method (DFE) for discretization of the advective term. A slope limiting procedure for DFE avoids numerical instabilities but creates very limited numerical dispersion for high Peclet numbers. An implicit finite differences scheme (FD) is used for the dispersive term.The unsaturated flow and transport model (Wamos-T) is applied to a variety of rigorous problems including transient flow, heterogeneous medium and abrupt variations of velocity in magnitude and direction due to time-varying boundary conditions. It produces accurate and mass-conservative solutions for a very large range of grid Peclet numbers. The Wamos-T model is a good and robust alternative for the simulation of mass transport in unsaturated domain.     

Application of continuous time random walk theory to nonequilibrium transport in soil

Continuous time random walk (CTRW) formulations have been demonstrated to provide a general and effective approach that quantifies the behavior of solute transport in heterogeneous media in field, laboratory, and numerical experiments. In this paper we first apply the CTRW approach to describe the sorbing solute transport in soils under chemical (or) and physical nonequilibrium conditions by curve-fitting. Results show that the theoretical solutions are in a good agreement with the experimental measurements. In case that CTRW parameters cannot be determined directly or easily, an alternative method is then proposed for estimating such parameters independently of the breakthrough curve data to be simulated. We conduct numerical experiments with artificial data sets generated by the HYDRUS-1D model for a wide range of pore water velocities (υ) and retardation factors (R) to investigate the relationship between CTRW parameters for a sorbing solute and these two quantities (υ, R) that can be directly measured in independent experiments. A series of best-fitting regression equations are then developed from the artificial data sets, which can be easily used as an estimation or prediction model to assess the transport of sorbing solutes under steady flow conditions through soil. Several literature data sets of pesticides are used to validate these relationships. The results show reasonable performance in most cases, thus indicating that our method could provide an alternative way to effectively predict sorbing solute transport in soils. While the regression relationships presented are obtained under certain flow and sorption conditions, the methodology of our study is general and may be extended to predict solute transport in soils under different flow and sorption conditions.     

One-dimensional simulation of solute transfer in saturated-unsaturated porous media using the discontinuous finite elements method

A one-dimensional transport model for simulating water flow and solute transport in homogeneous-heterogeneous, saturated-unsaturated porous media is presented. The model is composed of a combination of accurate numerical algorithms for solving the nonlinear Richard's and advection-dispersion equations (ADE). The mixed form of Richard's equation is solved using a standard finite element method (FEM) with primary variable switching. The transport equation is solved using operator splitting, with the discontinuous finite element method (DFE) for discretization of the advective term. A slope limiting procedure for DFE avoids numerical instabilities but creates very limited numerical dispersion for high Peclet numbers. An implicit finite differences scheme (FD) is used for the dispersive term. The unsaturated flow and transport model (Wamos-T) is applied to a variety of rigorous problems including transient flow, heterogeneous medium and abrupt variations of velocity in magnitude and direction due to time-varying boundary conditions. It produces accurate and mass-conservative solutions for a very large range of grid Peclet numbers. The Wamos-T model is a good and robust alternative for the simulation of mass transport in unsaturated domain.     

Evidence of one-dimensional scale-dependent fractional advection-dispersion

A semi-analytical inverse method and the corresponding program FADEMain for parameter estimation of the fractional advection-dispersion equation (FADE) were developed in this paper. We have analyzed Huang et al.'s [Huang, K., Toride, N., van Genuchten, M.Th., 1995. Experimental investigation of solute transport in large homogeneous and heterogeneous saturated soil columns. Trans. Porous Media 18, 283-302.] laboratory experimental data of conservative solute transport in 12.5-m long homogeneous and heterogeneous soil columns to test the non-Fickian dispersion theory of FADE. The dispersion coefficient was calculated by fitting the analytical solution of FADE to the measured data at different transport scales. We found that the dispersion coefficient increased exponentially with transport scale for the homogeneous column, whereas it increased with transport scale in a power law function for the heterogeneous column. The scale effect of the dispersion coefficient in the heterogeneous soil was much more significant comparing to that in the homogeneous soil. The increasing rate of dispersion coefficient versus transport distance was smaller for FADE than that for the advection-dispersion equation (ADE). Finite difference numerical approximations of the scale-dependent FADE were established to interpret the experimental results. The numerical solutions were found to be adequate for predicting scale-dependent transport in the homogeneous column, while the prediction for the heterogeneous column was less satisfactory.     

Development of the hybrid cells in series model to simulate ammonia nutrient pollutant transport along the Umgeni River

Discharge of organic waste results in high nutrient pollution of the water bodies which is a major menace to the environment. A high quantity of nutrients such as ammonia causes a reduction in the dissolved oxygen level and induces algal growth in the water bodies. Water quality models have been the tools to evaluate the rate at which streams can disperse the pollutants they receive. Many water quality models are flawed either because of their inadequacy to completely simulate the advection component of the pollutant transport, or because of the limited application of the models, due to inaccurate estimation of model parameters. The hybrid cell in series (HCIS) developed by Ghosh et al. (2004) has been able to overcome such difficulties associated with the mixing cell-based models. Thus, the current study focuses on developing an analytical solution for the pollutant transport of the ammonia concentration through the plug flow, the first and second well-mixed cells of the HCIS model. The HCIS model coupled with the first order kinetic equation for ammonia nutrient was developed to simulate the ammonia pollutant concentration in the water column. The ammonia concentration at various points along the river system was assessed by considering the effects of the transformation of ammonia to nitrite, the uptake of ammonia by the algae, the respiration rate of the algae and the input of benthic source to the ammonia concentration in the water column. The proposed model was tested using synthetic data, and the HCIS-NH₃ model simulations for spatial and temporal variation of ammonia pollutant transport were analysed. The simulated results of the HCIS-NH₃ model agreed with the Fickian-based advection-dispersion equation (ADE) for simulating ammonia concentration solved using an explicit finite difference scheme. The HCIS-NH₃ model also showed a good agreement with the observed data from the Umgeni River, except during rainy periods.     

Modeling of dispersion near roadways based on the vehicle-induced turbulence concept

A mathematical model is developed for dispersion near roadways by incorporating vehicle-induced turbulence (VIT) into Gaussian dispersion modeling using computational fluid dynamics (CFD). The model is based on the Gaussian plume equation in which roadway is regarded as a series of point sources. The Gaussian dispersion parameters are modified by simulation of the roadway using CFD in order to evaluate turbulent kinetic energy (TKE) as a measure of VIT. The model was evaluated against experimental carbon monoxide concentrations downwind of two major freeways reported in the literature. Good agreements were achieved between model results and the literature data. A significant difference was observed between the model results with and without considering VIT. The difference is rather high for data very close to the freeways. This model, after evaluation with additional data, may be used as a framework for predicting dispersion and deposition from any roadway for different traffic (vehicle type and speed) conditions.     

An efficient Lagrangian stochastic model of vertical dispersion in the convective boundary layer

We consider the one-dimensional case of vertical dispersion in the convective boundary layer (CBL) assuming that the turbulence field is stationary and horizontally homogeneous. The dispersion process is simulated by following Lagrangian trajectories of many independent tracer particles in the turbulent flow field, leading to a prediction of the mean concentration. The particle acceleration is determined using a stochastic differential equation, assuming that the joint evolution of the particle velocity and position is a Markov process. The equation consists of a deterministic term and a random term. While the formulation is standard, attention has been focused in recent years on various ways of calculating the deterministic term using the well-mixed condition incorporating the Fokker–Planck equation. Here we propose a simple parameterisation for the deterministic acceleration term by approximating it as a quadratic function of velocity. Such a function is shown to represent well the acceleration under moderate velocity skewness conditions observed in the CBL. The coefficients in the quadratic form are determined in terms of given turbulence statistics by directly integrating the Fokker–Planck equation. An advantage of this approach is that, unlike in existing Lagrangian stochastic models for the CBL, the use of the turbulence statistics up to the fourth order can be made without assuming any predefined form for the probability distribution function (PDF) of the velocity. The main strength of the model, however, lies in its simplicity and computational efficiency. The dispersion results obtained from the new model are compared with existing laboratory data as well as with those obtained from a more complex Lagrangian model in which the deterministic acceleration term is based on a bi-Gaussian velocity PDF. The comparison shows that the new model performs well.     

Fuel reformulation in the Metropolitan Area of Mexico City

Two models frequently used to simulate the dispersion of pollutants in the atmosphere have been compared. This is necessary because only a well-tested and well-calibrated simulation model can be a good representation of the reality of the dispersion of pollutants. The models evaluated (HYSPLIT_4 with its four variants and MEDIA) were run using as input parameters the same meteorological dataset (for 23-26 October 1994) from the French model ARPEGE. The following statistical criteria were compared: the space and time evolution of the pollutant cloud; the variation of statistical parameters in time and space; and the differences between the simulated and measured values of concentration in time for six different stations. The results emphasise the characteristics of the two models and their abilities in the framework of the air quality monitoring.     

Wind-tunnel study of concentration fields in street canyons

The paper presents results from a case study of gaseous pollutant dispersion in street canyons. Tracer-gas experiments were performed in a neutrally stratified wind tunnel. Vehicle emissions were simulated as line sources. Concentration profiles along building walls were measured. A two-dimensional street canyon was considered as the reference case. The influence of systematic parameter variations on the concentration field is studied and discussed. Building dimensions, upwind building configuration, wind direction and roof geometry were found to be important parameters. Data sets from the study may be used for evaluation of numerical models and for expert estimates of air quality in the urban environment     

Temperature modeling in activated sludge systems: a case study.

A model of temperature dynamics was developed as part of a general model of activated-sludge reactors. Transport of heat was described by the one-dimensional, advection-dispersion equation, with a source term based on a theoretical heat balance over the reactor. The model was compared to several reference models, including a tanks-in-series model and the dispersion model with heat components neglecting biochemical-energy inputs and other activated-sludge, heat-balance terms. All the models were tested under steady-state and dynamic conditions at a full-scale facility, the Rock Creek wastewater treatment plant in Hillsboro, Oregon, using meteorological data from a station located 16 km from the plant. The dispersion model and tanks-in-series model matched in situ temperature data with absolute-mean errors less than 0.1 degrees C. Neglecting biochemical-heat-energy inputs in the activated-sludge reactor underestimated temperatures by up to 0.5 degrees C. The biochemical-heat-energy inputs accounted for 30 to 40% of the total heat flux throughout the year.     

An experimental study of flow and transport in fractured slate

A series of laboratory tracer migration experiments in a single rock fracture have been performed, and the breakthrough curves have been interpreted using mathematical modelling. Discrepancies were observed between the experimental data and the predictions made using a simple advection-dispersion model. The potential reasons for these discrepancies have been investigated by applying more complex models: one model incorporates channelling of flow within the fracture, the other couples dispersion and advection in the fracture with rock-matrix diffusion. It is concluded that chanelling of flow can adequately explain the observed spreading behaviour; rock-matrix diffusion is not a significant mechanism influencing transport in these experiments.     

Migration of a water pulse through fractured porous media

Contaminant transport from waste-disposal sites is strongly affected by the presence of fractures and the degree of fracture matrix interaction. Characterization of potential contaminant plumes at such sites is difficult, both experimentally and numerically. Simulations of water flow through fractured rock were performed to examine the penetration depth of a large pulse of water entering such a system. Construction water traced with lithium bromide was released during the excavation of a tunnel at Yucca Mountain, Nevada, which is located in an unsaturated fractured tuff formation. Modeling of construction-water migration is qualitatively compared with bromide-to-chloride ratio (Br/Cl) data for pore-water salts extracted from drillcores. The influences of local heterogeneities in the fracture network and variations in hydrogeologic parameters were examined by sensitivity analyses and Monte Carlo simulations. The simulation results are qualitatively consistent with the observed Br/Cl signals, although these data may only indicate a minimum penetration depth, and water may have migrated farther through the fracture network.     

Biodegradation modelling of a dissolved gasoline plume applying independent laboratory and field parameters

Biodegradation of organic contaminants in groundwater is a microscale process which is often observed on scales of 100s of metres or larger. Unfortunately, there are no known equivalent parameters for characterizing the biodegradation process at the macroscale as there are, for example, in the case of hydrodynamic dispersion. Zero- and first-order degradation rates estimated at the laboratory scale by model fitting generally overpredict the rate of biodegradation when applied to the field scale because limited electron acceptor availability and microbial growth are not considered. On the other hand, field-estimated zero- and first-order rates are often not suitable for predicting plume development because they may oversimplify or neglect several key field scale processes, phenomena and characteristics. This study uses the numerical model BIO3D to link the laboratory and field scales by applying laboratory-derived Monod kinetic degradation parameters to simulate a dissolved gasoline field experiment at the Canadian Forces Base (CFB) Borden. All input parameters were derived from independent laboratory and field measurements or taken from the literature a priori to the simulations. The simulated results match the experimental results reasonably well without model calibration. A sensitivity analysis on the most uncertain input parameters showed only a minor influence on the simulation results. Furthermore, it is shown that the flow field, the amount of electron acceptor (oxygen) available, and the Monod kinetic parameters have a significant influence on the simulated results. It is concluded that laboratory-derived Monod kinetic parameters can adequately describe field scale degradation, provided all controlling factors are incorporated in the field scale model. These factors include advective–dispersive transport of multiple contaminants and electron acceptors and large-scale spatial heterogeneities.     

A two-region nonequilibrium model for solute transport in solution conduits in karstic aquifers

A two-region nonequilibrium model was used to calibrate initial solute-transport parameter estimates generated from tracer-breakthrough curves (TBCs) developed from tracer tests conducted in uni-axial solution conduits in karstic aquifers. Two-region nonequilibrium models account for partitioning of solute into mobile- and immobile-fluid regions to produce a more representative model fit to the strong tails associated with TBCs than do equilibrium models. The nonequilibrium model resulted in an increase in average flow velocities and a decrease in longitudinal dispersion coefficients over comparable estimates using an equilibrium model. Increases in velocity and decreases in dispersion were obtained at the expense of including parameters that describe solute partitioning and mass transfer rate for the mobile- and immobile-fluid regions. In addition, nonidentifiable sorption and mass transfer parameters for the immobile-fluid regions could only be described in terms of upper and lower bounds using readily determined identifiable ratios representing solute partitioning and system constraints based on known physical properties. The identifiable ratios and system constraints serve to minimize model nonuniqueness and renders the nonidentification problem trivial.