Chaotic Mixing in Three Dimensional Porous Media

Under steady flow conditions, the topological complexity inherent to all random 3D porous media imparts complicated flow and transport dynamics. It has been established that this complexity generates persistent chaotic advection via a three-dimensional (3D) fluid mechanical analogue of the baker's map which rapidly accelerates scalar mixing in the presence of molecular diffusion. Hence pore-scale fluid mixing is governed by the interplay between chaotic advection, molecular diffusion and the broad (power-law) distribution of fluid particle travel times which arise from the non-slip condition at pore walls. To understand and quantify mixing in 3D porous media, we consider these processes in a model 3D open porous network and develop a novel stretching continuous time random walk (CTRW) which provides analytic estimates of pore-scale mixing which compare well with direct numerical simulations. We find that chaotic advection inherent to 3D porous media imparts scalar mixing which scales exponentially with longitudinal advection, whereas the topological constraints associated with 2D porous media limits mixing to scale algebraically. These results decipher the role of wide transit time distributions and complex topologies on porous media mixing dynamics, and provide the building blocks for macroscopic models of dilution and mixing which resolve these mechanisms.Comment: 36 page

Non-Fickian dispersion in porous media : 2. Model validation from measurements at different scales

International audienceWe aim at testing and validating a mobile-immobile mass transfer model from a set of single-well injection withdrawal tracer tests in a heterogeneous porous aquifer. By varying the duration of the injection phase, different volumes of aquifer are investigated by the tracer. Hence, we focus the transport model validation not only on reproducing a single breakthrough curve (BTC) but also on the model's capacity to predict the amount of mixing as a function of the volume visited by the tracer. All the BTCs are strongly asymmetric, as expected when dispersion is controlled by diffusive mass transfers between the mobile water and the immobile water part of the porosity. However, the BTC cannot be modeled by a conventional mobile-immobile mass transfer model with a simple power law memory function. To account for that, we implement a continuous time random walk model in which the transition time distribution y (t), which is related to the excursion time probability of the tracer in the immobile domain, is a dual-slope power law distribution. The model best fits the BTC data set with a transitional regime controlled by y(t) t2 and an asymptotic regime characteristic of the conventional double-porosity model with y(t) t1.5 . This work emphasizes that high-resolution concentration measurement and multiple-scale tracer tests are required for assessing solute dispersion models in heterogeneous reservoirs and for subsequently obtaining reliable predictions

Enhanced reaction kinetics and reactive mixing scale dynamics in mixing fronts under shear flow for arbitrary Damk\"ohler numbers

Mixing fronts, where fluids of different chemical compositions mix with each other, are typically subjected to velocity gradients, ranging from the pore scale to the catchment scale due to permeability variations and flow line geometries. A common trait of these processes is that the mixing interface is strained by shear. Depending on the P\'eclet number $Pe$, which represents the ratio of the characteristic diffusion time to the characteristic advection time, and the Damk\"ohler number $Da$, which represents the ratio of the characteristic diffusion time to the characteristic reaction time, the local reaction rates can be strongly impacted by the dynamics of the mixing interface. This impact has been characterized mostly either in kinetics-limited or in mixing-limited conditions, that is, for either very low or very high $Da$. Here the coupling of shear flow and chemical reactivity is investigated for arbitrary Damk\"ohler numbers, for a bimolecular reaction and an initial interface with separated reactants. Approximate analytical expressions for the global production rate and reactive mixing scale are derived based on a reactive lamella approach that allows for a general coupling between stretching enhanced mixing and chemical reactions. While for $Pe<Da$, reaction kinetics and stretching effects are decoupled, a scenario which we name "weak stretching", for $Pe>Da$, we uncover a "strong stretching" scenario where new scaling laws emerge from the interplay between reaction kinetics, diffusion, and stretching. The analytical results are validated against numerical simulations. These findings shed light on the effect of flow heterogeneity on the enhancement of chemical reaction and the creation of spatially localized hotspots of reactivity for a broad range of systems ranging from kinetic limited to mixing limited situations

Hypermixing in linear shear flow

International audience[1] In this technical note we study mixing in a two‐dimensional linear shear flow. We derive analytical expressions for the concentration field for an arbitrary initial condition in an unbounded two‐dimensional shear flow. We focus on the solution for a point initial condition and study the evolution of (1) the second centered moments as a measure for the plume dispersion, (2) the dilution index as a measure of the mixing state, and (3) the scalar dissipation rate as a measure for the rate of mixing. It has previously been shown that the solute spreading grows with the cube of time and thus is hyperdispersive. Herein we demonstrate that the dilution index increases quadratically with time in contrast to a homogeneous medium, for which it increases linearly. Similarly, the scalar dissipation rate decays as t−3, while for a homogeneous medium it decreases more slowly as t−2. Mixing is much stronger than in a homogeneous medium, and therefore we term the observed behavior hypermixing

Mixing as a correlated aggregation process

Mixing describes the process by which scalars, such as solute concentration or fluid temperature, evolve from an initial heterogeneous state to uniformity under the stirring action of a fluid flow. Mixing occurs initially through the formation of scalar lamellae as a result of fluid stretching and later by their coalescence due to molecular diffusion. Owing to the linearity of the advection-diffusion equation, scalar coalescence can be envisioned as an aggregation process. While random aggregation models have been shown to capture scalar mixing across a range of turbulent flows, we demonstrate here that they are not accurate for most chaotic flows. In particular, we show that the spatial distribution of the number of lamellae in aggregates is highly correlated with their elongation and is also influenced by the fractal geometry that arises from the chaotic flow. The presence of correlations makes mixing less efficient than a completely random aggregation process because lamellae with similar elongations and scalar levels tend to remain isolated from each other. Based on these observations, we propose a correlated aggregation framework that captures the asymptotic mixing dynamics of chaotic flows and predicts the evolution of the scalar pdf based on the flow stretching statistics. We show that correlated aggregation is uniquely determined by a single exponent which quantifies the effective number of random aggregation events, and is dependent on the fractal dimension of the flow. These findings expand aggregation theories to a larger class of systems, which have relevance to various fundamental and applied mixing problems

Anomalous transport on regular fracture networks: Impact of conductivity heterogeneity and mixing at fracture intersections

We investigate transport on regular fracture networks that are characterized by heterogeneity in hydraulic conductivity. We discuss the impact of conductivity heterogeneity and mixing within fracture intersections on particle spreading. We show the emergence of non-Fickian transport due to the interplay between the network conductivity heterogeneity and the degree of mixing at nodes. Specifically, lack of mixing at fracture intersections leads to subdiffusive scaling of transverse spreading but has negligible impact on longitudinal spreading. An increase in network conductivity heterogeneity enhances both longitudinal and transverse spreading and leads to non-Fickian transport in longitudinal direction. Based on the observed Lagrangian velocity statistics, we develop an effective stochastic model that incorporates the interplay between Lagrangian velocity correlation and velocity distribution. The model is parameterized with a few physical parameters and is able to capture the full particle transition dynamics.United States. Dept. of Energy (Grant DE-SC0003907)MISTI (Hayashi Seed Fund

Scaling forms of particle densities for Lévy walks and strong anomalous diffusion

International audienceWe study the scaling behavior of particle densities for Lévy walks whose transition length r is coupled with the transition time t as |r| ∝ t α with an exponent α > 0. The transition-time distribution behaves as ψ(t) ∝ t −1−β with β > 0. For 1 q c. These results give insight into the possible origins of strong anomalous diffusion and anomalous behaviors in disordered systems in general