Solutal convection induced by dissolution

Abstract

The dissolution of minerals into water becomes significant in geomorphology, when the erosion rate is controlled by the hydrodynamics transport of the solute. Even in absence of an external flow, dissolution itself can induce a convection flow due to the action of gravity. Here we perform a study of the physics of solutal convection induced by dissolution. We simulate numerically the hydrodynamics and the solute transport, in a 2D geometry, corresponding to the case, where a soluble body is suddenly immersed in a quiescent solvent. We identify three regimes. At short timescale, a concentrated boundary layer grows by diffusion at the interface. After a finite onset time, the thickness and the density reach critical values which starts the destabilization of the boundary layer. Finally, the destabilization is such as we observe the emission of intermittent plumes. This last regime is quasi-stationnary: the structure of the boundary layer as well as the erosion rate fluctuate around constant values. Assuming that the destabilization of the boundary layer occurs at a specific value of the solutal Rayleigh number, we derive scaling laws both for fast and slow dissolution kinetics. Our simulations confirm this scenario by validating the scaling laws both for onset, and the quasi-stationary regime. We find a constant value of the Rayleigh number during the quasi-stationary regime showing that the structure of the boundary layer is well controlled by the solutal convection. Finally, we apply the scaling laws previously established to the case of real dissolving minerals. We predicts the typical dissolution rate in presence of solutal convection. Our results suggest that solutal convection could occur in more natural situations than expected. Even for minerals with a quite low saturation concentration, the erosion rate would increase as the dissolution would be controlled by the hydrodynamics.