Geology ReportsSearch

USGS · 70044274

Applying dispersive changes to Lagrangian particles in groundwater transport models

Abstract

Method-of-characteristics groundwater transport models require that changes in concentrations computed within an Eulerian framework to account for dispersion be transferred to moving particles used to simulate advective transport. A new algorithm was developed to accomplish this transfer between nodal values and advecting particles more precisely and realistically compared to currently used methods. The new method scales the changes and adjustments of particle concentrations relative to limiting bounds of concentration values determined from the population of adjacent nodal values. The method precludes unrealistic undershoot or overshoot for concentrations of individual particles. In the new method, if dispersion causes cell concentrations to decrease during a time step, those particles in the cell having the highest concentration will decrease the most, and those with the lowest concentration will decrease the least. The converse is true if dispersion is causing concentrations to increase. Furthermore, if the initial concentration on a particle is outside the range of the adjacent nodal values, it will automatically be adjusted in the direction of the acceptable range of values. The new method is inherently mass conservative.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Leonard F. Konikow. 2010-04-23. Applying dispersive changes to Lagrangian particles in groundwater transport models. https://doi.org/10.1007/s11242-010-9571-2

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

Determination of water retention in stratified porous materials

Predicted and measured water-retention values, θ ( ψ ), were compared for repacked, stratified core samples consisting of either a sand with a stone-bearing layer or a sand with a clay loam layer in various spatial orientations. Stratified core samples were packed in submersible pressure outflow cells, then water-retention measurements were performed between matric potentials, ψ , of 0 to -100 kPa. Predictions of θ ( ψ ) were based on a simple volume-averaging model using estimates of the relative fraction and θ ( ψ ) values of each textural component within a stratified sample. In general, predicted θ ( ψ ) curves resembled measured curves well, except at higher saturations in a sample consisting of a clay loam layer over a sand layer. In this case, the model averaged the air-entry of both materials, while the air-entry of the sample was controlled by the clay loam in contact with the cell's air-pressure inlet. In situ, avenues for air-entry generally exist around clay layers, so that the model should adequately predict air-entry for stratified formations regardless of spatial orientation of fine versus coarse layers. Agreement between measured and predicted volumetric water contents, θ , was variable though encouraging, with mean differences between measured and predicted θ values in the range of 10%. Differences in θ of this magnitude are expected due to variability in pore structure between samples, and do not indicate inherent problems with the volume averaging model. This suggets that explicit modeling of stratified formations through detailed characterization of the stratigraphy has the potential of yielding accurate θ ( ψ ) values. However, hydraulic-equilibration times were distinctly different for each variation in spatial orientation of textural layering, indicating that transient behavior during drainage in stratified formations is highly sensitive to the stratigraphic sequence of textural components, as well as the volume fraction of each textural component in a formation. This indicates that prolonged residence times of water, nutrients, and pollutants are likely within finer-textured layers, when ψ conditions have resulted in drainage of underlying coarser-textured strata.

Transport in Porous Media

Effects of clay dispersion on aquifer storage and recovery in coastal aquifers

Cyclic injection, storage, and withdrawal of freshwater in brackish aquifers is a form of aquifer storage and recovery (ASR) that can beneficially supplement water supplies in coastal areas. A 1970s field experiment in Norfolk, Virginia, showed that clay dispersion in the unconsolidated sedimentary aquifer occurred because of cation exchange on clay minerals as freshwater displaced brackish formation water. Migration of interstitial clay particles clogged pores, reduced permeability, and decreased recovery efficiency, but a calcium preflush was found to reduce clay dispersion and lead to a higher recovery efficiency. Column experiments were performed in this study to quantify the relations between permeability changes and clay mineralogy, clay content, and initial water salinity. The results of these experiments indicate that dispersion of montmorillonite clay is a primary contributor to formation damage. The reduction in permeability by clay dispersion may be expressed as a linear function of chloride content. Incorporating these simple functions into a radial, cross-sectional, variable-density, ground-water flow and transport model yielded a satisfactory simulation of the Norfolk field test - and represented an improvement over the model that ignored changes in permeability. This type of model offers a useful planning and design tool for ASR operations in coastal clastic aquifer systems.

Transport in Porous Media