Picture water droplets

Developing mathematical theories of the physical world: Open Access research on fluid dynamics from Strathclyde

Strathprints makes available Open Access scholarly outputs by Strathclyde's Department of Mathematics & Statistics, where continuum mechanics and industrial mathematics is a specialism. Such research seeks to understand fluid dynamics, among many other related areas such as liquid crystals and droplet evaporation.

The Department of Mathematics & Statistics also demonstrates expertise in population modelling & epidemiology, stochastic analysis, applied analysis and scientific computing. Access world leading mathematical and statistical Open Access research!

Explore all Strathclyde Open Access research...

Natural convection and the evolution of a reactive porous medium

Ritchie, Lindsey Thomson and Pritchard, David (2011) Natural convection and the evolution of a reactive porous medium. Journal of Fluid Mechanics, 673. pp. 286-317. ISSN 0022-1120

Final Published Version

Download (7MB) | Preview


We describe a mathematical model of buoyancy-driven flow and solute transport in a saturated porous medium, the porosity and permeability of which evolve through precipitation and dissolution as a mineral is lost or gained from the pore fluid. Imposing a vertically varying equilibrium solubility creates a density gradient which can drive convective circulation. We characterise the onset of convection using linear stability analysis, and explore the further development of the coupled reaction–convection system numerically. At low Rayleigh numbers, the effect of the reaction–permeability feedback is shown to be destabilising through a novel reaction–diffusion mechanism; at higher Rayleigh numbers, the precipitation and dissolution have a stabilising effect. Over longer time scales, reaction–permeability feedback triggers secondary instabilities in quasi-steady convective circulation, leading to rapid reversals in the direction of circulation. Over very long time scales, characteristic patterns of porosity emerge, including horizontal layering as well as the development of vertical chimneys of enhanced porosity. We discuss the implications of these findings for more comprehensive models of reactive convection in porous media.