Correspondence Between One- and Two-Equation Models for Solute Transport in Two-Region Heterogeneous Porous Media

In this work, we study the transient behavior of homogenized models for solute transport in two-region porous media. We focus on the following three models: (1) a time non-local, two-equation model (2eq-nlt). This model does not rely on time constraints and, therefore, is particularly useful in the short-time regime, when the timescale of interest (t) is smaller than the characteristic time (τ1) for the relaxation of the effective macroscale parameters (i.e., when t ≤ τ1); (2) a time local, two-equation model (2eq). This model can be adopted when (t) is significantly larger than (τ1) (i.e., when \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${t\gg\tau_{1}}$$\end{document}); and (3) a one-equation, time-asymptotic formulation (1eq∞). This model can be adopted when (t) is significantly larger than the timescale (τ2) associated with exchange processes between the two regions (i.e., when \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${t\gg\tau_{2}}$$\end{document}). In order to obtain insight into this transient behavior, we combine a theoretical approach based on the analysis of spatial moments with numerical and analytical results in several simple cases. The main result of this paper is to show that there is only a weak asymptotic convergence of the solution of (2eq) towards the solution of (1eq∞) in terms of standardized moments but, interestingly, not in terms of centered moments. The physical interpretation of this result is that deviations from the Fickian situation persist in the limit of long times but that the spreading of the solute is eventually dominating these higher order effects.