Abstract
The effects of a homogeneousheterogeneous reaction on steady micropolar fluid flow from a permeable stretching or shrinking sheet in a porous medium are numerically investigated in this paper. The model developed by Chaudhary and Merkin (Fluid Dyn. Res. 16:311333, 1995) for a homogeneousheterogeneous reaction in boundary layer flow with equal diffusivities for reactant and autocatalysis is used and extended in this study. The uniqueness of this problem lies in the fact that the solutions are possible for all values of the stretching parameter , while for (shrinking surface), solutions are possible only for a limited range of values. The effects of physical and fluid parameters such as the stretching parameter, micropolar parameter, permeability parameter, Schmidt number, strength of homogeneous and heterogeneous reaction parameter on the skin friction, velocity and concentration are analyzed, and these results are presented through graphs. The solute concentration at the surface is found to decrease with the strength of the homogeneous reaction, and to increase with heterogeneous reactions, the permeability parameter and stretching or shrinking parameters. The velocity at the surface was found to increase with the micropolar parameter.
Keywords:
micropolar fluid; permeable stretching; shrinking sheet; homogeneousheterogeneous reactions; porous medium1 Introduction
Micropolar fluids are fluids with internal structures in which coupling between the spin of each particle and the microscope velocity field is taken into account. They represent fluids consisting of rigid, randomly oriented or spherical particles suspended in a viscous medium, where the deformation of fluid particles is ignored. Micropolar fluid theory was introduced by Eringen [1] in order to describe physical systems, which do not satisfy the NavierStokes equations. The equations governing the micropolar fluid involve a spin vector and a microinertia tensor in addition to the velocity vector. The potential importance of micropolar fluids in industrial applications has motivated many researchers to extend the study in numerous ways to include various physical effects. The essence of the theory of micropolar fluid lies in particle suspension (Hudimoto and Tokuoka [2]), liquid crystals (Lockwood et al.[3]); animal blood (Ariman et al.[4]), exotic lubricants (Erigen [5]), etc. An excellent review of the various applications of micropolar fluid mechanics was presented by Ariman et al.[6].
Boundary layer flow over a stretching surface is important as it occurs in several engineering processes, for example, materials manufactured by extrusion. During the manufacturing process, a stretching sheet interacts with the ambient fluid both thermally and mechanically. The study of boundary layer flow caused by a stretching surface was initiated by Crane [7]. Recently, several works on the dynamic of the boundary layer flow over a stretching surface have appeared in literature (Dutta et al.[8], Hayat et al.[9], Ishak [10]). The effect of surface conditions on the micropolar flow driven by a porous stretching sheet was studied by Kelson and Desseaux [11]. Mohammadein and Gorla [12] examined the flow of micropolar fluids bounded by a stretching sheet with prescribed wall heat flux, viscous dissipation and internal heat generation. The effect of suction or injection at a stretching surface was studied by Erickson et al.[13] and Fox et al.[14]. The process of suction is used in many engineering activities such as thermal oil recovery, removal of reactants etc. Elbashbeshy and Bazid [15] studied the flow and heat transfer in a porous medium over a stretching surface. Bhargava et al.[16] investigated the flow of a mixed convection micropolar fluid driven by a porous stretching sheet with uniform section. Later, Bhargava et al.[17] studied the same flow of a micropolar flow over a nonlinear stretching sheet. Abel et al.[18] carried out a numerical study of hydromagnetic micropolar fluid flow due to horizontal/vertical stretching sheet using a shooting method. They highlighted a scientific approach for the choice of the missing initial values on which the convergence of the shooting method highly depends. Recently, Narayana and Sibanda [19] studied the effects of laminar flow of a nanoliquid film over an unsteady stretching sheet. Kameswaran et al.[20] studied hydromagnetic nanofluid flow due to a stretching or shrinking sheet with viscous dissipation and chemical reaction effects. Recently, Kameswaran et al.[21] studied homogeneousheterogeneous reactions in a nanofluid flow over a permeable stretching sheet.
Many chemically reacting systems involve both homogeneous and heterogeneous reactions, with examples occurring in combustion, catalysis and biochemical systems. The interaction between the homogeneous reactions in the bulk of fluid and heterogeneous reactions occurring on some catalytic surfaces is generally very complex, involving the production and consumption of reactant species at different rates both within the fluid and on the catalytic surfaces. A simple mathematical model for homogeneousheterogeneous reactions in stagnationpoint boundarylayer flow was initiated by Chaudhary and Merkin [22]. They modeled the homogeneous (bulk) reaction by isothermal cubic kinetics and the heterogeneous (surface) reaction was assumed to have firstorder kinetics. Later Chaudhary and Merkin [23] extended their previous work to include the effect of loss of the autocatalyst. They studied the numerical solution near the leading edge of a flat plate. A model for isothermal homogeneousheterogeneous reactions in boundary layer flow of a viscous fluid flow past a flat plate was studied by Merkin [24]. Ziabakhsh et al.[25] studied the problem of flow and diffusion of chemically reactive species over a nonlinearly stretching sheet immersed in a porous medium. Chambre and Acrivos [26] studied an isothermal chemical reaction on a catalytic in a laminar boundary layer flow. They found the actual surface concentration without introducing unnecessary assumptions related to the reaction mechanism. The effects of flow near the twodimensional stagnation point flow on an infinite permeable wall with a homogeneousheterogeneous reaction was studied by Khan and Pop [27]. They solved the governing nonlinear equations using the implicit finite difference method. It was observed that the mass transfer parameter considerably affects the flow characteristics. Khan and Pop [28] studied the effects of homogeneousheterogeneous reactions on the viscoelastic fluid toward a stretching sheet. They observed that the concentration at the surface decreased with an increase in the viscoelastic parameter.
The purpose of the present study is to analyze the influence of the permeability, the homogeneous and heterogeneous reaction on the micropolar fluid towards a stretching/shrinking sheet. We transformed the governing momentum and concentration equations into a system of ordinary differential equations using a similarity variable and then numerically solved the equations for some values of the governing parameters. To the best of authors knowledge, such study has not been reported earlier in the literature.
2 Mathematical formulation
Consider steady, incompressible twodimensional boundary layer flow of a micropolar fluid through a porous medium. The Cartesian coordinates x and y are taken along the surface and are normal to it, respectively, and u and v are the respective velocity components. The flow is generated due to stretching or shrinking of the sheet caused by the simultaneous application of two equal forces along the xaxis. Keeping the origin fixed, it is assumed that the surface is stretched/shrunk with a linear velocity , where is a constant with for a stretching sheet, for a shrinking sheet and for a static sheet. We consider a simple model for the interaction between a homogeneous (or bulk) reaction and a heterogeneous (or surface) reaction involving the two chemical species A and B in a boundary layer flow proposed by Merkin [24] and Chaudhary and Merkin [22,23] in the following form:
Here a and b are concentrations of chemical species A and B, and () are the rate constants. We also assume that both reaction processes are isothermal. It is also assumed that the ambient fluid moves with a velocity , where is a constant, in which there is a uniform concentration of reactant A and in which there is no autocatalyst B over a flat surface.
Under these assumptions, the governing equations lead to
where is the permeability of the porous medium, ϵ is the porosity, μ is the dynamic viscosity, is the effective dynamic viscosity, κ is the vortex viscosity, ρ is the fluid density, N is the microrotation, is the microinertia per unit mass and γ is the spin gradient viscosity defined as
where is the material or micropolar parameter. The corresponding boundary conditions are
where is the constant mass flux with for suction and for injection (blowing), respectively; n is a constant such that , where the case is called strong concentration (Guram and Smith [29]), indicates near the surface and represents concentrated particle flows in which the microelements close to the surface are unable to rotate (Jena and Mathur [30]). The case is used for the modeling of turbulent boundary layer flow (Peddieson [31]). The case indicates the vanishing of the antisymmetrical part of the stress tensor and denotes weak concentration (Ahmadi [32]), which is the case considered in the present study. We introduce the stream function ψ, where and , which satisfies the continuity equation identically. Using similarity variables of the form
the governing equations are written as the following system of ordinary differential equations:
where is the permeability parameter, is the Schmidt number, gives a measure of the strength of the homogeneous reaction, is the ratio of the diffusion coefficient and primes denote differentiation with respect to η. The boundary conditions (equation (2.9)) are written as
where is the mass transfer parameter with for suction, for injection and for an impermeable surface. is the stretching parameter where corresponds to a stretching surface and corresponds to a shrinking surface. measures the strength of the heterogeneous (surface) reaction and is the Reynolds number.
It is expected that the diffusion coefficients of chemical species A and B are of comparable size, which leads us to further assumption that the diffusion coefficients and are equal, i.e., (Chaudhary and Merkin [22]). This assumption leads to the following relation:
Thus, equations (2.13) and (2.14) reduce to
and the boundary conditions equation (2.16) take the form
The physical quantity of interest is the skin friction coefficient . It characterizes the surface drag. The shearing stress at the surface of the wall is given by
The skin friction coefficient is defined as
Using the similarity variables in the above equation, we obtain
where represents the local Reynolds number defined by . In the present paper, we consider only the case (suction) and (weak concentration).
3 Results and discussion
The system of ordinary differential equations (2.10), (2.11) and (2.17) along with the boundary conditions (2.15) and (2.19) are solved numerically for some values of λ, , , k, χ and Sc using Matlab bvp4c routine. We compare our results for a stretching sheet with those reported by Wang [33], Ishak et al.[34] and Rosali et al.[35] in Table 1; and for the shrinking sheet, we compare our results with those reported by Rosali et al.[35] in Table 2. These show an excellent agreement with our present results. This investigation confirms that the existence and uniqueness of solution depends on the stretching/shrinking sheet parameter. It is also found that the solutions of equations (2.10), (2.11) and (2.17) can be obtained for all values of λ, while for the case of a shrinking surface (), the governing equations have the solution only in the range of , where is a critical value of λ, which depends on the other parameter, and we have no solution for . A dual solution appears for the range . We have plotted the velocity profile, concentration profile and concentration at the surface for different parameters. In this problem, we have shown single as well as dual solution. As was mentioned by Merkin [36], Merrill et al.[37] and Rosali et al.[35], the first solution is stable and physically realizable, while the second solution is unstable. Although the second solution seems to be deprived of physical significance, it is interesting in nonlinear differential equation theory, since a similar equation may reappear in some other situations where the corresponding solution could have a more realistic meaning.
Table 1 . Comparison offor a stretching sheet obtained for different values ofλ, for fixed values of,and
Table 2 . Comparison offor a shrinking sheet obtained for different values ofλ, for fixed values of,and
From Table 1, it is clear that the skin friction is a decreasing function of λ. All values of the skin friction coefficient are positive for , while they are negative when . Physically, the negative values of the skin friction coefficient correspond to the surface exerting a drag force on the fluid and the opposite sign implies the inverse phenomenon. The skin friction coefficient is zero when regardless of the values of other parameters. This is because for , there is no shear stress at the surface as the surface and fluid move with the same velocity.
The effect of the stretching/shrinking parameter λ for is shown in Table 2. It is evident that initially the skin friction is an increasing function of λ, but it decreases after a certain value of λ. As mentioned earlier, the solution of the equations is possible only in the range of . It is evident that the value of depends on other parameters , s and n. The value of is approximate to −1.2465 when and , while introducing and , we get . It is also observed that the first solution is a decreasing function of λ, whereas the second solution is an increasing function of λ. It is interested to note that these both solutions coincide at .
The variation of the velocity and concentration profiles is plotted as a function of η for some values of λ in Figure 1. (i) For (stretching surface), the fluid velocity is becoming increasingly greater than the free stream. In this case the fluid velocity decreases with the value of η and converges at unity as per the condition. (ii) For (static surface), the fluid velocity initially is stationary, but with η value it increases in a nonlinear way. (iii) For (shrinking surface), the fluid velocity is initially negative, but it increases with η, and after a certain value of η, it becomes positive. For the concentration profile, all the curves are started from the origin and they increase nonlinearly with η to follow ‘S’ shape and finally reach unity according to the given condition. For and −1.3, the graphs contain a dual solution. From the boundary condition, it is clear that the velocity at the surface is equal to λ; and so with increase in λ, the momentum boundary layer thickness increases for the case of the first solution, while an opposite phenomenon appears for the second solution which concurs with the results reported by Bhattacharyya et al.[38]. It shows that for the first solution, the concentration boundary layer thickness increases with increase in λ, while it decreases for the second solution.
Figure 1 . Effect ofλon (a) and (b) for,,,,,,.
The effect of the micropolar parameter and suction parameter on the velocity and concentration profile is shown in Figure 2. A comparison is made for the Newtonian fluid () and the micropolar fluid (). It shows that momentum boundary layer thickness is more for the micropolar fluid rather than for the Newtonian one, i.e., the micropolar parameter increases the velocity of the fluid, which is similar to the results given by Ishak et al.[34]. A similar phenomenon is observed for concentration boundary thickness. The velocity of the fluid is reduced due to suction, and this leads to an increase in the solute concentration.
Figure 2 . Variation ofandson (a) and (b) for,,,,,.
The effect of the homogeneous and heterogeneous reaction on the concentration profile is shown in Figure 3. The effect of heterogeneous and homogeneous reactions separately is shown through Figures 3(a) and 3(b), respectively. We considered and so a dual solution appeared. It is evident that the concentration boundary layer of the reactants is increasing with η in both cases, and after a certain η value, they all coincide, i.e., after a certain η value, the homogeneous and heterogeneous reactions have no effect on the concentration of the reactants. This critical value of η () depends on the strength of the homogeneous reaction and increases with the value of K, but it does not depend on the strength of the heterogeneous reaction. A similar phenomenon is observed for the second solution. The graphs for the second solution with and 1 coincide. It is observed that the first solution is more stable and converges more easily than the second solution.
Figure 3 . Effect ofλon (a) and (b) for,,,,.
The concentration of the reactants depends on the Schmidt number (Sc) and heterogeneous reaction parameter. The variation of the concentration with K for different values of the Schmidt number and is shown in Figure 4. The Schmidt number is the ratio between a viscous diffusion rate and a molecular diffusion rate. For a fixed molecular diffusion rate, with increase in Schmidt number, the viscous diffusion rate increases, which reduces the velocity of the nanofluid and helps to increase the concentration of the nanofluid. Homogeneous reaction parameter and heterogeneous reaction parameter influence the reaction rate of the solute and help to reduce the concentration of the reactants, which concurs with the results reported by Khan and Pop [28].
Figure 4 . Concentration at the surface effect on (a) and (b) for,,,.
4 Conclusions
The present analysis investigates the effect of the homogeneous and heterogeneous reaction on the micropolar fluid flow through a porous medium past a porous stretching/shrinking sheet with suction. The momentum and concentration equations were transformed into a set of coupled nonlinear differential equations using similarity transformations and solved numerically by Matlab bvp4c package. We compared our results with those in the literature for some limiting case. A dual solution appeared for the shrinking sheet case. The effect of the dual solution is shown by tables and graphs. The momentum boundary layer thickness increased for the case of the first solution, while an opposite phenomenon appeared for the second solution and a similar phenomenon was observed for concentration profile. It was observed that the concentration at the surface decreased with the strengths of the homogeneous and heterogeneous reactions. The solute concentration, however, increased with the permeability and stretching/shrinking parameters. The velocity of the fluid and the concentration of the reactants at the surface increase with the stretching/shrinking parameter. Also, velocity increases due to the increase in micropolar parameter. The concentration of the reactants decreases with the strength of the homogeneous and heterogeneous reaction.
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
All the authors participated in the design of the study and helped to draft and proofread the manuscript. SS and PKK carried out the numerical computations.
Acknowledgements
The authors wish to thank the University of KwaZuluNatal for financial support.
References

Eringen, AC: Theory of micropolar fluids . J. Math. Mech.. 16, 1–18 (1966)

Hudimoto, B, Tokuoka, T: Twodimensional shear flows of linear micropolar fluids . Int. J. Eng. Sci.. 7, 515–522 (1969). Publisher Full Text

Lockwood, F, Benchaita, F, Friberg, S: Study of lyotropic liquid crystal in viscometric flow and elastohydrodynamics contact . Tribol. Trans.. 30, 539–548 (1987)

Ariman, T, Turk, MA, Sylvester, ND: On steady and pulsatile flow of blood . J. Appl. Mech.. 41, 1–7 (1974). Publisher Full Text

Eringen, AC: NonLocal Polar Field Theory, Academic Press, New York (1976)

Ariman, T, Turk, MA, Sylvester, ND: Microcontinuum fluid mechanics  a review . Int. J. Eng. Sci.. 12, 273–293 (1974). Publisher Full Text

Crane, LJ: Flow past a stretching plate . Z. Angew. Math. Phys.. 21, 645–647 (1970). Publisher Full Text

Dutta, BK, Roy, P, Gupta, AS: Temperature field in the flow over a stretching sheet with uniform heat flux . Int. Commun. Heat Mass Transf.. 12, 89–94 (1985). Publisher Full Text

Hayat, T, Javed, T, Abbas, Z: MHD flow of a micropolar fluid near a stagnation point towards a nonlinear stretching surface . Nonlinear Anal., Real World Appl.. 10(3), 1514–1526 (2009). Publisher Full Text

Ishak, A: Thermal boundary layer flow over a stretching sheet in a micropolar fluid with radiation effect . Meccanica. 45, 367–373 (2010). Publisher Full Text

Kelson, NA, Desseaux, A: Effect of surface condition on flow of micropolar fluid driven by a porous stretching sheet . Int. J. Eng. Sci.. 39, 1881–1897 (2001). Publisher Full Text

Mohammadein, AA, Gorla, RSR: Heat transfer in a micropolar fluid over a stretching sheet with viscous dissipation and internal heat generation . Int. J. Numer. Methods Heat Fluid Flow. 11, 50–58 (2001). Publisher Full Text

Erickson, LE, Fan, LT, Fox, VG: Heat and mass transfer on moving continuous flat plate with suction or injection . Ind. Eng. Chem. Fundam.. 5, 19–25 (1966). Publisher Full Text

Fox, VG, Erickson, LE, Fan, LT: Methods for solving the boundary layer equations for moving continuous flat surfaces with suction and injection . AIChE J.. 14, 726–736 (1968). Publisher Full Text

Elbashbeshy, EMA, Bazid, MAA: Heat transfer in a porous medium over a stretching surface with internal heat generation and suction or injection . Appl. Math. Comput.. 158, 799–807 (2004). Publisher Full Text

Bhargava, R, Kumar, L, Takhar, HS: Finite element solution of mixed convection micropolar fluid driven by a porous stretching sheet . Int. J. Eng. Sci.. 41, 2161–2178 (2003). Publisher Full Text

Bhargava, R, Sharma, S, Takhar, HS, Beg, OA, Bhargava, P: Numerical solutions for microplar transport phenomena over a nonlinear stretching sheet . Nonlinear Anal., Model. Control. 12, 45–63 (2007)

Abel, MS, Siddheshwar, PG, Mahesha, N: Numerical solution of the momentum and heat transfer equations for a hydromagnetic flow due to a stretching sheet of a nonuniform property micropolar liquid . Appl. Math. Comput.. 217, 5895–5909 (2011). Publisher Full Text

Narayana, M, Sibanda, P: Laminar flow of a nanoliquid film over an unsteady stretching sheet . Int. J. Heat Mass Transf.. 55, 7552–7560 (2012). Publisher Full Text

Kameswaran, PK, Narayana, M, Sibanda, P, Murthy, PVSN: Hydromagnetic nanofluid flow due to a stretching or shrinking sheet with viscous dissipation and chemical reaction effects . Int. J. Heat Mass Transf.. 55, 7587–7595 (2012). Publisher Full Text

Kameswaran, PK, Shaw, S, Sibanda, P, Murthy, PVSN: Homogeneousheterogeneous reactions in a nanofluid flow due to a porous stretching sheet . Int. J. Heat Mass Transf.. 57, 465–472 (2013). Publisher Full Text

Chaudhary, MA, Merkin, JH: A simple isothermal model for homogeneousheterogenous reactions in boundarylayer flow I. Equal diffusivities . Fluid Dyn. Res.. 16, 311–333 (1995). Publisher Full Text

Chaudhary, MA, Merkin, JH: A simple isothermal model for homogeneousheterogenous reactions in boundarylayer flow II. Unequal diffusivities . Fluid Dyn. Res.. 16, 335–359 (1995). Publisher Full Text

Merkin, JH: A model for isothermal homogeneousheterogeneous reactions in boundarylayer flow . Math. Comput. Model.. 24, 125–136 (1996). Publisher Full Text

Ziabakhsh, Z, Domairry, G, Bararnia, H, Babazadeh, H: Analytical solution of flow and diffusion of chemically reactive species over a nonlinearly stretching sheet immersed in a porous medium . J. Taiwan Inst. Chem. Eng.. 41, 22–28 (2010). Publisher Full Text

Chambre, PL, Acrivos, A: On chemical surfaces reactions in laminar boundary layer flows . J. Appl. Phys.. 27, 1322–1328 (1956). Publisher Full Text

Khan, WA, Pop, I: Flow near the twodimensional stagnation  point on an infinite permeable wall with a homogeneousheterogeneous reaction . Commun. Nonlinear Sci. Numer. Simul.. 15, 3435–3443 (2010). Publisher Full Text

Khan, WA, Pop, I: Effects of homogeneousheterogeneous reactions on the viscoelastic fluid toward a stretching sheet . J. Heat Transf.. 134, (2012) Article ID 064506

Guram, GS, Smith, AC: Stagnation flows of micropolar fluids with strong and weak interactions . Comput. Math. Appl.. 6, 213–233 (1980)

Jena, SK, Mathur, MN: Similarity solutions for laminar free convection flow of a thermomicropolar fluid past a nonisothermal vertical flat plate . Int. J. Eng. Sci.. 19, 1431–1439 (1981). Publisher Full Text

Peddieson, J: An application of the micropolar fluid model to the calculation of a turbulent shear flow . Int. J. Eng. Sci.. 10, 23–32 (1972). Publisher Full Text

Ahmadi, G: Selfsimilar solution of incompressible micropolar boundary layer flow over a semiinfinite plate . Int. J. Eng. Sci.. 14, 639–646 (1976). PubMed Abstract  Publisher Full Text

Wang, CY: Stagnation flow towards a shrinking sheet . Int. J. NonLinear Mech.. 43, 377–382 (2008). Publisher Full Text

Ishak, A, Lok, YY, Pop, I: Stagnation point flow over a shrinking sheet in a micropolar fluid . Chem. Eng. Commun.. 197, 1417–1427 (2010). Publisher Full Text

Rosali, H, Ishak, A, Pop, I: Micropolar fluid flow towards a stretching/shrinking sheet in a porous medium with suction . Int. Commun. Heat Mass Transf.. 39, 826–829 (2012). Publisher Full Text

Merkin, JH: On dual solutions occurring in mixed convection in a porous media . J. Eng. Math.. 20, 171–179 (1985)

Merrill, K, Beauchesne, M, Previte, J, Paullet, J, Weidman, P: Final steady flow near a stagnation point on a vertical surface in a porous medium . Int. J. Heat Mass Transf.. 49, 4681–4686 (2006). Publisher Full Text

Bhattacharyya, K: Dual solutions in boundary layer stagnationpoint flow and mass transfer with chemical reaction past a stretching/shrinking sheet . Int. Commun. Heat Mass Transf.. 38, 917–922 (2011). Publisher Full Text