Abstract
The modal expansion procedure has been used to analyze penetrative convection that arises when a thin unstable layer is embedded between two stable regions. The Boussinesq approximation is applied in which the effect of compressibility and stratification are neglected. Various calculations have been made, with one and two modes, for Rayleigh numbers ranging from the critical value to more than 105 times critical. The effect of decreasing the Prandtl number has also been investigated.
It is found that in the nonlinear regime, the convective motions penetrate substantially into the stable regions. The flux of kinetic energy plays a crucial role in such penetration, and its existence puts some requirements on the motions: in the single-mode case, they need to be three-dimensional. The extent of penetration amounts to about half of the thickness of the unstable layer on each side of it when the degree of instability and that of stability are comparable in the two domains; it increases as the stability of the outer region is lowered. The penetration depth appears to be independent of all other parameters defining the problem.