image thumbnail

Rayleigh Benard Convection

version (3.87 KB) by Suraj Shankar
Thermal gradient driven natural convection is simulated in a 2D domain


Updated 10 Sep 2012

View License

Natural convection due to thermal gradients are simulated in a 2D rectangular domain. The Navier-Stokes equations are solved by the pressure projection method on a staggered grid. The hyperbolic flux terms are discretized explicitly (CD, MacCormack and Richtmyer) while the diffusive terms are dealt with both explicitly and implicitly. The energy transport equation is explicitly discretized (CD) for the advective fluxes along with an option for implicit or explicit method for the conduction terms. The Pressure Poisson equation is solved implicitly. The top and bottom surfaces are isothermal while the sides are adiabatic. No-slip is enforced for velocity on all sides and homogenous boundary conditions are employed for pressure. Both velocity and temperature field are used for visualization purposes. For specific (critical) values of Pe,Gr and Re, Rayleigh Benard convective rolls are observed.

Cite As

Suraj Shankar (2022). Rayleigh Benard Convection (, MATLAB Central File Exchange. Retrieved .

MATLAB Release Compatibility
Created with R2010a
Compatible with any release
Platform Compatibility
Windows macOS Linux

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!