Computational modeling of biomagnetic micropolar blood flow and heat transfer in a two-dimensional non-Darcian porous medium

We study theoretically and computationally the incompressible, non-conducting, micropolar, biomagnetic (blood) flow and heat transfer through a two-dimensional square porous medium in an (x,y) coordinate system, bound by impermeable walls. The magnetic field acting on the fluid is generated by an electrical current flowing normal to the x-y plane, at a distance l beneath the base side of the square. The flow regime is affected by the magnetization B (0) and a linear relation is used to define the relationship between magnetization and magnetic field intensity. The steady governing equations for x-direction translational (linear) momentum, y-direction translational (linear) momentum, angular momentum (micro-rotation) and energy (heat) conservation are presented. The energy equation incorporates a special term designating the thermal power per unit volume due to the magnetocaloric effect. The governing equations are non-dimensionalized into a dimensionless (xi,eta) coordinate system using a set of similarity transformations. The resulting two point boundary value problem is shown to be represented by five dependent non-dimensional variables, f (xi) (velocity), f (eta) (velocity), g (micro-rotation), E (magnetic field intensity) and theta (temperature) with appropriate boundary conditions at the walls. The thermophysical parameters controlling the flow are the micropolar parameter (R), biomagnetic parameter (N (H) ), Darcy number (Da), Forchheimer (Fs), magnetic field strength parameter (Mn), Eckert number (Ec) and Prandtl number (Pr). Numerical solutions are obtained using the finite element method and also the finite difference method for Ec=2.476x10(-6) and Prandtl number Pr=20, which represent realistic biomagnetic hemodynamic and heat transfer scenarios. Temperatures are shown to be considerably increased with Mn values but depressed by a rise in biomagnetic parameter (N (H) ) and also a rise in micropolarity (R). Translational velocity components are found to decrease substantially with micropolarity (R), a trend consistent with Newtonian blood flows. Micro-rotation values are shown to increase considerably with a rise in R values but are reduced with a rise in biomagnetic parameter (N (H) ). Both translational velocities are boosted with a rise in Darcy number as is micro-rotation. Forchheimer number is also shown to decrease translational velocities but increase micro-rotation. Excellent agreement is demonstrated between both numerical solutions. The mathematical model finds applications in blood flow control devices, hemodynamics in porous biomaterials and also biomagnetic flows in highly perfused skeletal tissue.