Two-domain description of solute transport in heterogeneous porous media: Comparison between theoretical predictions and numerical experiments (original) (raw)

Solute transport in dual-permeability porous media

Water Resources Research, 2012

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.

Comparison of theory and experiment for solute transport in highly heterogeneous porous medium

Advances in Water Resources, 2007

In this work, the influence of non-equilibrium effects on solute transport in a weakly heterogeneous medium is discussed. Three macro-scale models (upscaled via the volume averaging technique) are investigated: (i) the two-equation non-equilibrium model, (ii) the one-equation asymptotic model and (iii) the one-equation local equilibrium model. The relevance of each of these models to the experimental system conditions (duration of the pulse injection, dispersivity values. . .) is analyzed. The numerical results predicted by these macroscale models are compared directly with the experimental data (breakthrough curves). Our results suggest that the preasymptotic zone (for which a non-Fickian model is required) increases as the solute input pulse time decreases. Beyond this limit, the asymptotic regime is recovered. A comparison with the results issued from the stochastic theory for this regime is performed. Results predicted by both approaches (volume averaging method and stochastic analysis) are found to be consistent.