Abstract
The paper discusses the effects of homogeneousheterogeneous reactions on stagnationpoint flow of a nanofluid over a stretching or shrinking sheet. The model presented describes mass transfer in copperwater and silverwater 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:
homogeneousheterogeneous reactions; nanofluids; volume fraction model; stretching or shrinking sheet; magnetic field strength1 Introduction
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 [1] was the first to consider the steady twodimensional 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.[2], Erickson et al.[3], Mucoglu and Chen [4], Grubka and Bobba [5], Karwe and Jaluria [6], Chen [7], AboEldahab and ElAziz [8], Salem and ElAziz [9], Ali [10], Ishak et al.[11]).
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.[12]. 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.[13] analyzed MHD stagnationpoint flow from a shrinking sheet. They found that dual solutions existed for small values of the magnetic parameter. The stagnationpoint flow over a stretching or shrinking sheet in a nanofluid was investigated by Bachok et al.[14]. They showed that adding nanoparticles to a base fluid increased the skin friction and heat transfer coefficients. Recently, Narayana and Sibanda [15] 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 Cuwater nanoliquid. However, the opposite appeared to be true in the case of an water nanofluid. Kameswaran et al.[16] 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 homogeneousheterogeneous reactions in the boundary layer flow of a viscous fluid past a flat plate was presented by Merkin [17]. He modeled the homogeneous reaction by a cubic autocatalysis process and the heterogeneous reaction by a firstorder process. Chaudhary and Merkin [18] analyzed homogeneousheterogeneous 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.[19] studied the diffusion of a chemically reactive species into a nonlinearly stretching sheet immersed in a porous medium. Chambre and Acrivos [20] studied isothermal chemical reactions on laminar boundary layer flow. The twodimensional stagnationpoint flow near an infinite permeable wall with a homogeneousheterogeneous reaction was studied by Khan and Pop [21], while Khan and Pop [22] and Bachok et al.[23] studied the effects of homogeneousheterogeneous reactions on fluid flow due to a stretching sheet. Recently, the effects of homogeneousheterogeneous reactions in nanofluid flow due to a porous stretching sheet were studied by Kameswaran et al.[24]. 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 Cuwater nanofluid.
This article presents a study of homogeneousheterogeneous 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 twodimensional steady boundary layer flow of an incompressible nanofluid over a stretching or shrinking sheet. A Cartesian coordinate system is used with the xaxis along the sheet and the yaxis 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 waterbased 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 homogeneousheterogeneous reaction model proposed by Chaudhary and Merkin [18] of the form
while on the catalyst surface we have the single, isothermal, firstorder 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 [26] 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
where , is the dimensionless stream function and .
The velocity components are given by
The concentrations of the chemical species A and B are represented as
where and are dimensionless concentrations. On using equations (8), (9), (11) and (12), equations (4)(7) transform to the following twopoint boundary value problem:
The nondimensional 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
where
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 [18]). In this case we have, from equations (16), (17) and (18),
Thus equations (14) and (15) reduce to
subject to the boundary conditions
It is quite straightforward to show that the skin friction coefficient , which characterizes the surface drag, is given by
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 Cuwater and Agwater 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 skinfriction coefficient were compared with the available published results of Jafar et al.[12], Bachok et al.[23] and Wang [28] in Tables 2, 3 and 4.
Table 2. Comparison offrom current results with Wang[28]and Jafar et al.[12]for various values ofλfor the stretching sheet whenand
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.[12] and Wang [28] 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.[23] in the absence of the particular physical parameter.
The effects of the Schmidt number, stretching, magnetic and chemical reaction parameters are shown in Figures 212.
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 Cuwater 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 Cuwater 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 [29].
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 Cuwater 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 Cuwater nanofluid. However, beyond , the opposite is true for the clear fluid and the Cuwater 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.[24].
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 Agwater than for Cuwater. 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 Cuwater 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 Agwater and Cuwater 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 wellknown 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 [30]. 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.
4 Conclusions
The effects of homogeneousheterogeneous 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 Cuwater and Agwater nanofluids. The boundary layer thickness for the first solution is always thinner than that for the second solution.
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
The work including proof reading was done by all the authors.
Acknowledgements
The authors are grateful to the University of KwaZuluNatal for financial support.
References

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

Grubka, LJ, Bobba, KM: Heat transfer characteristics of a continuous, stretching surface with variable temperature. J. Heat Transf.. 107, 248–250 (1985). PubMed Abstract  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

AboEldahab, EM, ElAziz, 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, ElAziz, 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 stagnationpoint flow towards a shrinking sheet. Int. J. Numer. Methods Heat Fluid Flow. 21, 61–72 (2011). Publisher Full Text

Bachok, N, Ishak, A, Pop, I: Stagnation point flow over a stretching/shrinking sheet in a nanofluid. Nanoscale Res. Lett.. 6, (2011) Article ID 623

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 homogeneousheterogeneous reactions in boundary layer flow. Math. Comput. Model.. 24, 125–136 (1996). Publisher Full Text

Chaudhary, MA, Merkin, JH: A simple isothermal model for homogeneousheterogeneous 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 twodimensional stagnationpoint 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

Bachok, N, Ishak, A, Pop, I: On the stagnationpoint flow towards a stretching sheet with homogeneousheterogeneous reactions effects. Commun. Nonlinear Sci. Numer. Simul.. 16, 4296–4302 (2011). 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

Tiwari, RK, Das, MK: Heat transfer augmentation in a twosided liddriven 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, AbuNada, 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. NonLinear Mech.. 43, 377–382 (2008). Publisher Full Text

Jat, RN, Chaudhary, S: MHD flow and heat transfer over a stretching sheet. Appl. Math. Sci.. 3, 1285–1294 (2009)

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