Numerical simulation of two-dimensional Rayleigh-Benard convection in an enclosure

In this Note, a numerical approach based on the finite volume method and a full multigrid acceleration is used, applied to the classical Rayleigh Benard convection problem. Fine grids corresponding to 256(2) nodes are used and Benchmark solutions are proposed for Rayleigh numbers ranging from 10(3) to 10(6). Some streamlines and isotherms are presented to analyze the natural convection flow patterns set up by the buoyancy force.