Magnetic Field Effect on Natural Convection Flow with Internal Heat Generation using Fast  –  Method

The magnetic field effect on laminar natural convection flow is investigated in a filled enclosure with internal heat generation using two-dimensional numerical simulation. The enclosure is heated by a uniform volumetric heat density and walls have constant temperature. A fixed magnetic field is applied to the enclosure. The dimensionless governing equations are solved numerically for the stream function, vorticity and temperature using finite difference method for various Rayleigh (Ra) and Hartmann (Ha) numbers in MATLAB software. The stream function equation is solved using fast Poisson's equation solver on a rectangular grid (POICALC function in MATLAB), voricity and temperature equations are solved using red-black Gauss-Seidel and bi-conjugate gradient stabilized (BiCGSTAB) methods respectively. The results show that the strength of the magnetic field has significant effects on the flow and temperature fields. For the square cavity, the maximum temperature reduces with increasing Ra number. It is also observed that at low Ra number, location of the maximum temperature is at the centre of the cavity and it shifts upwards with increase in Ra number. Circulation inside the enclosure and therefore the convection becomes stronger as the Ra number increases while the magnetic field suppresses the convective flow and the heat transfer rate. The ratio of the Lorentz force to the buoyancy force (Ha2/Ra) is as an index to compare the contribution of natural convection and magnetic field strength on heat transfer.