Eur. Phys. J. Plus (2021) 136: 885
https://doi.org/10.1140/epjp/s13360-021-01881-3

Regular Article

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

¹ Department of Applied Mathematics, Sri Padmavathi Mahila Visvavidyalayam, 517502, Tirupati, India
² Department of Mathematics, Sreenivasa Institute of Technology and Management Studies, 517127, Chittoor, India
³ Multi-Physical Engineering Sciences Group, Aeronautical/Mechanical Engineering, Salford University, School of Science, Engineering and Environment (SEE), M54WT, Manchester, UK
⁴ Department of Mathematics, School of Advanced Sciences, Vellore Institute of Technology, 6322014, Vellore, India

^c venkatadri.venki@gmail.com

Received: 2 May 2021
Accepted: 18 August 2021
Published online: 27 August 2021

Abstract

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 ($\omega -\psi$) 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³–10⁵), 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.