Source author record

Marco Dentz

Marco Dentz appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

9works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

9 published item(s)

preprint2022arXiv

Dispersion of motile bacteria in a porous medium

Understanding flow and transport of bacteria in porous media is crucial to technologies such as bioremediation, biomineralization or enhanced oil recovery. While physicochemical bacteria filtration is well-documented, recent studies showed that bacterial motility plays a key role in the transport process. Flow and transport experiments performed in microfluidic chips containing randomly placed obstacles confirmed that the distributions of non-motile particles stays compact, whereas for the motile strains, the distributions are characterized by both significant retention as well as fast downstream motion. For motile bacteria, the detailed microscopic study of individual bacteria trajectories reveals two salient features: (i) the emergence of an active retention process triggered by motility, (ii) enhancement of dispersion due to the exchange between fast flow channels and low flow regions in the vicinity of the solid grains. We propose a physical model based on a continuous time random walk approach. This approach accounts for bacteria dispersion via variable pore-scale flow velocities through a Markov model for equidistant particle speeds. Motility of bacteria is modeled by a two-rate trapping process that accounts for the motion towards and active trapping at the obstacles. This approach captures the forward tails observed for the distribution of bacteria displacements, and quantifies an enhanced hydrodynamic dispersion effect that originates in the interaction between flow at the pore-scale and bacterial motility. The model reproduces the experimental observations, and predicts bacteria dispersion and transport at the macroscale.

preprint2020arXiv

Heterogeneous Multi-Rate mass transfer models in OpenFOAM

We implement the Multi-Rate Mass Transfer (MRMT) model for mobile-immobile transport in porous media within the open-source finite volume library \textsc{OpenFOAM}\reg \citep{Foundation2014}. Unlike other codes available in the literature [Geiger, S., Dentz, M., Neuweiler, I., SPE Reservoir Characterisation and Simulation Conference and Exhibition (2011); Silva, O., Carrera, J., Dentz, M., Kumar, S., Alcolea, A., Willmann, M., Hydrology and Earth System Sciences 13, (2009)], we propose an implementation that can be applied to complex three-dimensional geometries and highly heterogeneous fields, where the parameters of the MRMT can arbitrarily vary in space. Furthermore, being built over the widely diffused OpenFOAM\reg library, it can be easily extended and included in other models, and run in parallel. We briefly describe the structure of the < multiContinuumModels > library that includes the formulation of the MRMT based on the works of [Haggerty, R., Gorelick, S.M., Water Resources Research 31, (1995)] and [F. Municchi and M. Icardi Phys. Rev. Research 2, 013041, (2020)]. The implementation is verified against benchmark solutions and tested on two- and three-dimensional random permeability fields. The role of various physical and numerical parameters, including the transfer rates, the heterogeneities, and the number of terms in the MRMT expansions is investigated. Finally, we illustrate the significant role played by heterogeneity in the mass transfer when permeability and porosity are represented using Gaussian random fields.

preprint2020arXiv

Probability density function (PDF) models for particle transport in porous media

Mathematical models based on probability density functions (PDF) have been extensively used in hydrology and subsurface flow problems, to describe the uncertainty in porous media properties (e.g., permeability modelled as random field). Recently, closer to the spirit of PDF models for turbulent flows, some approaches have used this statistical viewpoint also in pore-scale transport processes (fully resolved porous media models). When a concentration field is transported, by advection and diffusion, in a heterogeneous media, in fact, spatial PDFs can be defined to characterise local fluctuations and improve or better understand the closures performed by classical upscaling methods. In the study of hydrodynamical dispersion, for example, PDE-based PDF approach can replace expensive and noisy Lagrangian simulations (e.g. trajectories of drift-diffusion stochastic processes). In this work we derive a joint position-velocity Fokker-Planck equation to model the motion of particles undergoing advection and diffusion in in deterministic or stochastic heterogeneous velocity fields. After appropriate closure assumptions, this description can help deriving rigorously stochastic models for the statistics of Lagrangian velocities. This is very important to be able to characterise the dispersion properties and can, for example, inform velocity evolution processes in Continuous Time Random Walk (CTRW) dispersion models. The closure problem that arises when averaging the Fokker Planck equation shows also interesting similarities with the mixing problem and can be used to propose alternative closures for anomalous dispersion.

preprint2016arXiv

Continuous Time Random Walks for Non-Local Radial Solute Transport

This paper derives and analyzes continuous time random walk (CTRW) models in radial flow geometries for the quantification of non-local solute transport induced by heterogeneous flow distributions and by mobile-immobile mass transfer processes. To this end we derive a general CTRW framework in radial coordinates starting from the random walk equations for radial particle positions and times. The particle density, or solute concentration is governed by a non-local radial advection-dispersion equation (ADE). Unlike in CTRWs for uniform flow scenarios, particle transition times here depend on the radial particle position, which renders the CTRW non-stationary. As a consequence, the memory kernel characterizing the non-local ADE, is radially dependent. Based on this general formulation, we derive radial CTRW implementations that (i) emulate non-local radial transport due to heterogeneous advection, (ii) model multirate mass transfer (MRMT) between mobile and immobile continua, and (iii) quantify both heterogeneous advection in a mobile region and mass transfer between mobile and immobile regions. The expected solute breakthrough behavior is studied using numerical random walk particle tracking simulations. This behavior is analyzed by explicit analytical expressions for the asymptotic solute breakthrough curves. We observe clear power-law tails of the solute breakthrough for broad (power-law) distributions of particle transit times (heterogeneous advection) and particle trapping times (MRMT model). The combined model displays two distinct time regimes. An intermediate regime, in which the solute breakthrough is dominated by the particle transit times in the mobile zones, and a late time regime that is governed by the distribution of particle trapping times in immobile zones.

preprint2016arXiv

Continuous Time Random Walks for the Evolution of Lagrangian Velocities

We develop a continuous time random walk (CTRW) approach for the evolution of Lagrangian velocities in steady heterogeneous flows based on a stochastic relaxation process for the streamwise particle velocities. This approach describes persistence of velocities over a characteristic spatial scale, unlike classical random walk methods, which model persistence over a characteristic time scale. We first establish the relation between Eulerian and Lagrangian velocities for both equidistant and isochrone sampling along streamlines, under transient and stationary conditions. Based on this, we develop a space continuous CTRW approach for the spatial and temporal dynamics of Lagrangian velocities. While classical CTRW formulations have non-stationary Lagrangian velocity statistics, the proposed approach quantifies the evolution of the Lagrangian velocity statistics under both stationary and non-stationary conditions. We provide explicit expressions for the Lagrangian velocity statistics, and determine the behaviors of the mean particle velocity, velocity covariance and particle dispersion. We find strong Lagrangian correlation and anomalous dispersion for velocity distributions which are tailed toward low velocities as well as marked differences depending on the initial conditions. The developed CTRW approach predicts the Lagrangian particle dynamics from an arbitrary initial condition based on the Eulerian velocity distribution and a characteristic correlation scale.

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

Fluid Stretching as a Levy Process

We study the relation between flow structure and fluid deformation in steady two-dimensional random flows. Beyond the linear (shear flow) and exponential (chaotic flow) elongation paradigms, we find a broad spectrum of stretching behaviors, ranging from sub- to superlinear, which are dominated by intermittent shear events. We analyze these behaviors from first principles, which uncovers stretching as a result of the non-linear coupling between Lagrangian shear deformation and velocity fluctuations along streamlines. We derive explicit expressions for Lagrangian deformation and demonstrate that stretching obeys a coupled continous time random walk, which for broad distributions of flow velocities describes a Lévy walk for elongation. The derived model provides a direct link between the flow and deformation statistics, and a natural way to quantify the impact of intermittent shear events on the stretching behavior, which can have strong anomalous diffusive character.

preprint2016arXiv

Multi-rate mass transfer modeling of two-phase flow in highly heterogeneous fractured and porous media

We study modeling of two-phase flow in highly heterogeneous fractured and porous media. The flow behaviour is strongly influenced by mass transfer between a highly permeable (mobile) fracture domain and less permeable (immobile) matrix blocks. We quantify the effective two-phase flow behaviour using a multirate rate mass transfer (MRMT) approach. We discuss the range of applicability of the MRMT approach in terms of the pertinent viscous and capillary diffusion time scales. We scrutinize the linearization of capillary diffusion in the immobile regions, which allows for the formulation of MRMT in the form of a non-local single equation model. The global memory function, which encodes mass transfer between the mobile and the immobile regions, is at the center of this method. We propose two methods to estimate the global memory function for a fracture network with given fracture and matrix geometry. Both employ a scaling approach based on the known local memory function for a given immobile region. With the first method, the local memory function is calculated numerically, while the second one employs a parametric memory function in form of truncated power-law. The developed concepts are applied and tested for fracture networks of different complexity. We find that both physically based parameter estimation methods for the global memory function provide predictive MRMT approaches for the description of multiphase flow in highly heterogeneous porous media.

preprint2016arXiv

Pictures of Blockscale Transport: Effective versus Ensemble Dispersion and Its Uncertainty

Using a Lagrangian framework, we discuss different conceptualizations of the blockscale dispersion tensor. We distinguish effective and ensemble blockscale dispersion, which measure the impact of subscale velocity fluctuations on solute dispersion. Ensemble dispersion quantifies subscale velocity fluctuations between realizations, which overestimates the actual velocity variability. Effective dispersion on the other hand quantifies the actual blockscale velocity variability and thus reflects the impact of subscale velocity fluctuations on mixing and spreading. Based on these concepts, we quantify the impact of subscale velocity fluctuations on solute particle spreading and determine the governing equations for the coarse-grained concentration distributions. We develop analytical and semi-analytical expressions for the average and variance of the blockscale dispersion tensor in 3D flow fields as a function of the structural parameters characterizing the subsurface. Our results illustrate the relevance of the blockscale, the initial scale of the solute body and local-scale dispersion in controlling the uncertainty of the plume's dispersive behavior. The analysis performed in this work has implications in numerical modeling (i.e. grid design) and allows to quantify the uncertainty of the blockscale dispersion tensor.