ABSTRACT: Utilizing porous media in a new mathematical model to improve convective heat transfer characteristics in a variety of applications, such as radiation nuclear disposal storing, evaporation cooling, sieving, geological extraction, crude petroleum refining, and building heating and cooling, is becoming increasingly important. This study proposed a numerical analysis of the unsteady magnetohydrodynamic free convection flow of an exothermic fluid with Newtonian heating. This discovery reveals two types of solutions: steady state and unsteady state. After transforming the governing equation from dimensional form to dimensionless form, the steady state governing equation was solved by the Homotopy Perturbation Method. However, the implicit finite difference approach is used to solve the time-dependent governing equations numerically. The impact of various emerging parameters, namely the Hartmann number, Boit number, Darcy number, Navier slip parameter, and the Frank-Kamenetskii parameter, was discussed and graphically analyzed. During the computations and analysis, it was discovered that a minor rise in the Hartman number results in the Lorentz force, which streamlines the momentum barrier layer and hence slows the fluid flow. The fluid velocity, on the other hand, rose as the porous medium, thermal Biot number, slip parameter, and temperature field increased as the viscous reactive fluid parameter and Newtonian heating increased. The skin friction and Nusselt number were also examined and reported. By comparing the finding to an existing work, a great agreement was revealed.
