Abstract

This study aims to investigate the unsteady boundary-layer flow of a viscoelastic non-Newtonian fluid over a flat surface. The plate is suddenly jerked to move with uniform velocity in a uniform stream of non-Newtonian fluid. Purely analytic solution to governing nonlinear equation is obtained. The solution is highly accurate and valid for all values of the dimensionless time . Flow properties of the viscoelastic fluid are discussed through graphs.

1. Introduction

Boundary-layer flows over flat surfaces are fundamental for the understanding of aerodynamical properties of the flow such as wall friction and the dynamical drag. Due to the boundary-layer theory proposed by Prandtl in 1904 [1], it became possible to calculate the drag on a sphere accurately. Blasius [2] investigated the boundary-layer flow past a flat plate at zero incidence. The wall friction in Blasius flow was calculated to be approximately. The same two-dimensional boundary-layer flow was studied by Sakiadis [3] over a moving wall in a still fluid. Sakiadis reported the value of skin friction at the wall to be for his case. Sakiadis also predicted that the drag force for his flow problem was greater than the drag force for the Blasius flow [2]. The theoretical results of Sakiadis [3] were experimentally confirmed by Tsuo et al. [4]. In both studies [2, 3], either fluid or plate is assumed to be at rest. However, the most practical situation might be the one in which both the fluid and the plate are moving. Such a flow situation is more practical in studying the aerodynamical properties of the flow. Keeping this fact in mind, Klemp and Acrivos [5] considered two-dimensional boundary-layer flow of a viscous fluid over a flat surface moving in a stream of constant velocity. Hussaini et al. [6] determined similarity solution of the boundary layer equation with upstream moving wall. Later on in 2003, Fang [7, 8] studied the similarity solution of boundary-layer flow and heat transfer for steady case. In [7, 8], Fang studied both cases, namely, when the plate is moving in the direction of free stream or in opposite direction to the free stream. Recently, Mehmood et al. [9, 10] considered unsteady boundary-layer flow over an impulsively started moving plate in a free stream with parallel and antiparallel motion.

Boundary-layer flows over moving surfaces find important industrial applications such as in the manufacturing of food and paper, plastic sheet extrusion, application of coating of paints layers on surfaces, and many other activities (see for instance [11, 12]). In such industrial applications, the fluid is observed to have non-Newtonian character. Due to these important industrial applications, the study of non-Newtonian fluid flows over moving surface needs attention. In 1969, Fox et al. [13] investigated the non-Newtonian flow over a moving surface using the power-law model. The flow of a power-law fluid over a moving plate in a parallel free stream was investigated by Hassanien [14]. Further, Hassanien [15] investigated the boundary-layer flow over a moving plate in a free stream of second-grade fluid with heat transfer analysis.

In all the abovementioned studies for non-Newtonian fluids, the authors considered steady flow. However, there are many practical applications in which the flow is essentially unsteady. For unsteady problems governed by nonlinear equations, it is very difficult to obtain an analytic solution valid for all time. Perturbation methods are sometimes used to get analytic solution valid for small time [1618]. Currently, Liao [19, 20], Xu et al. [21], Cheng et al. [22], Wang [23], Abbasbandy et al. [24], Xu and Liao [25], Mehmood et al. [9, 10, 26, 27], and Munawar et al. [28] obtained purely analytic solution for unsteady flows by homotopy analysis method [29]. In [9, 10, 1928], the authors have shown their results to be valid uniformly for all time. Homotopy analysis method is a powerful analytic technique (introduced by Liao [29]) for highly nonlinear problems. The technique has been widely used by a number of researchers in nonlinear problems arising in science and engineering [3041].

In the present study, we extend the work of Mehmood and Ali [9] for the non-Newtonian case. The fluid considered is of second-grade type and the problem is solved by homotopy analysis method. The objective of this study is twofold: firstly to report a purely analytic solution to the considered unsteady problem and secondly to investigate the viscoelastic effects on the flow characteristics. The solution is highly accurate and is uniformly valid for all time in the whole spatial domain. The accuracy and convergence of present analytic solution are discussed in detail. The outlines of the paper are as follows.

Section 2 contains the mathematical formulation and the HAM solution of the problem. The issue of convergence and accuracy of HAM results is also discussed in Section 2. Section 3 consists of graphical representation of results and their discussion. Section 4 is reserved for concluding remarks.

2. Formulation of the Problem

Consider the unsteady incompressible flow of second-grade fluid over an infinite plate at . The fluid occupies the half space and at infinity it is assumed that fluid is flowing with the constant free-stream velocity . Initially, fluid is at rest; then at time , suddenly plate starts its motion with the constant velocity , where is the ratio of plate velocity to the free-stream velocity. For an incompressible homogeneous fluid of second-grade type, the Cauchy stress tensor is related to the deformation field through where is the pressure, is the identity tensor, is the dynamic viscosity, and are the viscoelastic parameters, and the kinematical tensors and are given by where is the material time derivative and is velocity gradient. This rheological model was first introduced by Rivlin and Ericksen [42].

Experimental data available for the large number of viscoelastic fluids suggests that, in order to satisfy the thermodynamical analysis, some restrictions must be put on the signs and magnitudes of the material parameters [43]

Under the above assumptions and conditions, the boundary layer equations governing the unsteady laminar flow of an incompressible viscoelastic fluid due to an impulsive motion of plate are given by where and are the components of the stress tensor given by

Using (6) and (7) in (5) and assuming the constant pressure, we have subject to the boundary conditions (when ) in which is the kinematic viscosity, is the fluid density, and and are the velocity components in the - and -directions, respectively. The initial condition (when ) is given by

Introducing the similarity transformations [9] the governing equation (8) readily transforms to with boundary conditions where is the dimensionless stream function, is the local Deborah number, and denotes differentiation with respect to . Notice that, for , the above equation reduces to that of the viscous fluid [9].

3. HAM Solution

To solve the problem (12) and (13) analytically, we use the well-known analytic technique homotopy analysis method. According to boundary conditions (13), it is convenient to express by the base functions in the following form: where are coefficients involved in the solution series. According to the solution expression (15) and boundary conditions (13), we choose the initial approximation and the linear operator satisfying the property where , and are constants and is the controlling parameter to be adjusted. From (12), we define the nonlinear operator

Taking as nonzero auxiliary parameter, we construct the zero-order deformation equation [29] subject to the boundary conditions where is the embedding parameter and obviously, when and , (20) has the solutions respectively. Thus, the variation of from to is continuous deformation of from to . Expanding in Taylor’s series with respect to , and using (22), we have where

Assume that is chosen so properly such that the series (24) is convergent at . Using (22) and (23), we can write

Differentiating the zero-order deformation equations (20) and (21) -times with respect to , at and then dividing by , we obtain the th-order deformation equation subject to the boundary conditions where The system of linear nonhomogeneous equations (27)-(28) can be solved up to high order of approximation with the help of the symbolic computation software such as Mathematica.

Let represent a special solution; then the general solution can be determined as follows: in which are constants of integration which can be determined with the help of boundary conditions (28) as under In this way, the complete solution for the present problem can be written in the form of an infinite series of functions; that is,

To prove that the series is an approximate solution of the system (12) and (13), it is necessary to show the convergence of the solution series (32). As mentioned by Liao [29], the convergence of the solution series strongly depends upon the auxiliary parameter once the initial guess and the linear operator have been selected. A rough estimate for the allowed regions of the values of , the so-called -curve, is important. We have plotted the -curve for our present problem in Figure 1. It is worth mentioning here that, in some problems with strong nonlinearity, the values of strongly depend upon the involved physical parameters (see for instance [44]). In Figure 1, it can be seen that the intervals of allowed values of are shrinking with increasing . Further, it can also be seen that, for higher values of , the interval of allowed values of shifts towards zero. However, to search a more appropriate value of , it is useful to calculate the residual errors. In Figure 2, we have plotted the residual error against for fixed values of the parameters involved. From Figure 2, it is clear that the error for is minimum at which is at the 20th order of approximation. In Figure 3, we have plotted the errors graph against the space variable also. It shows that in the boundary layer region the error fluctuates and dies out to be zero as one moves to the outer region. Though the error fluctuates in the boundary layer region, it remains in the acceptable limits. In order to prove the convergence of the solution series, it is recommended that the corrections to solution must become negligible in the succeeding orders of approximation. We apply the homotopy Padé approximation in order to accelerate the convergence of solution series. In Table 1, we report the HAM solution at different orders of Padé approximation. Clearly, there are no corrections up to four decimal places after the 8th order of Padé approximation. This proves the convergence and accuracy of HAM solution.

4. Graphical Illustration and Discussion of Results

To understand the physics of flow phenomenon, we have investigated the solution through graphs. In Figure 4, the longitudinal component of velocity is plotted for different values of the viscoelastic parameter . Clearly, strong non-Newtonian behavior of the fluid results in large skin friction at the solid wall. However, large values of help in reducing the thickness of the boundary layer. In Figure 5, the velocity function is plotted for different values of the parameter provided . The velocity increases at the plate which causes the skin friction at the plate to reduce. This is due to the fact that, for , the plate and the free stream progress in the same direction, but, for when the plate and free stream have opposite direction of progression, the skin friction increases considerably. Finally, the present HAM solution is shown to be uniformly valid for all time in the whole spatial domain in Figure 6. Figure 6 depicts that the velocity is very small at initial time and with the passage of time the flow develops and reaches its steady state due to the motion of the plate. From the figure, it can also be observed that the steady state is reached at (roughly). This shows that the vorticity diffusion takes place in a very short interval of time.

The effect of viscoelastic parameter on the boundary layer thickness is shown in Table 2. It is observed from the table that the thickness of the boundary layer, , has a declining behavior with an increase in . For the low value of , this drop in is less as compared with the high value of . Such kind of result can be expected because at low values of our governing equation behaves like that of Newtonian fluid [45]. At the low values of , shear stress produced in the flow is not so large so the boundary layer thickness increases. By increasing , the shear stress will grow more and more and results in the reduction of boundary layer thickness. We can observe such kind of results from the theory of polymeric liquid as well as experimental observations.

5. Concluding Remarks

In this study, we have considered the unsteady boundary-layer flow of a viscoelastic fluid over an impulsively started moving flat plate. Fluid at infinity was assumed to be flowing with a uniform free stream velocity. The governing nonlinear equations of an incompressible second-grade fluid are modelled using the similarity transformations. The resulting nonlinear problem is solved analytically with the help of HAM. The influence of various parameters of interest on the velocity profile is graphically illustrated. We find the following observations.(i)It is noticed that the boundary layer thickness decreases by increasing second-grade parameter .(ii)It is observed that as the time passes by, the unsteady velocity reaches the steady state showing that our analytic solution is valid for all the time throughout the spatial domain.(iii)It is observed that increasing (for ) reduces the drag at the plate whereas for the skin friction at the plate increases by increasing .

Conflict of Interests

The authors declare that they have no conflict of interests regarding the publication of this paper.