Abstract
The objective of the present study is to investigate the first-order chemical reaction and Soret and Dufour effects on an incompressible MHD combined free and forced convection heat and mass transfer of a micropolar fluid through a porous medium between two parallel plates. Assume that there are a periodic injection and suction at the lower and upper plates. The nonuniform temperature and concentration of the plates are assumed to be varying periodically with time. A suitable similarity transformation is used to reduce the governing partial differential equations into nonlinear ordinary differential equations and then solved numerically by the quasilinearization method. The fluid flow and heat and mass transfer characteristics for various parameters are analyzed in detail and shown in the form of graphs. It is observed that the concentration of the fluid decreases whereas the temperature of the fluid enhances with the increasing of chemical reaction and Soret and Dufour parameters.
1. Introduction
The flow through porous boundaries has many applications in science and technology such as water waves over a shallow beach, mechanics of the cochlea in the human ear, aerodynamic heating, flow of blood in the arteries, and petroleum industry. Several authors have studied theoretically the laminar flow in porous channels. Berman [1] considered the viscous fluid and analyzed the flow characteristics when it passed through the porous walls. Later the same problem for different permeability was studied by Terril and Shrestha [2]. The theory of micropolar fluids was introduced by Eringen [3] which are considered as an extension of generalized viscous fluids with microstructure. Examples for micropolar fluids include lubricants, colloidal suspensions, porous rocks, aerogels, polymer blends, and microemulsions. The same Berman problem with micropolar fluid was discussed by Sastry and Rama Mohan Rao [4]. The flow and heat transfer of micropolar fluid between two porous parallel plates was analyzed by Ojjela and Naresh Kumar [5]. Srinivasacharya et al. [6] obtained an analytical solution for the unsteady Stokes flow of micropolar fluid between two parallel plates. The effect of buoyancy parameter on flow and heat transfer of micropolar fluid between two vertical parallel plates was investigated by Maiti [7].
The study of MHD heat and mass transfer through porous boundaries has attracted many researchers in the recent past due to applications in engineering and science, such as oil exploration, boundary layer control, and MHD power generators. The steady incompressible free convection flow and heat transfer of an electrically conducting micropolar fluid in a vertical channel was studied by Bhargava et al. [8]. The laminar incompressible magnetohydrodynamic flow and heat transfer of micropolar fluid between porous disks was analyzed numerically by Ashraf and Wehgal [9]. Islam et al. [10] obtained a numerical solution for an incompressible unsteady magnetohydrodynamic flow through vertical porous medium. Nadeem et al. [11] discussed the unsteady MHD stagnation flow of a micropolar fluid through porous media. The effects of Hall and ion slip currents on micropolar fluid flow and heat and mass transfer in a porous medium between parallel plates with chemical reaction were considered by Ojjela and Naresh Kumar [12]. The MHD heat and mass transfer of micropolar fluid in a porous medium with chemical reaction and Hall and ion slip effects by considering variable viscosity and thermal diffusivity were investigated by Elgazery [13]. The mixed convection flow and heat transfer of an electrically conducting micropolar fluid over a vertical plate with Hall and ion slip effects was analyzed by Ayano [14].
When heat and mass transfer occurs simultaneously in a moving fluid, the energy flux caused by a concentration gradient is termed as diffusion thermoeffect, whereas mass fluxes can also be created by temperature gradients which is known as a thermal diffusion effect. These effects are studied as second-order phenomena and may have significant applications in areas like petrology, hydrology, and geosciences. The effect of thermophoresis on an unsteady natural convection flow and heat and mass transfer of micropolar fluid with Soret and Dufour effects was studied by Aurangzaib et al. [15]. Srinivasacharya and RamReddy [16] considered the problem of the steady MHD mixed convection heat and mass transfer of micropolar fluid through non-Darcy porous medium over a semi-infinite vertical plate with Soret and Dufour effects. Influence of the Soret and Dufour numbers on mixed convection flow and heat and mass transfer of non-Newtonian fluid in a porous medium over a vertical plate was analyzed by Mahdy [17]. Hayat and Nawaz [18] investigated analytically the effects of the Hall and ion slip on the mixed convection heat and mass transfer of second-grade fluid with Soret and Dufour effects. Rani and Kim [19] studied numerically the laminar flow of an incompressible viscous fluid past an isothermal vertical cylinder with Soret and Dufour effects. The effects of chemical reaction and Soret and Dufour on the mixed convection heat and mass transfer of viscous fluid over a stretching surface in the presence of thermal radiation were analyzed by Pal and Mondal [20]. Sharma et al. [21] studied the mixed convective flow, heat and mass transfer of viscous fluid in a porous medium past a radiative vertical plate with chemical reaction, and Soret and Dufour effects.
In the field of fluid mechanics many fluid flow problems are nonlinear boundary value problems. To solve these problems we can use a numerical technique, quasilinearization method which is a powerful technique having second-order convergence. Several authors (Lee and Fan [22], Hymavathi and Shanker [23], Huang [24], Motsa et al. [25], and Ojjela and Naresh Kumar [5, 12]) applied the quasilinearization method to solve the nonlinear boundary layer equations.
In the present study the effects of chemical reaction on two-dimensional mixed convection flow and heat transfer of an electrically conducting micropolar fluid in a porous medium between two parallel plates with Soret and Dufour have been considered. The reduced flow field equations are solved using the quasilinearization method. The effects of various parameters such as Hartmann number, inverse Darcy’s parameter, Schmidt number, Prandtl number, chemical reaction rate, Soret and Dufour numbers on the velocity components, microrotation, temperature distribution, and concentration are studied in detail and presented in the form of graphs.
2. Formulation of the Problem
Consider a two-dimensional laminar incompressible micropolar fluid flow through an elongated rectangular channel, as shown in Figure 1. Assume that the fluid is injected and aspirated periodically through the plates with injection velocity and suction velocity . Also the nonuniform temperature and concentration at the lower and upper plates are , and , , respectively. The region inside the parallel plates is subjected to porous medium and a constant external magnetic field of strength perpendicular to the -plane is considered.
The governing equations of the micropolar fluid flow and heat and mass transfer in the presence of buoyancy forces, magnetic field and in the absence of body forces, body couples are given by where is the buoyancy force and it is defined as .
Neglecting the displacement currents, the Maxwell equations and the generalized Ohm’s law arewhere , is induced magnetic field. Assume that the induced magnetic field is negligible compared to the applied magnetic field so that magnetic Reynolds number is small, the electric field is zero, and magnetic permeability is constant throughout the flow field.
The velocity and microrotation components areFollowing Ojjela and Naresh Kumar [5, 12] the velocity and microrotation components are The temperature and concentration distributions can be taken aswhere and , , , , , and are to be determined.
The boundary conditions for the velocity, microrotation, temperature, and concentration areSubstituting (8) and (9) in (2), (3), (4), and (5) then we getwhere the prime denotes the differentiation with respect to λ.
The dimensionless forms of temperature and concentration from (9) are The boundary conditions (10) in terms of , , , , , and areFor micropolar fluids the shear stress is given byThen the nondimensional shear stress at the lower and upper plates is For micropolar fluids, the couple stress is given byThen the nondimensional couple stress at lower and upper plates is
3. Solution of the Problem
The nonlinear equations (11) are converted into the following system of first-order differential equations by the substitutionThe boundary conditions in terms of , , , , , , , , , , , , areThe system of (18) is solved numerically subject to boundary conditions (19) using the quasilinearization method given by Bellman and Kalaba [26].
Let (, ) be an approximate current solution and let (, ) be an improved solution of (18). Using Taylor’s series expansion about the current solution by neglecting the second- and higher-order derivative terms, coupled first-order system (18) is linearized asTo solve for (, ), the solutions to seven separate initial value problems, denoted by , , , , , , and (which are the solutions of the homogeneous system corresponding to (20)) and (which is the particular solution of (20)), with the following initial conditions are obtained by using the 4th-order Runge-Kutta method:By using the principle of superposition, the general solution can be written aswhere , , , , , , and 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 Discussion
The system of nonlinear differential equations (18) subject to boundary conditions (19) is solved numerically by the quasilinearization method. The influences of various fluid and geometric parameters such as Soret number Sr, Dufour number Du, Hartmann number Ha, chemical reaction parameter Kr, Schmidt number Sc, Eckert number Ec, and Prandtl number Pr on nondimensional velocity components, microrotation, temperature distribution, and concentration are analyzed through graphs in the domain .
Figures 2 and 3 show the influence of Sr and Du on temperature and concentration. From these figures, it is evident that the temperature of the fluid increases whereas the concentration decreases with the increasing of Sr and Du. This is because of the difference between the temperatures of the fluid and surface as well as the difference between concentrations of the fluid and surface concentrations are increased with Sr and Du. Figure 4 describes the behavior of the temperature distribution and concentration for the various values of Kr. As Kr increases the temperature distribution of the fluid also increases, whereas the concentration decreases from the lower plate to the upper plate. It is clear that the increase in the Kr produces a decrease in the species concentration. This causes the concentration buoyancy effects to decrease as Kr increases. The effect of Ec on velocity components, microrotation, and temperature is presented in Figure 5. As Ec increases the radial velocity, microrotation and temperature are decreasing towards the upper plate. However, the axial velocity decreases towards the center of the channel and then increases. Since the Eckert number is the relation between the kinetic energy and enthalpy, as enthalpy increases, the temperature distribution decreases. Figure 6 elucidates the change in velocity components, microrotation, temperature distribution, and concentration for different values of Pr. It is observed that the axial velocity reaches highest value near the hot plate and the radial velocity, microrotation, and concentration decrease whereas the temperature increases with the increasing value of Pr. Physically, if Pr increases the thermal diffusivity decreases and this leads to the decrease in the heat transfer ability at the thermal boundary layer. Figure 7 displays the change in the velocity components, microrotation, temperature distribution, and concentration for several values of Ha. From this it is observed that when Ha increases, the temperature distribution also increases whereas the concentration decreases from the lower plate to upper plate and the axial velocity attains maximum value at the center of the plates. However, the radial velocity and microrotation increase towards the center of the plates and then decrease. This is due to the fact that the magnetic force retards the flow in both axial and radial directions. The variations in the velocity components, microrotation, temperature distribution, and concentration for different values of are shown in Figure 8. From these one can deduce that the temperature distribution is increasing with whereas the radial velocity, microrotation, and concentration are decreasing towards the upper plate and the axial velocity decreases towards the center of the plates and then increases because the resistance offered by the porosity of the medium is more than the resistance due to the magnetic lines of force.
(a)
(b)
(a)
(b)
(a)
(b)
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
(e)
(a)
(b)
(c)
(d)
(e)
(a)
(b)
(c)
(d)
(e)
5. Conclusions
The thermal diffusion and diffusion thermoeffects on combined free and forced convection magnetohydrodynamic flow of micropolar fluid in a porous medium between two parallel plates with chemical reaction are considered. The numerical solution of the transformed governing equations is obtained by the method of quasilinearization and the results are analyzed for various fluid and geometric parameters through graphs. From the results the following is concluded:(i)The influences of Sr and Du on temperature and concentration are similar.(ii)The temperature of the fluid is enhanced whereas the concentration of the fluid is decreased with the increasing of Ha and .(iii)Kr reduces the concentration and enhances the temperature of the fluid.(iv)Ec and Pr exhibit similar effects on the velocity components and microrotation whereas it is opposite in the case of temperature.
Nomenclature
| : | Injection suction ratio, |
| : | Time |
| : | Distance between two parallel plates |
| : | Injection velocity |
| : | Suction velocity |
| : | Fluid pressure |
| : | Velocity vector |
| : | Specific heat at constant temperature |
| : | Microrotation vector |
| : | Microrotation component |
| : | Eckert number, |
| : | Thermal conductivity |
| : | Viscosity parameter |
| : | Permeability of the medium |
| : | Chemical reaction rate |
| : | Velocity component in -direction |
| : | Velocity component in -direction |
| : | Prandtl number, |
| : | 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 |
| : | Hartmann number, |
| : | Inverse Darcy parameter, |
| : | Nondimensional viscosity parameter, |
| : | Nondimensional micropolar parameter, |
| : | Nondimensional micropolar parameter, |
| : | Temperature |
| : | Temperature of the lower plate |
| : | Temperature of the upper plate |
| : | Dimensionless temperature, |
| : | Concentration |
| : | Nondimensional concentration, |
| : | Concentration of the lower plate |
| : | Concentration of the upper plate |
| : | Mass diffusivity |
| : | Nondimensional chemical reaction parameter, |
| : | Schmidt number, |
| : | Thermal Grashof number, |
| : | Solutal Grashof number, |
| : | Sherwood number, |
| : | Mass transfer rate |
| : | Soret number, |
| : | Dufour number, |
| : | Mean temperature |
| : | Thermal diffusion ratio |
| : | Concentration susceptibility. |
| : | Dimensionless coordinate, |
| : | Gyro viscosity parameters |
| : | Dimensionless axial variable, |
| : | Kinematic viscosity |
| : | Fluid density |
| : | Fluid viscosity |
| : | Magnetic permeability |
| : | Electric conductivity |
| : | Nondimensional frequency parameter, . |
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The authors thank referees for their suggestions which have resulted in the improvement of the paper. Also one of the authors (N. Naresh Kumar) is grateful to the Defence Research and Development Organization, Government of India, for providing senior research fellowship.