Transport processes taking place in natural formations are often characterized by the spatial nonstationarity of flow field and the nonergodic condition of solute spreading. The former may originate from statistical inhomogeneity of porous media properties or from the influence of boundary conditions as well as from conditioning procedures on measured data, while the latter is generally due to the finite size of the solute source. In real-world applications these aspects sometimes dominate the transport phenomena, and known literature results based on statistical homogeneity and ergodicity assumptions cannot be properly applied. In this paper a method is proposed to handle different and concurrent causes of the flow field nonstationarity and to describe nonergodic transport of inert solutes by spatial moments in a domain of finite size. The goal is reached by expanding the steady state flow equation in Taylor series limited to the first-order and by the recursive application of finite element (FE) method. The unknowns are the piezometric head mean values and its derivatives in respect to the fluctuating porous media hydraulic conductivity. From the latter the velocity statistics are derived on the basis of Darcy's law according to the first-order approximation. By neglecting the pore-scale dispersion the spatial moments of a solute plume are a posteriori obtained by a consistent Lagrangian analysis starting from the knowledge of the velocity field covariance matrices without any restriction regarding the statistical homogeneity of flow and/or ergodicity conditions. The proposed method is here applied to investigate the influence of boundary conditions on the solute transport developing in a limited domain. The results obtained from different test cases give a velocity field behavior in agreement with literature findings and allow us to define the nonergodic transport in bounded domains. The comparison of the results with some ad hoc developed Monte Carlo simulations ensures the applicability of the proposed approach at least to the small heterogeneity cases here considered.
Nonstationary flow and non-egodic transport in random porous media
SALANDIN, PAOLO
2006
Abstract
Transport processes taking place in natural formations are often characterized by the spatial nonstationarity of flow field and the nonergodic condition of solute spreading. The former may originate from statistical inhomogeneity of porous media properties or from the influence of boundary conditions as well as from conditioning procedures on measured data, while the latter is generally due to the finite size of the solute source. In real-world applications these aspects sometimes dominate the transport phenomena, and known literature results based on statistical homogeneity and ergodicity assumptions cannot be properly applied. In this paper a method is proposed to handle different and concurrent causes of the flow field nonstationarity and to describe nonergodic transport of inert solutes by spatial moments in a domain of finite size. The goal is reached by expanding the steady state flow equation in Taylor series limited to the first-order and by the recursive application of finite element (FE) method. The unknowns are the piezometric head mean values and its derivatives in respect to the fluctuating porous media hydraulic conductivity. From the latter the velocity statistics are derived on the basis of Darcy's law according to the first-order approximation. By neglecting the pore-scale dispersion the spatial moments of a solute plume are a posteriori obtained by a consistent Lagrangian analysis starting from the knowledge of the velocity field covariance matrices without any restriction regarding the statistical homogeneity of flow and/or ergodicity conditions. The proposed method is here applied to investigate the influence of boundary conditions on the solute transport developing in a limited domain. The results obtained from different test cases give a velocity field behavior in agreement with literature findings and allow us to define the nonergodic transport in bounded domains. The comparison of the results with some ad hoc developed Monte Carlo simulations ensures the applicability of the proposed approach at least to the small heterogeneity cases here considered.Pubblicazioni consigliate
I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.