Numerical simulation of thermal management during natural convection in a porous triangular cavity containing air and hot obstacles

• Published in: European physical journal plus 2021-08, Vol.136 (8), p.885, Article 885
• Main Authors: chandanam, Veena, lakshmi, C. Venkata, Venkatadri, K., Bég, O. Anwar, Prasad, V. Ramachandra
• Format: Article
• Language: English
• Subjects: Applied and Technical Physics; Atomic; Barriers; Boundary conditions; Complex Systems; Condensed Matter Physics; Convection; Differential equations; Entropy; Finite difference method; Finite volume method; Free convection; Fuel cells; Heat transfer; Investigations; Isotherms; Magnetic fields; Magnetic materials; Magnetohydrodynamics; Materials processing; Mathematical and Computational Physics; Mathematical models; Molecular; Nanoparticles; Optical and Plasma Physics; Parameters; Partial differential equations; Physical properties; Physics; Physics and Astronomy; Porous materials; Prandtl number; Rayleigh number; Regular Article; Theoretical; Thermal management; Thermal simulation; Vorticity
• ISSN: 2190-5444; 2190-5444
• DOI: 10.1140/epjp/s13360-021-01881-3

A numerical study is presented of laminar viscous magnetohydrodynamic natural convection flow in a triangular-shaped porous enclosure filled with electrically conducting air and containing two hot obstacles. The mathematical model is formulated in terms of dimensional partial differential equations. The pressure gradient term is eliminated by using the vorticity–stream ( ω - ψ ) function approach. The emerging dimensionless governing equations are employed by the regular finite difference scheme along with thermofluidic boundary conditions. The efficiency of the obtained computed results for isotherms and streamlines is validated via comparison with previously published work. The impact of physical parameters on streamlines and temperature contours for an extensive range of Rayleigh number (Ra = 10 3 –10 5 ), Hartmann magnetohydrodynamic number (Ha = 5–30), Darcy parameter (Da = 0.0001–0.1) for fixed Prandtl number (Pr = 0.71) is considered. Numerical results are also presented for local and average Nusselt numbers along the hot base wall. Interesting features of the thermofluid behaviour are revealed. At lower Rayleigh number, the isotherms are generally parallel to the inclined wall and only distorted substantially near the obstacles at the left vertical adiabatic wall; however, with increasing Rayleigh number , this distortion is magnified in the core zone and simultaneously warmer zones expand towards the inclined cold wall. The simulations are relevant to magnetic materials processing and hybrid magnetic fuel cells.