Abstract
An analysis is carried out to study the heat transfer characteristics of steady twodimensional stagnationpoint flow of a copper (Cu)water nanofluid over a permeable stretching/shrinking sheet. The stretching/shrinking velocity and the ambient fluid velocity are assumed to vary linearly with the distance from the stagnationpoint. Results for the skin friction coefficient, local Nusselt number, velocity as well as the temperature profiles are presented for different values of the governing parameters. It is found that dual solutions exist for the shrinking case, while for the stretching case, the solution is unique. The results indicate that the inclusion of nanoparticles into the base fluid produces an increase in the skin friction coefficient and the heat transfer rate at the surface. Moreover, suction increases the surface shear stress and in consequence increases the heat transfer rate at the fluidsolid interface.
MSC: 34B15, 76D10.
Keywords:
nanofluid; stagnationpoint; stretching/shrinking sheet; suction/injection; heat transfer; dual solutions1 Introduction
Nanofluids are the suspension of metallic, nonmetallic or polymeric nanosized powders in base liquid which are employed to increase the heat transfer rate in various applications. The term nanofluid, first introduced by Choi [1], refers to the fluids with suspended nanoparticles. Most of the convectional heat transfer fluids such as water, ethylene glycol and mineral oils have low thermal conductivity and thus are inadequate to meet the requirements of today’s cooling rate. An innovative way of improving the thermal conductivities of such fluids is to suspend small solid particles in the base fluids to form slurries. An industrial application test was carried out by Liu et al.[2] and Ahuja [3], in which the effect of particle volumetric loading, size and flow rate on the slurry pressure drop and heat transfer behavior was investigated (Xuan and Li [4]). Experimental results by Eastman et al.[5] showed that an increase in thermal conductivity of approximately 60% is obtained for the nanofluid consisting of water and 5% volume fraction of CuO nanoparticles. The procedure for preparing a nanofluid is given in the paper by Xuan and Li [4].
Many of the publications on nanofluids are about understanding of their behaviors so that they can be utilized where straight heat transfer enhancement is paramount as in many industrial applications, nuclear reactors, transportation, electronics as well as biomedicine and food (see Ding et al.[6]). Nanofluid is a smart fluid, where the heat transfer capabilities can be reduced or enhanced at will. These fluids enhance thermal conductivity of the base fluid enormously, which is beyond the explanation of any existing theory. They are also very stable and have no additional problems, such as sedimentation, erosion, additional pressure drop and nonNewtonian behavior, due to the tiny size of nanoelements and the low volume fraction of nanoelements required for conductivity enhancement. Much attention has been paid in the past to this new type of composite material because of its enhanced properties and behavior associated with heat transfer, mass transfer, wetting and spreading as well as antimicrobial activities, and the number of publications related to nanofluids increases in an exponential manner. The enhanced thermal behavior of nanofluids could provide a basis for an enormous innovation for heat transfer intensification, which is of major importance to a number of industrial sectors including transportation, power generation, micromanufacturing, thermal therapy for cancer treatment, chemical and metallurgical sectors, as well as heating, cooling, ventilation and airconditioning. Nanofluids are also important for the production of nanostructured materials, for the engineering of complex fluids, as well as for cleaning oil from surfaces due to their excellent wetting and spreading behavior (Ding et al.[6]).
There are some nanofluid models available in the literature. Among the popular models are the model proposed by Buongiorno [7] and Tiwari and Das [8]. Buongiorno [7] noted that the nanoparticle absolute velocity can be viewed as the sum of the base fluid velocity and a relative velocity (that he calls the slip velocity). He considered in turn seven slip mechanisms: inertia, Brownian diffusion, thermophoresis, diffusiophoresis, Magnus effect, fluid drainage and gravity settling (Nield and Kuznetsov [9]). The nanofluid mathematical model proposed by Buongiorno [7] was very recently used by several researchers such as, among others, Nield and Kuznetsov [9,10], Kuznetsov and Nield [11,12], Khan and Pop [13], Khan and Aziz [14], Makinde and Aziz [15], Bachok et al.[16,17], etc. On the other hand, the Tiwari and Das model analyzes the behavior of nanofluids taking into account the solid volume fraction of the nanofluid. In the present paper, we study the flow and heat transfer characteristics near a stagnation region of a permeable stretching/shrinking sheet immersed in a Cuwater nanofluid using the Tiwari and Das model. It is worth mentioning that this model was recently employed in Refs. [1833], and the flow over a shrinking sheet was considered in Refs. [3444]. The velocity distribution of the twodimensional stagnation flow was first analyzed by Hiemenz (see White [45]) who discovered that this flow can be analyzed exactly by the NavierStokes equations. Homann (see White [45]) extended this problem to the axisymmetric stagnation flow and found that the solution differs a little from the plane flow, where the displacement and boundary layer thicknesses are slightly smaller and the wall shear stress is slightly larger. On the other hand, the temperature distributions of the Hiemenz and Homann flows were given by Goldstein [46] and Sibulkin [47], respectively. The governing partial differential equations are first transformed into a system of ordinary differential equations before being solved numerically. We study the effects of suction and injection at the boundary. Suction or injection of a fluid through the bounding surface, as, for example, in mass transfer cooling, can significantly change the flow field and, as a consequence, affect the heat transfer rate at the surface. In general, suction tends to increase the skin friction and heat transfer coefficients, whereas injection acts in the opposite manner (AlSanea [48]). Injection of fluid through a porous bounding heated or cooled wall is of general interest in practical problems involving film cooling, control of boundary layer, etc. This can lead to enhance heating (or cooling) of the system and can help to delay the transition from laminar flow (see Chaudhary and Merkin [49]). We mention to this end that studies of the boundary layer flows of a Newtonian (or regular) fluid past a permeable static or moving flat plate have been done by Merkin [50], Weidman et al.[51], Ishak et al.[52], Zheng et al.[53] and Zhu et al.[54,55], while Bachok et al.[32] have considered the boundary layers over a permeable moving surface in a nanofluid.
2 Mathematical formulation
Consider a stagnation flow of an incompressible nanofluid over a stretching/shrinking surface located at with a fixed stagnation point at . The stretching/shrinking velocity and the ambient fluid velocity are assumed to vary linearly from the stagnation point, i.e., and , where a and b are constant with . We note that and correspond to stretching and shrinking sheets, respectively. The simplified twodimensional equations governing the flow in the boundary layer of a steady, laminar, and incompressible nanofluid are (see Ahmad et al.[26])
subject to the boundary conditions
where u and v are the velocity components along the x and yaxes, respectively, is the mass transfer velocity, T is the temperature of the nanofluid, is the surface temperature, is the ambient temperature, is the viscosity of the nanofluid, is the thermal diffusivity of the nanofluid and is the density of the nanofluid, which are given by (Oztop and AbuNada [28])
Here, φ is the nanoparticle volume fraction, is the heat capacity of the nanofluid, is the thermal conductivity of the nanofluid, and are the thermal conductivities of the fluid and of the solid fractions, respectively, and and are the densities of the fluid and of the solid fractions, respectively. It should be mentioned that the use of the above expression for is restricted to spherical nanoparticles where it does not account for other shapes of nanoparticles (AbuNada [18]). Also, the viscosity of the nanofluid has been approximated by Brinkman [56] as viscosity of a base fluid containing dilute suspension of fine spherical particles.
The governing Eqs. (1)(3) subject to the boundary conditions (4) can be expressed in a simpler form by introducing the following transformation:
where η is the similarity variable and ψ is the stream function defined as and , which identically satisfies Eq. (1). Employing the similarity variables (6), Eqs. (2) and (3) reduce to the following ordinary differential equations:
subjected to the boundary conditions (4) which become
In the above equations, primes denote differentiation with respect to η, Pr is the Prandtl number, S is the suction/injection parameter and ε is the stretching/shrinking parameter defined respectively as
with for stretching and for shrinking.
The physical quantities of interest are the skin friction coefficient and the local Nusselt number , which are defined as
where the surface shear stress and the surface heat flux are given by
with and being the dynamic viscosity and thermal conductivity of the nanofluids, respectively. Using the similarity variables (6), we obtain
3 Numerical scheme
The nonlinear differential equations (7) and (8) along with the boundary conditions (9) form a twopoint boundary value problem (BVP) and are solved using a shooting method, by converting them into an initial value problem (IVP). This method is very well described in the recent papers by Bhattacharyya and Layek [57] and Bhattacharyya et al.[58]. In this method, we choose suitable finite values of η, say , which depend on the values of the parameters considered. First, the system of equations (7) and (8) is reduced to a firstorder system (by introducing new variables) as follows:
with the boundary conditions
Now we have a set of ‘partial’ initial conditions
As we notice, we do not have the values of and . To solve Eqs. (15) and (16) as an IVP, we need the values of and , i.e., and . We guess these values and apply the RungeKuttaFehlberg method, then see if this guess matches the boundary conditions at the very end. Varying the initial slopes gives rise to a set of profiles which suggest the trajectory of a projectile ‘shot’ from the initial point. That initial slope is sought which results in the trajectory ‘hitting’ the target, that is, the final value (Bailey et al.[59]).
To determine either the solution obtained is valid or not, it is necessary to check the velocity and the temperature profiles. The correct profiles must satisfy the boundary conditions at asymptotically. This procedure is repeated for other guessing values of and for the same values of parameters. If a different solution is obtained and the profiles satisfy the far field boundary conditions asymptotically but with different boundary layer thickness, then this solution is also a solution to the boundaryvalue problem (second solution).
4 Results and discussion
We have considered one type of nanoparticle, namely, copper (Cu), with water as the base fluid. The effects of the solid volume fraction of nanoparticles φ, the stretching/shrinking parameter ε and the suction/injection parameter S are analyzed. Following Oztop and AbuNada [28], AbuNada and Oztop [29] and Khanafer et al.[60], the value of the Prandtl number Pr is taken as 6.2 (for water) and the volume fraction of nanoparticles is from 0 to 0.2 () in which corresponds to the regular (Newtonian) fluid. The thermophysical properties of the base fluid (water) and the nanoparticles are given in Table 1. The numerical results are presented in Figures 16.
Figure 1 . Variation ofwithεfor different values ofSwhen.
Figure 2 . Variation ofwithεfor different values ofSwhenand.
Figure 3 . Variation of the skin friction coefficient withφfor different values ofSwhen.
Figure 4 . Variation of the local Nusselt number withφfor different values ofSwhenand.
Figure 5 . Velocity profiles for different values ofϕwhenand.
Figure 6 . Temperature profiles for different values ofϕwhen,and. Solid line for first solution and dash line for second solution.
Solving Eqs. (7) and (8) subject to the boundary conditions (9), it was found that dual solutions exist, which were obtained by setting different initial guesses for the missing values of and , where all profiles satisfy the far fields boundary conditions (9) asymptotically but with different shapes. The variations of and with ε are shown in Figures 1 and 2 for some values of the suction/injection parameter S. These figures show that there are regions of unique solutions for , dual (upper and lower branches) solutions for and no solutions for , where is the critical value of ε () beyond which Eqs. (7) and (8) have no solutions. Based on our computation, the critical values of ε, say , are −1.1098, −1.2465 and −1.4311 for , and , respectively. Thus, from this observation, the value of increases as S increases. Hence, suction widens the range of ε for which the solution exists, while injection acts in the opposite manner. Further, it should be mentioned that similar to other studies where dual solutions exist, we postulate that upper branch (the first) solutions of Eqs. (7) and (8) are stable and physically realizable, while the lower branch (the second) solutions are not. The procedure for showing this has been described by Merkin [50], Weidman et al.[51] and very recently by Postelnicu and Pop [61], so that we will not repeat it here.
Figures 3 and 4 illustrate the variations of the skin friction coefficient and the local Nusselt number given by Eqs. (13) and (14) with the nanoparticle volume fraction parameter φ for three different values of S with (stretching sheet). These figures show that these quantities increase almost linearly with φ. The presence of the nanoparticles in the fluids increases appreciably the effective thermal conductivity of the fluid and consequently enhances the heat transfer characteristics, as seen in Figure 4. Nanofluids have a distinctive characteristics, which is quite different from those of traditional solidliquid mixtures in which millimeter and/or micrometersized particles are involved. Such particles can clot equipment and can increase pressure drop due to settling effects. Moreover, they settle rapidly, creating substantial additional pressure (Khanafer et al.[60]). Figure 3 also shows that the effect of suction () is to increase the skin friction coefficient, and in consequence it increases the local Nusselt number, as presented in Figure 4.
The samples of velocity and temperature profiles for some values of parameters are presented in Figures 5 and 6. These profiles have essentially the same form as in the case of regular fluid (). The terms first solution and second solution refer to the curves shown in Figures 1 and 2, where the first solution has larger values of and compared to the second solution. Figures 5 and 6 show that the far field boundary conditions (9) are satisfied asymptotically, thus support the validity of the numerical results, besides supporting the existence of the dual solutions presented in Figures 1 and 2.
5 Conclusions
We have numerically studied the existence of dual similarity solutions in boundary layer flow over a stretching/shrinking sheet immersed in a nanofluid with suction and injection effects. Discussions were carried out for the effects of suction/injection parameter S, the nanoparticle volume fraction φ and the stretching/shrinking parameter ε on the skin friction coefficient and the local Nusselt number. It was found that dual solutions exist for the shrinking case, while for the stretching case, the solution is unique. The results indicate that the inclusion of nanoparticles into the base fluid produced an increase in the skin friction coefficient and the local Nusselt number. Moreover, these quantities increase with suction but decrease with injection.
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
The paper is the result of joint work of all authors who contributed equally to the final version of the paper. All authors read and approved the final manuscript.
Acknowledgements
The authors would like to thank the anonymous reviewers for their comments and suggestions which led to the improvement of this paper. The financial supports received from the Ministry of Higher Education, Malaysia (Project Code: FRGS/1/2012/SG04/UKM/01/1) and the Universiti Kebangsaan Malaysia (Project Code: DIP201231) are gratefully acknowledged.
References

Choi, US: Enhancing thermal conductivity of fluids with nanoparticles . ASME FED. 231, 99–103 (1995)

Liu, KV, Choi, US, Kasza, KE: Measurements of pressure drop and heat transfer in turbulent pipe flows of particulate slurries. Argonne National Laboratory Report, ANL8815, Argonne (1988)

Ahuja, AS: Augmentation of heat transfer in laminar flow of polystyrene suspensions . J. Appl. Phys.. 46, 3408–3425 (1975). Publisher Full Text

Xuan, Y, Li, Q: Heat transfer enhancement of nanofluids . Int. J. Heat Fluid Flow. 21, 58–64 (2000). Publisher Full Text

Eastman, JA, Choi, US, Li, S, Thompson, LJ, Lee, S: Enhanced thermal conductivity through the development of nanofluids . In: Komarneni S, Parker JC, Wollenberger HJ (eds.) Nanophase and Nanocomposite Materials II, pp. 3–11. MRS, Pittsburgh (1997)

Ding, Y, Chen, H, Wang, L, Yang, CY, He, Y, Yang, W, Lee, WP, Zhang, L, Huo, R: Heat transfer intensification using nanofluids . Kona. 25, 23–38 (2007)

Buongiorno, J: Convective transport in nanofluids . ASME J. Heat Transf.. 128, 240–250 (2006). 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

Nield, DA, Kuznetsov, AV: The ChengMinkowycz problem for natural convective boundarylayer flow in a porous medium saturated by a nanofluid . Int. J. Heat Mass Transf.. 52, 5792–5795 (2009). Publisher Full Text

Nield, DA, Kuznetsov, AV: The ChengMinkowycz problem for the doublediffusive natural convective boundary layer flow in a porous medium saturated by a nanofluid . Int. J. Heat Mass Transf.. 54, 374–378 (2011). Publisher Full Text

Kuznetsov, AV, Nield, DA: Natural convective boundarylayer flow of a nanofluid past a vertical plate . Int. J. Therm. Sci.. 49, 243–247 (2010). Publisher Full Text

Kuznetsov, AV, Nield, DA: Doublediffusive natural convective boundarylayer flow of a nanofluid past a vertical plate . Int. J. Therm. Sci.. 50, 712–717 (2011). Publisher Full Text

Khan, AV, Pop, I: Boundarylayer flow of a nanofluid past a stretching sheet . Int. J. Heat Mass Transf.. 53, 2477–2483 (2010). Publisher Full Text

Khan, WA, Aziz, A: Natural convection flow of a nanofluid over a vertical plate with uniform surface heat flux . Int. J. Therm. Sci.. 50, 1207–1214 (2011). Publisher Full Text

Makinde, OD, Aziz, A: Boundary layer flow of a nanofluid past a stretching sheet with a convective boundary condition . Int. J. Therm. Sci.. 50, 1326–1332 (2011). Publisher Full Text

Bachok, N, Ishak, A, Pop, I: Boundary layer flow of nanofluids over a moving surface in a flowing fluid . Int. J. Therm. Sci.. 49, 1663–1668 (2010). Publisher Full Text

Bachok, N, Ishak, A, Pop, I: Unsteady boundary layer flow of a nanofluid over a permeable stretching/shrinking sheet . Int. J. Heat Mass Transf.. 55, 2102–2109 (2012). Publisher Full Text

AbuNada, E: Application of nanofluids for heat transfer enhancement of separated flow encountered in a backward facing step . Int. J. Heat Fluid Flow. 29, 242–249 (2008). Publisher Full Text

Muthtamilselvan, M, Kandaswamy, P, Lee, J: Heat transfer enhancement of copperwater nanofluids in a liddriven enclosure . Commun. Nonlinear Sci. Numer. Simul.. 15, 1501–1510 (2010). Publisher Full Text

AbuNada, E, Oztop, HF: Effect of inclination angle on natural convection in enclosures filled with cuwater nanofluid . Int. J. Heat Fluid Flow. 30, 669–678 (2009). Publisher Full Text

Talebi, F, Houshang, A, Shahi, M: Numerical study of mixed convection flows in a square liddriven cavity utilizing nanofluid . Int. Commun. Heat Mass Transf.. 37, 79–90 (2010). Publisher Full Text

Bachok, N, Ishak, A, Pop, I: Flow and heat transfer over a rotating porous disk in a nanofluid . Physica B. 406, 1767–1772 (2011). Publisher Full Text

Bachok, N, Ishak, A, Nazar, R, Pop, I: Flow and heat transfer at a general three dimensional stagnation point flow in a nanofluid . Physica B. 405, 4914–4918 (2010). Publisher Full Text

Bachok, N, Ishak, A, Pop, I: Stagnationpoint flow over a stretching/shrinking sheet in a nanofluid . Nanoscale Res. Lett.. 6, 623 (2011). PubMed Abstract  BioMed Central Full Text  PubMed Central Full Text

Bachok, N, Ishak, A, Pop, I: Flow and heat transfer characteristics on a moving plate in a nanofluid . Int. J. Heat Mass Transf.. 55, 642–648 (2012). Publisher Full Text

Ahmad, S, Rohni, AM, Pop, I: Blasius and Sakiadis problems in nanofluids . Acta Mech.. 218, 195–204 (2011). Publisher Full Text

Yacob, NA, Ishak, A, Pop, I, Vajravelu, K: Boundary layer flow past a stretching/shrinking surface beneath an external uniform shear flow with a convective surface boundary condition in a nanofluid . Nanoscale Res. Lett.. 6, (2011) Article ID 314

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

AbuNada, E, Oztop, HF: Effect of inclination angle on natural convection in enclosures filled with Cuwater nanofluid . Int. J. Heat Fluid Flow. 30, 669–678 (2009). Publisher Full Text

Talebi, F, Houshang, A, Shahi, M: Numerical study of mixed convection flows in a square liddriven cavity utilizing nanofluid . Int. Commun. Heat Mass Transf.. 37, 79–90 (2010). Publisher Full Text

Rohni, AM, Ahmad, S, Pop, I: Boundary layer flow over a moving surface in a nanofluid beneath a uniform free stream . Int. J. Numer. Methods Heat Fluid Flow. 21, 828–846 (2011). Publisher Full Text

Bachok, N, Ishak, A, Pop, I: Boundary layer flow over a moving surface in nanofluid with suction or injection . Acta Mech. Sin.. 28, 34–40 (2012). Publisher Full Text

Yacob, NA, Ishak, A, Pop, I: FalknerSkan problem for a static or moving wedge in nanofluids . Int. J. Therm. Sci.. 50, 133–139 (2011). 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

Rosali, H, Ishak, A, Pop, I: Stagnation point flow and heat transfer over a stretching/shrinking sheet in a porous medium . Int. Commun. Heat Mass Transf.. 38, 1029–1032 (2011). Publisher Full Text

Bhattacharyya, K: Boundary layer flow and heat transfer over an exponentially shrinking sheet . Chin. Phys. Lett.. 28, (2011) Article ID 074701

Bhattacharyya, K: Dual solutions in unsteady stagnationpoint flow over a shrinking sheet . Chin. Phys. Lett.. 28, (2011) Article ID 084702

Bhattacharyya, K: Effects of radiation and heat source/sink on unsteady MHD boundary layer flow and heat transfer over a shrinking sheet with suction/injection . Front. Chem. Sci. Eng.. 5, 376–384 (2011). Publisher Full Text

Bhattacharyya, K, Pop, I: MHD boundary layer flow due to an exponentially shrinking sheet . Magnetohydrodynamics. 47, 337–344 (2011)

Bhattacharyya, K, Arif, MG, Pramanik, WA: MHD boundary layer stagnationpoint flow and mass transfer over a permeable shrinking sheet with suction/blowing and chemical reaction . Acta Tech.. 57, 1–15 (2012)

Bhattacharyya, K, Mukhopadhyay, S, Layek, GC, Pop, I: Effects of thermal radiation on micropolar fluid flow and heat transfer over a porous shrinking sheet . Int. J. Heat Mass Transf.. 55, 2945–2952 (2012). Publisher Full Text

Bhattacharyya, K, Vajravelu, K: Stagnationpoint flow and heat transfer over an exponentially shrinking sheet . Commun. Nonlinear Sci. Numer. Simul.. 17, 2728–2734 (2012). Publisher Full Text

Fang, T, Zhong, Y: Viscous flow over a shrinking sheet with an arbitrary surface velocity . Commun. Nonlinear Sci. Numer. Simul.. 15, 3768–3776 (2010). Publisher Full Text

Fang, T, Yao, S, Zhang, J, Aziz, A: Viscous flow over a shrinking sheet with a second order slip flow model . Commun. Nonlinear Sci. Numer. Simul.. 15, 1831–1842 (2010). Publisher Full Text

Goldstein, S: Modern Development in Fluid Dynamics, Oxford University Press, London (1938)

Sibulkin, M: Heat transfer near the forward stagnation point of a body of revolution . J. Aeronaut. Sci.. 19, 570–571 (1952)

AlSanea, SA: Mixed convection heat transfer along a continuously moving heated vertical plate with suction or injection . Int. J. Heat Mass Transf.. 47, 1445–1465 (2004). Publisher Full Text

Chaudhary, MA, Merkin, JH: The effect of blowing and suction on free convection boundary layers on vertical surface with prescribed heat flux . J. Eng. Math.. 27, 265–292 (1993). Publisher Full Text

Merkin, JH: A note on the similarity equations arising in free convection boundary layers with blowing and suction . Z. Angew. Math. Phys.. 45, 258–274 (1994). PubMed Abstract  Publisher Full Text

Weidman, PD, Kubitschek, DG, Davis, AMJ: The effect of transpiration on selfsimilar boundary layer flow over moving surfaces . Int. J. Eng. Sci.. 44, 730–737 (2006). Publisher Full Text

Ishak, A, Nazar, R, Pop, I: Boundary layer on a moving wall with suction and injection . Chin. Phys. Lett.. 24, 2274–2276 (2007). Publisher Full Text

Zheng, L, Wang, L, Zhang, X: Analytic solutions of unsteady boundary flow and heat transfer on a permeable stretching sheet with nonuniform heat source/sink . Commun. Nonlinear Sci. Numer. Simul.. 16, 731–740 (2011). Publisher Full Text

Zhu, J, Zheng, L, Zhang, X: Homotophy analysis method for hydromagnetic plane and axisymmetric stagnationpoint flow with velocity slip . World Acad. Sci., Eng. Technol.. 63, 151–155 (2010)

Zhu, J, Zheng, L, Zhang, X: The effect of the slip condition on the MHD stagnationpoint over a powerlaw stretching sheet . Appl. Math. Mech.. 31, 439–448 (2010). Publisher Full Text

Brinkman, HC: The viscosity of concentrated suspensions and solutions . J. Chem. Phys.. 20, 571–581 (1952). Publisher Full Text

Bhattacharyya, K, Layek, GC: Effects of suction/blowing on steady boundary layer stagnationpoint flow and heat transfer towards a shrinking sheet with thermal radiation . Int. J. Heat Mass Transf.. 54, 302–307 (2011). Publisher Full Text

Bhattacharyya, K, Mukhopadhyay, S, Layek, GC: Slip effects on boundary layer stagnationpoint flow and heat transfer towards a shrinking sheet . Int. J. Heat Mass Transf.. 54, 308–313 (2011). Publisher Full Text

Bailey, PB, Shampine, LF, Waltman, PE: Nonlinear Two Point Boundary Value Problems, Academic Press, New York (1968)

Khanafer, K, Vafai, K, Lightstone, M: Buoyancydriven heat transfer enhancement in a twodimensional enclosure utilizing nanofluids . Int. J. Heat Mass Transf.. 46, 3639–3653 (2003). Publisher Full Text

Postelnicu, A, Pop, I: FalknerSkan boundary layer flow of a powerlaw fluid past a stretching wedge . Appl. Math. Comput.. 217, 4359–4368 (2011). Publisher Full Text