Pore-scale modeling of dispersion in disordered porous media

We employ a direct pore-level model of incompressible flow that uses the modified moving particle semi-implicit (MMPS) method. The model is capable of simulating both unsteady- and steady-state flow directly in microtomography images of naturally-occurring porous media. We further develop this model to simulate solute transport in disordered porous media. The governing equations of flow and transport at the pore level, i.e., Navier-Stokes and convection-diffusion, are solved directly in the pore space mapped by microtomography techniques. Three naturally-occurring sandstones are studied in this work. We verify the accuracy of the model by comparing the computed longitudinal dispersion coefficients against the experimental data for a wide range of Peclet numbers, i.e., 5×10(-2)相似文献   

Choice of dispersion coefficients in reactive transport calculations on smoothed fields

Local dispersion dominates the mixing of compounds that are introduced separately into the subsurface and do not partition into any other than the aqueous phase. Thus, reactions between these compounds are controlled by dispersive mixing if they are not limited by kinetics. I quantify longitudinal dispersive mixing by the longitudinal effective dispersion coefficient of a conservative tracer introduced by a point-like injection [Water Resour. Res. 36 (12) (2000) 3591-3604]. In the upscaling of mixing-controlled reactive transport, I apply the mean velocity and the effective dispersion coefficient to the macroscopic transport calculations, whereas the reactive parameters on the macro-scale are identical to those on the local scale. The applicability of the approach is demonstrated for the transport of compounds undergoing a second-order irreversible bimolecular reaction. Ten realizations of a two-dimensional heterogeneous log-conductivity field are considered. Using the effective dispersion parameters, the overall mass balance is met in the ensemble average, whereas solute spreading is underestimated. I assess the lack of spreading by the difference between the expected macrodispersion and effective dispersion coefficients. I extend the approach to simulations on log-conductivity fields obtained by kriging of regularly spaced conductivity measurements. These fields contain the large-scale features of the true fields but do not resolve the small-scale variability. For the calculations on the kriged fields, the corresponding conditional covariance is substituted into the analytical expressions of effective dispersion, yielding a correction effective dispersion coefficient. The comparison between simulations on the fully resolved fields and on the kriged fields indicates that the approach is valid for wide plumes meeting the ergodicity condition. The high variability of mixing on small scales unresolved by kriging, however, leads to severe uncertainty when mixing-controlled reactions are predicted for narrow plumes.     

Transverse vertical dispersion in groundwater and the capillary fringe

Transverse dispersion is the most relevant process in mass transfer of contaminants across the capillary fringe (both directions), dilution of contaminants, and mixing of electron acceptors and electron donors in biodegrading groundwater plumes. This paper gives an overview on literature values of transverse vertical dispersivities alpha(tv) measured at different flow velocities and compares them to results from well-controlled laboratory-tank experiments on mass transfer of trichloroethene (TCE) across the capillary fringe. The measured values of transverse vertical dispersion in the capillary fringe region were larger than in fully saturated media, which is credited to enhanced tortuosity of the flow paths due to entrapped air within the capillary fringe. In all cases, the values observed for alpha(tv) were < 1 mm. The new measurements and the literature values indicate that alpha(tv) apparently declines with increasing flow velocity. The latter is attributed to incomplete diffusive mixing at the pore scale (pore throats). A simple conceptual model, based on the mean square displacement and the pore size accounting for only partial diffusive mixing at increasing flow velocities, shows very good agreement with measured and published data.     

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.     

A wind tunnel study of dense gas dispersion in a neutral boundary layer over a rough surface

Measurements of the vertical entrainment velocity into two-dimensional dense gas plumes over fully rough surfaces were carried out as part of a co-operative research programme with wind tunnel facilities in the USA. This paper presents results obtained for neutral boundary layer conditions in the EnFlo wind tunnel at the University of Surrey; a companion paper treats the stable boundary layer case. Entrainment velocities, W_E, were deduced from the streamwise development of the concentration field, non-dimensionalised with respect to the friction velocity in the undisturbed flow, u^*, and correlated with the plume Richardson number, Ri^*. Results for Richardson numbers in the range Ri^*<15 were found to be well fitted by the empirical expression: W_E/u^*=0.65/(1+0.2Ri^*). Flow visualisation studies showed layered plume structures with a sharp upper interface at higher Richardson numbers and in this regime turbulent motion below the interface became progressively more intermittent as Ri^* increased. Measured turbulence levels collapsed within such high Richardson number plumes and flow and dispersion were significantly affected by molecular processes. Up-welling above the source was observed when the emission speed exceeded the approach flow friction velocity, though there was no clear evidence that this affected plume behaviour away from the immediate vicinity of the source.     

Predictive modelling of dispersion controlled reactive plumes at the laboratory-scale

A model-based interpretation of laboratory-scale experimental data is presented. Hydrolysis experiments carried out using thin glass tanks filled with glass beads to construct a hypothetical and inert, homogeneous porous medium were analysed using a 2D numerical model. A new empirical formula, based upon results for non-reactive (tracer) experiments is used to calculate transversal dispersivity values for a range of grain sizes and any flow velocities. Combined with effective diffusion coefficients calculated from Stokes-Einstein type equations, plume lengths arising from mixing between two solutes can be predicted accurately using numerical modelling techniques. Moreover, pH and ion concentration profiles lateral to the direction of flow of the mixing species can be determined at any given point downstream, without the need for result fitting. In our case, this approach does not lead to overpredictions of lateral mixing, as previously reported when using parameters derived from non-reactive tracer experiments to describe reactive solute transport. The theory is based on the assumption of medium homogeneity.     

A new Lagrangian particle model for the simulation of dense gas dispersion

A Lagrangian stochastic model (MicroSpray), able to simulate the airborne dispersion in complex terrain and in presence of obstacles, was modified to simulate the dispersion of dense gas clouds. This is accomplished by taking into account the following processes: negative buoyancy, gravity spreading and the particle's reflection at the bottom computational boundary. Elevated and ground level sources, continuous and instantaneous emissions, time varying sources, plumes with initial momentum (horizontal, vertical or oblique in any direction), plumes without initial momentum are considered. MicroSpray is part of the model system MSS, which also includes the diagnostic MicroSwift model for the reconstruction of the 3-D wind field in presence of obstacles and orography. To evaluate the MSS ability to simulate the dispersion of heavy gases, its simulation performances are compared in detail to two field experiments (Thorney Island and Kit Fox) and to a chlorine railway accident (Macdona). Then, a comprehensive analysis considering several experiments of the Modelers Data Archive is presented. The statistical analysis on the overall available data reveals that the performance of the new MicroSpray version for dense gas releases is generally reliable. For instance, the agreement between concentration predictions and observations is within a factor of two in the 72% up to 99% of the occurrences for the case studies considered. The values of other performance measures, such as correlation coefficient, geometric mean bias and geometric variance, mostly set in the ranges indicated as good-model performances in the specialized literature.     

Estimation of local scale dispersion from local breakthrough curves during a tracer test in a heterogeneous aquifer: the Lagrangian approach

The local scale dispersion tensor, Dd, is a controlling parameter for the dilution of concentrations in a solute plume that is displaced by groundwater flow in a heterogeneous aquifer. In this paper, we estimate the local scale dispersion from time series or breakthrough curves, BTCs, of Br concentrations that were measured at several points in a fluvial aquifer during a natural gradient tracer test at Krauthausen. Locally measured BTCs were characterized by equivalent convection dispersion parameters: equivalent velocity, v(eq)(x) and expected equivalent dispersivity, [lambda(eq)(x)]. A Lagrangian framework was used to approximately predict these equivalent parameters in terms of the spatial covariance of log(e) transformed conductivity and the local scale dispersion coefficient. The approximate Lagrangian theory illustrates that [lambda(eq)(x)] increases with increasing travel distance and is much larger than the local scale dispersivity, lambda(d). A sensitivity analysis indicates that [lambda(eq)(x)] is predominantly determined by the transverse component of the local scale dispersion and by the correlation scale of the hydraulic conductivity in the transverse to flow direction whereas it is relatively insensitive to the longitudinal component of the local scale dispersion. By comparing predicted [lambda(eq)(x)] for a range of Dd values with [lambda(eq)(x)] obtained from locally measured BTCs, the transverse component of Dd, DdT, was estimated. The estimated transverse local scale dispersivity, lambda(dT) = DdT/U1 (U1 = mean advection velocity) is in the order of 10(1)-10(2) mm, which is relatively large but realistic for the fluvial gravel sediments at Krauthausen.     

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.     

Effects of sorption on transverse mixing in transient flows

Transverse mixing has been identified as a potentially limiting factor for natural attenuation of plumes originating from continuously emitting sources. Under steady-state flow conditions, dispersion is the only process leading to lateral mixing. This process is very slow and cannot explain the lateral spread of plumes observed in the field. When the flow direction fluctuates with time, transverse dispersion is slightly enhanced, but not very pronounced. Under these flow conditions, however, sorption can contribute to mixing into the mean transverse direction. If the reacting compounds differ in their strength of sorption, chromatographic mixing and separation alternate in time-periodic flows. For instantaneous sorption, the plumes may overlap within a stripe of fixed width. In contrast to sorption in local equilibrium, kinetic sorption contributes to mixing also for compounds with identical sorption strength. I derive an analytical expression for the equivalent transverse dispersion coefficient of a kinetically sorbing compound in a spatially uniform flow field undergoing sinusoidal fluctuations in time. This expression may be used for reactive transport calculations in an equivalent domain with constant flow. The effects are the strongest for compounds with a dimensionless partitioning coefficient of about unity, slow sorption kinetics, and slowly fluctuating velocities. For realistic parameters, kinetic sorption contributes to transverse mixing in the same range as heterogeneity.     

Model-based evaluation of controlled-release systems in the remediation of dissolved plumes in groundwater

Controlled-release, semi-passive reactive barrier systems have been recently developed as a long-term treatment option for controlling the spread of contaminant plumes in groundwater. This paper describes a new computer code, and applies it to study coupled processes of solute release, reaction, and mass transport in an in situ remediation scheme using the controlled release of potassium permanganate. Confidence with the modeling approach was developed by model verifications and simulating results of a pilot-scale test-cell experiment. Sensitivity analyses indicated the possibilities of treatment inefficiencies due to inability of transverse dispersion to mix the permanganate (MnO(4)(-)) within the zone of reaction, fluctuations in source strength due to variations in flow velocity, and the small length of treatment zone due to strong soil utilization of MnO(4)(-). Although problems associated with the fluctuating source strength and strong soil utilization can be addressed by optimizing the release rate, the inefficiency of transverse dispersion to create mixing could pose a serious limitation. Through a series of model simulations, a system of injection/withdrawal wells in a doublet arrangement was developed to facilitate lateral spreading and mixing of MnO(4)(-). A well-mixed, stable MnO(4)(-) zone with predetermined size (DxL=8m x 2m) and concentration ranges (1.5-20 mg l(-1)) was created by four 1-day injection/withdrawal pumping periods over 24 d. This type of mixing zone may persist for many years with periodic well mixing and replacements of exhausted controlled-release forms. Coupled use of the generalized code with field hydrologic data will help to optimize the design and operation of controlled-release systems in practice.     

Modelling plume rise and dispersion from pool fires

This paper describes an investigation into the behaviour of smoke plumes from pool fires, and the subsequent generation of empirical models to predict plume rise and dispersion from such a combustion source. Synchronous video records of plumes were taken from a series of small-scale (0.06–0.25m²) outdoor methanol/toluene pool fire experiments, and used to produce sets of images from which plume dimensions could be derived. Three models were used as a basis for the multiple regression analysis of the data set, in order to produce new equations for improved prediction. Actual plume observations from a large (20.7 m×14.2 m) aviation fuel pool fire were also used to test the predictions. The two theoretically based models were found to give a better representation of plume rise and dispersion than the empirical model based on measurements of small-scale fires. It is concluded that theoretical models tested on small-scale fires (heat output ≈70 kW) can be used to predict plume behaviour from much larger combustion sources (heat output ≈70 MW) under near neutral atmospheric conditions.     

Near-field dispersion modeling for regulatory applications

This paper evaluates the application of dispersion models to estimate near-field pollutant concentrations in two case studies. The Industrial Source Complex Short-Term Model (ISCST3) was evaluated with hexavalent chromium measurements collected within 100 m of two facilities in Barrio Logan, San Diego, CA. ISCST3 provided reasonable estimates for higher pollutant concentrations but underestimated lower concentrations. To understand the observed distribution of concentrations in Barrio Logan, a recently conducted tracer experiment was analyzed. The tracer, sulfur hexafluoride, was released at ambient temperature from an urban facility at the University of California at Riverside, and concentrations were measured within 20 m of the source. Modeling results indicated that Industrial Source Complex-Plume Rise Model Enhancement and American Meteorological Society/U.S. Environmental Protection Agency Regulatory Model-Plume Rise Model Enhancement overestimated high concentrations and underestimated low concentrations. A diagnostic study with a simple Gaussian dispersion model that incorporated site-specific meteorology was used to evaluate model results. This study found that incorporating lateral meandering for nonbuoyant urban plumes in Gaussian dispersion models could improve concentration estimates even when downwash is not considered. Incorporating a meandering component in ISCST3 resulted in improvements in estimating hexavalent chromium concentrations in Barrio Logan. Credible near-source concentration estimates depend on accurate characterization of emissions, onsite micrometeorology, and a method to account for lateral meandering in the near field.     

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.     

Experiments on vertical transverse mixing in a large-scale heterogeneous model aquifer

Vertical transverse mixing is known to be a controlling factor in natural attenuation of extended biodegradable plumes originating from continuously emitting sources. We perform conservative and reactive tracer tests in a quasi two-dimensional 14 m long sand box in order to quantify vertical mixing in heterogeneous media. The filling mimics natural sediments including a distribution of different hydro-facies, made of different sand mixtures, and micro-structures within the sand lenses. We quantify the concentration distribution of the conservative tracer by the analysis of digital images taken at steady state during the tracer-dye experiment. Heterogeneity causes plume meandering, leading to distorted concentration profiles. Without knowledge about the velocity distribution, it is not possible to determine meaningful vertical dispersion coefficients from the concentration profiles. Using the stream-line pattern resulting from an inverse model of previous experiments in the sand box, we can correct for the plume meandering. The resulting vertical dispersion coefficient is approximately approximately 4 x 10(-)(9) m(2)/s. We observe no distinct increase in the vertical dispersion coefficient with increasing travel distance, indicating that heterogeneity has hardly any impact on vertical transverse mixing. In the reactive tracer test, we continuously inject an alkaline solution over a certain height into the domain that is occupied otherwise by an acidic solution. The outline of the alkaline plume is visualized by adding a pH indicator into both solutions. From the height and length of the reactive plume, we estimate a transverse dispersion coefficient of approximately 3 x 10(-)(9) m(2)/s. Overall, the vertical transverse dispersion coefficients are less than an order of magnitude larger than pore diffusion coefficients and hardly increase due to heterogeneity. Thus, we conclude for the assessment of natural attenuation that reactive plumes might become very large if they are controlled by vertical dispersive mixing.     

Experimental Zn(II) retention in a sandy loam soil by very small columns

The use of column experiments, usually performed to better approximate field conditions, may provide information that is not available from batch experiments. In such experiments heavy metals are often adsorbed until saturation followed by desorption experiments. When the affinity of the metal to soil is high, the retention factor (R) could be greater than thousands and the duration of experiments can become impractically long. In order to use reasonable laboratory time, the flow rate should be increased or the column size decreased. The increase in flow rate produces undesirable kinetic and dispersion effects, so we used very small soil columns (pore volume = 0.31–0.70 ml) and relatively high flow rates (0.03–0.12 ml min⁻¹) in studies of Zn(II) adsorption and retention in soils. Conservative tracer flow column experiments under saturation conditions were carried out to determine flow parameters for different flow rates. Column pore volume (V_p), Peclet numbers (Pe) and longitudinal dispersion coefficients (D_L) were determined from breakthrough curves. The effect of type of electrolyte and ionic strength on the Zn(II) retention onto soil was determined. The influence of flow rate and bed height on the retention coefficient and on the mass transfer zone was also studied. The effect of different influent Zn(II) concentrations on the R values obtained was analyzed. Freundlich parameters from column experiments were compared with batch ones. The leaching efficiency of different electrolytes, salts of weak organic acids and EDTA was also studied.     

Atmospheric dispersion under non-homogeneous and unsteady conditions

Accidental releases of toxic material can give rise to serious air pollution situations. The behaviour of such releases is frequently quite different from that of normal chimney plumes and is highly unpredictable since their extent and occurrence cannot be known in advance. Thus one has to resort to computational simulation procedures for reasonable estimates of pollutant dispersion parameters in the atmosphere. In the present study efforts are made to develop such simulation procedures. Computer software is developed which can be obtained free of cost from the authors to simulate pollutant dispersal from continuous/short-term, steady/fluctuating, surface/elevated sources in stationary/non-stationary and homogeneous/non-homogeneous atmosphere using basic Gaussian theory. Validity and limitations of developed software are demonstrated through comparison of results obtained using different programs and field measurements. In general, the predicted results appear to be quite logical and the computer software prepared may prove to be a valuable aid in estimating the extent of air pollution.     

A high-resolution non-invasive approach to quantify oxygen transport across the capillary fringe and within the underlying groundwater

Oxygen transport across the capillary fringe is relevant for many biogeochemical processes. We present a non-invasive technique, based on optode technology, to measure high-resolution concentration profiles of oxygen across the unsaturated/saturated interface. By conducting a series of quasi two-dimensional flow-through laboratory experiments, we show that vertical hydrodynamic dispersion in the water-saturated part of the capillary fringe is the process limiting the mass transfer of oxygen. A number of experimental conditions were tested in order to investigate the influence of grain size and horizontal flow velocity on transverse vertical dispersion in the capillary fringe. In the same setup, analogous experiments were simultaneously carried out in the fully water-saturated zone, therefore allowing a direct comparison with oxygen transfer across the capillary fringe. The outcomes of the experiments under various conditions show that oxygen transport in the two zones of interest (i.e., the unsaturated/saturated interface and the saturated zone) is characterized by very similar transverse dispersion coefficients. An influence of the capillary fringe morphology on oxygen transport has not been observed. These results may be explained by the narrow grain size distribution used in the experiments, leading to a steep decline in water saturation at the unsaturated/saturated interface and to the absence of trapped gas in this transition zone. We also modeled flow (applying the van Genuchten and the Brooks-Corey relationships) and two-dimensional transport across the capillary fringe, obtaining simulated profiles of equivalent aqueous oxygen concentration that were in good agreement with the observations.     

A controlled field experiment on groundwater contamination by a multicomponent DNAPL: creation of the emplaced-source and overview of dissolved plume development

A unique field experiment has been undertaken at the CFB Borden research site to investigate the development of dissolved chlorinated solvent plumes from a residual dense non-aqueous phase liquid (DNAPL) source. The "emplaced-source" tracer test methodology involved a controlled emplacement of a block-shaped source of sand containing chlorinated solvents below the water table. The gradual dissolution of this residual DNAPL solvent source under natural aquifer conditions caused dissolved solvent plumes of trichloromethane (TCM), trichloroethene (TCE) and perchloroethene (PCE) to continuously develop down gradient. Source dissolution and 3-D plume development were successfully monitored via 173 multilevel samplers over a 475-day tracer test period prior to site remediation research being initiated. Detailed groundwater level and hydraulic conductivity data were collected. Development of plumes with concentrations spanning 1-700,000 micrograms/1 is described and key processes controlling their migration identified. Plumes were observed to be narrow due to the weakness of transverse dispersion processes and long due to advection and significant longitudinal dispersion, very limited sorptive retardation and negligible, if any, attenuation due to biodegradation or abiotic reaction. TCM was shown to be essentially conservative, TCE very nearly conservative and PCE, consistent with its greater hydrophobicity, more retarded yet having a greater mobility than observed in previous Borden field tests. The absence of biodegradation was ascribed to the prevailing aerobic conditions and lack of any additional biodegradable carbon substrates. The transient groundwater flow regime caused significant transverse lateral plume movement, plume asymmetry and was likely responsible for most of the, albeit limited, transverse horizontal plume spreading. In agreement with the widespread incidence of extensive TCE and PCE plumes throughout the industrialized world, the experiment indicates such solvent plumes are likely to be highly mobile and persistent, at least in aquifers that are aerobic and have low sorption potential (low foc content).     

Overview of Petroleum Environmental Research Forum (PERF) dense gas dispersion modeling project

This paper provides a background for and an overview of the results of a comprehensive study of transport and dispersion of dense gas plumes over rough surfaces typical of industrial sites. The Petroleum Environmental Research Forum (PERF) 93-16 project involved model development and evaluations using observations from three wind tunnels and from the Kit Fox field experiment. Detailed discussions of the results of the research are given in the other papers in this special issue. The wind tunnel experiments produced data showing that the resulting best-fit vertical entrainment formula was close to (i.e., within about 30%) the vertical entrainment formulas already in use by current models, which were derived primarily from observations over smooth surfaces. Observations from the Kit Fox field experiment demonstrated the validity of the entrainment curves derived from the wind tunnel data. The Kit Fox data were also used to evaluate algorithms for along-wind dispersion and cloud advection speeds for short-duration releases typical of an industrial site, and to evaluate the HEGADAS dense gas dispersion model.