Open AccessLocal Solute Sinks and Sources Cause Erroneous Dispersion Fluxes in Transport Simulations with the Convection–Dispersion Equation
Author(s) -
Peters Andre,
Iden Sascha C.,
Durner Wolfgang
Publication year - 2019
Publication title -
vadose zone journal
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.036
H-Index - 81
ISSN - 1539-1663
DOI - 10.2136/vzj2019.06.0064
Subject(s) - dispersion (optics) , flux (metallurgy) , molecular diffusion , convection , diffusion , porous medium , mechanics , groundwater , evaporation , water transport , chemistry , vadose zone , flow (mathematics) , thermodynamics , water flow , soil science , environmental science , soil water , porosity , geology , physics , geotechnical engineering , optics , metric (unit) , operations management , organic chemistry , economics
Core Ideas The CDE predicts non‐physical dispersive solute flux against water flow direction. This wrongly directed dispersive solute flux can be as high as the convective flux. This problem occurs around zones with an increase or decrease in concentration. Particular cases are vaporization of water, root water uptake, and biodegradation. The convection–dispersion equation (CDE) is the most widely used model for simulating the transport of dissolved substances in porous media. The dispersion term in the CDE lumps molecular diffusion and hydromechanical dispersion into an effective diffusive solute flux. This is possible by describing hydrodynamic dispersion with Fick's first law of diffusion. We critically analyzed this concept for specific water flow situations where the solute concentration is locally increased by processes like root water uptake or water evaporation. The local accumulation of solutes in these situations leads to high concentration gradients and a dispersive solute flux component opposite to the direction of the water flux. This is physically wrong because it assumes that molecules or ions are moving against the flow direction by dispersion. The aim of this study was to investigate the magnitude of the resulting error by means of numerical modeling. We simulated solute transport from a groundwater table to a bare soil surface during steady‐state evaporation using the HYDRUS‐1D code. The simulations showed that in the region where dissolved substances accumulate due to the transition from liquid water to vapor, the resulting incorrect dispersive flux against the mean transport direction can reach the same order of magnitude as the convective solute flux. Under such conditions, application of the CDE is questionable.