Sensitivity Analysis of Pulsatile Hydromagnetic Biofluid Flow and Heat Transfer with Non Linear Darcy-Forchheimer Drag


1 Department of General Studies (Mathematics group), Jubail University College, Jubail, Saudi Arabia

2 Department of Mathematics, Indian Institute of Technology, Roorkee - 247667, India

3 Department of Education in Science and Mathematics, Regional Institute of Education, Bhubneswar (NCERT), India

4 Department of Mathematics, DAV University College, Jalandhar, India


In the present paper we examine the pulsatile hydromagnetic flow and heat transfer of a non-Newtonian biofluid through a saturated non-Darcian porous medium channel. The upper plate of the channel is heated and the lower plate is cooled. The Nakamura-Sawada rheological model is employed which provides a higher yield stress than the Casson model. A Darcy-Forchheimer porous medium drag force model is incorporated to simulate blood vessel blockage with deposits in the cardiovascular system. Viscous heating is also included in the energy equation. The governing conservation equations for mass, momentum and energy equation are transformed into a system of nonlinear, coupled ordinary differential equations and these are solved numerically using finite element method. The effect of other important parameters such as magnetohydrodynamic parameter (Nm), Reynolds number (Re), Eckert number (Ec), Darcian parameter (), Forchheimer parameter (NF) and Prandtl number on velocity and temperature profiles are studied graphically. Spatial-temporal velocity and temperature profile visualizations are also presented. Numerical results shows that normalized fluid velocity (U) increases throughout the channel (-1 < Y < 1) with an increase in Reynolds number, Darcian parameter, steady pressure gradient parameter and rheological parameter; conversely velocity is decreased with the increase in magnetic parameter and Forchheimer quadratic drag parameter. Higher Eckert number (Ec = 3) is also found to have a considerable effect on temperature () profile. Finite difference numerical computations are also compared with the finite element solutions to verify efficiency and accuracy.