In partially miscible two-layer systems within a gravity field, buoyancy-driven convective motions can appear when one phase dissolves with a finite solubility into the other one. We investigate the influence of chemical reactions on such convective dissolution by a linear stability analysis of a reaction-diffusion-convection model. We show theoretically that a chemical reaction can either enhance or decrease the onset time of the convection, depending on the type of density profile building up in time in the reactive solution. We classify the stabilizing and destabilizing scenarios in a parameter space spanned by the solutal Rayleigh numbers. As an example, we experimentally demonstrate the possibility to enhance the convective dissolution of gaseous ${\mathrm{CO}}_2$ in aqueous solutions by a classical acid-base reaction.

The dynamics of $A+B\to C$ fronts in horizontal solution layers can be influenced by buoyancy-driven convection as soon as the densities of $A$, $B$, and $C$ are not all identical. Such convective motions can lead to front propagation even in the case of equal diffusion coefficients and initial concentration of reactants for which reaction-diffusion (RD) scalings predict a nonmoving front. We show theoretically that the dynamics in the presence of convection can in that case be predicted solely on the basis of the knowledge of the one-dimensional RD density profile across the front.

When autocatalytic chemical fronts propagate in thin layers of solution in contact with air, they can induce capillary flows due to surface tension gradients across the front (Marangoni flows). We investigate here such an interplay between autocatalytic reactions, diffusion, and Marangoni effects with a theoretical model coupling the incompressible Navier-Stokes equations to a conservation equation for the autocatalytic product concentration in the absence of gravity and for isothermal conditions. The boundary condition at the open liquid/air interface takes the surface activity of this product into account and introduces the solutal Marangoni number M representing the intensity of the coupling between hydrodynamics and reaction-diffusion processes. Positive and negative Marangoni numbers correspond, respectively, to the cases where the product decreases or increases surface tension behind the front. We show that, in both cases, such coupled systems reach an asymptotic dynamics characterized by a steady fluid vortex traveling at a constant speed with the front and deforming it, with, however, an asymmetry between the results for positive and negative M. A parametric study shows that increased propagation speed, front deformation, and possible transient oscillating dynamics occur when the absolute value of M is increased.

The convective dissolution of carbon dioxide (CO₂) in salted water is theoretically studied to determine how parameters such as CO₂ pressure, salt concentration, and temperature impact the short-time characteristics of the buoyancy-driven instability. On the basis of a parameter-free dimensionless model, we perform a linear stability analysis of the time-dependent concentration profiles of CO₂ diffusing into the aqueous solution. We explicit the procedure to transform the predicted dimensionless growth rate and wavelength of the convective pattern into dimensional ones for typical laboratory-scale experiments in conditions close to room temperature and atmospheric pressure. This allows to investigate the implicit influence of the experimental parameters on the characteristic length and time scales of the instability. We predict that increasing CO₂ pressure, or decreasing salt concentration or temperature destabilizes the system with regard to convection, leading to a faster dissolution of CO₂ into salted water.