Numerical solutions of stability and convective flow in an infinite ho
rizontal water layer, including density inversion, have been obtained
using a finite element code. The evolution of the temperature field an
d flow pattern near the onset of convection are studied in detail. It
is known that natural convection develops primarily in the lower unsta
bly stratified layer. Of interest is the penetration of the convection
rolls into the upper stably stratified layer and concurrent liquid en
trainment as a function of the increasing Rayleigh number at different
aspect ratios. Individual convection rolls may grow and expand before
splitting up into two roll cells. It is shown that changing the aspec
t ratio influences critical Rayleigh number, flow symmetry, flow patte
rn, and transitions between flow patterns. Numerical results on heatin
g from above or from below, agree well with available results in the l
iterature. A correlation to predict critical Rayleigh numbers is given
for the case of heating from above.