We study penetrative convection of a fluid confined between two horizontal plates, the temperatures of which are such that a temperature of maximum density lies between them. The range of Rayleigh numbers studied is $Ra = \left[10^6, 10^8 \right]$ and the Prandtl numbers are $Pr = 1$ and $11.6$. An evolution equation for the growth of the convecting region is obtained through an integral energy balance. We identify a new non-dimensional parameter, $��$, which is the ratio of temperature difference between the stable and unstable regions of the flow; larger values of $��$ denote increased stability of the upper stable layer. We study the effects of $��$ on the flow field using well-resolved lattice Boltzmann simulations, and show that the characteristics of the flow depend sensitively upon it. For the range $��= \left[0.01, 4\right]$, we find that for a fixed $Ra$ the Nusselt number, $Nu$, increases with decreasing $��$. We also investigate the effects of $��$ on the vertical variation of convective heat flux and the Brunt-V��is��l�� frequency. Our results clearly indicate that in the limit $��\rightarrow 0$ the problem reduces to that of the classical Rayleigh-B��nard convection.

12 pages, 19 figures