Convection in a compressible fluid with infinite Prandtl number

Abstract

An approximate set of equations is derived for a compressible liquid of infinite Prandtl number. These are referred to as the anelastic-liquid equations. The approximation requires the product of absolute temperature and volume coefficient of thermal expansion to be small compared to one. A single parameter defined as the ratio of the depth of the convecting layer,d, to the temperature scale height of the liquid,HT, governs the importance of the non-Boussinesq effects of compressibility, viscous dissipation, variable adiabatic temperature gradients and non-hydrostatic pressure gradients. Whend/HT[Lt ] 1 the Boussinesq equations result, but whend/HTisO(1) the non-Boussinesq terms become important. Using a time-dependent numerical model, the anelastic-liquid equations are solved in two dimensions and a systematic investigation of compressible convection is presented in whichd/HTis varied from 0·1 to 1·5. Both marginal stability and finite-amplitude convection are studied. Ford/HT[les ] 1·0 the effect of density variations is primarily geometric; descending parcels of liquid contract and ascending parcels expand, resulting in an increase in vorticity with depth. Whend/HT> 1·0 the density stratification significantly stabilizes the lower regions of the marginal state solutions. At all values ofd/HT[ges ] 0·25, an adiabatic temperature gradient proportional to temperature has a noticeable stabilizing effect on the lower regions. Ford/HT[ges ] 0·5, marginal solutions are completely stabilized at the bottom of the layer and penetrative convection occurs for a finite range of supercritical Rayleigh numbers. In the finite-amplitude solutions adiabatic heating and cooling produces an isentropic central region. Viscous dissipation acts to redistribute buoyancy sources and intense frictional heating influences flow solutions locally in a time-dependent manner. The ratio of the total viscous heating in the convecting system, ϕ, to the heat flux across the upper surface,Fu, has an upper limit equal tod/HT. This limit is achieved at high Rayleigh numbers, when heating is entirely from below, and, for sufficiently large values ofd/HT, Φ/Fuis greater than 1·00.