The paper discusses the effects of homogeneous-heterogeneous reactions on stagnation-point flow of a nanofluid over a stretching or shrinking sheet. The model presented describes mass transfer in copper-water and silver-water nanofluids. The governing system of equations is solved numerically, and the study shows that dual solutions exist for certain suction/injection, stretching/shrinking and magnetic parameter values. Comparison of the numerical results is made with previously published results for special cases.
Keywords:homogeneous-heterogeneous reactions; nanofluids; volume fraction model; stretching or shrinking sheet; magnetic field strength
Problems involving fluid flow over stretching or shrinking surfaces can be found in many manufacturing processes such as in polymer extrusion, wire and fiber coating, foodstuff processing, etc. Crane  was the first to consider the steady two-dimensional flow of a Newtonian fluid driven by a stretching elastic flat sheet which moved in its own plane with velocity varying linearly with the distance from a fixed point. This study was subsequently extended by many authors to explore various aspects of heat transfer in a fluid surrounding a stretching sheet (Tsou et al., Erickson et al., Mucoglu and Chen , Grubka and Bobba , Karwe and Jaluria , Chen , Abo-Eldahab and El-Aziz , Salem and El-Aziz , Ali , Ishak et al.).
The magnetohydrodynamic effect has important engineering applications in electrical motors. Heat transfer over a stretching or shrinking sheet subject to an external magnetic field, viscous dissipation and joule effects was studied by Jafar et al.. They observed that the flow and heat transfer characteristics for a shrinking sheet were quite different from those of a stretching sheet. Lok et al. analyzed MHD stagnation-point flow from a shrinking sheet. They found that dual solutions existed for small values of the magnetic parameter. The stagnation-point flow over a stretching or shrinking sheet in a nanofluid was investigated by Bachok et al.. They showed that adding nanoparticles to a base fluid increased the skin friction and heat transfer coefficients. Recently, Narayana and Sibanda  investigated the laminar flow of a nanoliquid film over an unsteady stretching sheet. They noticed that the effect of an increase in the nanoparticle volume fraction was to reduce the axial velocity and free stream velocity in the case of a Cu-water nanoliquid. However, the opposite appeared to be true in the case of an -water nanofluid. Kameswaran et al. studied the effects of viscous dissipation and a chemical reaction in hydromagnetic nanofluid flow due to a stretching or shrinking sheet. They found that the velocity profiles decreased with an increase in the nanoparticle volume fraction, while the opposite was true in the case of temperature and concentration profiles. The study further showed that liquids with nanoparticle suspensions were better suited for effective cooling of the stretching sheet due to their enhanced conductivity and thermal properties.
Most chemically reacting systems involve both homogeneous and heterogeneous reactions (combustion, catalysis and biochemical systems). The simple combustion model helps us to understand the combustion phenomenon in many complex engineering applications such as in aircraft and rocket engines. A model for isothermal homogeneous-heterogeneous reactions in the boundary layer flow of a viscous fluid past a flat plate was presented by Merkin . He modeled the homogeneous reaction by a cubic autocatalysis process and the heterogeneous reaction by a first-order process. Chaudhary and Merkin  analyzed homogeneous-heterogeneous reactions in boundary layer flow. They presented a numerical solution of the boundary layer equations near the leading edge of a flat plate. Ziabakhsh et al. studied the diffusion of a chemically reactive species into a nonlinearly stretching sheet immersed in a porous medium. Chambre and Acrivos  studied isothermal chemical reactions on laminar boundary layer flow. The two-dimensional stagnation-point flow near an infinite permeable wall with a homogeneous-heterogeneous reaction was studied by Khan and Pop , while Khan and Pop  and Bachok et al. studied the effects of homogeneous-heterogeneous reactions on fluid flow due to a stretching sheet. Recently, the effects of homogeneous-heterogeneous reactions in nanofluid flow due to a porous stretching sheet were studied by Kameswaran et al.. They found that the concentration at the surface decreased with the strength of the heterogeneous reaction. In the case of a shrinking sheet, they showed that the velocity profiles decreased with increasing nanoparticle volume fraction in the case of a Cu-water nanofluid.
This article presents a study of homogeneous-heterogeneous reactions on MHD nanofluid stagnation point flow due to a stretching or shrinking sheet. The transformed nonlinear conservation equations are solved numerically.
2 Mathematical formulation
Consider a two-dimensional steady boundary layer flow of an incompressible nanofluid over a stretching or shrinking sheet. A Cartesian co-ordinate system is used with the x-axis along the sheet and the y-axis normal to the sheet. The flow configuration and the coordinate system are shown in Figure 1. The velocity of the outer flow is of the form and the velocity of the stretching or shrinking sheet is , where and are constants. The fluid is a water-based nanofluid containing copper (Cu) or silver (Ag) nanoparticles. The base fluid and the nanoparticles are in thermal equilibrium with no slip occurring between them. We assume the simple homogeneous-heterogeneous reaction model proposed by Chaudhary and Merkin  of the form
while on the catalyst surface we have the single, isothermal, first-order reaction
where a and b are the concentrations of the chemical species A and B, and and are the rate constants. We assume that both reaction processes are isothermal. Under these assumptions, the boundary layer equations governing the flow can be written in the dimensional form [18,25]
The boundary conditions for equations (3)-(6) are given in the form
where u, v are the velocity components in the x and y directions, respectively, and are the respective diffusion species coefficients of A and B, is a positive constant. The effective dynamic viscosity of the nanofluid was given by Brinkman  as
where ϕ is the solid volume fraction of nanoparticles. The effective density of the nanofluids is given as
Figure 1. Physical model and coordinate system.
Here, the subscripts nf, f and s represent the thermophysical properties of the nanofluid, the base fluid and nanoparticles, respectively.
The continuity equation (3) is satisfied by introducing a stream function ψ such that
The velocity components are given by
The concentrations of the chemical species A and B are represented as
The non-dimensional parameters in equations (13)-(18) are the magnetic parameter M, the Schmidt number Sc, the measure of the strength of the homogeneous reaction K, the ratio of diffusion coefficients δ, the mass transfer parameter , with for suction and for injection, the measure of the strength of the heterogeneous reaction , the Reynolds number Re and is the stretching parameter, with for stretching and for shrinking, respectively. These parameters are respectively defined as
In most applications, we expect the diffusion coefficients of chemical species A and B to be of a comparable size. This leads us to making a further assumption that the diffusion coefficients and are equal, i.e., (Chaudhary and Merkin ). In this case we have, from equations (16), (17) and (18),
Thus equations (14) and (15) reduce to
subject to the boundary conditions
3 Results and discussion
The system of ordinary differential equations (13) and (23) with boundary conditions (16) and (24) was solved numerically using Matlab bvp4c routine. We considered Cu-water and Ag-water nanofluids. The thermophysical properties of the nanofluids used in this paper are given in Table 1.
In order to determine the accuracy of our numerical results, the present results for the skin-friction coefficient were compared with the available published results of Jafar et al., Bachok et al. and Wang  in Tables 2, 3 and 4.
Tables 2 and 3 give the coefficient for different parameter values. Table 2 gives a comparison of the present results with those obtained by Jafar et al. and Wang  when , for different values of the stretching parameter. We observe that for increasing λ, the present results are in good agreement with results in the literature.
Table 4 gives the values of for selected λ when , . We note here that with decreasing λ the first solution decreases while the second solution increases. These results are in good agreement with the results obtained by Bachok et al. in the absence of the particular physical parameter.
Figure 2. Effects ofλon velocity, when,,,,.
Figure 2 shows the effects of both stretching and shrinking on the velocity profiles in the case of a Cu-water nanofluid. We observe that in both cases, the velocity profiles increase with the parameter λ. Further, we note that for a shrinking sheet, the velocity in the case of a Cu-water nanofluid is larger than that of a clear fluid. The opposite is, however, true for the case of a stretching sheet. The momentum boundary layer thickness decreases as λ increases and the flow has an inverted boundary layer structure when . The findings in the case of a clear fluid are similar to the results obtained by Jat and Chaudhary .
Figures 3(a) and 3(b) illustrate the effect of the magnetic parameter, nanoparticle volume fraction and stretching or shrinking parameters on the velocity profiles. We note that for both stretching and shrinking sheets, the fluid velocity increases with λ and M. Furthermore, increasing the value of M also causes thinning of the boundary layer. This implies an increase in the velocity gradient . Thus the magnetic field enhances the fluid motion in the boundary layer in the case of a clear fluid. The same trend is observed in the case of a Cu-water nanofluid. We also observe in the case of injection that, with increasing magnetic parameter, the increment in the momentum boundary layer is more significant than in the case of suction.
Figure 3. Effects ofMon velocity when (a),,,and (b),,,.
Figures 4(a) and 4(b) illustrate the effects of the stretching or shrinking parameter and volume fraction on the solute concentration when , , , , . We observe that the concentration profiles increase with the stretching parameter. We note also that the concentration increases as λ varies from to for both the clear fluid and the Cu-water nanofluid. However, beyond , the opposite is true for the clear fluid and the Cu-water nanofluid.
Figure 4. Variation of concentration for different values ofλ, when (a),,,,and (b),,,,.
The variation of for different values of K and is shown in Figures 5 and 6, respectively. From Figure 5 we observe that the concentration at the surface decreases as the strength of the heterogeneous reaction increases. This is simply explained by the fact that the strength of the chemical reaction depends on the concentration. On the other hand, from Figure 6, we found that decreases with increasing K and . These findings are similar to the results reported by Kameswaran et al..
Figure 5. Effects ofandKon concentration, when,,,,.
Figure 6. Effects ofKandon concentration, when,,,.
Figures 7 and 8 show the effects of K and λ on the concentration when the other parameters are fixed. We note, as expected, that the wall concentration decreases as the strength of the homogeneous reaction increases. The level of decrease is, however, more significant in the case of Ag-water than for Cu-water. Figure 8 shows the influence of stretching on the concentration profiles. It can be seen that the concentration decreases when the sheet is shrunk and increases with stretching.
The effects of suction/injection parameter on the velocity and concentration profiles are presented in Figure 9 for a Cu-water nanofluid. These profiles are qualitatively similar to the profiles obtained in the case of a clear fluid. Here we note dual solutions, i.e., the first and second solutions. The velocity and concentration profiles decrease with increasing in the case of the first solution. The opposite is, however, true in the case of the second solution. The far field boundary conditions are asymptotically satisfied, thus supporting the validity of the numerical solutions and the existence of dual solutions. In Figure 9(b), the concentration increases more rapidly with increasing suction/injection in the case of the second solution.
Figure 9. Effects ofon (a) velocity and (b) concentration profiles, when,,,,and.
Figure 10 shows the effect of the magnetic parameter on and for both Ag-water and Cu-water nanofluids. The velocity and concentration profiles are higher for Cu particles compared to Ag nanoparticles in the absence of the magnetic parameter. Increasing the value of M causes thinning of the boundary layer. The magnetic field enhances the fluid motion at the boundary, and for clear fluids these results are similar to well-known results in the literature.
Figure 10. Effects of magnetic parameterMon (a) velocity and (b) concentration profiles, when,,,,and.
The variation of the velocity and concentration profiles with the stretching/shrinking parameter is shown in Figure 11. The dual velocity and concentration profiles decrease with an increase in the magnitude of λ in the case of the first solution and increases in the case of the second solution. It is to be noted that momentum boundary layer thickness for the second solution is thicker than for the first solution. For the case of a clear fluid, the results are similar to those obtained by Bhattacharyya . We also note that the velocity gradient at the surface increases with λ, which is consistent with the results predicted from the computation of the skin friction coefficient.
Figure 11. Effects ofλon (a) velocity and (b) concentration profiles, when,,and.
The variation of the reduced skin friction coefficient with λ is shown in Figure 12(a). The values of are positive when and negative when . Physically, a positive implies that the fluid exerts a drag force on the plate and a negative sign implies the opposite. It is evident that dual solutions of equations (13) and (23) subject to the boundary conditions (16) and (24) exist when . There is a critical value of for which the first and second solutions meet. This critical value depends on the values of the other embedded parameters, and we found, for instance, that when , when and when . These results show that increases with M. However, for , the solution is not unique, there being two solutions for each λ. Figure 12(b) shows that increases with the magnetic parameter M.
Figure 12. Effects ofλon (a)and (b), when,and.
The effects of homogeneous-heterogeneous reactions in MHD nanofluid flow due to a stretching or shrinking sheet have been studied. The transformed governing nonlinear differential equations have been solved numerically. Dual solutions for the velocity and concentration distributions have been obtained for some values of the stretching/shrinking, suction/injection and magnetic parameters. The effects of physical and fluid parameters on the velocity, concentration and skin friction have been analyzed. It was observed that for both the cases of shrinking and stretching sheets, the fluid velocity increased with the magnetic parameter. The concentration at the surface decreased as the strength of heterogeneous reactions increased for both Cu-water and Ag-water nanofluids. The boundary layer thickness for the first solution is always thinner than that for the second solution.
The authors declare that they have no competing interests.
The work including proof reading was done by all the authors.
The authors are grateful to the University of KwaZulu-Natal for financial support.
Crane, LJ: Flow past a stretching plate. Z. Angew. Math. Phys.. 21, 645–647 (1970). Publisher Full Text
Tsou, FK, Sparrow, EM, Goldstein, RJ: Flow and heat transfer in the boundary layer on a continuous moving surface. Int. J. Heat Mass Transf.. 10, 219–235 (1967). Publisher Full Text
Erickson, LE, Fan, LT, Fox, VG: Heat and mass transfer on a moving continuous flat plate with suction or injection. Ind. Eng. Chem. Fundam.. 5, 19–25 (1966). Publisher Full Text
Mucoglu, A, Chen, TS: Mixed convection on inclined surfaces. J. Heat Transf.. 101, 422–426 (1979). Publisher Full Text
Karwe, MV, Jaluria, Y: Fluid flow and mixed convection transport from a moving plate in rolling and extrusion processes. J. Heat Transf.. 110, 655–661 (1988). Publisher Full Text
Chen, CH: Laminar mixed convection adjacent to vertical, continuously stretching sheets. Heat Mass Transf.. 33, 471–476 (1998). Publisher Full Text
Abo-Eldahab, EM, El-Aziz, MA: Blowing/suction effect on hydromagnetic heat transfer by mixed convection from an inclined continuously stretching surface with internal heat generation/absorption. Int. J. Therm. Sci.. 43, 709–719 (2004). Publisher Full Text
Salem, AM, El-Aziz, MA: Effect of Hall currents and chemical reaction on hydromagnetic flow of a stretching vertical surface with internal heat generation/absorption. Appl. Math. Model.. 32, 1236–1254 (2008). Publisher Full Text
Ali, ME: Heat transfer characteristics of a continuous stretching surface. Wärme- Stoffübertrag.. 29, 227–234 (1994). PubMed Abstract
Ishak, A, Nazar, R, Pop, I: Boundary layer flow and heat transfer over an unsteady stretching vertical surface. Meccanica. 44, 369–375 (2009). Publisher Full Text
Jafar, K, Nazar, R, Ishak, A, Pop, I: MHD flow and heat transfer over stretching/shrinking sheets with external magnetic field, viscous dissipation and joule effects. Can. J. Chem. Eng.. 90, 1336–1346 (2012). Publisher Full Text
Lok, YY, Ishak, A, Pop, I: MHD stagnation-point flow towards a shrinking sheet. Int. J. Numer. Methods Heat Fluid Flow. 21, 61–72 (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
Merkin, JH: A model for isothermal homogeneous-heterogeneous reactions in boundary layer flow. Math. Comput. Model.. 24, 125–136 (1996). Publisher Full Text
Chaudhary, MA, Merkin, JH: A simple isothermal model for homogeneous-heterogeneous reactions in boundary layer flow. I. Equal diffusivities. Fluid Dyn. Res.. 16, 311–333 (1995). 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 surface reactions in laminar boundary layer flows. J. Appl. Phys.. 27, 1322–1328 (1956). Publisher Full Text
Khan, WA, Pop, I: Flow near the two-dimensional stagnation-point on an infinite permeable wall with a homogeneous-heterogeneous reaction. Commun. Nonlinear Sci. Numer. Simul.. 15, 3435–3443 (2010). Publisher Full Text
Bachok, N, Ishak, A, Pop, I: On the stagnation-point flow towards a stretching sheet with homogeneous-heterogeneous reactions effects. Commun. Nonlinear Sci. Numer. Simul.. 16, 4296–4302 (2011). Publisher Full Text
Kameswaran, PK, Shaw, S, Sibanda, P, Murthy, PVSN: Homogeneous-heterogeneous reactions in a nanofluid flow due to a porous stretching sheet. Int. J. Heat Mass Transf.. 57, 465–472 (2013). Publisher Full Text
Tiwari, RK, Das, MK: Heat transfer augmentation in a two-sided lid-driven differentially heated square cavity utilizing nanofluids. Int. J. Heat Mass Transf.. 50, 2002–2018 (2007). Publisher Full Text
Brinkman, HC: The viscosity of concentrated suspensions and solutions. J. Chem. Phys.. 20, 571–581 (1952). Publisher Full Text
Oztop, HF, Abu-Nada, E: Numerical study of natural convection in partially heated rectangular enclosures filled with nanofluids. Int. J. Heat Fluid Flow. 29, 1326–1336 (2008). Publisher Full Text
Wang, CY: Stagnation flow towards a shrinking sheet. Int. J. Non-Linear Mech.. 43, 377–382 (2008). Publisher Full Text
Bhattacharyya, K: Dual solutions in boundary layer stagnation-point flow and mass transfer with chemical reaction past a stretching/shrinking sheet. Int. Commun. Heat Mass Transf.. 38, 917–922 (2011). Publisher Full Text