Influence of Chemical Reaction on Heat and Mass Transfer Flow of a Micropolar Fluid over a Permeable Channel with Radiation and Heat Generation

The effects of chemical reaction on heat andmass transfer flowof amicropolar fluid in a permeable channel with heat generation and thermal radiation is studied.The Rosseland approximations are used to describe the radiative heat flux in the energy equation.The model contains nonlinear coupled partial differential equations which have been transformed into ordinary differential equation by using the similarity variables. The relevant nonlinear equations have been solved by Runge-Kutta-Fehlberg fourth fifth-order method with shooting technique. The physical significance of interesting parameters on the flow and heat transfer characteristics as well as the local skin friction coefficient, wall couple stress, and the heat transfer rate are thoroughly examined.


Introduction
The micropolar fluid theory is the one of the most important non-Newtonian fluid models described by Eringen [1].This theory shows microrotation effects as well as microinertia and has many applications such as polymer fluids, liquid crystals, animal bloods, unusual lubricants, colloidal and suspension solutions, colloidal fluids, liquid crystals, and polymeric suspension.The extensive reviews of the micropolar fluid theory and its applications can be found in Eringen [2] and Lukaszewicz [3].The unsteady mixed convection flow of a micropolar fluid from a vertical surface in the presence of viscous dissipation and the buoyancy force has been studied by El-Aziz [4].Ashraf et al. [5] investigated the micropolar fluids flow through a porous channel.Sheikholeslami et al. [6] proposed the flow and heat transfer of micropolar fluid in a permeable channel.Prakash and Muthtamilselvan [7] investigated radiation effect on MHD micropolar fluid flow in porous vertical channel.Darvishi et al. [8] analyzed numerically the micropolar fluid flow from a porous channel.Sherief et al. [9] studied the motion along axis of a circular cylindrical pore in a micropolar fluid of a slip spherical particle.Mosayebidorcheh [10] analyzed the flow of micropolar fluid over a porous channel with changing walls.Recently, many authors [11][12][13][14] have studied the micropolar fluid flow for different fluid properties over different geometries.However, the effects of chemical reaction on micropolar fluid flow over a permeable channel in the presence of radiation and heat generation have not been considered in the above investigations.
Combined heat and mass transfer flows in the presence of chemical reaction have numerous applications in engineering and geophysics such as drying, geothermal reservoirs, dehydration at the surface of a water body, drying of porous solids, geothermal pool, thermal insulation, enhanced oil recovery, cooling of nuclear reactors, fibrous insulation, evaporation at the surface of a water body, pollution studies, cooling the polymer production and manufacturing of ceramics, energy transfer in a wet cooling tower, the flow in a desert, and oxidation and synthesis materials.Mohamed and Abo-Dahab [15] studied the effects of chemical reaction and thermal radiation on hydromagnetic free convection heat and mass transfer for a micropolar fluid bounded by a semi-infinite vertical porous plate in the presence of heat generation.In recent years, many researchers have studied and reported the effect of first-order chemical reaction [16][17][18][19][20][21][22][23][24][25].
Influence of thermal radiation on flow and heat transfer study has become more important industrially.The heat transfer and temperature profile of a micropolar fluid over different geometries can be affected significantly at high temperature.Bhattacharyya et al. [26] considered thermal radiation effect on micropolar fluid flow and heat transfer over a porous shrinking sheet.Hussain et al. [27] analyzed radiation effects on the thermal boundary layer flow of a micropolar fluid towards a permeable stretching sheet.Oahimire and Olajuwon [28] investigated the influence of Hall current and thermal radiation on heat and mass transfer of a chemically reacting MHD flow of a micropolar fluid through a porous medium.Mabood et al. [29] studied effects of nonuniform heat source/sink and Soret on MHD non-Darcian convective flow past a stretching sheet in a micropolar fluid with radiation.
The effect of heat generation on heat transfer is an important issue in view of various physical problems.Ziabakhsh et al. [30] analyzed the micropolar fluid flow with heat generation.Singh and Kumar [31] considered the melting effect in stagnation-point flow of micropolar fluid towards a stretching/shrinking surface.Bakr [32] investigated the effects of chemical reaction and heat source magnetoconvection and mass transfer flow of a micropolar fluid in a rotating frame of reference.The heat generation/absorption effects on MHD flow and heat transfer of micropolar fluid through a stretching surface have been proposed by Mahmoud and Waheed [33].Abbasi et al. [34] examined the flow of Maxwell nanofluid in the presence of heat generation/absorption.Mliki et al. [35] investigated the influence of nanoparticle Brownian motion and heat generation/absorption over linear/sinusoidally heated cavity in the presence of magnetohydrodynamic natural convection.Sheikholeslami and Ganji [36] studied three-dimensional heat and mass transfer flow of nanofluid over a rotating system.Thermal radiation effects on mixed convection flow and heat transfer of a micropolar fluid through an unsteady stretching surface with heat generation/absorption are presented by Singh and Kumar [37].
Motivated by the above studies and applications, the present work explores effects of chemical reaction on heat and mass transfer flow of a micropolar fluid over a permeable channel in the presence of radiation and heat generation.The equations of continuity, momentum, angular momentum, energy, and concentration have been reduced to a system of nonlinear ordinary differential equations by similarity transforms which are solved by Runge-Kutta-Fehlberg method with shooting technique.It is expected that the results obtained from present paper will provide important information to the audience.To the best of our knowledge, such type of study is not investigated before in the scientific literature.

Mathematical Formulation
The heat and mass transfer flow of a micropolar fluid in a permeable channel with chemical reaction is considered in the present work.The thermal radiation and heat source are incorporated in the energy equation.The graphical model of Micropolar fluid flow between permeable walls the problem has been given along with flow configuration and coordinate system in Figure 1.The assumptions of the problem in detail can be found in [6].The governing equations of boundary layer are given in the following form: where  and V indicate the velocity components in the  and  directions, respectively,  is the fluid density,  is the dynamic viscosity,  is the material parameter,  is the angular or microrotation velocity,  is the fluid pressure,  is the microinertia density,   = ( + /2) is the microrotation viscosity,  is the fluid temperature,   is the specific heat at constant pressure,  is the fluid concentration,  is the thermal conductivity,   is the radiative heat flux,  0 is the heat generation coefficient,  * is the molecular diffusivity, and  0 is the chemical reaction rate coefficient.Using Rosseland's approximation for radiation, we obtain where  is the Stefan-Boltzmann constant and  0 is the absorption coefficient.We consider that the temperature variation within the flow is such that  4 may be expanded in a Taylor's series.Expanding  4 about  ∞ and neglecting higher order terms, we get The appropriate boundary conditions for the flow are where boundary parameters  (0 ≤  ≤ 1) indicate the degree to which the microelements are free to rotate near the channel walls.The case when  = 0 is called strong concentration of microelements, which implies  = 0 near the wall surface.This represents concentrated particle flow where the microelements close to the wall surface are unable to rotate.In the case when  = 1/2, this indicates vanishing of the antisymmetric part of the stress tensor and denoted weak concentration of microelement and  = 1 is used for the modeling of turbulent boundary layer flow.In this paper, the authors considered  = 1/2 for which the governing equations can be reduced to the classical problem of steady boundary layer flow of a viscous incompressible fluid near the channel wall.Equations ( 2), ( 3), ( 4), (6), and ( 8) can be transformed into a set of nonlinear ordinary differential equations by using the following similarity transformations: where  2 =  1 −  and  2 =  1 − , with  and  as constants.The stream function  is defined as The coupled system of transformed nonlinear ordinary differential equations is Boundary conditions in nondimensional form are where  1 = / is the coupling number,  2 = V  /ℎ 2 is the spin-gradient viscosity parameter,  3 = /ℎ 2 is the micropolar material constant, Ec = V 2 0 /ℎ 3    is the local Eckert number,  =  0 ℎ/  V 0 is the heat generation parameter, Pe ℎ = Pr Re and Pe  = Sc Re are the Peclet numbers for the diffusion of heat and the diffusion of mass, Re = (V 0 /)ℎ is the Reynolds number,  = 3 0 /16 3 ∞ is the thermal radiation parameter, Gr =   ℎ 4 / 2 is the Grashof number, Pr =   / is the Prandtl number, Sc = / * is the Schmidt number, and  =  0 ℎ/V 0 is the chemical reaction parameter.
The other parameters of physical interest are the local Nusselt Nu  and Sherwood Sh  numbers, which are defined as follows: , where   and   are the local heat flux and mass flux, respectively, which are defined as Now using ( 10) and ( 15) in ( 14), we get Nu  = −  (−1) ,

Method of Solution
In this present paper, Runge-Kutta-Fehlberg fourth fifthorder method has been employed to solve the system of nonlinear ordinary differential equations (12) with the boundary conditions given by ( 13) for different values of governing parameters.The RKF45 method has a procedure to determine if the appropriate step size ℎ is being used.The formula of fifth-order Runge-Kutta-Fehlberg method can be defined as follows: The computation of the error can be achieved by subtracting the fifth-order from the fourth-order method: If the error goes beyond a specified antechamber, the results can be recalculated using a smaller step size.The approach to computing the new step size is shown as follows: The numerical computations were carried out with Δ = 0.01.The variation of the dimensionless velocity, microrotation, temperature, and concentration is ensured to be less than 10 −6 between any two successive iterations for the convergence criterion.

Results and Discussion
In order to study, the effects of various governing physical parameters on the flow, heat, and mass transfer numerical computations are carried out for 0.
, and 0 ≤  ≤ 10 while the Eckert number Ec = 0.01 is fixed.A critical analysis with previously published work is done in Table 1 and these results are found to be in very good agreement.
The effects of the Reynolds number Re on velocity   () and microrotation () are shown in Figures 2 and 3.It is seen from Figure 2 that velocity profile   () decreases near the lower channel wall, while it increases near the upper channel wall when the Reynolds number Re increases.It is clear from Figure 3 that an increase in the magnitude of Re leads to decrease in microrotation profile ().The microrotation profile () for the coupling number  1 , the spin-gradient viscosity parameter  2 , and the micropolar material constant  3 are displayed in Figures 4-6.It is noticed from Figure 4 that the microrotation profile increases with  1 at lower channel wall for 0 ≤  < −0.65, but reverse trend is found for the case  > −0.65.It is depicted from Figure 5 that the value of microrotation profile () increases as the value of  2 increases.It is evident from Figure 6 that the values of microrotation are lower for higher values of  3 .
Figures 7-9 show the temperature distribution () with collective variation in heat source parameter , thermal radiation parameter , and Peclet number Pe ℎ .Figure 7 exhibits that temperature () considerably increases with an increase in . Figure 8 indicates that an increase in  induces a decrease in temperature ().It is noticed from Figure 9 that a rise in Pe ℎ causes rapid increase in temperature.
The effects of the Peclet number Pe  and chemical reaction parameter  on the concentration profile () are shown in Figures 10 and 11.Figures 10 and 11 indicate that concentration decreases with increasing values of Pe  and .
The variation of the Nusselt number (dimensionless heat transfer rate at the surface) is displayed for different parameters in Figures 12-14.Figure 12 depicts the behavior of the Nusselt's number against Reynolds number Re with various values of heat source parameter .It is clear that, with the increasing of , the heat transfer rate decreases.Figure 13 enlightens the variation of the heat transfer rates with Reynolds number Re for various values of thermal radiation parameter .It is evident from Figure 13 that Nusselt number monotonically increases with a rise in  values.Figure 14

Conclusions
The present paper deals with numerical analysis of chemical reaction effects on heat and mass transfer flow of a micropolar fluid over a permeable channel in the presence of radiation and heat generation.The system of nonlinear partial differential equations was converted to a system of ordinary differential equations and then is solved numerically using the Runge-Kutta-Fehlberg method along with shooting method.From the above discussion, the important results are summarized as follows: velocity and Reynolds number are proportional to each other at upper channel wall.
(ii) Microrotation decreases with increase in the value of coupling number, micropolar material constant, and Reynolds number, but it increases with increase in the value of spin-gradient viscosity parameter.
(iii) Temperature increases with heat generation parameter and Pe ℎ and Pe  are the Peclet numbers for the diffusion of heat and the diffusion of mass, and temperature decreases with thermal radiation parameter.
(iv) Concentration is inversely proportional to the Peclet number for the diffusion of mass and chemical reaction parameter.

Figure 10 :
Figure 10: Concentration profile () for various values of Pe  .
is noted that heat transfer rate decreases with a rise in Pe ℎ values.Also, a negligible change in heat transfer rate with increasing values of Reynolds number Re is seen from Figures 12-14.Figures 15 and 16 displayed the variation of the Sherwood number or mass transfer rate for different parameters.Figures 15 and 16 depict the behavior of the mass transfer rate against Reynolds number Re with different values of Peclet number Pe  and chemical reaction parameter .The mass transfer rate rises with increasing values of the Peclet number Pe  and chemical reaction parameter .It is also clear from these figures that there is no change in mass transfer rate with rise in Reynolds number Re values.