Bi-dimensional plume generated by the convective dissolution of an extended source of CO 2

Niloy De (1) , Patrice Meunier (2) , Yves Méheust (3) , François Nadal (1)

1 Loughborough University
2 IRPHE - Institut de Recherche sur les Phénomènes Hors Equilibre
3 GR - Géosciences Rennes

The dynamics of a bi-dimensional solutal plume generated by the convective dissolution of an extended source of carbon dioxide (CO2) is studied from both experimental and numerical standpoints. In the experiments, the CO2 is injected into a 1-mm-thick water filled Hele-Shaw cell at pressures ranging between 1 and 5 bars. The plume is visualized using a pH-sensitive dye, and the velocity field is measured by means of standard particle image velocimetry (PIV). Concentration and velocity fields are scrutinized over one order of magnitude in the Rayleigh number (Ra), and for two different values of the Darcy number (Da). In order to extend the explored range of these dimensionless numbers, experiments are complemented by a comprehenvive set of numerical simulations. At low Darcy numbers (Da < 1), the morphology of the numerical plume is fairly close to that of the analytical solution previously derived by Wooding [J. Fluid Mech. 15, 527 (1963)] for the case of a line source in a porous medium. At larger Darcy numbers (Da > 1), the numerical solution departs on several aspects from this solution and exhibits closer similarity to the analytical solution derived by Spalding and Cruddace [Int. J. Heat Mass Transfer. 3, 55 (1961)] for the case of a line source in a viscous fluid. As the introduction of an additional length scale (the lateral size of the source) in the problem allows for the definition of a length-based Rayleigh number Ra, the respective amplitudes of the plume velocity w, plume width sigma, and dimensionless total injection rate (i.e., the Nusselt number Nu) can be explicitly obtained as functions of Ra and Da. The scaling laws are reasonably close to those obtained from simple mass conservation considerations: (i) w similar to Ra, sigma similar to Ra-1/2, and Nu similar to Ra-1/2 for Ra Da < 1, and (ii) w similar to (Ra/Da)(1/2), sigma similar to (Ra/Da)(-1/4), and Nu similar to (Ra/Da)(1/4) for Ra Da > 1. Although the concentration field measurements are marred by a large uncertainty due to the logarithmic nature of the pH-concentration relationship, the experimental data obtained from the PIV measurements show good agreement with the numerical results.

Keywords: AXISYMMETRICAL CREEPING PLUME; SCARBON-DIOXIDE; PLANAR BOUNDARY; LARGE RAYLEIGH; STABILITY; ONSET; MEDIA

Domains: Fluid mechanics [physics.class-ph]; Hydrology

Niloy De, Patrice Meunier, Yves Méheust, François Nadal. Bi-dimensional plume generated by the convective dissolution of an extended source of CO 2. Physical Review Fluids, 2021, 6 (6), pp.063503. ⟨10.1103/PhysRevFluids.6.063503⟩. ⟨insu-03281162⟩