Pietro de Anna

Mixing and reaction kinetics in porous media: an experimental pore scale quantification.

We propose a new experimental set up to characterize mixing and reactive transport in porous media with a high spatial resolution at the pore scale. The analogous porous medium consists of a Hele-Shaw cell containing a single layer of cylindrical solid grains built by soft lithography. On the one hand, the measurement of the local, intrapore, conservative concentration field is done using a fluorescent tracer. On the other hand, considering a fast bimolecular reaction A + B → C occurring as A displaces B, we quantify the rate of product formation from the spatially resolved measurement of the pore scale reaction rate, using a chemiluminescent reaction. The setup provides a dynamical measurement of the local concentration field over 3 orders of magnitude and allows investigating a wide range of Péclet and Damköhler numbers by varying the flow rate within the cell and the local reaction rate. We use it to study the kinetics of the reaction front between A and B. While the advection-dispersion (Fickian) theory, applied at the continuum scale, predicts a scaling of the cumulative mass of product C as MC ∝ √t, the experiments exhibit two distinct regimes in which the produced mass MC evolves faster than the Fickian behavior. In both regimes the front rate of product formation is controlled by the geometry of the mixing interface between the reactants. Initially, the invading solute is organized in stretched lamellae and the reaction is limited by mass transfer across the lamella boundaries. At longer times the front evolves into a second regime where lamellae coalesce and form a mixing zone whose temporal evolution controls the rate of product formation. In this second regime, the produced mass of C is directly proportional to the volume of the mixing zone defined from conservative species. This interesting property is indeed verified from a comparison of the reactive and conservative data. Hence, for both regimes, the direct measurement of the spatial distribution of the pore scale reaction rate and conservative component concentration is shown to be crucial to understanding the departure from the Fickian scaling as well as quantifying the basic mechanisms that govern the mixing and reaction dynamics at the pore scale.

Pore-scale intermittent velocity structure underpinning anomalous transport through 3-D porous media

We study the nature of non-Fickian particle transport in 3-D porous media by simulating fluid flow in the intricate pore space of real rock. We solve the full Navier-Stokes equations at the same resolution as the 3-D micro-CT (computed tomography) image of the rock sample and simulate particle transport along the streamlines of the velocity field. We find that transport at the pore scale is markedly anomalous: longitudinal spreading is superdiffusive, while transverse spreading is subdiffusive. We demonstrate that this anomalous behavior originates from the intermittent structure of the velocity field at the pore scale, which in turn emanates from the interplay between velocity heterogeneity and velocity correlation. Finally, we propose a continuous time random walk model that honors this intermittent structure at the pore scale and captures the anomalous 3-D transport behavior at the macroscale.

The filamentary structure of mixing fronts and its control on reaction kinetics in porous media flows

The mixing dynamics resulting from the combined action of diffusion, dispersion, and advective stretching of a reaction front in heterogeneous flows leads to reaction kinetics that can differ by orders of magnitude from those measured in well-mixed batch reactors. The reactive fluid invading a porous medium develops a filamentary or lamellar front structure. Fluid deformation leads to an increase of the front length by stretching and consequently a decrease of its width by compression. This advective front deformation, which sharpens concentration gradients across the interface, is in competition with diffusion, which tends to increase the interface width and thus smooth concentration gradients. The lamella scale dynamics eventually develop into a collective behavior through diffusive coalescence, which leads to a disperse interface whose width is controlled by advective dispersion. We derive a new approach that quantifies the impact of these filament scale processes on the global mixing and reaction kinetics. The proposed reactive filament model, based on the elementary processes of stretching, coalescence, and fluid particle dispersion, provides a new framework for predicting reaction front kinetics in heterogeneous flows.

Anomalous kinetics in diffusion limited reactions linked to non-Gaussian concentration probability distribution function

We investigate anomalous reaction kinetics related to segregation in the one-dimensional reaction-diffusion system A + B → C. It is well known that spatial fluctuations in the species concentrations cause a breakdown of the mean-field behavior at low concentration values. The scaling of the average concentration with time changes from the mean-field t(-1) to the anomalous t(-1/4) behavior. Using a stochastic modeling approach, the reaction-diffusion system can be fully characterized by the multi-point probability distribution function (PDF) of the species concentrations. Its evolution is governed by a Fokker-Planck equation with moving boundaries, which are determined by the positivity of the species concentrations. The concentration PDF is in general non-Gaussian. As long as the concentration fluctuations are small compared to the mean, the PDF can be approximated by a Gaussian distribution. This behavior breaks down in the fluctuation dominated regime, for which anomalous reaction kinetics are observed. We show that the transition from mean field to anomalous reaction kinetics is intimately linked to the evolution of the concentration PDF from a Gaussian to non-Gaussian shape. This establishes a direct relationship between anomalous reaction kinetics, incomplete mixing and the non-Gaussian nature of the concentration PDF.

Spatial model of autocatalytic reactions.

Biological cells with all of their surface structure and complex interior stripped away are essentially vesicles--membranes composed of lipid bilayers which form closed sacs. Vesicles are thought to be relevant as models of primitive protocells, and they could have provided the ideal environment for prebiotic reactions to occur. In this paper, we investigate the stochastic dynamics of a set of autocatalytic reactions, within a spatially bounded domain, so as to mimic a primordial cell. The discreteness of the constituents of the autocatalytic reactions gives rise to large sustained oscillations even when the number of constituents is quite large. These oscillations are spatiotemporal in nature, unlike those found in previous studies, which consisted only of temporal oscillations. We speculate that these oscillations may have a role in seeding membrane instabilities which lead to vesicle division. In this way synchronization could be achieved between protocell growth and the reproduction rate of the constituents (the protogenetic material) in simple protocells.

Interface evolution during radial miscible viscous fingering.

We study experimentally the miscible radial displacement of a more viscous fluid by a less viscous one in a horizontal Hele-Shaw cell. For the range of tested injection rates and viscosity ratios we observe two regimes for the evolution of the fluid-fluid interface. At early times the interface length increases linearly with time, which is typical of the Saffman-Taylor instability for this radial configuration. However, as time increases, the interface growth slows down and scales as ∼t(1/2), as one expects in a stable displacement, indicating that the overall flow instability has shut down. Surprisingly, the crossover time between these two regimes decreases with increasing injection rate. We propose a theoretical model that is consistent with our experimental results, explains the origin of this second regime, and predicts the scaling of the crossover time with injection rate and the mobility ratio. The key determinant of the observed scalings is the competition between advection and diffusion time scales at the displacement front, suggesting that our analysis can be applied to other interfacial-evolution problems such as the Rayleigh-Bénard-Darcy instability.

Pore-scale mechanisms for the enhancement of mixing in unsaturated porous media and implications for chemical reactions: MIXING IN UNSATURATED POROUS MEDIA

Porous media in which different fluid phases coexist are common in nature (e.g., vadose zone and gas-oil reservoirs). In partially saturated porous media, the intricate spatial distributions of the wetting and nonwetting phases causes their flow to be focused onto preferential paths. Using a novel 2-D experimental setup allowing pore-scale measurement of concentration fields in a controlled unsaturated flow, we highlight mechanisms by which mixing of an invading fluid with the resident fluid is significantly enhanced when decreasing saturation. The mean scalar dissipation rate is observed to decrease slowly in time, while under saturated conditions it decays rapidly. This slow decrease is due to sustained longitudinal solute fingering, which causes concentration gradients to remain predominantly transverse to the average flow. Consequently, the effective reactivity is found to be much larger than under saturated conditions. These results provide new insights into the role that multiphase flows play on mixing/reaction in porous media.