Numerical study of magneto-convection inside an enclosure using enhanced stream function-vorticity formulation

The buoyancy-driven magneto-hydrodynamic ow is investigated using a two dimensional numerical simulation in a filled square enclosure, which is heated by a uniform volumetric heat density. the top and bottom walls of the enclosure are adiabatic and the side walls have constant temperature. a fixed magnetic field is applied in the direction normal to the left side wall of the enclosure. the stationary dimensionless governing equations are solved numerically for the stream function, vorticity and temperature, with finite difference method, using matlab software. the prandtl (pr) number of the uid is assumed to be 0.733; the rayleigh (ra) number is made to vary from 10^4 to 10^7 and the hartmann (ha) number is between 0 and 10^3. the stream function equation is solved using fast poisson's equation solver on a rectangular grid (poicalc function in matlab) vorticity, and temperature equations are solved using the red-black gauss-seidel and biconjugate gradient stabilized (bicgstab) methods, respectively. the proposed method is fast, and there is no need for the variables under-relaxation. it is interesting to know that convergence of this method is superior to other segregated algorithms. the ratio of the lorentz force to the buoyancy force (ha2/ra) is introduced as an index to compare the contribution of natural convection with the magnetic field strength on heat transfer.