A New Numerical Solution of Maxwell Fluid over a Shrinking Sheet in the Region of a Stagnation Point

The mathematical model for the incompressible two-dimensional stagnation flow of a Maxwell fluid towards a shrinking sheet is proposed. The developed equations are used to discuss the problem of being two dimensional in the region of stagnation point over a shrinking sheet. The nonlinear partial differential equations are transformed to ordinary differential equations by firsttaking boundary-layer approximations into account and then using the similarity transformations. The obtained equations are then solved by using a successive linearisation method. The influence of the pertinent fluid parameters on the velocity is discussed through the help of graphs.


Introduction
Many engineering fluid mechanical problems have been solved using this boundary-layer theory rendering results which compare well with experimental observations-at least as far as Newtonian fluids are concerned. In spite of the success of this theory for Newtonian fluids, an extension of the theory to non-Newtonian fluids has turned out to be a rather formidable task 1-3 . The main difficulty in reaching to a general boundary-layer theory for non-Newtonian fluids lies obviously in the diversity of these fluids in their constitutive behavior, simultaneous viscous, and elastic properties such that differentiating between those effects which arise as a result of a fluid's shear-dependent viscosity from those which are attributable to the fluid's elasticity becomes virtually impossible. But some mathematical models have been proposed to fit well with the experimental observations. The simplest model for the rheological effects of viscoelastic fluids is the Maxwell model where the dimensionless relaxation time is small. However, in some more concentrated polymeric fluids the Maxwell model is also used for large dimensionless relaxation time. Some recent investigations which deal with the flows of Maxwell fluids are given in 4-8 . The viscous flows due to a stretching sheet have attracted the attention of many researchers. This is due to their several applications in polymer processing industries, environmental pollution, biological process, aerodynamic extrusion of plastic sheets, glass fiber production of the boundary layer along a liquid film and condensation process, the cooling, and/or drying of paper and textiles. In 1966, Erickson et al. 9 investigated the problem of heat and mass transfer on a moving continuous plate with suction and injection. In view of above applications, Crane 10 extended Sakiadis 11, 12 study of boundary layer flow for stretching sheet. The literature on the stretching sheet topic is quite extensive and hence cannot be described in detail. However, some most recent works of eminent researchers regarding the flow over a stretching sheet may be mentioned in the studies 13-20 . Literature survey indicates that no attention has been given to the shrinking flow 21-25 for Maxwell fluid for two-dimensional stagnation flow and axisymmetric stagnation flow towards an axisymmetric shrinking over a shrinking sheet. There are few situations like rising shrinking balloon, and so forth, the standard stretching phenomena is not useful therefore, the shrinking phenomena are used. Wang 26 has discussed the stagnation flow towards a shrinking sheet. According to him, solution does not exist for a shrinking sheet in an otherwise still fluid since vorticity could not be confined in a boundary layer. However, in the presence of stagnation flow to contain the vorticity, similarity solution may exist. Mention may be made to the works of 27-33 . The objective of this paper is two-fold: first, to formulate the Maxwell fluid for two-dimensional stagnation flow towards a shrinking sheet; second, to calculate the numerical solution of transformed nonlinear ordinary differential equations via the successive linearisation method. The numerical results for boundary layer Maxwell fluid in the stagnation flow induced by shrinking sheet by means of SLM are yet not available in the literature.

Mathematical Modeling
For the steady, incompressible stagnation flow equations for Maxwell fluid towards a shrinking sheet are of the following form: where u, v, and w are the velocity components along x, y, and z-axis, U ax is the free stream velocity, ν is the kinematic viscosity, and λ is the relaxation time.
The boundary conditions for the present problem are given by where a is the strength of the stagnation flow, b is stretching rate shrinking is b < 0 , and c is the location of the stretching origin.

Mathematical Problems in Engineering 3
We define the following similarity variables 26 : Here primes denote differentiation with respect to η, and β λc is the Deborah number. The boundary conditions are where α b/a, α > 0, and α < 0 correspond to stretching and shrinking sheets, while α 0 correspond to planar stagnation flow towards a stationary sheet. Moreover, when α 1, the flow will be with no boundary layer.

Solution Method
In this section the implementation of the successive linearisation method SLM is solving the governing boundary value problems. The SLM see, e.g., 34-38 for details is based on transforming the governing nonlinear boundary value problem into an iterative scheme made up of linear differential equations which are subsequently solved using analytical or numerical methods wherever possible. Details of the formulae for implementing the SLM in one dimensional equations are given in 36 .
In applying the SLM on 2.5 , we set where f i are unknown functions that are obtained by iteratively solving the linearized version of the governing equations assuming that f i 0 ≤ m ≤ i − 1 are known from previous iterations. The algorithm starts with an initial approximation f 0 η which is chosen to satisfy the boundary conditions 2.7 . A suitable initial guess in this example is To investigate the SLM application on 2.5 , we write the equation as

3.4
Substituting 3.1 in 2.5 and linearising the resulting equation, we obtain the following governing equation: subject to the boundary conditions

3.7
We solve 3.5 using the Chebyshev spectral collocation method. To allow for numerical implementation of the spectral method, the physical region 0, ∞ is truncated to 0, L , where L is chosen to be sufficiently large. The truncated region is further transformed to the space −1, 1 using the transformation As with any other numerical approximation method some sort of discretization is introduced to the interval −1, 1 . To this end, we choose the Gauss-Lobatto collocation points see, e.g., 39-42 to define the nodes in −1, 1 as In general, the derivative of order p for the function f ξ can be expressed by The entries of D can be computed in different ways see, e.g., 39-42 . In this work we use the method proposed by Trefethen 42 in the cheb.m Matlab m-file. Thus, applying the spectral method, with derivative matrices on linearised equations 3.5 and 3.6 leads to the following linear matrix system: with the boundary conditions After imposing the boundary conditions 3.13 , the solutions for f i are obtained by iteratively solving the system 3.12 starting from the initial approximation 3.2 . Once the solution for f η has been found from solving 3.13 , the Chebyshev pseudospectral method is directly applied on 2.6 which is now linear. This gives the following: with the boundary conditions h ξ N 1, h ξ 0 0, 3.17

Results and Discussion
In this section we give results obtained by the successive linearisation method for selected values of important parameters of the Maxwell fluid over a shrinking sheet in the region of a stagnation point. To check the accuracy of the SLM employed in the current study, a comparison of the initial values is made with published results from literature. In particular, comparison is made with the skin friction coefficient f 0 results and results for h 0 of Wang 26 for the special case when β 0. Tables 1 and 2 give the results for the friction coefficient f 0 and h 0 for α > 0 and α < 0, respectively. It can be seen from these tables that the present SLM results are in excellent agreement with those of Wang 26 . We also observe that SLM results converge rapidly with full convergence to the numerical results of 26 reached after only 3 or 4 iterations. This confirms the validity of the SLM method for solving the problem under study. For some negative values of α it was found that multiple solutions exist. To generate the multiple solution different initial approximations were used as first guesses of the SLM iteration scheme. The multiple solutions for f 0 and h 0 are resolved when α −1.15 and illustrated in Table 2 for the case β 0. Again, good agreement with the results of 26 is observed. The presence of dual solutions has also been observed in related investigations stagnation flow in the presence of shrinking, see, e.g., 26-33 . The effect of the shrinking parameter α α < 0 of f η and h η is depicted on Figure 1. We observe that shrinking of the sheet decelerates the fluid flow. Figure 2 shows the effect of the stretching parameter α α > 0 on f η and h η . We observe that the boundary layer thickness is decreased when values of α increase. Physically, this can be explained as follows: for fixed values of a corresponding to the stretching of the sheet, an increase in b implies an increase in straining motion near the stagnation region resulting in the increased acceleration of the external stream, and this leads to the thinning of the boundary layer.

Mathematical Problems in Engineering
The influence of the Deborah number β, for both the shrinking α < 0 and stretching α > 0 , can be seen from Figures 3 and 4. As shown in these figures, the velocity decreases with the increase of the Deborah number β or the relaxation time.

Concluding Remarks
The two-dimensional equations for Maxwell fluid are derived in this paper. The boundary layer in the region of stagnation point over a two-dimensional shrinking sheet is discussed for Maxwell fluid. The similarity transformations are used to transform the partial differential equations to ordinary ones, and hence a numerical solution is obtained using the SLM. The novel technique was compared against previous studies and excellent agreement was observed. The study also confirmed the existence of a dual solution. We also observed the relatively rapid convergence of the SLM. The results are presented graphically and the effects of the parameters are discussed.