Upscaling of Anomalous Pore-Scale Dispersion

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.

Pore-Scale Mixing and the Evolution of Hydrodynamic Dispersion in Porous Media

We study the interplay of pore-scale mixing and network-scale advection through heterogeneous porous media, and its role for the evolution and asymptotic behavior of hydrodynamic dispersion. In a Lagrangian framework, we identify three fundamental mechanisms of pore-scale mixing that determine large scale particle motion, namely, the smoothing of intrapore velocity contrasts, the increase of the tortuosity of particle paths, and the setting of a maximum time for particle transitions. Based on these mechanisms, we derive a theory that predicts anomalous and normal hydrodynamic dispersion in terms of the characteristic pore length, Eulerian velocity distribution, and Péclet number.

A review of spatial Markov models for predicting pre-asymptotic and anomalous transport in porous and fractured media

Heterogeneity across a broad range of scales in geologic porous media often manifests in observations of non-Fickian or anomalous transport. While traditional anomalous transport models can successfully make predictions in certain geological systems, increasing evidence suggests that assumptions relating to independent and identically distributed increments constrain where and when they can be reliably applied. A relatively novel model, the Spatial Markov model (SMM), relaxes the assumption of independence. The SMM belongs to the family of correlated continuous time random walks and has shown promise across a wide range of transport problems relevant to natural porous media. It has been successfully used to model conservative as well as more recently reactive transport in highly complex flows ranging from pore scales to much larger scales of interest in geology and subsurface hydrology. In this review paper we summarize its original development and provide a comprehensive review of its advances and applications as well as lay out a vision for its future development.

Solute transport in porous media studied by lattice Boltzmann simulations at pore scale and x-ray tomography experiments

With the aid of nondestructive microfocus x-ray computed tomography (CT), we performed three-dimensional (3D) tracer dispersion experiments on randomly unconsolidated packed beds. Plumes of nonreactive sodium iodide solution were point injected into a sodium chloride solvent as a tracer for the evaluation of the dispersion process. The asymptotic dispersion coefficient was obtainable within the experimental scale and was summarized over Péclet numbers from 11.7 to ∼860. Then, the lattice Boltzmann method and moment propagation method were used to elucidate the mechanisms embedded in the dispersion phenomenon. The methods were rigorously verified against the classical Taylor dispersion problem and extended to simulate fluid flow and tracer dispersion in high-resolution 3D digital porous structures from CT. The method of moments, Lagrangian velocity correction function, and dilution index were thoroughly analyzed to evaluate the dispersion behaviors. Numerical simulations revealed ballistic and superdiffusive regimes at the transient times, whereas asymptotic dispersion behaviors appear at longer characteristic times. Besides, the observed transient times unanimously persist over convective length scales of around 12 particles transversely and 16 particles longitudinally. The estimated dispersion coefficients from simulation are in consistence with the experimental result. Furthermore, the simulation also enabled the identification of regimes, including diffusive, power law, and mechanical dispersion. Thus, the proposed experimental and computational schemes are of practical means to study dispersion behaviors by direct pore scale imaging and modeling.