1
|
Bauer D, Zaafouri Z, Batôt G, Coasne B. From Transient to Stationary Transport in Porous Networks under Various Adsorption Conditions and Kinetics. J Phys Chem B 2022; 126:6125-6135. [PMID: 35973006 DOI: 10.1021/acs.jpcb.2c02769] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We investigate the interplay between adsorption and transport in a two-dimensional porous medium by means of an extended Lattice Boltzmann technique within the Two-Relaxation-Time framework. We focus on two canonical adsorption thermodynamics and kinetics formalisms: (1) the Henry model in which the adsorbed amount scales linearly with the free adsorbate concentration and (2) the Langmuir model that accounts for surface saturation upon adsorption. We simulate transport of adsorbing and nonadsorbing particles to investigate the effect of the adsorption/desorption ratio k, initial free adsorbate concentration c0, surface saturation Γ∞, and Peclet numbers Pe on their dispersion behavior. In all cases, despite marked differences between the different adsorption models, the three following transport regimes are observed: diffusion-dominated regime, transient regime and Gaussian or nearly Gaussian dispersion regime. On the one hand, at short times, the intermediate transient regime strongly depends on the system's parameters with the shape of the concentration field at a given time being dependent on the amount of particles adsorbed shortly after injection. On the other hand, at longer times, the influence of the initial condition attenuates as particles sample sufficiently the adsorbed and nonadsorbed states. Once such dynamical equilibrium is reached, transport becomes Gaussian (i.e., normal) or nearly Gaussian in the asymptotic regime. Interestingly, the characteristic time scale to reach equilibrium, which varies drastically with the system's parameters, can be much longer than the actual simulation time. In practice, such results reflect many experimental situations such as in water treatment where dispersion is found to remain anomalous (non-Gaussian), even if transport is considered over long macroscopic times.
Collapse
Affiliation(s)
- Daniela Bauer
- IFP Energies Nouvelles, 1 & 4 Av. Bois Préau, 92852 Rueil Malmaison, France
| | - Zaineb Zaafouri
- IFP Energies Nouvelles, 1 & 4 Av. Bois Préau, 92852 Rueil Malmaison, France
| | - Guillaume Batôt
- IFP Energies Nouvelles, 1 & 4 Av. Bois Préau, 92852 Rueil Malmaison, France
| | - Benoit Coasne
- Université Grenoble Alpes, CNRS, LIPhy, 38000 Grenoble, France
| |
Collapse
|
2
|
Alexandre A, Mangeat M, Guérin T, Dean DS. How Stickiness Can Speed Up Diffusion in Confined Systems. PHYSICAL REVIEW LETTERS 2022; 128:210601. [PMID: 35687439 DOI: 10.1103/physrevlett.128.210601] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/22/2021] [Accepted: 04/15/2022] [Indexed: 06/15/2023]
Abstract
The paradigmatic model for heterogeneous media used in diffusion studies is built from reflecting obstacles and surfaces. It is well known that the crowding effect produced by these reflecting surfaces slows the dispersion of Brownian tracers. Here, using a general adsorption desorption model with surface diffusion, we show analytically that making surfaces or obstacles attractive can accelerate dispersion. In particular, we show that this enhancement of diffusion can exist even when the surface diffusion constant is smaller than that in the bulk. Even more remarkably, this enhancement effect occurs when the effective diffusion constant, when restricted to surfaces only, is lower than the effective diffusivity with purely reflecting boundaries. We give analytical formulas for this intriguing effect in periodic arrays of spheres as well as undulating microchannels. Our results are confirmed by numerical calculations and Monte Carlo simulations.
Collapse
Affiliation(s)
- A Alexandre
- Laboratoire Ondes et matière d'Aquitaine, CNRS/University of Bordeaux, F-33400 Talence, France
| | - M Mangeat
- Center for Biophysics and Department for Theoretical Physics, Saarland University, D-66123 Saarbrücken, Germany
| | - T Guérin
- Laboratoire Ondes et matière d'Aquitaine, CNRS/University of Bordeaux, F-33400 Talence, France
| | - D S Dean
- Laboratoire Ondes et matière d'Aquitaine, CNRS/University of Bordeaux, F-33400 Talence, France
- Team MONC, INRIA Bordeaux Sud Ouest, CNRS UMR 5251, Bordeaux INP, University Bordeaux, F-33400 Talence, France
| |
Collapse
|
3
|
Zaafouri Z, Batôt G, Nieto-Draghi C, Rotenberg B, Bauer D, Coasne B. Lattice Boltzmann method for adsorption under stationary and transient conditions: Interplay between transport and adsorption kinetics in porous media. Phys Rev E 2021; 104:015314. [PMID: 34412268 DOI: 10.1103/physreve.104.015314] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/06/2021] [Accepted: 07/01/2021] [Indexed: 11/07/2022]
Abstract
A numerical method based on the Lattice Boltzmann formalism is presented to capture the effect of adsorption kinetics on transport in porous media. Through the use of a general adsorption operator, canonical models such as Henry and Langmuir adsorption as well as more complex adsorption mechanisms involving collective behavior with lateral interactions and surface aggregation can be investigated using this versatile model. By extending the description of adsorption phenomena to kinetic regimes with any underlying adsorption model, this effective technique allows assessing the coupled dynamics resulting from advection, diffusion, and adsorption in pores not only in stationary conditions but also under transient conditions (i.e., in regimes where the adsorbed amount evolves with time due to diffusion and advection). As illustrated in this paper, the development of such an approach provides a simple tool to determine the reciprocal effect of molecular flow and dispersion on adsorption kinetics. In this context, the use of a Lattice Boltzmann-based approach is important as it allows considering porous media of any morphology and topology. Beyond fundamental implications, this efficient method allows treating real engineering conditions such as pollutant dispersion or surfactant injection in a flowing liquid in soils and porous rocks.
Collapse
Affiliation(s)
- Zaineb Zaafouri
- IFP Energies Nouvelles, 1 & 4 Av. Bois Préau, 92852 Rueil Malmaison, France and Université Grenoble Alpes, CNRS, LIPhy, 38000 Grenoble, France
| | - Guillaume Batôt
- IFP Energies Nouvelles, 1 & 4 Av. Bois Préau, 92852 Rueil Malmaison, France
| | | | | | - Daniela Bauer
- IFP Energies Nouvelles, 1 & 4 Av. Bois Préau, 92852 Rueil Malmaison, France
| | - Benoit Coasne
- Université Grenoble Alpes, CNRS, LIPhy, 38000 Grenoble, France
| |
Collapse
|
4
|
Sherman T, Paster A, Porta G, Bolster D. A spatial Markov model for upscaling transport of adsorbing-desorbing solutes. JOURNAL OF CONTAMINANT HYDROLOGY 2019; 222:31-40. [PMID: 30803748 DOI: 10.1016/j.jconhyd.2019.02.003] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/01/2018] [Revised: 01/31/2019] [Accepted: 02/04/2019] [Indexed: 06/09/2023]
Abstract
The Spatial Markov Model (SMM) is an upscaled model with a strong track record in predicting upscaled behavior of conservative solute transport across hydrologic systems. Here we propose an SMM that can account for reactive linear adsorption and desorption processes and test it on a simple benchmark problem: flow and transport through an idealized periodic wavy channel. The methodology is built using trajectories that are obtained from a single high resolution random walk simulation of conservative transport across one periodic element. Our approach encodes information about where a particle starts at the inlet, where it leaves at the outlet, how long it takes to cross the domain and one additional piece of information, the number of times a particle strikes the boundary, with the objective of predicting large scale transport with arbitrary linear adsorption and desorption rates. Our benchmark problem demonstrates that predictions made with our proposed SMM agree favorably with results from direct numerical simulations, which resolve the full transport problem.
Collapse
Affiliation(s)
- Thomas Sherman
- Dept. of Civil and Environmental Engineering and Earth Sciences, University of Notre Dame, IN, USA
| | | | - Giovanni Porta
- Dipartimento di Ingegneria Civile e Ambientale, Politecnico di Milano, Piazza L. Da Vinci, 32, 20133 Milano, Italy
| | - Diogo Bolster
- Dept. of Civil and Environmental Engineering and Earth Sciences, University of Notre Dame, IN, USA.
| |
Collapse
|
5
|
Asta A, Levesque M, Rotenberg B. Moment propagation method for the dynamics of charged adsorbing/desorbing species at solid-liquid interfaces. Mol Phys 2018. [DOI: 10.1080/00268976.2018.1461944] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
Affiliation(s)
- Adelchi Asta
- Sorbonne Université, CNRS, Physicochimie des électrolytes et nanosystèmes interfaciaux, UMR PHENIX , Paris, France
| | - Maximilien Levesque
- PASTEUR, Département de chimie, École normale supérieure, PSL University, Sorbonne Université, CNRS , 75005 Paris, France
| | - Benjamin Rotenberg
- Sorbonne Université, CNRS, Physicochimie des électrolytes et nanosystèmes interfaciaux, UMR PHENIX , Paris, France
| |
Collapse
|
6
|
Vanson JM, Boutin A, Klotz M, Coudert FX. Transport and adsorption under liquid flow: the role of pore geometry. SOFT MATTER 2017; 13:875-885. [PMID: 28074205 DOI: 10.1039/c6sm02414a] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
We study here the interplay between transport and adsorption in porous systems with complex geometries under fluid flow. Using a lattice Boltzmann scheme extended to take into account the adsorption at solid/fluid interfaces, we investigate the influence of pore geometry and internal surface roughness on the efficiency of fluid flow and the adsorption of molecular species inside the pore space. We show how the occurrence of roughness on pore walls acts effectively as a modification of the solid/fluid boundary conditions, introducing slippage at the interface. We then compare three common pore geometries, namely honeycomb pores, inverse opal, and materials produced by spinodal decomposition. Finally, we quantify the influence of those three geometries on fluid transport and tracer adsorption. This opens perspectives for the optimization of materials' geometries for applications in dynamic adsorption under fluid flow.
Collapse
Affiliation(s)
- Jean-Mathieu Vanson
- École Normale Supérieure, PSL Research University, Département de Chimie, Sorbonne Universités-UPMC Univ Paris 06, CNRS UMR 8640 PASTEUR, 24 rue Lhomond, 75005 Paris, France. and Laboratoire de Synthèse et Fonctionnalisation des Céramiques, UMR 3080 Saint Gobain CREE/CNRS, 550 Avenue Alphonse Jauffret, 84306 Cavaillon, France
| | - Anne Boutin
- École Normale Supérieure, PSL Research University, Département de Chimie, Sorbonne Universités-UPMC Univ Paris 06, CNRS UMR 8640 PASTEUR, 24 rue Lhomond, 75005 Paris, France.
| | - Michaela Klotz
- Laboratoire de Synthèse et Fonctionnalisation des Céramiques, UMR 3080 Saint Gobain CREE/CNRS, 550 Avenue Alphonse Jauffret, 84306 Cavaillon, France
| | - François-Xavier Coudert
- Chimie ParisTech, PSL Research University, CNRS, Institut de Recherche de Chimie Paris, 75005 Paris, France.
| |
Collapse
|
7
|
Chahibi Y, Pierobon M, Akyildiz IF. Pharmacokinetic Modeling and Biodistribution Estimation Through the Molecular Communication Paradigm. IEEE Trans Biomed Eng 2015; 62:2410-20. [DOI: 10.1109/tbme.2015.2430011] [Citation(s) in RCA: 38] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
|
8
|
Vanson JM, Coudert FX, Rotenberg B, Levesque M, Tardivat C, Klotz M, Boutin A. Unexpected coupling between flow and adsorption in porous media. SOFT MATTER 2015; 11:6125-33. [PMID: 26139013 DOI: 10.1039/c5sm01348h] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/22/2023]
Abstract
We study the interplay between transport and adsorption in porous systems under a fluid flow, based on a lattice Boltzmann scheme extended to account for adsorption. We performed simulations on well-controlled geometries with slit and grooved pores, investigating the influence of adsorption and flow on dispersion coefficient and adsorbed density. In particular, we present a counterintuitive effect where fluid flow induces heterogeneity in the adsorbate, displacing the adsorption equilibrium towards downstream adsorption sites in grooves. We also present an improvement of the adsorption-extended lattice Boltzmann scheme by introducing the possibility for saturating Langmuir-like adsorption, while earlier work focused on linear adsorption phenomena. We then highlight the impact of this change in situations of high concentration of adsorbate.
Collapse
Affiliation(s)
- Jean-Mathieu Vanson
- École Normale Supérieure, PSL Research University, Département de Chimie, Sorbonne Universités - UPMC Univ Paris 06, CNRS UMR 8640 PASTEUR, 24 rue Lhomond, 75005 Paris, France.
| | | | | | | | | | | | | |
Collapse
|
9
|
Berezhkovskii AM, Dagdug L, Bezrukov SM. Biased diffusion in three-dimensional comb-like structures. J Chem Phys 2015; 142:134101. [PMID: 25854222 PMCID: PMC4385101 DOI: 10.1063/1.4916310] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2014] [Accepted: 03/12/2015] [Indexed: 11/14/2022] Open
Abstract
In this paper, we study biased diffusion of point Brownian particles in a three-dimensional comb-like structure formed by a main cylindrical tube with identical periodic cylindrical dead ends. It is assumed that the dead ends are thin cylinders whose radius is much smaller than both the radius of the main tube and the distance between neighboring dead ends. It is also assumed that in the main tube, the particle, in addition to its regular diffusion, moves with a uniform constant drift velocity. For such a system, we develop a formalism that allows us to derive analytical expressions for the Laplace transforms of the first two moments of the particle displacement along the main tube axis. Inverting these Laplace transforms numerically, one can find the time dependences of the two moments for arbitrary values of both the drift velocity and the dead-end length, including the limiting case of infinitely long dead ends, where the unbiased diffusion becomes anomalous at sufficiently long times. The expressions for the Laplace transforms are used to find the effective drift velocity and diffusivity of the particle as functions of its drift velocity in the main tube and the tube geometric parameters. As might be expected from common-sense arguments, the effective drift velocity monotonically decreases from the initial drift velocity to zero as the dead-end length increases from zero to infinity. The effective diffusivity is a more complex, non-monotonic function of the dead-end length. As this length increases from zero to infinity, the effective diffusivity first decreases, reaches a minimum, and then increases approaching a plateau value which is proportional to the square of the particle drift velocity in the main tube.
Collapse
Affiliation(s)
- Alexander M Berezhkovskii
- Program in Physical Biology, Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health, Bethesda, Maryland 20892, USA
| | - Leonardo Dagdug
- Program in Physical Biology, Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health, Bethesda, Maryland 20892, USA
| | - Sergey M Bezrukov
- Program in Physical Biology, Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health, Bethesda, Maryland 20892, USA
| |
Collapse
|
10
|
Dagdug L, Berezhkovskii AM, Skvortsov AT. Aris-Taylor dispersion in tubes with dead ends. J Chem Phys 2014; 141:024705. [PMID: 25028036 DOI: 10.1063/1.4885854] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
This paper deals with transport of point Brownian particles in a cylindrical tube with dead ends in the presence of laminar flow of viscous fluid in the cylindrical part of the tube (Poiseuille flow). It is assumed that the dead ends are identical and are formed by spherical cavities connected to the cylindrical part of the tube by narrow necks. The focus is on the effective velocity and diffusivity of the particles as functions of the mean flow velocity and geometric parameter of the tube. Entering a dead end, the particle interrupts its propagation along the tube axis. Later it returns, and the axial motion continues. From the axial propagation point of view, the particle entry into a dead end and its successive return to the flow is equivalent to the particle reversible binding to the tube wall. The effect of reversible binding on the transport parameters has been previously studied assuming that the particle survival probability in the bound state decays as a single exponential. However, this is not the case when the particle enters a dead end, since escape from the dead end is a non-Markovian process. Our analysis of the problem consists of two steps: First, we derive expressions for the effective transport parameters in the general case of non-Markovian binding. Second, we find the effective velocity and diffusivity by substituting into these expressions known results for the moments of the particle lifetime in the dead end [L. Dagdug, A. M. Berezhkovskii, Yu. A. Makhnovskii, and V. Yu. Zitserman, J. Chem. Phys. 127, 224712 (2007)]. To check the accuracy of our theory, we compare its predictions with the values of the effective velocity and diffusivity obtained from Brownian dynamics simulations. The comparison shows excellent agreement between the theoretical predictions and numerical results.
Collapse
Affiliation(s)
- Leonardo Dagdug
- Physics Department, Universidad Autonoma Metropolitana-Iztapalapa, 09340 Mexico DF, Mexico
| | - Alexander M Berezhkovskii
- Mathematical and Statistical Computing Laboratory, Division of Computational Bioscience, Center for Information Technology, National Institutes of Health, Bethesda, Maryland 20892, USA
| | - Alexei T Skvortsov
- Maritime Division, Defence Science and Technology Organisation, Fishermans Bend, Vic 3207, Australia
| |
Collapse
|
11
|
Hlushkou D, Gritti F, Guiochon G, Seidel-Morgenstern A, Tallarek U. Effect of Adsorption on Solute Dispersion: A Microscopic Stochastic Approach. Anal Chem 2014; 86:4463-70. [DOI: 10.1021/ac500309p] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Affiliation(s)
- Dzmitry Hlushkou
- Department
of Chemistry, Philipps-Universität Marburg, Hans-Meerwein-Strasse, 35032 Marburg, Germany
- Max-Planck-Institut für Dynamik komplexer technischer Systeme, Sandtorstrasse 1, 39106 Magdeburg, Germany
| | - Fabrice Gritti
- Department
of Chemistry, The University of Tennessee, Knoxville, Tennessee 37996−1600, United States
| | - Georges Guiochon
- Department
of Chemistry, The University of Tennessee, Knoxville, Tennessee 37996−1600, United States
| | - Andreas Seidel-Morgenstern
- Max-Planck-Institut für Dynamik komplexer technischer Systeme, Sandtorstrasse 1, 39106 Magdeburg, Germany
| | - Ulrich Tallarek
- Department
of Chemistry, Philipps-Universität Marburg, Hans-Meerwein-Strasse, 35032 Marburg, Germany
| |
Collapse
|
12
|
Berezhkovskii AM, Skvortsov AT. Aris-Taylor dispersion with drift and diffusion of particles on the tube wall. J Chem Phys 2014; 139:084101. [PMID: 24006968 DOI: 10.1063/1.4818733] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
A laminar stationary flow of viscous fluid in a cylindrical tube enhances the rate of diffusion of Brownian particles along the tube axis. This so-called Aris-Taylor dispersion is due to the fact that cumulative times, spent by a diffusing particle in layers of the fluid moving with different velocities, are random variables which depend on the realization of the particle stochastic trajectory in the radial direction. Conceptually similar increase of the diffusivity occurs when the particle randomly jumps between two states with different drift velocities. Here we develop a theory that contains both phenomena as special limiting cases. It is assumed (i) that the particle in the flow can reversibly bind to the tube wall, where it moves with a given drift velocity and diffusivity, and (ii) that the radial and longitudinal diffusivities of the particle in the flow may be different. We derive analytical expressions for the effective drift velocity and diffusivity of the particle, which show how these quantities depend on the geometric and kinetic parameters of the model.
Collapse
Affiliation(s)
- Alexander M Berezhkovskii
- Mathematical and Statistical Computing Laboratory, Division of Computational Bioscience, Center for Information Technology, National Institutes of Health, Bethesda, Maryland 20892, USA
| | | |
Collapse
|
13
|
Levesque M, Duvail M, Pagonabarraga I, Frenkel D, Rotenberg B. Accounting for adsorption and desorption in lattice Boltzmann simulations. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2013; 88:013308. [PMID: 23944584 DOI: 10.1103/physreve.88.013308] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/29/2013] [Indexed: 06/02/2023]
Abstract
We report a Lattice-Boltzmann scheme that accounts for adsorption and desorption in the calculation of mesoscale dynamical properties of tracers in media of arbitrary complexity. Lattice Boltzmann simulations made it possible to solve numerically the coupled Navier-Stokes equations of fluid dynamics and Nernst-Planck equations of electrokinetics in complex, heterogeneous media. With the moment propagation scheme, it became possible to extract the effective diffusion and dispersion coefficients of tracers, or solutes, of any charge, e.g., in porous media. Nevertheless, the dynamical properties of tracers depend on the tracer-surface affinity, which is not purely electrostatic and also includes a species-specific contribution. In order to capture this important feature, we introduce specific adsorption and desorption processes in a lattice Boltzmann scheme through a modified moment propagation algorithm, in which tracers may adsorb and desorb from surfaces through kinetic reaction rates. The method is validated on exact results for pure diffusion and diffusion-advection in Poiseuille flows in a simple geometry. We finally illustrate the importance of taking such processes into account in the time-dependent diffusion coefficient in a more complex porous medium.
Collapse
|
14
|
Levesque M, Bénichou O, Rotenberg B. Molecular diffusion between walls with adsorption and desorption. J Chem Phys 2013; 138:034107. [DOI: 10.1063/1.4775742] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|