Anomalous transport and chaotic advection in homogeneous porous media

Chaotic Transport of Solutes in Unsaturated Porous Media

Unsaturated porous media, characterized by the combined presence of several immiscible fluid phases in the pore space, are highly relevant systems in nature, because they control the fate of contaminants and the availability of nutrients in the subsoil. However, a full understanding of the mechanisms controlling solute mixing in such systems is still missing. In particular, the role of saturation in the development of chaotic solute mixing has remained unexplored. Using three-dimensional numerical simulations of flow and transport at the pore scale, built upon X-ray tomograms of a porous medium at different degrees of liquid (wetting)-phase saturation, we show the occurrence of chaotic dynamics in both the deformation of the solute plume, as characterized by computed chaos metrics (Lyapunov exponents), and the mixing of the injected solute. Our results show an enhancement of these chaotic dynamics at lower saturation and their occurrence even under diffusion-relevant conditions over the medium's length, also being strengthened by larger flow velocities. These findings highlight the dominant role of the pore-scale spatial heterogeneity of the system, enhanced by the presence of an immiscible phase (e.g., air), on the mixing efficiency. This represents a stepping stone for the assessment of mixing and reactions in unsaturated porous media.

Mixing in Porous Media: Concepts and Approaches Across Scales

This review provides an overview of concepts and approaches for the quantification of passive, non-reactive solute mixing in steady uniform porous media flows across scales. Mixing in porous media is the result of the interaction of spatial velocity fluctuations and diffusion or local-scale dispersion, which may lead to the homogenization of an initially segregated system. Velocity fluctuations are induced by spatial medium heterogeneities at the pore, Darcy or regional scales. Thus, mixing in porous media is a multiscale process, which depends on the medium structure and flow conditions. In the first part of the review, we discuss the interrelated processes of stirring, dispersion and mixing, and review approaches to quantify them that apply across scales. This implies concepts of hydrodynamic dispersion, approaches to quantify mixing state and mixing dynamics in terms of concentration statistics, and approaches to quantify the mechanisms of mixing. We review the characterization of stirring in terms of fluid deformation and folding and its relation with hydrodynamic dispersion. The integration of these dynamics to quantify the mechanisms of mixing is discussed in terms of lamellar mixing models. In the second part of this review, we discuss these concepts and approaches for the characterization of mixing in Poiseuille flow, and in porous media flows at the pore, Darcy and regional scales. Due to the fundamental nature of the mechanisms and processes of mixing, the concepts and approaches discussed in this review underpin the quantitative analysis of mixing phenomena in porous media flow systems in general.

Modeling Mixing in Stratified Heterogeneous Media: The Role of Water Velocity Discretization in Phase Space Formulation

Modeling solute transport in heterogeneous porous media faces two challenges: scale dependence of dispersion and reproducing mixing separately from spreading. Both are crucial since real applications may require km scales whereas reactions, often controlled by mixing, may occur at the pore scale. Methods have been developed in response to these challenges, but none has satisfactorily characterized both processes. In this paper, we propose a formulation based on the Water Mixing Approach extended to account for velocity variability. Velocity is taken as an independent variable, so that concentration depends on time, space and velocity. Therefore, we term the formulation the Multi-Advective Water Mixing Approach. A new mixing term between velocity classes emerges in this formulation. We test it on Poiseuille’s stratified flow using the Water Parcel method. Results show high accuracy of the formulation in both dispersion and mixing. Moreover, the mixing process exhibits Markovianity in space even though it is modeled in time.

A spatial Markov model for upscaling transport of adsorbing-desorbing solutes

The Spatial Markov Model (SMM) is an upscaled model with a strong track record in predicting upscaled behavior of conservative solute transport across hydrologic systems. Here we propose an SMM that can account for reactive linear adsorption and desorption processes and test it on a simple benchmark problem: flow and transport through an idealized periodic wavy channel. The methodology is built using trajectories that are obtained from a single high resolution random walk simulation of conservative transport across one periodic element. Our approach encodes information about where a particle starts at the inlet, where it leaves at the outlet, how long it takes to cross the domain and one additional piece of information, the number of times a particle strikes the boundary, with the objective of predicting large scale transport with arbitrary linear adsorption and desorption rates. Our benchmark problem demonstrates that predictions made with our proposed SMM agree favorably with results from direct numerical simulations, which resolve the full transport problem.

Recirculation zones induce non-Fickian transport in three-dimensional periodic porous media

In this work, the influence of pore space geometry on solute transport in porous media is investigated performing computational fluid dynamics pore-scale simulations of fluid flow and solute transport. The three-dimensional periodic domains are obtained from three different pore structure configurations, namely, face-centered-cubic (fcc), body-centered-cubic (bcc), and sphere-in-cube (sic) arrangements of spherical grains. Although transport simulations are performed with media having the same grain size and the same porosity (in fcc and bcc configurations), the resulting breakthrough curves present noteworthy differences, such as enhanced tailing. The cause of such differences is ascribed to the presence of recirculation zones, even at low Reynolds numbers. Various methods to readily identify recirculation zones and quantify their magnitude using pore-scale data are proposed. The information gained from this analysis is then used to define macroscale models able to provide an appropriate description of the observed anomalous transport. A mass transfer model is applied to estimate relevant macroscale parameters (hydrodynamic dispersion above all) and their spatial variation in the medium; a functional relation describing the spatial variation of such macroscale parameters is then proposed.

Solvable continuous-time random walk model of the motion of tracer particles through porous media

We consider the continuous-time random walk (CTRW) model of tracer motion in porous medium flows based on the experimentally determined distributions of pore velocity and pore size reported by Holzner et al. [M. Holzner et al., Phys. Rev. E 92, 013015 (2015)PLEEE81539-375510.1103/PhysRevE.92.013015]. The particle's passing through one channel is modeled as one step of the walk. The step (channel) length is random and the walker's velocity at consecutive steps of the walk is conserved with finite probability, mimicking that at the turning point there could be no abrupt change of velocity. We provide the Laplace transform of the characteristic function of the walker's position and reductions for different cases of independence of the CTRW's step duration τ, length l, and velocity v. We solve our model with independent l and v. The model incorporates different forms of the tail of the probability density of small velocities that vary with the model parameter α. Depending on that parameter, all types of anomalous diffusion can hold, from super- to subdiffusion. In a finite interval of α, ballistic behavior with logarithmic corrections holds, which was observed in a previously introduced CTRW model with independent l and τ. Universality of tracer diffusion in the porous medium is considered.

Experimental investigation of transverse mixing in porous media under helical flow conditions

Plume dilution and transverse mixing can be considerably enhanced by helical flow occurring in three-dimensional heterogeneous anisotropic porous media. In this study, we perform tracer experiments in a fully three-dimensional flow-through chamber to investigate the effects of helical flow on plume spiraling and deformation, as well as on its dilution. Porous media were packed in angled stripes of materials with different grain sizes to create blocks with macroscopically anisotropic hydraulic conductivity, which caused helical flows. Steady-state transport experiments were carried out by continuously injecting dye tracers at different inlet ports. High-resolution measurements of concentration and flow rates were performed at 49 outlet ports. These measurements allowed quantifying the spreading and dilution of the solute plumes at the outlet cross section. Direct evidence of plume spiraling and visual proof of helical flow was obtained by freezing and slicing the porous media at different cross sections and observing the dye-tracer distribution. We simulated flow and transport to interpret our experimental observations and investigate the effects of helical flow on mixing-controlled reactive transport. The simulation results were evaluated using metrics of reactive mixing such as the critical dilution index and the length of continuously injected steady-state plumes. The results show considerable reaction enhancement, quantified by the remarkable decrease of reactive plume lengths (up to four times) in helical flows compared to analogous scenarios in uniform flows.

Effect of the vertical component of diffusion on passive scalar transport in an isolated vortex model

On the basis of the ellipsoidal vortex model and a Monte-Carlo-type diffusion simulation, we examine the flux and ensuing distribution of passive fluid particles through the boundary of an idealized geophysical vortex. Our focus is on features that the horizontal and vertical diffusion components introduce into the fluid particle transport. We examine the concurrent effect of both components, and we compare it with the only horizontal diffusion impact. We analyze the ellipsoid vortex model in two cases: (i) the steady state when the ellipsoid is motionless, i.e., there is no variation in its axes' lengths, and consequently the exterior fluid is not being stirred; (ii) the perturbed case when the ellipsoid rotates periodically, varying it axes' lengths, which results in the appearance of stirred fluid outside the ellipsoid. Influenced by diffusion, a fluid particle is now permitted to move to another vertical horizon, thus there is an increased possibility that the particle will eventually be located in the exterior stirred region rather than in the ellipsoid vortex with the regular dynamics. This is because the area of the horizontal section of the ellipsoid vortex decreases with depth, but the region of stirred exterior fluid extends significantly deeper. Numerical calculations show that factoring in the vertical component of diffusion significantly affects scalar spreading in the horizontal plane in the perturbed case, while in the steady state the vertical component of diffusion only induces dispersion linear growth according to a Gaussian distribution.