Source author record

Yves Méheust

Yves Méheust appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

physics.flu-dyn physics.geo-ph cond-mat.soft physics.comp-ph

Catalog footprint

What is connected

6works

4topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2021arXiv

Convective dissolution of CO$_2$ in 2D and 3D porous media: the impact of hydrodynamic dispersion

Convective dissolution is the process by which CO$_2$ injected in deep geological formations dissolves into the aqueous phase, which allows storing it perennially by gravity. The process results from buoyancy-coupled Darcy flow and solute transport. Proper theoretical modeling of the process should consider in the transport equation a diffusive term accounting for hydrodynamics (or, mechanical) dispersion, with an effective diffusion coefficient that is proportional to the local interstitial velocity. A few two-dimensional (2D) numerical studies, and three-dimensional (3D) experimental investigations, have investigated the impact of hydrodynamic dispersion on convection dynamics, with contradictory conclusions. Here, we investigate systematically the impact of the dispersion strength $S$ (relative to molecular diffusion), and of the anisotropy $α$ of its tensor, on convective dissolution in 2D and 3D geometries. We use a new numerical model and analyze the solute fingers' number density (FND), penetration depth and maximum velocity; the onset time of convection; the dissolution flux in the quasi-constant flux regime; the mean concentration of the dissolved CO2; and the scalar dissipation rate. The efficiency of convective dissolution over long times is observed to be mostly controlled by the onset time of convection. For most natural porous media ($α= 0.1$), the onset time is found to increase as a function of $S$, in agreement with previous experimental findings and in stark contrast to previous numerical findings. However, if $α$ is sufficiently large this behavior is reversed. Furthermore, results in 3D are fully consistent with the 2D results on all accounts, except that in 3D the onset time is slightly smaller, the dissolution flux in the quasi-constant flux regime is slightly larger, and the dependence of the FND on the dispersion parameters is impacted by $Ra$.

preprint2021arXiv

Sharp transition to strongly anomalous transport in unsaturated porous media

The simultaneous presence of liquid and gas in porous media increases flow heterogeneity compared to saturated flows. However, so far the impact of saturation on flow statistics and transport dynamics remained unclear. Here, we develop a theoretical framework that captures the impact of flow reorganization on the statistics of pore-scale fluid velocities, due to the presence of gas in the pore space, which leads to the development of a highly-structured flow field. Preferential flow is distributed spatially through the denoted backbone and flow recirculation occurs in dead-end regions branching from it. This induces a marked change in the scaling of the velocity PDF compared to the saturated case, and a sharp transition to strongly anomalous transport. We develop a transport model based on the continuous time random walk theory that successfully predicts advective transport dynamics for all saturation degrees. Our results provide a new modelling framework linking phase heterogeneity to flow heterogeneity and to transport in unsaturated media.

preprint2016arXiv

Enhanced reaction kinetics and reactive mixing scale dynamics in mixing fronts under shear flow for arbitrary Damköhler 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éclet number $Pe$, which represents the ratio of the characteristic diffusion time to the characteristic advection time, and the Damköhler 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öhler 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.

preprint2016arXiv

Reactive strip method for mixing and reaction in two dimensions

A numerical method to efficiently solve for mixing and reaction of scalars in a two-dimensional flow field at large Péclet numbers but otherwise arbitrary Damköhler numbers is reported. We consider a strip of one reactant in a pool of another reactant, both of which are advected with the known velocity field. We first establish that the system evolution for such a system under certain conditions is described by a locally one-dimensional reaction-diffusion problem. The approximation of a locally one-dimensional dynamics is true for cases where the strip thickness is smaller than the local radius of curvature and also when the strip thickness is smaller than the distance between adjacent strips. We first demonstrate the method for the transport of a conservative scalar under a linear shear flow, point vortex and a chaotic sine flow. We then proceed to consider the situation with a simple bimolecular reaction between two reactants to yield a single product. The methodology presented herewith essentially generalizes nontrivially the Diffusive Strip Method developed by Meunier and Villermaux (J. Fluid Mech. 662, 134-172 (2010)) to address passive scalar transport, to the generalized situation with multiple reacting species. In essence, the reduction of dimensionality of the problem, which renders the 2D problem 1D, allows one to efficiently model reactive transport under high Péclet numbers which are otherwise prohibitively difficult to resolve from classical finite difference or finite element based methods.

preprint2013arXiv

Structure-dependent mobility of a dry aqueous foam flowing along two parallel channels

The velocity of a two-dimensional aqueous foam has been measured as it flows through two parallel channels, at a constant overall volumetric flow rate. The flux distribution between the two channels is studied as a function of the ratio of their widths. A peculiar dependence of the velocity ratio on the width ratio is observed when the foam structure in the narrower channel is either single staircase or bamboo. In particular, discontinuities in the velocity ratios are observed at the transitions between double and single staircase and between single staircase and bamboo. A theoretical model accounting for the viscous dissipation at the solid wall and the capillary pressure across a film pinned at the channel outlet predicts the observed non-monotonic evolution of the velocity ratio as a function of the width ratio. It also predicts quantitatively the intermittent temporal evolution of the velocity in the narrower channel when it is so narrow that film pinning at its outlet repeatedly brings the flow to a near stop.

preprint2012arXiv

Biphasic flow: structure and upscaling, consequences on macroscopic transport properties

In disordered porous media, two-phase flow of immiscible fluids (biphasic flow) is organized in patterns that sometimes exhibit fractal geometries over a range of length scales, depending on the capillary, gravitational and viscous forces at play. These forces, as well as the boundary conditions, also determine whether the flow leads to the appearance of fingering pathways, i.e., unstable flow, or not. We present here a short review of these aspects, focusing on drainage and summarizing when these flows are expected to be stable or not, what fractal dimensions can be expected, and in which range of scales. We base our review on experimental studies performed in two-dimensional Hele-Shaw cells, or addressing three dimensional porous media by use of several imaging techniques. We first present configurations in which solely capillary forces and gravity play a role. Next, we review configurations in which capillarity and viscosity are the main forces at play. Eventually, we examine how the microscopic geometry of the fluid clusters affects the macroscopic transport properties. An example of such an upscaling is illustrated in detail: For air invasion in a mono-layer glass-bead cell, the fractal dimension of the flow structures and the associated scale-ranges, are shown to depend on the displacement velocity. This controls the relationship between saturation and the pressure difference between the two phases at the macroscopic scale. We provide in this case expressions for dynamic capillary pressure and residual fluid phase saturations.