Abstract

This paper investigates the hydromagnetic boundary layer flow and heat transfer of a non-Newtonian Casson fluid in the neighborhood of a stagnation point over a stretching surface in the presence of velocity and thermal slips at the boundary. The governing partial differential equations are transformed into nonlinear ordinary differential equations using similarity transformations. The analytic solutions are developed by a homotopy analysis method (HAM). The results pertaining to the present study indicate that the flow and temperature fields are significantly influenced by Casson parameter (), the magnetic parameter , the velocity slip parameter , and the thermal slip parameter . An increase in the velocity slip parameter causes decrease in the flow velocity, while an increase in the value of the thermal slip parameter causes increase in the temperature of the fluid. It is also observed that the velocity at a point decreases with increase in .

1. Introduction

The problems of flow and heat transfer in the boundary layer adjacent to a continuous moving surface have received great attention during the last decades owing to the abundance of practical applications in chemical and manufacturing processes, such as polymer extrusion, continuous casting of metals, glass fibre production, hot rolling of paper, and wire drawing. Sakiadis [1] was the first, among others, to investigate the flow behavior on continuous solid surface. Thereafter, numerous investigations were made on the flow and heat transfer over a stretching surface in different directions [28].

All the previous researchers restricted their analyses to flow and heat transfer for the Newtonian fluid. In recent years, it has been observed that a number of industrial fluids such as molten plastics, polymeric liquids, blood, food stuff, and slurries exhibit non-Newtonian fluid behavior. Different types of non-Newtonian fluids are viscoelastic fluid, couple stress fluid, micropolar fluid, power-law fluid, Casson fluid, and many others. Rajagopal et al. [9] and Siddappa and Abel [10] studied the flow of a viscoelastic fluid over a linear stretching sheet. Troy et al. [11], Lawrence and Rao [12], and McLeod and Rajagopal [13] discussed the problem of uniqueness/nonuniqueness of the flow of a non-Newtonian viscoelastic fluid over a stretching sheet. Rajagopal et al. [9] analyzed the solutions for the flow of viscoelastic fluid over a stretching sheet. This study was further generalized to investigate the flow of short memory fluid of type Walter's liquid B by several authors, such as Andersson [14], Rollins and Vajravelu [15], and Abel and Veena [16].

Although different models are proposed to explain the behavior of non-Newtonian fluids, the most important non-Newtonian fluid possessing a yield value is the Casson fluid. This fluid has significant applications in polymer processing industries and biomechanics. We can define a casson fluid as a shear thinning fluid which is assumed to have an infinite viscosity at zero rate of shear. Casson's constitutive equations are found to describe accurately the flow curves of suspensions of pigments in lithographic varnishes used for preparation of printing inks and silicon suspensions [17]. Various experiments performed on blood with varying haematocrits, anticoagulants, temperatures, and so forth strongly suggest the behavior of blood as a casson fluid [1820]. In particular, casson fluid model describes the flow characteristics of blood more accurately at low shear rates and when it flows through small vessels [21].

Fluid flow in microelectromechanical systems (MEMS) has become an interesting topic because in microscale dimensions, the fluid flow behavior deviates significantly from the traditional no-slip flow. Under the microscale dimensions, the fluid motion still obeys the Navier-Stokes equations but under slip velocity boundary conditions. For large scale problems with low density, the fluid can be modeled as a rarefied gas, and rarefied gas flows with slip boundary conditions are often countered in the microscale devices. The nonadherence of the fluid to a solid boundary, known as velocity slip, is a phenomenon observed in certain circumstances. Partial slips occur for fluids with particulate such as emulsions, suspensions, foams, and polymer solutions. Fluids exhibiting slip are important in technological applications such as in the polishing of artificial heart valves and internal cavities. With a slip at the wall boundary, the flow behavior and the shear stress in the fluid are quite different from those in the no-slip flows. The slip flows in different flow configurations were studied in recent years (see [2227]). The MHD flow under slip conditions over a permeable shrinking surface was solved analytically by Fang et al. [28], and they reported that the velocity slip at the shrinking surface greatly affects the velocity distribution and drag forces on the wall. Bhattacharyya et al. [29] studied the effects of partial slip on the boundary layer stagnation-point flow and heat transfer towards a shrinking surface. The slip effect on MHD boundary layer flow over a flat plate was also considered by Bhattacharyya et al. [30]. Ariel et al. [31] considered the flow of a viscoelastic fluid over a stretching sheet with partial slip, and Ariel [32] also studied the slip effect on stagnation-point flow of an elasticoviscous fluid over a wall. Fang et al. [33] obtained the exact analytic solution of MHD flow under slip condition over a permeable stretching sheet.

However, in the literature, articles are not available on the MHD boundary layer stagnation-point flow of a casson fluid over a stretching sheet. The present paper aims to study the slip effect on MHD boundary layer flow and heat transfer of a casson fluid over a stretching sheet using the homotopy analysis method (HAM) which is successively applied to various interesting problems [3441]. The self-similar equations are solved analytically by HAM for a range of values of the physical parameters. The results obtained are then compared with those from the available literature for some particular values of the physical parameters, and it is found that they are in a good agreement. The results pertaining to the present study indicate that the flow and the temperature fields are greatly influenced by the velocity and thermal slip parameters.

2. Flow Analysis

Consider the steady magnetohydrodynamic (MHD) flow of a non-Newtonian Casson fluid near the stagnation point over a stretching surface coinciding with the plane , the flow being confined to . Two equal and opposite forces are applied along the -axis so that the surface is stretched keeping the origin fixed. A uniform magnetic field of strength is applied in a direction normal to the surface. It is assumed that the velocity distribution far from the surface is given by , where is a constant. The stretching surface has a uniform temperature and the temperature far away from the surface is . The flow configuration is shown in Figure 1.

The rheological equation of state for an isotropic and incompressible flow of a Casson fluid can be written as, see Nakamura and Sawada [42], where is the rate-of-strain tensor, is the Casson's coefficient of viscosity, is the product of the component of deformation rate with itself, is the critical value of the product of the component of the rate-of-strain tensor with itself, is the yield stress of the fluid, and are the velocity components.

The MHD equations for this problem in the boundary layer near the stretching surface are where and are the velocity components along the and directions, respectively, is the density of the fluid, is the electrical conductivity, is the uniform magnetic field along -axis, is the thermal diffusivity, is the specific heat at constant pressure, is the kinematic viscosity, and is the Casson parameter. In writing (4), we have neglected the induced magnetic field since the magnetic Reynolds number for the flow is assumed to be very small. This assumption is justified for flow of electrically conducting fluids such as liquid metals, for example, mercury, liquid sodium, and so forth (Shercliff [43]).

The slip boundary conditions for the problem are where is the velocity slip factor, is the thermal slip factor, and , are constants. For and , no-slip condition can be recovered.

The continuity equation can be satisfied by introducing a stream function such that The momentum and energy equations can be transformed into the corresponding nonlinear ordinary differential equations by the following transformations: where is the independent similarity variable and where , is a positive constant. Using these relations, we get the transformed nonlinear ordinary differential equations from (4) and (5) as where is the dimensionless velocity ratio parameter, is the magnetic parameter, is the Prandtl number, is the Eckert number, and a prime denotes differentiation with respect to the similarity variable . The transformed boundary conditions for and are where is the dimensionless velocity slip parameters is the dimensionless thermal slip parameter. It is to be noted that when the non-Newtonian (Casson) parameter , the Casson flow problem reduces to viscous flow problem.

The physical quantities of interest are the skin friction coefficient and the local Nusselt number , which are defined as where is the shear stress along the stretching surface and is the heat flux from the stretching surface, which are given by Hence using (8) we get where is the local Reynolds number.

3. Analytical Solutions for and

For the explicit analytical solution of (9)–(11) by Homotopy Analysis Method (HAM), the velocity and the temperature distributions can be expressed by the set of base functions in the form: where and are coefficients. Then from (16) and the boundary conditions (11), it is straightforward to choose as our initial approximations for and , respectively. We choose the auxiliary linear operators as with the properties where are constants. Let be an embedding parameter and , denote the nonzero auxiliary linear operators, then the zeroth order deformation equations are subject to the boundary conditions where the nonlinear operators and are defined as When , we have from (20) that When , the zeroth order deformation equations (20)–(22) are equivalent to the original equations (9)–(11) so that we have So as the embedding parameter increases from to , and vary from their initial approximations and to their exact solutions and , respectively. Expanding and in Taylor series with respect to the embedding parameter , we can write where Note that the zero order deformation equations (20) contain nonzero auxiliary parameters and , respectively. Thus and are dependent upon these parameters. Assume that and are so chosen that the series (26) and (27) are convergent at . Hence we have from (25)–(27), Differentiating the zeroth-order deformation equations (20)–(22) times with respect to and then dividing by and finally setting , we have the th order deformation equations as where The general solutions of (30)–(32) are where and denote the special solutions of (30) and (31) and the integral constants are determined by the following boundary conditions: The linear th order deformation equations (30) and (31) can be solved one after the other by using the symbolic computation software MATHEMATICA.

4. Convergence of the HAM Solution

Equation (29) gives the analytical solution of the problem in series form. As pointed out by Liao [35], the convergence region and rate of convergence of the solution series given by HAM depend upon auxiliary parameters and . Hence these auxiliary parameters provide us with a convenient way to adjust and control the convergence region and rate of the solution series. To select appropriate values of these parameters, we display the so called and curves at th order approximations. Figure 2 reveals that there exists horizontal line segment for curves in and for curve in . The obtained series solutions given by (29) converge in the whole region of when . Here we have employed the homotopy-Pade approximation (see Mahapatra et al. [40]) instead of traditional-Pade approximation for the rapid convergence of the series solutions. Table 1 is made to show the convergence of the solutions for and using homotopy-Pade approximation.

5. Results and Discussion

Equations (9) and (10) subject to the boundary conditions (11) are solved analytically for some values of the governing parameters , , , , , , and . In order to verify the accuracy of the present method, we have compared our results with those of Mahapatra and Gupta [4] and Ishak et al. [6] for the skin friction coefficient for different values of (in the case of ,   and ) in Table 2. The comparison shows good agreement for each value of . Therefore, the present results obtained are accurate enough.

This section describes the influence of some important physical parameters on the velocity and thermal profiles. For this purpose, Figures 311 are displayed. Attention is focussed on variations of the velocity slip parameter , the magnetic parameter , the thermal slip parameter , the casson parameter , the velocity ratio parameter , and the Prandtl number on the velocity and the temperature distributions.

Figure 3 shows the influence of on the velocity component . It is observed that is a decreasing function of the slip parameter . Physically this is explained as follows: when slip occurs (for nonzero value of ), the fluid velocity near the sheet is no longer equal to the sheet stretching velocity; that is, a velocity slip exists. With the increase in , such a slip velocity increases. Furthermore, increasing the value of will decrease the flow velocity because under the slip condition, the pulling of the stretching sheet can be only partly transmitted to the fluid. The boundary layer thickness also decreases as the slip parameter increases. The effect of the magnetic field parameter on velocity is depicted in Figure 4. From this figure, it is observed that the transverse magnetic field contributes to the reduction in the velocity profile and boundary later thickness. This is evident from the fact that applied transverse magnetic field produces a body force, to be precise, the Lorentz force, which opposes the motion. The resistance offered to the flow is responsible in decreasing the fluid velocity. The influence of the Casson parameter on the velocity profile is displayed in Figure 5. The figure reveals that as increases, the velocity and the boundary layer thickness decrease. Hence, it is quite obvious that the magnitude of the velocity is greater in casson fluid when compared with the viscous fluid. The effect of the velocity ratio parameter on the velocity field is shown in Figure 6. The figure indicates that fluid velocity increases with the increase of .

To see the variations of , , , (Prandtl number), and (Eckert number) on the temperature profile , Figures 711 are plotted. The effect of the thermal slip parameter on temperature distribution is displayed in Figure 7. As the thermal slip increases, less heat is transferred from the sheet to the fluid and consequently the temperature decreases. Figure 8 depicts the effect of on temperature . It is observed that temperature and the thermal boundary layer thickness increase with increasing . Figure 9 exhibits the temperature profile for different values of the velocity ratio parameter . The figure shows that the temperature and the thermal boundary layer thickness decrease with an increase in . Hence, stronger free stream velocity causes a reduction in the temperature and the thermal boundary layer thickness.

Figure 10 demonstrates the effect of the Prandtl number on the temperature profile . The temperature and the thermal boundary layer thickness decrease with the increase of . Physically this is explained as follows. An increase in the Prandtl number means an increase of fluid viscosity, which causes a decrease in the temperature distribution. An enhancement in the Eckert number results in an increase in temperature , and this observation can easily be visualized in Figure 11.

Tables 3 and 4 show the variation of the reduced skin friction coefficient and the coefficient of the reduced Nusselt number for different values of the physical parameters , , , and considering other parameters fixed. From Table 3, it is observed that is an increasing function of the dimensionless parameters and . Also decreases with the increase of as long as and increases with the increase of for . Table 4 reveals that is a decreasing function of and while an increasing function of .

6. Conclusion

The present work considers the hydromagnetic steady flow and heat transfer of a non-Newtonian casson fluid in the neighborhood of a stagnation point over a stretching surface with partial slip at the boundary. The governing equations are formulated and transformed into a set of ordinary differential equations by similarity transformation. The resulting equations are solved analytically using homotopy analysis method. The main points of this study are as follows.(i)The velocity decreases with increasing values of , , and . (ii)The temperature decreases with increasing values of , , and , and increases with increasing values of and . (iii)The magnitude of velocity is greater in case of Casson fluid when compared with the viscous fluid.(iv)Increase of the velocity slip parameter causes decrease in the flow velocity, and the same qualitative result holds for the thermal slip parameter on the temperature. (v)The present results in limiting cases () are found in excellent agreement with those of Mahapatra and Gupta [4] and Ishak et al. [6].

Acknowledgment

The author thanks the referees for their valuable comments and suggestions which enabled and improved presentation of the paper.