The impact of geochemistry on convective mixing in a gravitationally unstable diffusive boundary layer in porous media: CO2 storage in saline aquifers

2011 ◽  
Vol 673 ◽  
pp. 480-512 ◽  
Author(s):  
KARIM GHESMAT ◽  
HASSAN HASSANZADEH ◽  
JALAL ABEDI

The storage of carbon dioxide and acid gases in deep geological formations is considered a promising option for mitigation of greenhouse gas emissions. An understanding of the primary mechanisms such as convective mixing and geochemistry that affect the long-term geostorage process in deep saline aquifers is of prime importance. First, a linear stability analysis of an unstable diffusive boundary layer in porous media is presented, where the instability occurs due to a density difference between the carbon dioxide saturated brine and the resident brine. The impact of geochemical reactions on the stability of the boundary layer is examined. The equations are linearised, and the obtained system of eigenvalue problems is solved numerically. The linear stability results have revealed that geochemistry stabilises the boundary layer as reaction consumes the dissolved carbon dioxide and makes the density profile, as the source of instability, more uniform. A detailed physical discussion is also presented with an examination of vorticity and concentration eigenfunctions and streamlines' contours to reveal how the geochemical reaction may affect the hydrodynamics of the process. We also investigate the effects of the Rayleigh number and the diffusion time on the stability of a boundary layer coupled with geochemical reactions. Nonlinear direct numerical simulations are also presented, in which the evolution of density-driven instabilities for different reaction rates is discussed. The development of instability is precisely studied for various scenarios. The results indicate that the boundary layer will be more stable for systems with a higher rate of reaction. However, our quantitative analyses show that more carbon dioxide may be removed from the supercritical free phase as the measured flux at the boundary is always higher for flow systems coupled with stronger geochemical reactions.

2017 ◽  
Vol 19 (1) ◽  
pp. 644-655 ◽  
Author(s):  
Parama Ghoshal ◽  
Min Chan Kim ◽  
Silvana S. S. Cardoso

We quantify the destabilising effect of a first-order chemical reaction on the fingering instability of a diffusive boundary layer in a porous medium.


2015 ◽  
Vol 764 ◽  
pp. 296-315 ◽  
Author(s):  
X. Fu ◽  
L. Cueto-Felgueroso ◽  
D. Bolster ◽  
R. Juanes

AbstractMotivated by the process of $\text{CO}_{2}$ convective mixing in porous media, here we study the formation of rock-dissolution patterns that arise from geochemical reactions during Rayleigh–Bénard–Darcy convection. Under the assumption of instantaneous chemical equilibrium, we adopt a formulation of the local reaction rate as a function of scalar dissipation rate, a measure that depends solely on flow and transport, and chemical speciation, which is a measure that depends only on the equilibrium thermodynamics of the chemical system. We use high-resolution simulations to examine the interplay between the density-driven hydrodynamic instability and the rock dissolution reactions, and analyse the impact of geochemical reactions on the macroscopic mass exchange rate. We find that dissolution of carbonate rock initiates in regions of locally high mixing, but that the geochemical reaction shuts down significantly earlier than shutdown of convective mixing. This early shutdown feature reflects the important role that chemical speciation plays in this hydrodynamics–reaction coupled process. Finally, we extend our analysis to three dimensions and explore the morphology of dissolution patterns in three dimensions.


2019 ◽  
Vol 126 ◽  
pp. 40-54 ◽  
Author(s):  
Seyed Mostafa Jafari Raad ◽  
Hamid Emami-Meybodi ◽  
Hassan Hassanzadeh

1985 ◽  
Vol 160 ◽  
pp. 181-210 ◽  
Author(s):  
R. H. Nilson

Countercurrent flow may be induced by opposing buoyancy forces associated with compositional gradients and thermal gradients within a fluid. The occurrence and structure of such flows is investigated by solving the double-diffusive boundary-layer equations for steady laminar convection along a vertical wall of finite height. Non-similar solutions are derived using the method of matched asymptotic expansions, under the restriction that the Lewis and Prandtl numbers are both large. Two sets of asymptotic solutions are constructed, assuming dominance of one or the other of the buoyancy forces. The two sets overlap in the central region of the parameter space; each set matches up with neighbouring unidirectional similarity solutions at the respective borderlines of incipient counterflow.Interaction between the buoyancy mechanisms is controlled by their relative strength R and their relative diffusivity Le. Flow in the outer thermal boundary layer deviates from single-diffusive thermal convection, depending upon the magnitude of the parameter RLe. Flow in the inner compositional boundary layer deviates from single-diffusive compositional convection, depending upon the magnitude of $RLe^{\frac{1}{3}}$.


Sign in / Sign up

Export Citation Format

Share Document