Finite element and network electrical simulation of rotating magnetofluid flow in nonlinear porous media with inclined magnetic field and hall currents


O. Bég Anwar, S. Rawat, J. Zueco, L. Osmond, R. S. R. Gorla




A mathematical model is presented for viscous hydromagnetic flow through a hybrid non-Darcy porous media rotating generator. The system is simulated as steady, incompressible flow through a nonlinear porous regime intercalated between parallel plates of the generator in a rotating frame of reference in the presence of a strong, inclined magnetic field A pressure gradient term is included which is a function of the longitudinal coordinate. The general equations for rotating viscous magnetohydrodynamic flow are presented and neglecting convective acceleration effects, the two-dimensional viscous flow equations are derived incorporating current density components, porous media drag effects, Lorentz drag force components and Hall current effects. Using an appropriate group of dimensionless variables, the momentum equations for primary and secondary flow are rendered non-dimensional and shown to be controlled by six physical parameters-Hartmann number (Ha), Hall current parameter (Nh), Darcy number (Da), Forchheimer number (Fs), Ekman number (Ek) and dimensionless pressure gradient parameter (Np), in addition to one geometric parameter- the orientation of the applied magnetic field. Several special cases are extracted from the general model, including the non-porous case studied earlier by Ghosh and Pop (2006). A numerical solution is presented to the nonlinear coupled ordinary differential equations using both the Network Simulation Method and Finite Element Method, achieving excellent agreement. Additionally very good agreement is also obtained with the earlier analytical solutions of Ghosh and Pop (2006). for selected Ha, Ek and Nh values. We examine in detail the effects of magnetic field, rotation, Hall current, bulk porous matrix drag, second order porous impedance, pressure gradient and magnetic field inclination on primary and secondary velocity distributions and also frictional shear stresses at the plates. Primary velocity is seen to decrease with an increase in Hall current parameter (Nh) with the converse observed for the secondary velocity.