Analysis of an Euler implicit‐mixed finite element scheme for reactive solute transport in porous media

In this article, we analyze an Euler implicit‐mixed finite element scheme for a porous media solute transport model. The transporting flux is not assumed given, but obtained by solving numerically the Richards equation, a model for subsurface fluid flow. We prove the convergence of the scheme by estimating the error in terms of the discretization parameters. In doing so we take into account the numerical error occurring in the approximation of the fluid flow. The article, is concluded by numerical experiments, which are in good agreement with the theoretical estimates. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

Error estimates for an Euler implicit-mixed finite element scheme for reactive transport in saturated/unsaturated soil

We present convergence results for a fully discrete scheme based on the mixed finite element (MFE) method and an one-step Euler implicit (EI) method for simulating reactive solute transport in saturated/unsaturated soil. The results considered the low regularity of the solution of the degenerate parabolic equation describing the water flow in porous media. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Large Scale Mixing for Immiscible Displacement in Heterogeneous Porous Media

We derive a large scale mixing parameter for a displacement process of one fluid by another immiscible one in a two-dimensional heterogeneous porous medium. The mixing of the displacing fluid saturation due to the heterogeneities of the permeabilities is captured by a dispersive flux term in the large scale homogeneous flow equation. By making use of the stochastic approach we develop a definition of the dispersion coefficient and apply a Eulerian perturbation theory to determine explicit results to second order in the fluctuations of the total velocity. We apply this method to a uniform flow configuration as well as to a radial one. The dispersion coefficient is found to depend on the mean total velocity and can therefore be time varying. The results are compared to numerical multi-realization calculations. We found that the use of single phase flow stochastics cannot capture all phenomena observed in the numerical simulations.     

Source-like solution for radial imbibition into a homogeneous semi-infinite porous medium

We describe the imbibition process from a point source into a homogeneous semi-infinite porous material. When body forces are negligible, the advance of the wetting front is driven by capillary pressure and resisted by viscous forces. With the assumption that the wetting front assumes a hemispherical shape, our analytical results show that the absorbed volume flow rate is approximately constant with respect to time, and that the radius of the wetting evolves in time as r ≈ t(1/3). This cube-root law for the long-time dynamics is confirmed by experiments using a packed cell of glass microspheres with average diameter of 42 μm. This result complements the classical one-dimensional imbibition result where the imbibition length l ≈ t(1/2), and studies in axisymmetric porous cones with small opening angles where l ≈ t(1/4) at long times.