Abstract

This paper presents an incompressible two-dimensional heat and mass transfer of an electrically conducting micropolar fluid flow in a porous medium between two parallel plates with chemical reaction, Hall and ion slip effects. Let there be periodic injection or suction at the lower and upper plates and the nonuniform temperature and concentration at the plates are varying periodically with time. The flow field equations are reduced to nonlinear ordinary differential equations using similarity transformations and then solved numerically by quasilinearization technique. The profiles of velocity components, microrotation, temperature distribution and concentration are studied for different values of fluid and geometric parameters such as Hartmann number, Hall and ion slip parameters, inverse Darcy parameter, Prandtl number, Schmidt number, and chemical reaction rate and shown in the form of graphs.

1. Introduction

The flow and heat transfer of a fluid through porous channels is of great importance in both engineering and science. Applications of these were found in different areas such as oil exploration, geothermal energy extractions, the boundary layer control, extrusion of polymer fluids, solidification of liquid crystals, cooling of a metallic plate in a bath, MHD power generators, and suspension solutions. Many authors investigated the two-dimensional incompressible fluid flow through porous channels theoretically. Terrill and Shrestha [1] considered the laminar incompressible flow of an electrically conducting viscous fluid through a porous channel and gave a solution for a large suction Reynolds number and Hartmann number. Eringen [2, 3] initiated the theory of micropolar fluids and this theory constitutes a subclass of microfluids. Attia [4] considered the problem on MHD flow of a dusty fluid in a circular pipe with Hall and ion slip and obtained a series solution for reduced governing equations. A perturbation solution was obtained for heat and mass transfer in a rotating vertical channel with Hall current by Ahmed and Zueco [5]. Eldabe and Ouaf [6] investigated the problem on MHD flow and heat and mass transfer of a micropolar fluid over a stretching surface with ohmic heating. Magyari et al. [7] have studied Stokes’ first problem for micropolar fluid and solved the problem analytically by Laplace transforms and numerically by Valkó-Abate procedure. Bhattacharyya et al. [8] studied the effects of chemical reaction on the boundary layer flow of viscous fluid and a numerical solution obtained by shooting method. The problem of steady incompressible MHD flow of a micropolar fluid between concentric porous cylinders was discussed by Srinivasacharya and Shiferaw [9] who obtained a numerical solution by quasilinearization technique. The steady micropolar fluid flow between two vertical porous parallel plates in the presence of magnetic field was considered by Bhargava et al. [10] and the governing differential equations have been solved using the quasilinearization method. Abdulaziz et al. [11] examined the MHD fully developed natural convection flow of micropolar fluid between parallel vertical plates and the reduced governing nonlinear equations are solved by using HAM method. Das [12] studied the effects of thermal radiation and chemical reaction on incompressible flow of electrically conducting micropolar fluid over an inclined plate. Pal et al. [13] investigated the problem of oscillatory mixed convection-radiation of a micropolar fluid in a rotating system with Hall current and chemical reaction effects and obtained a solution using perturbation method. Rahman and Al-Lawatia [14] examined micropolar fluid flow over a stretching sheet with variable concentration and solved numerically using the shooting method. Bakr [15] considered the steady as well as unsteady MHD micropolar fluid with constant heat source and chemical reaction effect in a rotating frame of reference and solved by using the perturbation method. Patil and Kulkarni [16] studied the free convective flow of a polar fluid in a porous medium in the presence of an internal heat generation with chemical reaction effect. The magnetomicropolar fluid flow and heat and mass transfer in a porous medium with Hall and ion slip effects were discussed by Motsa and Shateyi [17] who obtained a numerical solution using the successive linearization method. Olajuwon [18] studied the heat and mass transfer of an electrically conducting micropolar fluid flow over a horizontal flat plate in the presence of transverse magnetic field with the chemical reaction and Hall effects and the problem is solved by using Runge-Kutta method. Ariman and Cakmak [19] analyzed the three basic viscous flows of micropolar fluid between parallel plates and solved analytically. Modather et al. [20] examined an incompressible MHD flow and heat and mass transfer of a micropolar fluid over a vertical plate in a porous medium. Kim [21] investigated an incompressible micropolar fluid flow past a semi-infinite plate in a porous medium and obtained an analytical solution by perturbation method. Reddy Gorla et al. [22] obtained a numerical solution for the steady boundary layer flow and heat transfer of a micropolar fluid over a flat plate. The MHD flow of viscoelastic fluid over a porous stretching sheet was analyzed by Hymavathi and Shanker [23] and a numerical solution was obtained by using the quasilinearization method. Bhatnagar et al. [24] discussed the steady incompressible laminar flow of viscoelastic fluid through a porous cylindrical annulus and the reduced governing equations are solved numerically using quasilinearization method.

This paper deals with the unsteady two-dimensional incompressible MHD flow and heat transfer of a micropolar fluid in a porous medium between parallel plates with chemical reaction, Hall and ion slip effects. The flow is generated due to periodic suction or injection at the plates and the reduced flow field equations are solved numerically using the quasilinearization technique. The effects of various parameters on velocity components, temperature distribution, microrotation, and concentration are studied and shown graphically.

2. Formulation of the Problem

Consider an incompressible electrically conducting micropolar fluid flow through a porous medium between two parallel plates at and (Figure 8). The lower and upper plates are subjected to periodic injection and suction of the forms Real () and Real (), respectively. Assume that the temperature and concentration at the lower and upper plates are , and , , respectively.

The equations governing the flow and heat and mass transfer in the presence of a magnetic field and in the absence of body forces and body couples are given by Shercliff [25]: The coefficients , , , ,   in the above equations are related by the inequalities Neglecting the displacement currents, the Maxwell equations and the generalized Ohm’s law are where , is induced magnetic field, is the Hall parameter, is the ion slip parameter, and is magnetic permeability.

For mixed suction case, that is, , following Terril and Shreshta [1], we take the velocity, microrotation, temperature distribution, and concentration as where where , is the average entrance velocity, , and , , , , , and are functions of to be determined.

The boundary conditions of the velocity components, microrotation, temperature, and concentration are

Substituting (9) into (2), (3), (4), and (5) and then comparing the real parts on both sides, we get where prime denotes the differentiation with respect to and .

From (9), the dimensionless forms of temperature and concentration are where is the Eckert number and is the dimensionless axial variable.

The boundary conditions (10) in terms of , and are

3. Solution of the Problem

The nonlinear differential equations (11) are converted into the system of first order differential equations by the following substitution: The boundary conditions (13) in terms of are The system of equations (14) is solved numerically subject to the boundary conditions (15) using the quasilinearization method given by Bellman and Kalaba [26].

Let be an approximate current solution and be an improved solution of (14). Using Taylor’s series expansion about the current solution by neglecting the second and higher order derivative terms, the coupled first order system (14) is linearized as

To solve for , the solution to seven separate initial value problems, denoted by , , , , , , (which are the solutions of the homogeneous system corresponding to (16), and (which is the particular solution of (16), with the following initial conditions, is obtained by using the 4th order Runge-Kutta method: By using the principle of superposition, the general solution can be written as where are the unknown constants and are determined by considering the boundary conditions at . This solution is then compared with solution at the previous step and further iteration is performed if the convergence has not been achieved.

4. Results and Discussions

To understand the flow characteristics in a better way, the numerical results of the axial velocity , radial velocity , microrotation , temperature distribution , and concentration are calculated correct to six places of decimal for various values of nondimensional chemical reaction parameter Kr, ion slip parameter , Hall parameter , inverse Darcy parameter , Hartmann number , Prandtl number , and Schmidt number in the domain .

The effect of Ha on velocity components, microrotation, and temperature is presented in Figures 1(a) to 1(d). It is observed that as Ha increases the axial velocity also increases towards the center of the plates and then decreases because of the Lorentz force which tends to resist the flow and the radial velocity, microrotation, and temperature distribution are increasing from the lower plate to the upper plate. Figures 2(a) to 2(d) give the effect of on velocity components, microrotation, and temperature distribution. It can be analyzed that as increases the axial velocity also increases in the center of the channel, whereas the temperature decreases towards the upper plate. However, the radial velocity and microrotation are decreasing up to the center of the channel and then increase. This is because of the fact that the decreases the resistive force imposed by magnetic field. The effect of on velocity components, microrotation, and temperature distribution has been presented in Figures 3(a) to 3(d), respectively, and these profiles follow the same trend of . The effect of on velocity components, microrotation, and temperature distribution is shown in Figures 4(a) to 4(d). It can be observed that as increases the radial velocity and microrotation are increasing up to the center of the channel and then decrease towards the upper plate and the axial velocity gives the maximum effect at the center of the plates. However, the temperature increases from lower plate to upper plate. The effects of Kr and Sc on concentration are shown in Figures 5 and 6, respectively. From Figure 5 it is understood that as Kr increases the concentration decreases because the chemical reaction increases the rate of interfacial mass transfer. It is depicted from Figure 6 that the concentration is decreasing as Sc increases and this causes the concentration buoyancy effect. Figure 7 shows the effect of Pr on temperature. It is noticed that as Pr increases the temperature distribution also increases towards the upper plate and this is due to the fact that as Pr increases the thickness of the boundary layer decreases and this gives the larger temperature values.

5. Conclusions

In the present paper, the effects of chemical reaction and Hall and ion slip on unsteady two-dimensional laminar flow of an electrically conducting micropolar fluid through a porous medium between parallel plates with heat and mass transfer are considered. The reduced nonlinear ordinary differential equations are solved by using the quasilinearization method. The results are presented in the form of graphs for various values of fluid and geometric parameters and from these the following is concluded.(i)The chemical reaction rate reduces the concentration while the Prandtl number enhances the temperature.(ii)The effects of Hall parameter and ion slip parameter on velocity components, microrotation, and temperature are similar.(iii)The Hartmann number and inverse Darcy’s parameter have the same result for velocity components, microrotation, and temperature.

Nomenclature

:Distance between parallel plates
:Suction velocity
:Injection velocity
:Injection suction ratio,
:Fluid pressure
:Velocity vector
:Specific heat at constant temperature
:Microrotation vector
:Microrotation component
:Eckert number,
:Thermal conductivity
:Viscosity parameter
:Permeability parameter
:Velocity component in -direction
:Velocity component in -direction
:Average entrance velocity
Pr:Prandtl number,
Re:Suction Reynolds number,
:Gyration parameter
:Current density
:Nondimensional gyration parameter,
:Total magnetic field
:Induced magnetic field
:Magnetic flux density
:Rate of deformation tensor
:Electric field
Ha:Hartmann number,
:Inverse Darcy parameter,
:Nondimensional viscosity parameter,
:Nondimensional micropolar parameter,
:Nondimensional micropolar parameter,
:Temperature
:Temperature at the lower plate
:Temperature at the upper plate
:Dimensionless temperature,
:Dimensionless concentration,
:Concentration
:Mass diffusivity
:Chemical reaction rate
:Nondimensional chemical reaction parameter,
Sc:Modified Schmidt number, .
Greek Letters
Dimensionless coordinate,
:Gyroviscosity parameter
:Dimensionless
:Fluid density
:Fluid viscosity
:Magnetic permeability
:Electric conductivity
:Nondimensional frequency parameter,
:Ion slip parameter
:Hall parameter
:Hall and ion slip parameters, .

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors are thankful to the referees and Professor J. V. Ramanamurthy, NIT Warangal, for their valuable suggestions which led to the improvement of the paper.