Complex Boundary Value Problems of Nonlinear Differential Equations 2014
View this Special IssueResearch Article  Open Access
Different Approximations to the Solution of UpperConvected Maxwell Fluid over a Porous Stretching Plate
Abstract
In the present paper, we consider an incompressible magnetohydrodynamic flow of twodimensional upperconvected Maxwell fluid over a porous stretching plate with suction and injection. The nonlinear partial differential equations are reduced to an ordinary differential equation by the similarity transformations and taking into account the boundary layer approximations. This equation is solved approximately by means of the optimal homotopy asymptotic method (OHAM). This approach is highly efficient and it controls the convergence of the approximate solutions. Different approximations to the solution are given, showing the exceptionally good agreement between the analytical and numerical solutions of the nonlinear problem. OHAM is very efficient in practice, ensuring a very rapid convergence of the solutions after only one iteration even though it does not need small or large parameters in the governing equation.
1. Introduction
The flow of nonNewtonian fluids has been analyzed by numerous researchers, because, in practical applications, nonNewtonian fluids are more appropriate than Newtonian fluids. For such fluids, the NavierStokes theory becomes inadequate and obviously there are several constitutive equations that have been proposed to various types of nonNewtonian fluids. Examples of the flow of nonNewtonian fluids occur in a large variety of applications: synthetic fibres, drilling muds, food processing, plastic polymers, and so on. There are numerous constitutive equations that have been proposed to various types of nonNewtonian fluids; a majority of nonNewtonian fluids models are concerned with the fluids of grade two or three. One type of fluids in which the relaxation type phenomena can be considered is known as Maxwell model. Some investigations in this field are made by the mathematicians, engineers, physicians, and computer scientists. Sakiadis [1, 2] first studied various aspects of the stretching problem: the flow due to a semiinfinite horizontally moving wall in an ambient fluid. PhanThien [3] considered the plane and axisymmetric stagnation flows in a Maxwell fluid, using the shooting and boundary element method, like Zheng et al. [4]. Sadeghy et al. [5] considered the problem of hydrodynamic Sakiadis flow of an upperconvected Maxwell fluid over a rigid plate moving steadily in an otherwise quiescent fluid. Homotopy analysis method is used by Hayat et al. [6] to solve nonlinear differential equation of the upperconvected Maxwell fluid over a porous stretching plate. Also, the same method is considered by Hayat and Sajid [7], Abbas et al. [8, 9], and Hayat et al. [10]. Dual solutions in mixed convection flow near stagnation point on a vertical porous plate have been presented by Ishak et al. [11]. Sahoo [12] investigated the effects of partial slip in the MHD flow and mass transfer of an electrically conducted second grade fluid past an axisymmetric stretching sheet.
Taking into considerations these analyses, the objective of the present paper is to propose an accurate procedure to nonlinear differential equation of the magnetohydrodynamic flow of an upperconvected Maxwell fluid over a porous stretching plate using OHAM. A version of the OHAM is applied in this study to derive highly accurate analytical expressions of the solutions. Our procedure does not depend upon any small or large parameters, contradistinguishing from other known methods in literature. The main advantage of this approach is the control of the convergence of approximate solutions in a very rigorous way. A very good agreement was found between our approximate solutions and numerical solutions, which proves that our method is very efficient and accurate.
2. Equation of Motion
If we consider the steady incompressible twodimensional flow of an upperconvected Maxwell fluid over a porous stretching plate, then the constitutive equation for the Cauchy stress tensor is given by such that the extra tensor satisfies where is the viscosity, is the relaxation time, and the first RivlinEricksen tensor is defined as follows: in which is the material time derivative and is the velocity.
For the magnetohydrodynamic Sakiadis flow, by imposing a uniform magnetic field along the direction and neglecting the induced magnetic field, the equations which govern the steady flow can be written as where , are the velocity components, is the pressure, is the density, is the electrical conductivity, and , , , are the components of the extra tensor .
Using the boundary layer approximations [5, 13] the flow in the absence of the pressure gradient is governed by (4) and where is the kinematic viscosity of the fluid.
The relevant initial/boundary conditions for the flow problem are in which is the stretching rate and is the suction velocity and is the injection velocity.
Introducing the stream function such that then (4) is identically satisfied and (7) becomes
Introducing the similarity transformations equation (9) becomes
Substituting (12) into (10), it obtains the governing equation in the following form: where and .
Now, the initial/boundary conditions (8) are with , where corresponds to suction velocity and for injection velocity.
In what follows, the nonlinear differential equation (13) with initial/boundary conditions (14) can be solved analytically using OHAM.
3. Basic Ideas of the Optimal Homotopy Asymptotic Method
Equation (13) with initial/boundary conditions (14) can be written in a more general form as follows: where is a given nonlinear differential operator depending on the unknown function , subjected to the initial/boundary condition
Let be an initial approximation of and an arbitrary linear operator such as
It should be emphasized that this linear operator is not unique.
If denotes an embedding parameter and is an analytic function, then we construct a homotopy [14â€“18]: with the properties where is an arbitrary auxiliary convergencecontrol function depending on variable and on a number of arbitrary parameters unknown now and will be determined later.
Let us consider the function in the form
By substituting (21) into equation obtained by means of homotopy (18), and, then, equating the coefficients of and , we obtain
From (23), we obtain the governing equation of given by (17) and the governing equation of ; that is, where we find the following expression for the nonlinear operator: where the functions and , , are known and depend on the function and also on the nonlinear operator, being a known integer number.
In this way, taking into account (20), from (21), for , we obtain the firstorder approximate solution which becomes
It should be emphasized that and are governed by the linear equations (17) and (24), respectively, with boundary conditions that come from the original problem. It is known that the general solution of nonhomogeneous linear equation (24) is equal to the sum of general solution of the corresponding homogeneous equation and of some particular solutions of the nonhomogeneous equation. However, the particular solutions are readily selected only in the exceptional cases.
In what follows we do not solve (24), but, from the theory of differential equations, taking into considerations the method of variation of parameters, Cauchy method, method of influence function, the operator method [19], and so on is more convenient to consider the unknown function , in the form where within expression of appear linear combinations of some functions , some terms which are given by the corresponding homogeneous equation, and the unknown parameters , . In the sum, appear an arbitrary number of terms.
For instance, if , then we can choose . Similarly, if , then we can choose . In the case when , we can choose or . We have large freedom to choose the value of . We cannot demand to be solutions of (24) but , given by (26) with given by (27), are the solutions of (15). This is an underlying idea of our method. The convergence of the approximate solution given by (26) depends upon the auxiliary functions , . There are many possibilities to choose these functions . We try to choose so that within (27) the term is of the same shape with the term given by (25). The firstorder approximate solution also depends on the parameters , . The values of these parameters can be optimally identified via various methods, such as the leastsquare method, the Galerkin method, the collocation method, and the Ritz method. The first option should be minimizing the square residual error: where the residual is given by
The unknown parameters can be identified from the conditions:
With these parameters known (called optimal convergencecontrol parameters), the firstorder approximate solution given by (26) is well determined.
It should be emphasized that our procedure contains the auxiliary functions , , , which provides us with a simple way to adjust and control the convergence of the approximate solutions. It is very important to properly choose these functions which appear in the construction of the firstorder approximation.
4. Different Approximations to the Solution of UpperConvected Maxwell Fluid with OHAM
In what follows we apply our procedure to obtain approximate solutions of (13) and (14). For this purpose, we choose the linear operator of the following form:
We mention that the linear operator is not unique. Also, we have freedom to choose where is an unknown positive parameter and will be determined later. The initial approximation can be obtained from (17) with initial/boundary conditions
Equation (17) with the linear operators (31) and (32) has the following solution: while (17) with the linear operators (33) has the following solution:
In what follows we consider only the linear operator given by (31) and the initial approximation given by (35).
The nonlinear operator corresponding to nonlinear differential equation (13) is defined by
By substituting (35) into (37), it holds that
Comparing (25) and (38), one can get
The function given by (27) becomes where we have the freedom to choose a lot of possibilities for the unknown functions as follows:
Substituting (41) into (40), we have
The firstorder approximate solution given by (26) is obtained from (35) and (42):
In this way, we can find other solutions:
It is clear that we can obtain many other solutions.
5. Numerical Results
We illustrate the accuracy of our procedure for different values of the coefficients , , and . Also, we represent graphically the behavior of the functions and and we compare results obtained through our procedure with numerical results.
First, we consider , , and which corresponds to suction velocity. The optimal convergencecontrol parameters are determined by means of the leastsquare method and are as follows: and so on.
The firstorder approximate solutions obtained by means OHAM, for the above cases ()â€“(), become as follows:
In the last case, we consider , , and for injection velocity. The optimal convergencecontrol parameters are as follows: and so on.
The firstorder approximate solutions can be written for only three cases in the following forms:
In Figures 1 and 2, a comparison between the firstorder approximate solutions and , respectively, for the suction velocity in the case () and numerical results is plotted, while, in Figures 3 and 4, a comparison between the firstorder approximate solutions and , respectively, for the injection velocity in the case () is plotted.
In Tables 1, 2, 3, 4, and 5, we present a comparison between the firstorder approximate solutions (46)â€“(62), respectively, with numerical results for some values of variable .





In Tables 6, 7, and 8, we present the relative error between the firstorder approximate solutions (46)â€“(62) and numerical results; for some values of variable .