1
|
Tanaka R, Almarcha C, Nagatsu Y, Méheust Y, Meunier P. Chemically Enhanced Convective Dissolution of CO_{2}. PHYSICAL REVIEW LETTERS 2024; 132:084002. [PMID: 38457725 DOI: 10.1103/physrevlett.132.084002] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/04/2023] [Accepted: 01/10/2024] [Indexed: 03/10/2024]
Abstract
Convective dissolution, one of the main mechanisms for geological storage of CO_{2}, occurs when supercritical or gas CO_{2} dissolves partially into an aqueous solution, thus triggering downward convection of the denser CO_{2}-enriched liquid. Chemical reaction in the liquid can greatly enhance the process. Here, experimental measurements of convective flow inside a cylinder filled with a sodium hydroxide (NaOH) solution show that the plume's velocity can be increased tenfold as compared to a situation with no NaOH. This tremendous effect is predicted by a model with no adjusting parameters.
Collapse
Affiliation(s)
- R Tanaka
- Department of Chemical Engineering, Tokyo University of Agriculture and Technology, 2-24-16 Naka-cho, Koganei, Tokyo 184-8588, Japan
| | - C Almarcha
- Aix Marseille Université, CNRS, Centrale Marseille, IRPHE UMR 7342, 13384 Marseille, France
| | - Y Nagatsu
- Department of Chemical Engineering, Tokyo University of Agriculture and Technology, 2-24-16 Naka-cho, Koganei, Tokyo 184-8588, Japan
| | - Y Méheust
- Université Rennes, CNRS, Géosciences Rennes, UMR 6118, 35000 Rennes, France
| | - P Meunier
- Aix Marseille Université, CNRS, Centrale Marseille, IRPHE UMR 7342, 13384 Marseille, France
| |
Collapse
|
2
|
Kabbadj S, Rongy L, De Wit A. Effect of variable solubility on reactive dissolution in partially miscible systems. Phys Rev E 2023; 107:065109. [PMID: 37464620 DOI: 10.1103/physreve.107.065109] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2023] [Accepted: 06/07/2023] [Indexed: 07/20/2023]
Abstract
When two partially miscible systems are put in contact, one phase, A, can dissolve into the other one with a given solubility. Chemical reactions in the host phase can impact this dissolution by consuming A and by generating products that impact the solubility of A. Here, we study theoretically the optimal conditions for transfer of a reactant A in a host phase containing a species B when a bimolecular A + B → C reaction generates a product C that linearly decreases the solubility of A. We have quantified numerically the influence of this variable solubility on the reaction-diffusion (RD) concentration profiles of all species in the host phase, on the temporal evolution of the position of the reaction front, and on the flux of A through the interface. We have also computed the analytical asymptotic concentration profiles, solutions at long times of the RD governing equations. For a fixed negative effect of C on the solubility of A, an increase in the initial concentration of reactant B or an increase in the diffusion rate of species B and C results in a larger flux of A and hence a larger amount of A dissolved in the host solution at a given time. However, when the influence of C on the solubility increases, the mass transfer decreases. Our results help understand to what extent a chemical reaction can optimize the reactive transfer of a solute to a host phase with application to, among other things, the geological sequestration of carbon dioxide in an aquifer.
Collapse
Affiliation(s)
- S Kabbadj
- Nonlinear Physical Chemistry Unit, Université libre de Bruxelles, CP231, Boulevard du Triomphe, 1050 Bruxelles, Belgium
| | - L Rongy
- Nonlinear Physical Chemistry Unit, Université libre de Bruxelles, CP231, Boulevard du Triomphe, 1050 Bruxelles, Belgium
| | - A De Wit
- Nonlinear Physical Chemistry Unit, Université libre de Bruxelles, CP231, Boulevard du Triomphe, 1050 Bruxelles, Belgium
| |
Collapse
|
3
|
Bratsun D, Mizev A, Utochkin V, Nekrasov S, Shmyrova A. Nonlinear development of convective patterns driven by a neutralization reaction in immiscible two-layer systems. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2023; 381:20220178. [PMID: 36842984 DOI: 10.1098/rsta.2022.0178] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/22/2022] [Accepted: 11/14/2022] [Indexed: 06/18/2023]
Abstract
This article provides the results of a theoretical and experimental study of buoyancy-driven instabilities triggered by a neutralization reaction in an immiscible two-layer system placed in a vertical Hele-Shaw cell. Flow patterns are predicted by a reaction-induced buoyancy number [Formula: see text], which we define as the ratio of densities of the reaction zone and the lower layer. In experiments, we observed the development of cellular convection ([Formula: see text]), the fingering process with an aligned line of fingertips at a slightly denser reaction zone ([Formula: see text]) and the typical Rayleigh-Taylor convection for [Formula: see text]. A mathematical model includes a set of reaction-diffusion-convection equations written in the Hele-Shaw approximation. The model's novelty is that it accounts for the water produced during the reaction, a commonly neglected effect. The persisting regularity of the fingering during the collapse of the reaction zone is explained by the dynamic release of water, which compensates for the heavy fluid falling and stabilizes the pattern. Finally, we present a stability map on the plane of the initial concentrations of solutions. Good agreement between the experimental data and theoretical results is observed. This article is part of the theme issue 'New trends in pattern formation and nonlinear dynamics of extended systems'.
Collapse
Affiliation(s)
- Dmitry Bratsun
- Applied Physics Department, Perm National Research Polytechnic University, Perm 614990, Russia
| | - Alexey Mizev
- Applied Physics Department, Perm National Research Polytechnic University, Perm 614990, Russia
- Hydrodynamic Stability Lab, Institute of Continuous Media Mechanics, Perm 614013, Russia
| | - Vladimir Utochkin
- Applied Physics Department, Perm National Research Polytechnic University, Perm 614990, Russia
| | - Svyatoslav Nekrasov
- Applied Physics Department, Perm National Research Polytechnic University, Perm 614990, Russia
| | - Anastasia Shmyrova
- Applied Physics Department, Perm National Research Polytechnic University, Perm 614990, Russia
- Hydrodynamic Stability Lab, Institute of Continuous Media Mechanics, Perm 614013, Russia
| |
Collapse
|
4
|
Jotkar M, Rongy L, De Wit A. Chemically-driven convective dissolution. Phys Chem Chem Phys 2019; 21:19054-19064. [PMID: 31468054 DOI: 10.1039/c9cp03044a] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
When a solute A dissolves in a host phase with a given solubility, the resulting density stratification is stable towards convection if the density profile increases monotonically along the gravity field. We theoretically and numerically study the convective destabilization by reaction of this dissolution when A reacts with a solute B present in the host phase to produce C via an A + B→C type of reaction. In this reactive case, composition changes can give rise to non-monotonic density profiles with a local maximum. A convective instability can then be triggered locally in the zone where the denser product overlies the less dense bulk solution. First, we perform a linear stability analysis to identify the critical conditions for this reaction-driven convective instability. Second, we perform nonlinear simulations and compare the critical values of the control parameters for the onset of convection in these simulations with those predicted by linear stability analysis. We further show that the asymptotic dissolution flux of A can be increased in the convective regime by increasing the difference ΔRCB = RC-RB between the Rayleigh numbers of the product C and reactant B above a critical value and by increasing the ratio β = B0/A0 between the initial concentration B0 of reactant B and the solubility A0 of A. Our results indicate that chemical reactions can not only initiate convective mixing but can also give rise to large dissolution fluxes, which is advantageous for various geological applications.
Collapse
Affiliation(s)
- M Jotkar
- Université libre de Bruxelles (ULB), Faculté des Sciences, Nonlinear Physical Chemistry Unit, C.P. 231, 1050 Brussels, Belgium.
| | - L Rongy
- Université libre de Bruxelles (ULB), Faculté des Sciences, Nonlinear Physical Chemistry Unit, C.P. 231, 1050 Brussels, Belgium.
| | - A De Wit
- Université libre de Bruxelles (ULB), Faculté des Sciences, Nonlinear Physical Chemistry Unit, C.P. 231, 1050 Brussels, Belgium.
| |
Collapse
|
5
|
Jotkar M, De Wit A, Rongy L. Enhanced convective dissolution due to an A + B → C reaction: control of the non-linear dynamics via solutal density contributions. Phys Chem Chem Phys 2019; 21:6432-6442. [PMID: 30839024 DOI: 10.1039/c8cp07642a] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
Chemical reactions can have a significant impact on convective dissolution in partially miscible stratifications in porous media and are able to enhance the asymptotic flux with respect to the non-reactive case. We numerically study such reactive convective dissolution when the dissolving species A increases the density of the host phase upon dissolution and reacts with a reactant B present in the host phase to produce C by an A + B → C reaction. Upon varying the difference ΔRCB = RC-RB between the Rayleigh numbers of the product C and the reactant B, we identify four regimes with distinct dynamics when the diffusion coefficients are the same. When ΔRCB < 0, the density profiles are non-monotonic and the non-linear dynamics are seen to depend on the relative values of the density at the interface and the initial density of the host phase. For ΔRCB > 0, the monotonic density profiles are destabilizing with respect to the non-reactive case above a certain critical value ΔRcr. We analyze quantitatively the influence of varying ΔRCB and the ratio β = B0/A0 of the initial concentration of B and the solubility of A on the asymptotic steady flux, the wavelength of the fingers and the position of the reaction front. In the context of CO2 geological sequestration, understanding how such reactions can enhance the dissolution flux is crucial for improving the efficiency and safety of the process.
Collapse
Affiliation(s)
- M Jotkar
- Université libre de Bruxelles (ULB), Faculté des Sciences, Nonlinear Physical Chemistry Unit, C.P. 231, 1050 Brussels, Belgium.
| | - A De Wit
- Université libre de Bruxelles (ULB), Faculté des Sciences, Nonlinear Physical Chemistry Unit, C.P. 231, 1050 Brussels, Belgium.
| | - L Rongy
- Université libre de Bruxelles (ULB), Faculté des Sciences, Nonlinear Physical Chemistry Unit, C.P. 231, 1050 Brussels, Belgium.
| |
Collapse
|
6
|
Abstract
When a solute A dissolves into a host fluid containing a reactant B, an A + B → C reaction can influence the convection developing because of unstable density gradients in the gravity field. When A increases density and all three chemical species A, B and C diffuse at the same rate, the reactive case can lead to two different types of density profiles, i.e., a monotonically decreasing one from the interface to the bulk and a non-monotonic profile with a minimum. We study numerically here the nonlinear reactive convective dissolution dynamics in the more general case where the three solutes can diffuse at different rates. We show that differential diffusion can add new dynamic effects like the simultaneous presence of two different convection zones in the host phase when a non-monotonic profile with both a minimum and a maximum develops. Double diffusive instabilities can moreover affect the morphology of the convective fingers. Analysis of the mixing zone, the reaction rate, the total amount of stored A and the dissolution flux further shows that varying the diffusion coefficients of the various species has a quantitative effect on convection.
Collapse
|
7
|
Ghoshal P, Cardoso SSS. Reactive convective-dissolution in a porous medium: stability and nonlinear dynamics. Phys Chem Chem Phys 2018; 20:21617-21628. [PMID: 30101260 DOI: 10.1039/c8cp03064b] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We investigate the effects of a dissolution reaction, A(aq) + B(s) → C(aq), on the gravitational instability and nonlinear dynamic behaviour of a diffusive boundary layer in a porous medium. Our linear stability and numerical results reveal that, unexpectedly, even when the density contribution of the soluble product C is smaller than that of the dissolved solute A, the chemical reaction can destabilize the layer and accelerate the onset of convection. However, for a very light product, the reaction stabilizes the layer. We show that these widely disparate characteristics of the reactive-diffusive layer are outcomes of the nonlinear competition between two reaction effects, the destabilizing sharpening of the solute concentration gradient and associated increase in the solute diffusive flux, and the stabilizing replacement of the solute by a less dense product near the interface.
Collapse
Affiliation(s)
- Parama Ghoshal
- Department of Chemical Engineering and Biotechnology, University of Cambridge, Cambridge CB2 3RA, UK.
| | | |
Collapse
|