Abstract
In this paper the thin film flow of an Oldroyd 6-constant fluid on a vertically moving belt is investigated. The basic equation of a non-Newtonian fluid in a container with a wide moving belt which passes through the container moving vertically upward with constant velocity, is reduced to an ordinary nonlinear differential equation. This equation is solved approximately by means of the Optimal Homotopy Asymptotic Method (OHAM). The solutions take into account the behavior of Newtonian and non-Newtonian fluids. Our procedure intended for solving nonlinear problems does not need small parameters in the equation and provides a convenient way to control the convergence of the approximate solutions.
1 Introduction
It is known that in practical applications the behavior of viscoelastic fluids cannot be represented by that of Newtonian fluids. The inadequacy of the classical Navier-Stokes model to describe the fluids which do not obey the Newtonian postulate that the stress tensor is multiple of the shear strain, has led to the development of several models of non-Newtonian fluids. During the past few decades, various constitutive equations have been proposed for these non-Newtonian fluids. The differential equations that describe non-Newtonian fluids are in general nonlinear and highly complicated, but there is not a single constitutive equation available in the literature by which one can study the flow behavior of viscoelastic fluids. Thin film flows have attracted the attention of numerous researchers, due to their wide applications in engineering and science. Such applications include for example microchip production, the linings of mammalian lungs and so on.
Last few years, many researchers have focused on different problems related to such non-Newtonian fluids. The Oldroyd fluid has acquired a special status amongst the many fluids of the rate type, as it included as special cases the classical Newtonian fluid and the Maxwell fluid [1]. This fluid takes into account the properties of stress relaxation and retardation. The two-dimensional steady flow of an incompressible Maxwell fluid in a corner formed by two planes, one of which is sliding past the other at a certain angle, was first investigated by Strauss [2]. Hancock et al have investigated the effects of inertial forces by constructing regular perturbation series for the stream function in [3]. Rajagopal and Bhatnagar [4] presented two solutions for the flow of an Oldroyd-B fluid: the flow past an infinite porous plate and the longitudinal and torsional oscillations of an infinitely long rod of finite radius. The two-dimensional steady slow flow of an Oldroyd 6-constant fluid between intersecting planes, one of which is fixed and the other moving, has been analyzed by Baris [5]. The effects of the non-Newtonian parameter on the flow pattern are carefully delimited. There is, unlike the case of Newtonian fluid, a secondary flow near the corner. Such flow are of considerable practical interest as they locally appear in a cylinder with a moving piston, or at the edge of a blade used to scrape up liquid from a surface.
Very recently, Hayat et al [6] studied some steady unidimensional flow of an Oldroyd 8-constant MHD fluid in bounded domain. The MHD solutions for a Newtonian fluid as well as those corresponding to the Oldroyd 3 and 6-constant fluids, a Maxwell fluid and a second grade one, appear as limiting cases of the obtained solutions. Hayat et al [7] and Wang et al [8] considered Couette and Poiseuille flows of an Oldroyd 6-constant fluid with magnetic field. Sajid et al [9] investigated the problem of wire coating by withdrawal from a bath of MHD Oldroyd 8-constant fluid, using homotopy analysis method, and Khan et al [10] also discussed an Oldroyd 8-constant MHD fluid but between coaxial cylinders subject to partial slip at the boundaries. The Oldroyd 6-constant fluid and fluid-slip on the three nonlinear boundary value problems are studied by Hayat et al in [11]. Ellahi et al [12] investigated an Oldroyd 8-constant fluid over a suddenly moved plate and Sajid and Hayat [13] obtained an exact solution for thin film flow of an Oldroyd 8-constant fluid. Using Fourier sine and Laplace transforms, Vieru et al [14], established exact solutions for the unsteady flow of an incompressible generalized Oldroyd-B fluid due to an infinite constantly accelerating plate. Hayat et al [15] obtained the exact solution of a thin film flow of an Oldroyd 6-constant fluid and Ellahi et al [16] for flows of an Oldroyd 8-constant fluid with nonlinear slip conditions. The energetic balance for the motion of an Oldroyd-B fluid due to an impulsively moved plate is presented by Fetecău et al in [17]. 3D flow of a generalized Oldroyd-B fluid due to a constant pressure gradient between two side walls perpendicular to a plate is presented by Zheng et al in [18]. Shah et al [19] used the optimal homotopy asymptotic method in the study of the wire coating in a pressure type die with the bath of Oldroyd 8-constant fluid with pressure gradient. The nonlinear differential equations of the steady thin film flow of an Oldroyd 8-constant fluid is solved by Siddiqui et. al.[20] applying the variational iteration method and Adomian decomposition method. In [21], the slip-effect on the rotating flows on an Oldroyd-B fluid in a porous space is examinated by Hayat et al. The fluid is permeated by a transverse magnetic field. The governing equations include the modified Darcy’s law for an Oldroyd-B fluid.
Analytical solutions to nonlinear differential equations play an important role in the study of the flow of an Oldroyd different types fluids, but is difficult to find these solutions in the presence of strong nonlinearity. Many new approaches have been proposed to find and develop approximate solutions of nonlinear differential equations. Perturbation methods have been applied to determine approximate solutions to weakly nonlinear problems [22]. But the use of perturbation theory in many problems is invalid for parameters beyond a certain specified range. Other procedures have been proposed such as: the Adomian decomposition method [23], some linearization methods [24],[25], various modified Lindsedt-Poincaré methods [26], the optimal homotopy perturbation method [27],[28] and so on.
In this paper we consider flows of the Oldroyd 6-constant fluid over a moving belt. A version of the optimal homotopy asymptotic method (OHAM) is applied in this study to derive highly accurate analytical expressions of solutions. Our procedure does not depend upon any small or large parameters, contradistinguishing from other known methods. 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 Governing equations
The Cauchy stress T is:
where −pI - is the indeterminate part of the stress due to the constraint of incompressibility and the extra stress tensor S is defined by
where λ1, λ2, λ3, λ4, λ5, λ6, λ7, μ are material constants, A1 - is the first Rivlin-Ericksen tensor defined by
where L is the spatial velocity gradient, L = gradV, V being the velocity vector. The contravariant derivative D/Dt in terms of the material derivative d/dt is defined by
It should be noted that when λi = 0, i = 1, . . . 7, the model reduces to the classical linearly viscous model of a Newtonian fluid. When λ3 = λ4 = λ5 = λ6 = λ7 = 0 it reduces to a 3-constant model of an Oldroyd-B fluid, while for λ2 = λ3 = λ4 = λ5 = λ6 = λ7 = 0 it reduces to a Maxwell model. For λ1 = λ3 = λ5 = λ6 = λ7 = 0 the model describes a second grade fluid and for λ5 = λ6 = λ7 = 0 it reduces to a Johnson-Segalman model. For λ6 = λ7 = 0 the model reduces to the Oldroyd 6-constant.
Neglecting thermal effects and body forces, the equations of motion of an incompressible fluid are:
where ρ is the density and t is the time.
If the extra tensor and velocity are respectively:
From Eq. (7), the continuity equation (5) is satisfied identically and the momentum equation (6) becomes:
where
and modified pressure p̂ becomes:
From Eqs. (9) and (8) we obtain after non-dimensionalization:
where α1 and α2 are given in terms of the nondimensional material constants and m is a nondimensional gravity parameter. The boundary conditions are
3 Basic ideas of optimal homotopy asymptotic method
Eq. (14) with boundary conditions (15) can be written in a more general form
where N is a given nonlinear differential operator depending on the unknown function Φ(x), subject to the boundary conditions
Let Φ0(x) be an initial approximation of Φ(x) and L an arbitrary linear operator such as
where f(x) is an arbitrary continuous function.
We remark that this operator L is not unique.
If p ∈ [0, 1] denotes an embedding parameter and F is a function, then we propose to construct a homotopy [32],[29],[30],[31]:
with the properties
where H(x, Ci) ≠ 0 is an arbitrary auxiliary convergence-control function depending on variable x and on a number of auxiliary parameters C1, C2, …, Cm which will be defined later.
Let the function F be in the form
By substituting Eq. (22) into equation obtained by means of the homotopy (19)
and equating the coefficients of like powers of p, we obtain the governing equation of Φ0(x) given by Eq. (18) and the governing equation of Φ1(x, Ci), Φ2(x, Ci) and so on. If the series (22) is convergent at p = 1, one has
But in particular we consider only the first-order approximate solution
and the homotopy (19) in the form
Equating only the coefficients of p0 and p1 into Eq. (26), we obtain the governing equation of Φ0(x) given by Eq. (18) and the governing equation of Φ1(x, Ci) i.e.
It should be emphasized that Φ0(x) and Φ1(x, Ci) are governed by the linear Eqs. (18) and (27), respectively with boundary conditions that come from the original problem, which can be easily solved. The convergence of the approximate solution (25)depends upon the auxiliary convergence-control function H(x, Ci). There are many possible forms of the function H(x, Ci). Basically the shape of H(x, Ci) must follow the terms appearing in the Eq. (27). Therefore, we try to choose H(x, Ci) so that in Eq. (27), the product H(x, Ci)N(Φ0(x)) be of the same shape with N(Φ0(x)).
Now, substituting Eq. (25) into Eq. (16) results in the following residual
At this moment, the first-order approximate solution given by Eq. (25) depends on the parameters C1, C2, …, Cm and these parameters can be optimally identified via various methods, such as the least square method, the Galerkin method, the Kantorowich method, the collocation method or by minimizing the square residual error
where a and b are two values depending on the given problem. The unknown parameters C1, C2, …, Cm can be identified from the condition
With these parameters known (namely convergence-control parameters), the first-order approximate solution (25) is well-determinate.
4 Application of OHAM to thin film flow of an Oldroyd 6-constant fluid
In what follows we apply our procedure to obtain approximate solutions of Eqs. (14) and (15). For this purpose, we choose the linear operator in two cases.
4.1 Case 1.
In the first case we suppose that the linear operator has the form
where prime denotes derivative with respect to x, and Φ(x) = v(x) then the nonlinear operator becomes
The initial conditions (15) become
From (33) we obtain two subcases:
– 4.1.a
In the first subcase we consider the following initial conditions
The Eq. (18) can be written as
Its solution is
The second function v1(x) is obtained from Eq. (27):
where
Substituting Eq. (36) into (38), the nonlinear operator N becomes
where
The expression (39) is a polynomial one and the auxiliary function H1(x, Ci) from Eq. (37) is chosen such that the product H1N be of the same form [32]. Therefore we choose the auxiliary function H1(x, Ci) in the form
where C1, C2, …, C5 are unknown parameters.
There are many possible forms of the function H1(x, Ci) which appears into Eq. (37).
For example we can consider the following possibilities:
or
and so on.
Inserting Eqs. (41) and (39) into Eq. (37) we obtain equation:
By solving Eq. (42) we obtain
The first-order approximate solution of Eqs. (14) and (34) is obtained from Eq. (24):
where v0(x) and v1(x) are given by Eqs. (36) and (43) respectively.
The residual given by Eq. (28) becomes
with v̄ obtained from Eq. (44).
– 4.1.b
In the second subcase, from Eq. (33) we have
The Eq. (18) becomes
with the solution
For the expression (48) of v0(x), the nonlinear operator N becomes
where
Eq. (27) for the function v1(x) is
where the auxiliary function H1(x, Ci) is chosen in the form
Substituting Eqs. (49) and (52) into Eq. (51), and thereafter solving this equation, we obtain the solution
The first-order approximate solution of Eqs. (14) and (46) is
where v0(x) and v1(x) are given by Eqs. (48) and (53) respectively.
4.2 Case 2.
In the second case, the linear operator can be chosen in the form:
where λ is an unknown parameter, and the nonlinear operator N is given by Eq. (32).
From Eq. (15)2 or Eq. (33) we can consider the following two subcases:
– 4.2.a
In the first subcase, the initial conditions are
The Eq. (18) becomes
The solution of Eq. (57) has the form
Substituting Eq. (58) into Eq. (32), the nonlinear operator N(v0) becomes:
where
The second function v1(x) is obtained from Eq. (27):
The expression (59) of the nonlinear operator N(v0) is an exponential one and therefore, the auxiliary function H2(x, Ci) can be chosen in the form
where C1, C2, and C3 are unknown parameters. Also we can use the following expressions for the auxiliary function:
or
and so on.
Now, substituting Eqs. (62) and (59) into Eq. (61), we obtain the equation
From Eq. (63) we can obtain the solution v1(x), such that the first-order approximate solution of Eqs. (14) and (56) in this subcase, can be written in the form
where v0(x) is given by Eq. (58) and
In the last subcase, the initial conditions are
The Eq. (18) has the form
and has the solution
By means of Eq. (68), the nonlinear operator N(v0) becomes:
where
For the equation in v1(x) we choose the auxiliary function H2(x, Ci) in the form
The Eq. (27) can be written, by means of Eqs. (71) and (69) as
The first-order approximate solution in the last subcase has the form
where
5 Numerical results and discussions
We illustrate the accuracy of our procedure for different values of the coefficients α1, α2 and m and for every the auxiliary function H1 and H2. Also we represent graphically the variation of velocity profile of the belt problem for an Oldroyd 6-constant fluid. In all cases and subcases, the values of the parameters C1, C2, … and λ are determined using the least square method and the Wolfram Mathematica 6.0 software. In this way the parameters Ci,, i = 1, 2, . . . (namely the convergence-control parameters) are optimally determined and the first-order approximate solutions are known for different values of the coefficients which appear in the Eq. (14).
If we take into consideration, the solution given by Eq. (44), for α1 = 0.5, α2 = 0.5, m = 1 and v′(1) = 0 (Newtonian case), it is interesting to remark that from Eq. (40) results that K1 = K2 = 0 and therefore v1(x) = 0. It follows that,
In this example we consider α1 = 0.5, α2 = 0.75, m = 1 and v′(1) = 0 (non-Newtonian case). The values of the convergence-control parameters for the case 4.1.a are
The first-order approximate solution (44) can be written in the form
In Table 1 we present a comparison between the first-order approximate solution given by Eq. (75) and numerical results for some values of variable x and the corresponding relative errors.
For α1 = 0.5, α2 = 1, m = 1 and v′(1) = 0 we obtain for Eq. (44)
and the first-order approximate solution in the form:
In Table 2 a comparison between the first-order approximate solution given by Eq. (76) and numerical results, and the corresponding relative errors, is presented.
If we consider α1 = 0.5, α2 = 0.75, m = 1.5 and v′(1) = 0, then the convergence-control parameters for the case 4.1.a are
and therefore the first-order approximate solution has the form:
This approximate solution is compared in Table 3 with the numerical results
In the case α1 = 0.5, α2 = 1, m = 1.5 and v′(1) = 0, the convergence-control parameters for the case 4.1.a are
The first-order approximate solution becomes:
The OHAM solution (78) is compared in Table 4 with the numerical results.
For α1 = −0.5, α2 = 0.5, m = 1 and
The first-order approximate solution given by Eq. (54) can be written as
In Table 5 we present a comparison between the approximate solution (79) obtained by means of OHAM, and numerical results, and the corresponding relative errors.
x | vnumeric | v̄OHAM, Eq. (75) | relative error = |vnumeric − v̅OHAM| |
---|---|---|---|
0 | 1. | 1. | 0. |
0.1 | 0.8866085725 | 0.8866821774 | 7.3 · 10−5 |
0.2 | 0.7875621671 | 0.7875887398 | 2.6 · 10−5 |
0.3 | 0.7025338925 | 0.7025310249 | 2.8 · 10−6 |
0.4 | 0.6308917264 | 0.6308667405 | 2.4 · 10−5 |
0.5 | 0.5717826705 | 0.5717915984 | 8.9 · 10−6 |
0.6 | 0.5245086113 | 0.5245288249 | 2.0 · 10−5 |
0.7 | 0.4884393924 | 0.4884348043 | 4.5 · 10−6 |
0.8 | 0.4630600752 | 0.4630345952 | 2.5 · 10−5 |
0.9 | 0.4480006226 | 0.4479965346 | 4.0 · 10−6 |
1 | 0.4430235508 | 0.4430506264 | 2.7 · 10−5 |
x | vnumeric | v̅OHAM, Eq. (76) | relative error = |vnumeric − v̅OHAM| |
---|---|---|---|
0 | 1. | 1. | 0. |
0.1 | 0.8571788783 | 0.8572229943 | 4.4 · 10−5 |
0.2 | 0.7361604138 | 0.7361798440 | 1.9 · 10−5 |
0.3 | 0.6360415408 | 0.6360346979 | 6.8 · 10−6 |
0.4 | 0.5548519523 | 0.5548435245 | 8.4 · 10−6 |
0.5 | 0.4903083144 | 0.4903091700 | 8.5 · 10−7 |
0.6 | 0.4402488104 | 0.4402657870 | 1.6 · 10−5 |
0.7 | 0.4029480153 | 0.4029501998 | 2.1 · 10−6 |
0.8 | 0.3771255723 | 0.3771041676 | 2.1 · 10−5 |
0.9 | 0.3619470473 | 0.3619379098 | 9.1 · 10−6 |
1 | 0.3569430437 | 0.3569716548 | 2.8 · 10−5 |
In this case we consider α1 = −0.5, α2 = 0.75, m = 1 and
The first-order approximate solution given by Eq. (54) is
In Table 6 we present a comparison between the present solution (80) and numerical results and corresponding relative errors.
x | vnumeric | v̄OHAM, Eq. (77) | relative error = |vnumeric − v̄OHAM| |
---|---|---|---|
0 | 1. | 1. | 0. |
0.1 | 0.8122030427 | 0.8122723528 | 6.9 · 10−5 |
0.2 | 0.6480524030 | 0.6480757122 | 2.3 · 10−5 |
0.3 | 0.5075106818 | 0.5074963533 | 1.4 · 10−5 |
0.4 | 0.3899587013 | 0.3899662747 | 7.5 · 10−6 |
0.5 | 0.2943579028 | 0.2943935300 | 3.5 · 10−5 |
0.6 | 0.2193544932 | 0.2193300304 | 2.4 · 10−5 |
0.7 | 0.1631889587 | 0.1631613755 | 2.7 · 10−5 |
0.8 | 0.1243031560 | 0.1243139043 | 1.0 · 10−5 |
0.9 | 0.1014691098 | 0.1014847907 | 1.5 · 10−5 |
1. | 0.0939598983 | 0.0939116440 | 4.8 · 10−5 |
x | vnumeric | v̅OHAM, Eq. (78) | relative error = |vnumeric − v̄OHAM| |
---|---|---|---|
0 | 1. | 0.9999999999 | 2.2 · 10−16 |
0.1 | 0.7497063820 | 0.7498345352 | 1.2 · 10−4 |
0.2 | 0.5327963111 | 0.5327805533 | 1.5 · 10−5 |
0.3 | 0.3498945869 | 0.3498605615 | 3.4 · 10−5 |
0.4 | 0.2010814564 | 0.2012282330 | 1.4 · 10−4 |
0.5 | 0.0852885962 | 0.0853862885 | 9.7 · 10−5 |
0.6 | −0.0006296037 | −0.0008156336 | 1.8 · 10−4 |
0.7 | −0.0615295648 | −0.0616797274 | 1.5 · 10−4 |
0.8 | −0.1020045446 | −0.1018253040 | 1.7 · 10−4 |
0.9 | −0.1251782938 | −0.1250838453 | 9.4 · 10−5 |
1. | −0.1327530483 | −0.133080361 | 3.2 · 10−4 |
For α1 = −0.5, α2 = 1, m = 1 and
and the first-order approximate solution given by Eq. (54) is
In Table 7 we present a comparison between the solution (81) and numerical results, and corresponding relative errors.
In the following Tables, we present only comparisons between the present solutions obtained by means of OHAM and numerical results and corresponding relative errors.
In Figs. 1-23 the velocity v̅ is plotted against the horizontal distance x. Figs 1 and 2 show the variation of velocity field for Newtonian fluid, and Figs 3-6 show the variation of velocity field for non-Newtonian case, for different values of the parameter α2 and α1 = 0.5. It is observed that the speed of belt decreases with the increase of m. Figs 7-12 show the variation of the velocity field for the value α1 = −0.5 and different values of the parameter α2. In these cases the speed of belt increases with the increase of m for different values of the auxiliary function. In Figs 13-16 we present the variation of velocity field for α1 = 0.5 and different values of the parameters m and α2.
Here we observe that for any values of the parameter m, the speed of belt increases with the increase of α2. On the other hand, we remark from Figs 13-16 that the speed of belt decreases for any values of α1 > 0. But, from Figs 17-20 we deduce that the speed of belt increases for any values of α1 > 0. We observe in these last cases like Figs 13-16 that the speed of belt increases with the increase of α2. In Figs 21-23 present graphically the variation of velocity profile with increasing parameter α1. It is observed that the speed of belt decreases with the increase of α1. In all cases, we deduce that the speed of belt decreases for α1 = 0.5 and speed of belt increases for α1 = −0.5 and for any values of m.
x | vnumeric | v̄OHAM, Eq. (79) | relative error = |vnumeric − v̅OHAM| |
---|---|---|---|
0 | 1. | 1. | 0. |
0.1 | 1.2229062807 | 1.2229013971 | 4.8 · 10−6 |
0.2 | 1.4377929759 | 1.4377807450 | 1.2 · 10−5 |
0.3 | 1.6445808407 | 1.6445867997 | 5.9 · 10−6 |
0.4 | 1.8432159625 | 1.8432255865 | 9.6 · 10−6 |
0.5 | 2.0336046251 | 2.0336047334 | 1.0 · 10−7 |
0.6 | 2.2156306053 | 2.2156224054 | 8.1 · 10−6 |
0.7 | 2.3891355914 | 2.3891330225 | 2.5 · 10−6 |
0.8 | 2.5539053793 | 2.5539142068 | 8.8 · 10−6 |
0.9 | 2.7096483135 | 2.7096510721 | 2.7 · 10−6 |
1. | 2.8559539636 | 2.8559439901 | 9.9 · 10−6 |
x | vnumeric | v̄OHAM, Eq. (80) | relative error = |vnumeric − v̅OHAM| |
---|---|---|---|
0 | 1. | 1. | 0. |
0.1 | 1.2518351629 | 1.2518238888 | 1.1 · 10−5 |
0.2 | 1.4920691495 | 1.4920531530 | 1.5 · 10−5 |
0.3 | 1.7207387643 | 1.72074398947 | 5.2 · 10−6 |
0.4 | 1.9379031842 | 1.93792002203 | 1.6 · 10−5 |
0.5 | 2.1436214696 | 2.14362291135 | 1.4 · 10−6 |
0.6 | 2.3379068070 | 2.33789447422 | 1.2 · 10−5 |
0.7 | 2.5207347633 | 2.52072993082 | 4.8 · 10−6 |
0.8 | 2.6920220811 | 2.69203422220 | 1.2 · 10−5 |
0.9 | 2.8516000667 | 2.85160438618 | 4.3 · 10−6 |
1. | 2.9991636824 | 2.99914870737 | 1.4 · 10−5 |
x | vnumeric | v̄OHAM, Eq. (81) | relative error = |vnumeric − v̅OHAM| |
---|---|---|---|
0 | 1. | 1. | 0. |
0.1 | 1.2839691268 | 1.2839541027 | 1.5 · 10−5 |
0.2 | 1.5522080417 | 1.5521893463 | 1.8 · 10−5 |
0.3 | 1.8048997685 | 1.8049032366 | 3.4 · 10−6 |
0.4 | 2.0422481654 | 2.0422731488 | 2.4 · 10−5 |
0.5 | 2.2645158331 | 2.2645193174 | 3.4 · 10−6 |
0.6 | 2.4718989261 | 2.4718821917 | 1.6 · 10−5 |
0.7 | 2.6645698769 | 2.6645620846 | 7.7 · 10−6 |
0.8 | 2.8426469875 | 2.8426623139 | 1.5 · 10−5 |
0.9 | 3.0061618853 | 3.0061679939 | 6.1 · 10−6 |
1 | 3.1549980906 | 3.1549778194 | 2.0 · 10−5 |
x | vnumeric | v̄OHAM | relative error = |vnumeric − v̅OHAM| |
---|---|---|---|
0. | 1. | 1. | 0. |
0.1 | 1.2609509091 | 1.2607564771 | 1.9 · 10−4 |
0.2 | 1.5099315977 | 1.5092509182 | 6.8 · 10−4 |
0.3 | 1.7469050126 | 1.7471433465 | 2.3 · 10−4 |
0.4 | 1.9718544689 | 1.9726191668 | 7.6 · 10−4 |
0.5 | 2.1847229322 | 2.1849754601 | 2.5 · 10−4 |
0.6 | 2.3853998014 | 2.3850003567 | 3.9 · 10−4 |
0.7 | 2.5737023327 | 2.5735062125 | 1.9 · 10−4 |
0.8 | 2.7493402910 | 2.7498662597 | 5.2 · 10−4 |
0.9 | 2.9118599070 | 2.9120716519 | 2.1 · 10−4 |
1 | 3.0605243186 | 3.0595422315 | 9.8 · 10−4 |
x | vnumeric | v̄OHAM | relative error = |vnumeric − v̄OHAM| |
---|---|---|---|
0 | 1. | 1. | 0. |
0.1 | 1.3085076286 | 1.3086293331 | 1.2 · 10−4 |
0.2 | 1.5989200672 | 1.5982295694 | 6.9 · 10−4 |
0.3 | 1.8714226404 | 1.8715831764 | 1.6 · 10−4 |
0.4 | 2.1262167929 | 2.1270250021 | 8.0 · 10−4 |
0.5 | 2.3635521486 | 2.3638995934 | 3.4 · 10−4 |
0.6 | 2.5835932960 | 2.5832535309 | 3.3 · 10−4 |
0.7 | 2.7864530758 | 2.7863211230 | 1.3 · 10−4 |
0.8 | 2.9721424944 | 2.9728245480 | 6.8 · 10−4 |
0.9 | 3.1405027991 | 3.1408732113 | 3.7 · 10−4 |
1 | 3.2910606175 | 3.2900917909 | 9.6 · 10−4 |
x | vnumeric | v̄OHAM | relative error = |vnumeric − v̄OHAM| |
---|---|---|---|
0 | 1. | 1. | 0. |
0.1 | 1.3620021521 | 1.3619942603 | 7.8 · 10−6 |
0.2 | 1.6989187493 | 1.6989013673 | 1.7 · 10−5 |
0.3 | 2.0111542818 | 2.0111557494 | 1.4 · 10−6 |
0.4 | 2.2991388016 | 2.2991755901 | 3.6 · 10−5 |
0.5 | 2.5634672715 | 2.5634733579 | 6.0 · 10−6 |
0.6 | 2.8046481809 | 2.8046272723 | 2.0 · 10−5 |
0.7 | 3.0231881146 | 3.0231849910 | 3.1 · 10−6 |
0.8 | 3.2195266334 | 3.2195597898 | 3.3 · 10−5 |
0.9 | 3.3939541762 | 3.3939648992 | 1.0 · 10−5 |
1 | 3.5464381070 | 3.5464081807 | 2.9 · 10−5 |
x | vnumeric | v̄OHAM | relative error = |vnumeric − v̄OHAM| |
---|---|---|---|
0 | 1. | 1. | 0. |
0.1 | 0.9053448275 | 0.9052579961 | 8.6 · 10−5 |
0.2 | 0.8203448275 | 0.8203616740 | 1.6 · 10−5 |
0.3 | 0.7453448275 | 0.7453758909 | 3.1 · 10−5 |
0.4 | 0.6803448275 | 0.6803496596 | 4.8 · 10−6 |
0.5 | 0.6253448275 | 0.6253174083 | 2.7 · 10−5 |
0.6 | 0.5803448275 | 0.5803001231 | 4.4 · 10−5 |
0.7 | 0.5453448275 | 0.5453063866 | 3.8 · 10−5 |
0.8 | 0.5203448275 | 0.5203333245 | 1.1 · 10−5 |
0.9 | 0.5053448275 | 0.5053674704 | 2.2 · 10−5 |
1 | 0.5003448275 | 0.5003855577 | 4.0 · 10−5 |
x | vnumeric | v̄OHAM | relative error = |vnumeric − v̄OHAM| |
---|---|---|---|
0 | 1. | 1. | 0. |
0.1 | 0.8866085725 | 0.8863361236 | 2.7 · 10−4 |
0.2 | 0.7875621671 | 0.7875077487 | 5.4 · 10−5 |
0.3 | 0.7025338925 | 0.7026438847 | 1.0 · 10−4 |
0.4 | 0.6308917264 | 0.6309536994 | 6.1 · 10−5 |
0.5 | 0.5717826705 | 0.5717301604 | 5.2 · 10−5 |
0.6 | 0.5245086113 | 0.5243482669 | 1.6 · 10−4 |
0.7 | 0.4884393924 | 0.4882588929 | 1.8 · 10−4 |
0.8 | 0.4630600752 | 0.4629799268 | 8.0 · 10−5 |
0.9 | 0.4480006226 | 0.4480861921 | 8.5 · 10−5 |
1 | 0.4430235508 | 0.4431992337 | 1.7 · 10−4 |
x | vnumeric | v̄OHAM | relative error = |vnumeric − v̄OHAM| |
---|---|---|---|
0 | 1. | 1.0000000000 | 6.6 · 10−16 |
0.1 | 0.8571788783 | 0.8574415758 | 2.6 · 10−4 |
0.2 | 0.7361604138 | 0.7371670017 | 1.0 · 10−3 |
0.3 | 0.6360415408 | 0.6363653381 | 3.2 · 10−4 |
0.4 | 0.5548519523 | 0.5540906717 | 7.6 · 10−4 |
0.5 | 0.4903083144 | 0.4886618248 | 1.6 · 10−3 |
0.6 | 0.4402488104 | 0.4383723210 | 1.8 · 10−3 |
0.7 | 0.4029480153 | 0.4016994772 | 1.2 · 10−3 |
0.8 | 0.3771255723 | 0.3772291464 | 1.0 · 10−4 |
0.9 | 0.3619470473 | 0.3635699289 | 1.6 · 10−3 |
1 | 0.3569430437 | 0.3593259550 | 2.3 · 10−3 |
x | vnumeric | v̄OHAM | relative error = |vnumeric − v̄OHAM| |
---|---|---|---|
0 | 1. | 1.0000000018 | 1.8 · 10−9 |
0.1 | 0.8577489626 | 0.8577105235 | 3.8 · 10−5 |
0.2 | 0.7302489626 | 0.7302905954 | 4.1 · 10−5 |
0.3 | 0.6177489626 | 0.6177941197 | 4.5 · 10−5 |
0.4 | 0.5202489626 | 0.5202636998 | 1.4 · 10−5 |
0.5 | 0.4377489626 | 0.4377305945 | 1.8 · 10−5 |
0.6 | 0.3702489626 | 0.3702146522 | 3.4 · 10−5 |
0.7 | 0.3177489626 | 0.3177242586 | 2.4 · 10−5 |
0.8 | 0.2802489626 | 0.2802562722 | 7.3 · 10−6 |
0.9 | 0.2577489626 | 0.2577959802 | 4.7 · 10−5 |
1 | 0.2502489626 | 0.2503170184 | 6.8 · 10−5 |
x | vnumeric | v̄OHAM | relative error = |vnumeric − v̄OHAM| |
---|---|---|---|
0 | 1 | 1. | 0. |
0.1 | 0.8122030427 | 0.8121518633 | 5.1 · 10−5 |
0.2 | 0.6480524030 | 0.6479781411 | 7.4 · 10−5 |
0.3 | 0.5075106818 | 0.5075931624 | 8.2 · 10−5 |
0.4 | 0.3899587013 | 0.3902155536 | 2.5 · 10−4 |
0.5 | 0.2943579028 | 0.2945933839 | 2.3 · 10−4 |
0.6 | 0.2193544932 | 0.2193140489 | 4.0 · 10−5 |
0.7 | 0.1631889587 | 0.1629677570 | 2.2 · 10−4 |
0.8 | 0.1243031560 | 0.1242157441 | 8.7 · 10−5 |
0.9 | 0.1014691098 | 0.1018089431 | 3.3 · 10−4 |
1 | 0.0939598983 | 0.0945851040 | 6.2 · 10−4 |
x | vnumeric | v̄OHAM | relative error = |vnumeric − v̄OHAM| |
---|---|---|---|
0 | 1. | 0.9999999999 | 3.5 · 10−15 |
0.1 | 0.7497063820 | 0.7498561224 | 1.4 · 10−4 |
0.2 | 0.5327963111 | 0.5319377563 | 8.5 · 10−4 |
0.3 | 0.3498945869 | 0.3497527318 | 1.4 · 10−4 |
0.4 | 0.2010814564 | 0.2022779050 | 1.1 · 10−3 |
0.5 | 0.0852885962 | 0.0867149009 | 1.4 · 10−3 |
0.6 | −0.0006296037 | −0.0003046742 | 3.2 · 10−4 |
0.7 | −0.0615295648 | −0.0622411985 | 7.1 · 10−4 |
0.8 | −0.1020045446 | −0.1024854604 | 4.8 · 10−4 |
0.9 | −0.1251782938 | −0.1243191394 | 8.5 · 10−4 |
1 | −0.1327530483 | −0.1309114642 | 1.8 · 10−3 |
x | vnumeric | v̄OHAM | relative error = |vnumeric – v̄OHAM| |
---|---|---|---|
0 | 1. | 1.0000000000 | 3.5 · 10−15 |
0.1 | 1.2229384382 | 1.2230063532 | 6.7 · 10−5 |
0.2 | 1.4378247530 | 1.4378525565 | 2.7 · 10−5 |
0.3 | 1.6446154175 | 1.6445994145 | 1.6 · 10−5 |
0.4 | 1.8432480259 | 1.8432271985 | 2.0 · 10−5 |
0.5 | 2.0336362803 | 2.0336453830 | 9.1 · 10−6 |
0.6 | 2.2156629647 | 2.2157024590 | 3.9 · 10−5 |
0.7 | 2.3891700476 | 2.3891957459 | 2.5 · 10−5 |
0.8 | 2.5539441869 | 2.5538810605 | 6.3 · 10−5 |
0.9 | 2.7096939733 | 2.7094820881 | 2.1 · 10−4 |
1 | 2.8560115627 | 2.8556993037 | 3.1 · 10−4 |
x | vnumeric | v̄OHAM | relative error = |vnumeric – v̄OHAM| |
---|---|---|---|
0 | 1. | 1. | 0. |
0.1 | 1.2518656079 | 1.2518515888 | 1.4 · 10−5 |
0.2 | 1.4920877143 | 1.4920910613 | 3.3 · 10−6 |
0.3 | 1.7207480054 | 1.7207531923 | 5.1 · 10−6 |
0.4 | 1.9379133831 | 1.9379068594 | 6.5 · 10−6 |
0.5 | 2.1436299399 | 2.1436163389 | 1.3 · 10−5 |
0.6 | 2.3379146593 | 2.3379104050 | 4.2 · 10−6 |
0.7 | 2.5207436438 | 2.5207591446 | 1.5 · 10−5 |
0.8 | 2.6920346069 | 2.6920579596 | 2.3 · 10−5 |
0.9 | 2.8516191935 | 2.8516179705 | 1.2 · 10−6 |
1 | 2.9991957021 | 2.9991619152 | 3.3 · 10−5 |
x | vnumeric | v̄OHAM | relative error = |vnumeric – v̄OHAM| |
---|---|---|---|
0 | 1. | 1.0000000000 | 7.1 · 10−15 |
0.1 | 1.2839954362 | 1.2839898927 | 5.5 · 10−6 |
0.2 | 1.5522090538 | 1.5522083786 | 6.7 · 10−7 |
0.3 | 1.8048758830 | 1.8048779484 | 2.0 · 10−6 |
0.4 | 2.0422291469 | 2.0422291573 | 1.0 · 10−8 |
0.5 | 2.2644931580 | 2.2644898117 | 3.3 · 10−6 |
0.6 | 2.4718738134 | 2.4718705633 | 3.2 · 10−6 |
0.7 | 2.6645442766 | 2.6645460387 | 1.7 · 10−6 |
0.8 | 2.8426240012 | 2.8426304402 | 6.4 · 10−6 |
0.9 | 3.0061452466 | 3.0061463212 | 1.0 · 10−6 |
1 | 3.1549957250 | 3.1549849565 | 1.0 · 10−5 |
x | vnumeric | v̄OHAM | relative error = |vnumeric – v̄OHAM| |
---|---|---|---|
0 | 1. | 1. | 0. |
0.1 | 1.2609808120 | 1.2603664398 | 6.1 · 10−4 |
0.2 | 1.5099516324 | 1.5093655735 | 5.8 · 10−4 |
0.3 | 1.7469198626 | 1.7467254400 | 1.9 · 10−4 |
0.4 | 1.9718652899 | 1.9721741750 | 3.0 · 10−4 |
0.5 | 2.1847302867 | 2.1854400932 | 7.0 · 10−4 |
0.6 | 2.3854057847 | 2.3862516582 | 8.4 · 10−4 |
0.7 | 2.5737099736 | 2.5743374079 | 6.2 · 10−4 |
0.8 | 2.7493549252 | 2.7494260892 | 7.1 · 10−5 |
0.9 | 2.9118887501 | 2.9112465530 | 6.4 · 10−4 |
1 | 3.0605838497 | 3.0595277920 | 1.0 · 10−3 |
x | vnumeric | v̄OHAM | relative error = |vnumeric – v̄OHAM| |
---|---|---|---|
0 | 1. | 0.9999999999 | 2.8 · 10−14 |
0.1 | 1.3085289178 | 1.3084934573 | 3.5 · 10−5 |
0.2 | 1.5989107009 | 1.5989174796 | 6.7 · 10−6 |
0.3 | 1.8713842220 | 1.8713975406 | 1.3 · 10−5 |
0.4 | 2.1261798486 | 2.1261654650 | 1.4 · 10−5 |
0.5 | 2.3635084974 | 2.3634750640 | 3.3 · 10−5 |
0.6 | 2.5835451658 | 2.5835323438 | 1.2 · 10−5 |
0.7 | 2.7864041509 | 2.7864403839 | 3.6 · 10−5 |
0.8 | 2.9720990828 | 2.9721581992 | 5.9 · 10−5 |
0.9 | 3.1404738992 | 3.1404724101 | 1.4 · 10−6 |
1 | 3.2910682398 | 3.2909802656 | 8.7 · 10−5 |
x | vnumeric | v̄OHAM | relative error = |vnumeric – v̄OHAM| |
---|---|---|---|
0 | 1. | 0.9999999999 | 2.2 · 10−13 |
0.1 | 1.3620158757 | 1.3619756005 | 4.0 · 10−5 |
0.2 | 1.6988832705 | 1.6988902831 | 7.0 · 10−6 |
0.3 | 2.0110642417 | 2.0110793966 | 1.5 · 10−5 |
0.4 | 2.2990536573 | 2.2990383227 | 1.5 · 10−5 |
0.5 | 2.5633704925 | 2.5633333487 | 3.7 · 10−5 |
0.6 | 2.8045422031 | 2.8045269744 | 1.5 · 10−5 |
0.7 | 3.0230772590 | 3.0231164379 | 3.9 · 10−5 |
0.8 | 3.2194185417 | 3.2194842324 | 6.5 · 10−5 |
0.9 | 3.3938604328 | 3.3938593947 | 1.0 · 10−6 |
1 | 3.5463864691 | 3.5462883840 | 9.8 · 10−5 |
6 Conclusions
In the present work we propose Optimal Homotopy Asymptotic Method to obtain approximate analytical solutions for nonlinear differential equation of thin film flow of an Oldroyd 6-constant fluid. The validity of our procedure was demonstrated choosing the linear operator in two cases and different expressions of the convergence-control functions. Some representative examples are given and very good agreement was found between the approximate analytical results and numerical simulation results. The proposed procedure is valid even if the nonlinear differential equation does not contain any small or large parameter. We examine quantitative effect of the parameters α1, α2 and m and the relative errors of approximate solutions in comparison with numerical results. OHAM is an approach proposed for the first time by Marinca and Herisanu in the study of the thin film flow [29], then is used in many other studies as: the steady flow of a fourth-grade fluid past a porous plate [30], the MHD Jeffery-Hamel flow [31] and so on [32], such that for achieving a very accurate solution, our procedure ensured a very rapid convergence after only one iteration. Instead of infinite series, the OHAM searches for only a few terms and does not need recurrence formula. The parameters which appear in the composition of the auxiliary functions and in the linear operator are optimally identified via various methods by mathematically rigorous point of view. A large number of parameters in the auxiliary functions lead to more accurate results.
In all cases presented in this paper, for different values of the parameters α1, α2 and m, we obtain an excellent agreement of the first-order approximate solutions. Also, the exact solutions of Newtonian fluid is obtained. It is worth mentioning that the proposed method is straightforward, concise and can be applied to other nonlinear problems.
References
1 J.G. Oldroyd, Proc. Roy. Soc. Lond. A 200, 523 (1950).10.1098/rspa.1950.0035Search in Google Scholar
2 K. Strauss, Rheol. Acta 14, 1058 (1975).10.1007/BF01515901Search in Google Scholar
3 C. Hancock, E. Lewis, N.K. Moffat, J. Fluid Mech. 112, 315 (1981).Search in Google Scholar
4 K. R. Rajagopal, R.K. Bhatnagar, Acta Mechanica 113, 233 (1995).10.1007/BF01212645Search in Google Scholar
5 S. Baris, Acta Mechanica 147, 125 (2001).10.1007/BF01182357Search in Google Scholar
6 T. Hayat, M. Khan, S. Asghar, Acta Mechanica 168, 213 (2004).10.1007/s00707-004-0085-2Search in Google Scholar
7 T. Hayat, M. Khan, M. Ayub, J. Math. Anal. Appl. 298, 225 (2004).Search in Google Scholar
8 Y. Wang, T. Hayat, K. Hutter, Int. J. Non-Linear Mech. 40, 49 (2005).Search in Google Scholar
9 M. Sajid, A.M. Siddiqui, T. Hayat, Int. J. Eng. Sci. 45, 381 (2007).Search in Google Scholar
10 M. Khan, T. Hayat, M. Ayub, Computers and Mathematics with Applications 53, 1088 (2007).10.1016/j.camwa.2006.05.029Search in Google Scholar
11 T. Hayat, M. Khan, M. Ayub, J. Comput. Appl. Math. 202, 402 (2007).Search in Google Scholar
12 R. Ellahi, T. Hayat, T. Javed, S. Asghar, Mathematical and Computer Modelling 48, 1191 (2008).10.1016/j.mcm.2007.12.017Search in Google Scholar
13 M. Sajid, T. Hayat, Physics Letters A 372, 1827 (2008).10.1016/j.physleta.2007.10.069Search in Google Scholar
14 D. Vieru, C. Fetecău, C. Fetecău, Appl. Math. Comput. 201, 834 (2008).Search in Google Scholar
15 T. Hayat, R. Ellahi, F.M. Mahomed, Commun. Nonl. Sci. Numer. Simulat. 14, 133 (2009).Search in Google Scholar
16 R. Ellahi, T. Hayat, F.M. Mahomed, A. Zeeshan, Commun. Nonl. Sci. Numer. Simulat. 15, 322 (2010).Search in Google Scholar
17 C. Fetecău, T. Hayat, J. Zierep, M. Sajid, Nonlinear Analysis: Real World Applications 12, 1 (2011).10.1016/j.nonrwa.2009.12.009Search in Google Scholar
18 L. Zheng, Z. Guo, X. Zhang, Nonlinear Analysis: Real World Applications 12, 3499 (2011).10.1016/j.nonrwa.2011.06.010Search in Google Scholar
19 R.A. Shah, S. Islam, A.M. Siddiqui, T. Haroon, Computers and Mathematics with Applications 63, 695 (2012).10.1016/j.camwa.2011.11.033Search in Google Scholar
20 A.M. Siddiqui, A.A. Farooq, T. Haroon, M.A. Rana, Can. J. Phys. 92, 1196 (2014).Search in Google Scholar
21 T. Hayat, S. Zaib, S. Asghar, A.A. Hendi, Appl. Math. Mech. - Engl. Ed. 34(4), 411 (2012).10.1007/s10483-012-1560-7Search in Google Scholar
22 A. Nayfeh, Problems in Perturbation (Wiley New York, 1985).Search in Google Scholar
23 G. Adomian, J. Math. Anal. Appl. 135, 501 (1988).Search in Google Scholar
24 A. Beléndez, C. Pascual, C. Neipp, T. Beléndez, A. Hernández, Int. J. Non-Linear Numer. Simulat. 9, 9 (2001).Search in Google Scholar
25 J.I. Ramos, Appl. Math. Comput. 196, 483 (2008).Search in Google Scholar
26 J.N. He, Int. J. Non-Linear Mechanics 37, 309 (2002).10.1016/S0020-7462(00)00116-5Search in Google Scholar
27 N. Herişanu, V. Marinca, Z. Naturforsch 670, 509 (2012).10.1155/2012/465023Search in Google Scholar
28 V. Marinca, N. Herişanu, In: T.E. Simos et. al. (Ed.), Numerical Analysis and Appl. Math. ICNAAM, AIP Conf. Proc., Sep. 19-25, 2012, Kos, Greece, (Amer. Inst. Physics, New York 2012).Search in Google Scholar
29 V. Marinca, N. Herişanu, I. Nemeş, Cent. Eur. J. Phys. 6, 648 (2008).Search in Google Scholar
30 V. Marinca, N. Herişanu, C. Bota, B. Marinca, Applied Mathematics Letters 22, 245 (2009).10.1016/j.aml.2008.03.019Search in Google Scholar
31 V. Marinca, N. Herişanu, Mathematical Problems in Engineering, Article ID 169056 (2011).10.1155/2011/169056Search in Google Scholar
32 V. Marinca, N. Herişanu, Nonlinear Dynamical Systems in Engineering - Some Approximate Approaches (Springer Verlag, Heidelberg, 2011).Search in Google Scholar
© 2016 Remus-Daniel Ene, published by De Gruyter Open
This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 3.0 License.