Investigation of the Effects of Some Physical Parameters and Hall Current on MHD Fluid Flow with Heat Flux over a Porous Medium

Show more

1. Introduction

The problems of fluid flow continue to attract the attention of so many researchers in engineering science, applied mathematics and applied geophysics. Heat and mass transfer is used in aerodynamic, extrusion of plastic sheets and other engineering processes which include chemical engineering and petroleum engineering [1]. Flows through porous media are very much prevalent in nature and therefore the study of flows through porous medium has become important to many scientists. In view of the above [2], studied MHD free and forced convection flow of Newtonian fluid through a porous medium in an infinite vertical plate in the presence of thermal radiation, heat transfer and surface temperature oscillation. [3] examined heat and mass transfer in MHD viscous elastic fluid flow through a porous medium over a stretching sheet with chemical reaction. In the last six years many investigations dealing with heat flow problem over a vertical porous plate with variable suction, soret effects, heat absorption, dissipation and have been investigated by so many researchers. Hall current effects on unsteady magnetohydrodynamics flow between streaching sheet and an oscillating porous upper parallel plate with constant suction [4]. Mass Transfer and Free Convection Effects on the Flow Past an Accelerated Vertical Plate with Variable Suction or Injection, was illustrated by [5]. [6] studied radiation effects on heat and mass transfer over a vertical plate with Newtonian fluid. Hall current effects on unsteady Mhd fluid flow with radiative heat flux and heat source over a porous medium, was analytically examined by [7]. [8] conducted Analytical study of MHD radiation-convection with surface temperature oscillation and secondary flow effects. [9] have investigated the effects of hall current on the steady hydromagnetic flow in an incompressible viscous fluid. [2] studied steady and unsteady magneto hydrodynamic (MHD) viscous, incompressible free and forced convective flow of an electrically conducting Newtonian fluid through a porous medium in the presence of appreciable thermal radiation heat transfer and surface temperature oscillation taking hall current into account. In addition [10] analyzed MHD free convection and mass transfer flows with hall current, viscous dissipation joul heating and thermal diffusion. [11] analyzed MHD free convection and mass transfer flow in a porous media with simultaneous rotating fluid. [12] had discussed the MHD unsteady free convective Walter’s memory flow with constant suction and heat sink. [13] investigated radiation effects on heat and mass transfer over an exponentially accelerated infinite vertical plate with chemical reaction. Combined influence of radiation absorption and hall current effects on MHD double diffusive free convective flow past a stretching sheet, was demonstrated [14]. [15] analyzed hall current effects on MHD free convective heat and mass transfer flow past an oscillating vertical porous plate in the presence of radiation. Effects of radiation, chemical reaction and porosity of the medium on MHD flow past an oscillating plate with heat and mass transfer was examined by [16]. [17] studied Hall current effects on unsteady magnetohydrodynamics flow between stretching sheet and an oscillating porous upper parallel plate with constant suction. The effect of slip conditions and Hall current on unsteady MHD flow of a viscoelastic fluid past on infinite vertical porous plate through porous medium was studied by [18], the research investigated that primary velocity significantly decreases with the increases of Hall current parameter. Considering [18], in this research the effects of Hall current on MHD fluid flow are investigated to include radiative heat flux and concentration over a porous medium.

2. Problem Formulation

Consider the flow of incompressible memory fluid in an infinite plane with heat and mass transfer, under the influence of an induced magnetic field and constant suction. The *x*-axis is taken along the plane in the upward direction and a straight line perpendicular to that of the *y*-axis. All fluid properties are assumed constant. Since the fluid is conducting, the magnetic Reynolds number is much less than unity and hence the induced magnetic field is not neglected.

The equations governing the flow under Boussineqs approximation are:

Continuity equation:

$\frac{\partial v}{\partial y}=0$

Momentum equation

$\begin{array}{l}\frac{\partial {u}^{\prime}}{\partial t}+{\upsilon}_{0}\frac{\partial {u}^{\prime}}{\partial {y}^{\prime}}=\upsilon \frac{{\partial}^{2}{u}^{\prime}}{\partial {y}^{2}}-\upsilon \frac{{u}^{\prime}}{{K}^{*}}+g{\beta}^{*}\left({{C}^{\prime}}_{0}-{{C}^{\prime}}_{d}\right)+g{\beta}^{*}\left({{T}^{\prime}}_{0}-{{T}^{\prime}}_{d}\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}-\frac{\sigma {\beta}_{0}^{2}}{\rho \left(1+{m}^{2}\right)}\left(u+mw\right)-{K}_{1}\left\{\frac{{\partial}^{3}{u}^{\prime}}{\partial t\partial {y}^{2}}\right\}\end{array}$ (1)

$\frac{\partial {w}^{\prime}}{\partial {t}^{\prime}}+{\upsilon}_{0}\frac{\partial {w}^{\prime}}{\partial {y}^{\prime}}=\upsilon \frac{{\partial}^{2}{w}^{\prime}}{\partial {{y}^{\prime}}^{2}}-\upsilon \frac{{w}^{\prime}}{{K}^{*}}-\frac{\sigma {\beta}_{0}^{2}}{\rho \left(1+{m}^{2}\right)}\left(w-mu\right)-{K}_{1}\left\{\frac{{\partial}^{3}{w}^{\prime}}{\partial t\partial {{y}^{\prime}}^{2}}\right\}$ (2)

Energy equation:

$\frac{\partial {T}^{\prime}}{\partial {t}^{\prime}}+{\upsilon}_{0}\frac{\partial {T}^{\prime}}{\partial {y}^{\prime}}=\frac{k}{\rho cp}\frac{{\partial}^{2}{T}^{\prime}}{\partial {{y}^{\prime}}^{2}}-\frac{1}{\rho cp}\frac{\partial {q}_{r}}{\partial y}+\frac{1}{\rho cp}{\left(\frac{\partial {u}^{\prime}}{\partial {y}^{\prime}}\right)}^{2}$ (3)

Concentration equation:

$\frac{\partial {C}^{\prime}}{\partial {t}^{\prime}}+{\upsilon}_{0}\frac{\partial {C}^{\prime}}{\partial {y}^{\prime}}=D\frac{{\partial}^{2}{C}^{\prime}}{\partial {{y}^{\prime}}^{2}}-K{C}^{\prime}$ (4)

The initial boundary conditions are:

$\begin{array}{l}u=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}w=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}T={T}_{d}+\left({T}_{0}-{T}_{d}\right)\epsilon {\text{e}}^{i\omega t},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}C={C}_{d}+\left({C}_{0}-{C}_{d}\right)\epsilon {\text{e}}^{\iota \omega t}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{at}\text{\hspace{0.17em}}\text{\hspace{0.17em}}y=0\\ u=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}w=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}T=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}C=0\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{at}\text{\hspace{0.17em}}\text{\hspace{0.17em}}y=d\end{array}\}$ (5)

where *u* is the velocity of the fluid in the *x*-direction and *v* in the *y*-direction, *T* is the temperature of the fluid, *C* is the concentration of the fluid, *g* is the acceleration due to gravity,
${\beta}^{*}$ are the kinematic viscosity, *K* is the thermal conductivity and *Cp* is the specific heat capacity of the fluid at constant pressure. *t* is the time,
$\sigma $ is the electrical conductivity of the fluid and
${\mu}_{e}$ is the magnetic permeability. *T*_{0} is the temperature of the plate and *T _{d}* is the temperature of the fluid far away from plate.

Introducing the following non-dimensional quantities,

$\begin{array}{l}{{w}^{\prime}}_{1}=\frac{U}{e}{w}_{1},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\partial {{w}^{\prime}}_{1}=\frac{U}{e}\text{d}{w}_{1},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\partial}^{2}{{w}^{\prime}}_{1}=\frac{U}{e}{\text{d}}^{2}{w}_{1},\\ y=\frac{{y}^{\prime}}{d},t=\frac{t{\upsilon}_{0}}{\upsilon}\Rightarrow \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{t}^{\prime}=\frac{\upsilon t}{{\upsilon}_{0}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\partial {t}^{\prime}=\frac{\upsilon}{{\upsilon}_{0}}\text{d}t,\\ \theta =\frac{{T}^{\prime}-{{T}^{\prime}}_{d}}{{{T}^{\prime}}_{0}-{{T}^{\prime}}_{d}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}M=\frac{\sigma {\beta}_{0}^{2}\upsilon}{\rho {\upsilon}_{0}^{2}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{K}^{*}=\frac{{K}_{0}{\upsilon}_{0}^{2}}{{\upsilon}^{2}}\\ C=\frac{{C}^{\prime}-{{C}^{\prime}}_{d}}{{{C}^{\prime}}_{0}-{{C}^{\prime}}_{d}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}Gc=\frac{\upsilon {\beta}^{*}\upsilon \left({{C}^{\prime}}_{\omega}-{{C}^{\prime}}_{\infty}\right)}{U{\upsilon}_{0}^{2}}\end{array}\}$ (6)

The term $\frac{\partial {q}_{r}}{\partial y}$ represents the radiative heat flux. By using Rosseland approximation, the radiation heat flux ${q}_{r}=-\frac{4{\sigma}^{\ast}\partial {T}^{4}}{3{a}_{R}\partial y}$, where and Stephen and

Boltzmann constants and mean absorption coefficient respectively. We assume that the temperature difference within the flow is such that may be expanded in a Taylor’s series.

Hence, expanding ${T}^{4}$ about ${T}_{d}$ and neglecting higher order terms, we get

${{T}^{\prime}}^{4}={{T}^{\prime}}_{d}^{4}+4{T}_{d}^{3}{T}^{\ast}=4{T}_{d}^{3}{T}^{\prime}-3{{T}^{\prime}}_{d}^{4}$

We assume the following solutions:

$\begin{array}{l}{u}^{\prime}={{u}^{\prime}}_{0}+\epsilon {{u}^{\prime}}_{1}{\text{e}}^{i\omega t}\\ {w}^{\prime}={{w}^{\prime}}_{0}+\epsilon {{w}^{\prime}}_{1}{\text{e}}^{i\omega t}\\ {\theta}^{\prime}={{\theta}^{\prime}}_{0}+\epsilon {{\theta}^{\prime}}_{1}{\text{e}}^{i\omega t}\\ {C}^{\prime}={{C}^{\prime}}_{0}+\epsilon {{C}^{\prime}}_{1}{\text{e}}^{i\omega t}\end{array}\}$ (7)

$\frac{{\partial}^{2}{u}_{0}}{\partial {y}^{2}}+\frac{\partial {u}_{0}}{\partial y}-L{u}_{0}-J{w}_{0}+Gc{C}_{0}=0$ (8)

where $L=\left(\frac{1}{{K}_{s}}+\frac{M}{1+{m}^{2}}\right)$, $J=\left(\frac{Mm}{1+{m}^{2}}\right)$

$\frac{{\text{d}}^{2}{w}_{0}}{\text{d}{y}^{2}}-\frac{\text{d}{w}_{0}}{\text{d}y}-L{w}_{0}+J{u}_{0}=0$ (9)

Combining (8) and (9) using complex variable method, we have,

$\frac{{\text{d}}^{2}F}{\text{d}{y}^{2}}-\frac{\text{d}F}{\text{d}y}-FL-FJ=-Gc{C}_{0}$ (10)

where $F=\left({u}_{0}+i{w}_{0}\right)$ and $i=\sqrt{-1}$

$\frac{{\text{d}}^{2}F}{\text{d}{y}^{2}}-\frac{\text{d}F}{\text{d}y}-{P}_{1}F=-Gc{C}_{0}$ (11)

where ${P}_{\text{1}}=\left(L+J\right)$

$\Rightarrow \left(1-i\omega {K}_{1}\right)\frac{{\text{d}}^{2}{u}_{1}}{\text{d}{y}^{2}}-\frac{\text{d}{u}_{1}}{\text{d}y}-Ln{u}_{1}+J{w}_{1}=-Gc{{C}^{\prime}}_{1}$ (12)

where $Ln=\left(L+i\omega \right)$

$\Rightarrow \left(1-i\omega {K}_{1}\right)\frac{{\text{d}}^{2}{w}_{1}}{\text{d}{y}^{2}}-\frac{\text{d}{w}_{1}}{\text{d}y}-\left(L+i\omega \right){w}_{1}+J{u}_{1}=0$

$\Rightarrow \left(1-i\omega {K}_{1}\right)\frac{{\text{d}}^{2}{w}_{1}}{\text{d}{y}^{2}}-\frac{\text{d}{w}_{1}}{\text{d}y}-{L}_{0}{w}_{1}+J{u}_{1}=0$ (13)

Using the method of complex variable, combining (12) and (13) to have,

${P}_{2}\frac{{\text{d}}^{2}H}{\text{d}{y}^{2}}-\frac{\text{d}H}{\text{d}y}-{P}_{3}H=-Gc{C}_{1}$ (14)

where $H=\left({u}_{1}+i{w}_{1}\right)$, ${P}_{2}=\left(1-i\omega {K}_{1}\right)$, $L=\left({L}_{0}+i\omega \right)$.

The new boundary conditions are,

$\begin{array}{l}{{u}^{\prime}}_{0}={{u}^{\prime}}_{1}=0\text{\hspace{0.17em}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{{w}^{\prime}}_{0}={{w}^{\prime}}_{1}=0\text{\hspace{0.17em}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{{\theta}^{\prime}}_{0}={{\theta}^{\prime}}_{1}=1\text{\hspace{0.17em}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{{C}^{\prime}}_{0}={{C}^{\prime}}_{1}=1\text{\hspace{0.17em}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{at}\text{\hspace{0.17em}}\text{\hspace{0.17em}}y=0\\ {{u}^{\prime}}_{0}={{u}^{\prime}}_{1}=0\text{\hspace{0.17em}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{{w}^{\prime}}_{0}={{w}^{\prime}}_{1}=0\text{\hspace{0.17em}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{{\theta}^{\prime}}_{0}={{\theta}^{\prime}}_{1}=0\text{\hspace{0.17em}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{{C}^{\prime}}_{0}={{C}^{\prime}}_{1}=0\text{\hspace{0.17em}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{at}\text{\hspace{0.17em}}\text{\hspace{0.17em}}y=1\end{array}\}$ (15)

3. Method of Solution

To solve for mass diffusion, therefore we assume:

Concentration to be

$\begin{array}{l}C\left(y,t\right)={C}_{0}+{C}_{1}\epsilon \mathrm{exp}\left(i\omega t\right)\text{\hspace{0.17em}}\\ C\left(y,t\right)={A}_{1}\mathrm{exp}\left({m}_{1}y\right)+{A}_{2}\mathrm{exp}\left({m}_{2}y\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}+\left({A}_{3}\mathrm{exp}\left({m}_{3}y\right)+{A}_{4}\mathrm{exp}\left({m}_{4}y\right)\right)\epsilon \mathrm{exp}\left(i\omega t\right)\end{array}$ (16)

To solve for the momentum equation then

$F\left(h\right)={A}_{5}\mathrm{exp}\left({m}_{5}y\right)+{A}_{6}\mathrm{exp}\left({m}_{6}y\right)$ (17)

$\frac{{\text{d}}^{2}H}{\text{d}{y}^{2}}-{P}_{4}\frac{\text{d}H}{\text{d}y}-{P}_{5}H=0$ (18)

$H\left(y\right)={A}_{7}\mathrm{exp}\left({m}_{7}y\right)+{A}_{8}\mathrm{exp}\left({m}_{8}y\right)+{D}_{3}\mathrm{exp}\left({m}_{3}y\right)+{D}_{4}\mathrm{exp}\left({m}_{4}y\right)$ (19)

Therefore, the solution of primary velocity is assumed to be

$F\left({u}_{0}+{u}_{1}\right)={u}_{0}+{u}_{1}\epsilon {\text{e}}^{i\omega t}$

$\begin{array}{l}F\left({u}_{0}+{u}_{1}\right)\\ ={A}_{5}\mathrm{exp}\left({m}_{5}y\right)+{A}_{6}\mathrm{exp}\left({m}_{6}y\right)+{D}_{1}\mathrm{exp}\left({m}_{1}y\right)+{D}_{2}\mathrm{exp}\left({m}_{2}y\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}+{A}_{7}\mathrm{exp}\left({m}_{7}y\right)+{A}_{8}\mathrm{exp}\left({m}_{8}y\right)+{D}_{3}\mathrm{exp}\left({m}_{3}y\right)+{D}_{4}\mathrm{exp}\left({m}_{4}y\right)\epsilon {\text{e}}^{i\omega t}\end{array}$ (20)

Then, the secondary velocity also is as follows:

$H\left({w}_{0}+{w}_{1}\right)={w}_{0}+{w}_{1}\epsilon {\text{e}}^{i\omega t}$

$\begin{array}{l}H\left({w}_{0}+{w}_{1}\right)\\ =i\{{A}_{5}\mathrm{exp}\left({m}_{5}y\right)+{A}_{6}\mathrm{exp}\left({m}_{6}y\right)+{D}_{1}\mathrm{exp}\left({m}_{1}y\right)+{D}_{2}\mathrm{exp}\left({m}_{2}y\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+{A}_{7}\mathrm{exp}\left({m}_{7}y\right)+{A}_{8}\mathrm{exp}\left({m}_{8}y\right)+{D}_{3}\mathrm{exp}\left({m}_{3}y\right)+{D}_{4}\mathrm{exp}\left({m}_{4}y\right)\epsilon {\text{e}}^{i\omega t}\}\end{array}$ (21)

From the energy equation we have,

$\frac{\partial {T}^{\prime}}{\partial {t}^{\prime}}+{\upsilon}_{0}\frac{\partial {T}^{\prime}}{\partial {y}^{\prime}}=\frac{k}{\rho cp}\frac{{\partial}^{2}{T}^{\prime}}{\partial {{y}^{\prime}}^{2}}-\frac{1}{\rho cp}\frac{\partial {q}_{r}}{\partial y}+\frac{1}{\rho cp}{\left(\frac{\partial {u}^{\prime}}{\partial {y}^{\prime}}\right)}^{2}$

Substituting the value of $\partial {T}^{\prime}$ yields

$\frac{\partial \theta}{\partial t}+\frac{\partial \theta}{\partial y}=\frac{1}{Pr}\frac{{\partial}^{2}\theta}{\partial {y}^{2}}+\frac{4R{\partial}^{2}\theta}{3Pr\partial {y}^{2}}+Ec{\left(\frac{\partial {u}^{\prime}}{\partial {y}^{\prime}}\right)}^{2}$

where $Ec=\frac{\upsilon}{\rho {C}_{p}\left({{T}^{\prime}}_{0}-{{T}^{\prime}}_{d}\right)}$

$\frac{\partial \theta}{\partial t}+\frac{\partial \theta}{\partial y}=\frac{E}{Pr}\frac{{\partial}^{2}\theta}{\partial {y}^{2}}+Ec{\left(\frac{\partial {u}^{\prime}}{\partial {y}^{\prime}}\right)}^{2}$

where $E=\left(1+\frac{4R}{3}\right)$

$\frac{\partial \theta}{\partial t}+\frac{\partial \theta}{\partial y}=\frac{E}{Pr}\frac{{\partial}^{2}\theta}{\partial {y}^{2}}+Ec{\left(\frac{\partial {u}^{\prime}}{\partial {y}^{\prime}}\right)}^{2}$ (22)

We also assumed solution to be

$\begin{array}{l}\theta ={\theta}_{0}+{\theta}_{1}\epsilon {\text{e}}^{i\omega t}\\ u={u}_{0}+{u}_{1}\epsilon {\text{e}}^{i\omega t}\end{array}\}$

Therefore, the solution of the energy equation is as follows:

$\begin{array}{c}\theta =\{A+B\mathrm{exp}\left(Pr{E}_{2}y\right)+{D}_{11}\mathrm{exp}\left(2{m}_{5}y\right)+{D}_{12}\mathrm{exp}\left(2{m}_{6}y\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+{D}_{13}\mathrm{exp}\left(2{m}_{1}y\right)+{D}_{14}\mathrm{exp}\left(2{m}_{2}y\right)+{D}_{15}\mathrm{exp}\left({m}_{5}+{m}_{6}\right)y\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+{D}_{16}\mathrm{exp}\left({m}_{5}+{m}_{1}\right)y+{D}_{17}\mathrm{exp}\left({m}_{5}+{m}_{2}\right)y+{D}_{18}\mathrm{exp}\left({m}_{6}+{m}_{1}\right)y\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+{D}_{19}\mathrm{exp}\left({m}_{6}+{m}_{2}\right)y+{D}_{20}\mathrm{exp}\left({m}_{1}+{m}_{2}\right)y\}\end{array}$

$\begin{array}{l}\text{\hspace{0.05em}}+\{{A}_{9}\mathrm{exp}\left({m}_{9}y\right)+{A}_{10}\mathrm{exp}\left({m}_{10}y\right)+{D}_{21}\mathrm{exp}\left({m}_{5}+{m}_{7}\right)y\\ \text{\hspace{0.05em}}+{D}_{22}\mathrm{exp}\left({m}_{5}+{m}_{8}\right)y+{D}_{23}\mathrm{exp}\left({m}_{5}+{m}_{3}\right)y+{D}_{24}\mathrm{exp}\left({m}_{5}+{m}_{4}\right)y\\ \text{\hspace{0.05em}}+{D}_{25}\mathrm{exp}\left({m}_{6}+{m}_{7}\right)y+{D}_{26}\mathrm{exp}\left({m}_{6}+{m}_{8}\right)y+{D}_{27}\mathrm{exp}\left({m}_{6}+{m}_{3}\right)y\\ \text{\hspace{0.05em}}+{D}_{28}\mathrm{exp}\left({m}_{6}+{m}_{4}\right)y+{D}_{29}\mathrm{exp}\left({m}_{1}+{m}_{7}\right)y+{D}_{30}\mathrm{exp}\left({m}_{1}+{m}_{8}\right)y\end{array}$

$\begin{array}{l}\text{\hspace{0.05em}}+{D}_{31}\mathrm{exp}\left({m}_{1}+{m}_{3}\right)y+{D}_{32}\mathrm{exp}\left({m}_{1}+{m}_{4}\right)y+{D}_{33}\mathrm{exp}\left({m}_{2}+{m}_{7}\right)y\\ \text{\hspace{0.05em}}+{D}_{34}\mathrm{exp}\left({m}_{2}+{m}_{8}\right)y+{D}_{35}\mathrm{exp}\left({m}_{2}+{m}_{3}\right)y\\ \text{\hspace{0.05em}}+{D}_{36}\mathrm{exp}\left({m}_{2}+{m}_{4}\right)y\}\epsilon \mathrm{exp}\left(i\omega t\right)\end{array}$ (23)

4. Results and Discussion

4.1. Velocity Profiles

Figures 1-10 represent velocity profile for the flow.

4.2. Concentration Profiles

Figure 11 and Figure 12 represent concentration profiles for the flow.

4.3. Temperature Profiles

Figure 13 and Figure 14 represent temperature profiles for different values of Prandtle number (*Pr*) and radiation parameter (*R*) of the flow.

Figure 1. Primary velocity profiles for different values of mass Grashof number (*Gc*).This indicates that the primary velocity attends its peak at *y* = 0.5, in the same vein primary velocity increase with increase in the values of *Gc*.

Figure 2. Secondary velocity profiles for different values of Thermal Grashof number (*Gc*). The least value was attended at *y* = 0.45, consequently the velocity of the fluid increase with increase in values of *Gc*.

Figure 3. Primary velocity profiles for different values of viscoelastic parameter (*K*_{1}). In this flow primary velocity decrease with increasing values of *K*_{1}. This is in agreement with general belief that viscoelastic causes drag in a flow of a fluid.

Figure 4. Secondary velocity profiles for different values of *K*_{1}, this demonstrated the same behavior as in primary velocity.

Figure 5. Primary velocity profiles for different values of Schmidt number (*Sc*). As *Sc* increases the primary velocity decreases.

Figure 6. Secondary velocity profiles for different values of *Sc* indicated the same trend as in the case of primary velocity.

Figure 7. Primary velocity profiles for different values of hall current parameter (*m*). When *m* = 0 the velocity reached its maximum at 0.19 units, *y* = 0.45, as the values of m increases the velocity reduced to as low as 0.08 units. This means that hall current has adverse effects on flow of this fluid.

Figure 8. Secondary velocity profiles for different values of *m*. when *m* = 0 the flow was low, it became lowest as the values of m increased.

Figure 9. Primary velocity profiles for different values of magnetic field (*M*). The effects of *M* in this flow are insignificant, even though it reduces the flow of the fluid.

Figure 10. Primary velocity profiles for different values of magnetic field (*M*).

Figure 11. Concentration profiles for different values of chemical reaction parameter (*Kr*).

Figure 12. Concentration profiles for different values of Schmidt number (*Sc*).

Figure 13. Temperature profiles for different values of Prantle number (*Pr*). The graph demonstrated uniformity where *Pr* = 0.71 and 1.0, while the last two values did not show significant effect on temperature.

Figure 14. Temperature profiles for different values of radiation parameter (*R*). The effects of radiation parameter are insignificant as demonstrated in the figure above.

5. Conclusion

Investigation of effects of some physical parameters and hall current on MHD fluid flow with heat flux over a porous medium is studied by transforming the governing partial differential equations into ordinary differential equations which are then solved using perturbation techniques. The result of the flow variables indicates that the fluid temperature is reduced by increasing Prandtl number (*Pr*) and radiation parameter (*R*). Concentration is reduced with increase in Schmidt number (*Sc*) and chemical reaction parameter (*K*). The primary velocity decrease with increasing prandtl (*Pr*), radiation parameter and hall-current while the opposite trend is observed in secondary velocity. The primary velocity increases with increase in mass Grashoof number (*Gr*) and thermal Grasshoof number (*Gc*) also the reverse is the case in secondary velocity. The primary velocity decreases with increase in *M*, *s* and *Sc*.

References

[1] Alam, M.S., Rahaman, M.M. and Satar, M.A. (2008) Effects of Variable Suction and Thermophoresis on Steady MHD Combined Free-Forced Convective Heat and Mass Transfer Flow over a Semi Infinite Permeable Inclined Plate in the Presence of Thermal Radiation. International Journal of Thermal Science, 47, 758-765.

https://doi.org/10.1016/j.ijthermalsci.2007.06.006

[2] Veera Krishna, M. and Sujatha, P. (2011) MHD Free and Forced Convection Flow of Newtonian Fluid through a Porous Medium in an Infinite Vertical Plate in Presence of Thermal Radiation Heat Transfer and Surface Temperature Oscillation. International Journal of Physics and Mathematical Science, 1, 33-52.

[3] Saleh, M.A., Mohammed, A.A.B. and Mahmoud, S.E. (2010) Heat and Mass Transfer in MHD Viscous Elastic Fluid Flow through a Porous Medium over a Stretching Sheet with Chemical Reaction, Journal of Applied Mathematics, 1, 446-455.

https://doi.org/10.4236/am.2010.16059

[4] Malraju, C.R., Narravula, A.R. and Sibyale, V.K. (2011) Hallcurrent Effects on Unsteady Magnetohydrodynamics Flow between Stretching Sheet and an Oscillating Porous Upper Parallel Plate with Constant Suction. Thermal Science, 15, 527-536.

https://doi.org/10.2298/TSCI1102527C

[5] Kafousias, N.G. and Raptis, A.A. (1981) Mass Transfer and Free Convection Effects on the Flow Past an Accelerated Vertical Plate with Variable Suction or Injection. International Journal of Engineering Science, 26, 11-22.

[6] Uwanta, I.J. and Ahmed, A. (2012) Radiation Effects on Heat and Mass Transfer over a Vertical Plate with Newtonian Fluid. International Journal of Engineering Research and Technology (IJERT), 1, 1-15.

[7] Ahmed, A., Uwanta, I.J. and Sarki, M.N. (2015) Hall Current Effects on Unsteady MHD Fluid Flow with Radiative Heat Flux and Heat Source over a Porous Medium. British Journal of Mathematics and Computer Science, 6, 233-246.

https://doi.org/10.9734/BJMCS/2015/13849

[8] Anwerbeg, O. and Ghosh, S.K. (2010) Analytical Study of Magneto Hydrodynamic Radiation-Convection with Surface Temperature Oscillation and Secondary Flow Effects. International Journal of Applied Mathematics and Mechanics, 6, 1-22.

[9] Gupta, A.S. (1975) Hydro Magnetic Flow past a Porous Flat Plate with Hall Effects. Acta Mechanica, 22, 281-287.

https://doi.org/10.1007/BF01170681

[10] Singh, A.K. and Singh, N.P. (2003) Heat and Mass Transfer in MHD Flow of a Viscous Fluid Past a Vertical Plate under Oscillatory Suction Velocity. Indian Journal of Pure and Applied Mathematics, 34, 429-442.

[11] Ferdows, M., Kainoota, K. and Crepeau, J.C. (2008) MHD Free Convection and Mass Transfer Flow in a Porous Media with Simultaneous Rotating Fluid. International Journal of Dynamics of Fluids, 4, 69-82.

[12] Ramana Murthy, M.V., et al. (2007) MHD Unsteady Free Convective Walter’s Memory Flow with Constant Suction and Heat Sink. AAPN Journal of Engineering and Applied Science, 2, 12-16.

[13] Ahmed, A., Sarki, M.N. and Ahmad, M. (2014) Radiation Effects on Heat and Mass Transfer over an Exponentially Accelerated Infinite Vertical Plate with Chemical Reaction. Proceeding of the International Multiconference of Engineers and Computer Scientists, Vol. II, Hong Kong, 12-14 March 2014, 914-918.

[14] Sreedevi, R.R., Parasoda, D.R.V. and Chamkha, A.J. (2016) Combined Influence of Radiation Absorption and Hall Current Effects on MHD Double Diffusive Free Convective Flow past a Stretching Sheet. Ain Shams Engineering Journal, 7, 383-397.

https://doi.org/10.1016/j.asej.2015.11.024

[15] Babu, D. and Raju, G.S.S. (2017) Hall Current Effects on MHD Free Convective Heat and Mass Transfer Flow past an Oscillating Vertical Porous Plate in the Presence of Radiation, Thermal and Mass Diffusion with Chemical Reaction. International Journal of Mathematics and Computer Application Research (IJMCAR), 7, 1-22.

[16] Rajput, U.S. and Kumar, G. (2017) Effects of Radiation, Chemical Reaction and Porosity of the Medium on MHD Flow past an Oscillating Plate with Heat and Mass Transfer. Aryabhatta Journal of Mathematics and Informatics, 9, 230-241.

[17] Changal, M.R., Narravula, A.R. and Silbiya, V.K. (2011) Hall Current Effects on Unsteady Magnetohydrodynamics Flow between Stretching Sheet and an Oscillating Porous Upper Parallel Plate with Constant Suction. Thermal Science, 15, 527-536.

https://doi.org/10.2298/TSCI1102527C

[18] Kumar, R. and Chand, K. (2011) Effect of Slip Conditions and Hall Current on Unsteady MHD Flow of a Viscoelastic Fluid Past on Infinite Vertical Porous Plate through Porous Medium. International Journal of Engineering Science and Technology, 3, 3124-3133.