This article addresses the phenomenon of boundary layer flow and heat transfer in respect of the motion of second grade viscous fluid over an unsteady stretchable surface. Optimal variational iteration method (OVIM) is employed to solve the governing differential equation. OVIM is the modified version of variational iteration method (VIM). The convergence of the obtained solution as well as the effect of pertinent parameters on the velocity components are discussed, tabulated, and graphed. Behaviour of coefficient of skin-friction is also tabulated and explained for different parameters. Comparison of residual errors sought via use of VIM and OVIM is also presented.
I. INTRODUCTION
Study of boundary layer flow along with the phenomenon of heat transfer over stretching sheet fascinated many researchers for last two or three decades due to its manifold applications in engineering disciplines, like biomedical, metallurgy and chemical engineering. The accoutrements inured for surgical purposes also make use of this phenomenon. Specifically, we encounter with such type of flow while working in aerodynamics, inflation of balloons, drafting of plastic films etc. Concept of boundary layer flow over a stretching surface inaugurated in his innovative work by Sakiadis.28,29 After the introduction of such type of flows with incomparable utilities, many renowned mathematicians and engineers investigated different rheological based on boundary layer flow, see Refs. 1, 7–9, 18, and 32.
As all physical phenomenons are time varying, so its natural to discuss unsteady flows, in this respect the hydrodynamic viscous thin fluid film flow over a sheet stretching was studied by Wang,45 and by many others. Soon after,2,3,6,19,45 made extensions to Wang’s analysis from different perspectives. Wang43 and Wang et al.44 also scrutinised the time–dependent flows of viscous as well as power law fluids over stretchable surface, they adopted homotopy analysis method to find the solution.4 provided an analytically similar solution for flow over a flat plate convective surface. Moreover, Makinde and Aziz20 elongated the work contained in Ref. 4 by considering the flow over vertical stretchable sheet. Shahzad and Ali30,31 acquiesced an approximate series solution for flow of non-Newtonian fluid by constant vertical stretching of sheet. Bachok et al.17 inquired unsteady boundary layer flow and transfer of heat over a surface that was stretched in vertical. Sheikholeslami et al.35,37 made use of Adomian decomposition method. Further more, study of viscous fluid under the effects of magnetic field with the help of optimal homotopy asymptotic method has been carried out in Refs. 34 and 39. The homotopy perturbation and homotopy analysis methods have been used to investigate the flow of condensation film on tilted revolving disk and flow across parallel plates.33,38 Sheikholeslami and Ganji36 developed a systematic differential transform method for transfer of heat in conjunction with Brownian effects across parallel plates.
For second grade fluids, a number of attempts already have been made to determine analytic and closed form solutions of velocity field and stream functions against various geometrical and physical arrangements. Coscia and Galdi10 studied classical problems related to existence, stability, and uniqueness of solutions. Moreover, Dunn and Fosdick11 mainly put focus on stability pertaining to flow problems. Ting41 found solutions for unsteady second grade fluid flowing in bounded regions. Hayat et al.13 discussed some unsteady solutions while considering stream functions of particular forms. Further more, Siddiqui and Kaloni40 and Aziz et al.5 provided similar results by varying the form of stream function. Motivated by above mentioned contributions, in this paper we shall investigate the unsteady flow of a fluid film flowing in 2-dimensions over a surface with constant stretching. Precisely the current work has some significant differences with the prior work done so far. We deal with a no–linear problem that arose here analytically with the help of OVIM. The idea of OVIM for initial value problems was introduced by.12,15,21,42 The OVIM for boundary value problems is an efficient technique developed by Tauseef et al.22 OVIM for boundary value problems is a modified form of VIM. In this technique, all initial and boundary conditions are consumed before constructing an iterative algorithm. Thus we construct an iterative algorithm which is free from undetermined coefficients. In most of the previous iterative algorithms for boundary value problems, each iteration contains the undetermined coefficients. Consuming all the boundary Conditions in the developed algorithm reduces the computational volume and enhances convergence. In OVIM we insert an auxiliary parameter h in the algorithm that facilitates fast convergence of the approximate solution in the given domain of the problem. The minimum of the norm 2 of the residual error is used to find the optimal choice of the parameter which is not in the VIM, see Refs. 12, 22, 23, and 26.
II. MATHEMATICAL FORMULATION
Consider an unsteady incompressible second grade fluid flowing in xy-plane in the form of thin film having uniform width h(t) that is placed on a horizontal sheet. The fluid is set into motion by stretching the sheet constantly along x-axis, and the y-axis is perpendicular to the direction of flow. By keeping in mind the geometry of the problem, we shall consider velocity field of the form
The constitutive equation for a second grade incompressible homogenous fluid27 is
where T is the Cauchy stress tensor, p is the pressure, I is the identity tensor, μ the dynamic viscosity, α1 and α2 are the material constants. The first two Rivlin-Ericksen tensors are
where is time derivative and superscript † denotes usual matrix transpose. Note that the Clausius-Duhemin equality (i.e. μ ≥ 0, α1 ≥ 0, α1 + α2 = 0.) is also satisfied. For α1 = α2 = 0, the equation(3) gives the constitutive equation of viscous fluid. Using Equations (1), (2) and (3) lead to governing equations for the unsteady boundary layer flow that are given by
where ν and ρ denote the kinematic viscosity and density of the fluid respectively. We assume the corresponding boundary conditions are
where c1 and c2 are positive real numbers and 1 − c2t > 0. Consider the dimensionless transformations19
where η is the dimensionless local similarity variable and ψ(x, y) is the stream function defined by and ν is the kinematic viscosity of the second grade fluid. The continuity equation (4) is automatically satisfied and from (2)–(6), one can write
Here denotes the unsteadiness parameter, is the parameter of second grade with no dimension, the prime indicates the derivative with reference to η. Moreover, χ is the film thickness with no dimension and used to denote the value of η at the free surface so that (5) gives, see Ref. 19
Observe that h(t) get decreased monotonically with the increase in time and χ remains constant which depends only upon S, see Ref. 19. The detailed discussion on such cases can be found in the Refs. 2, 3, 6, 19, 43, and 44.
The shear stress τ on the surface of the thin liquid film sheet is
and the local skin-friction coefficient or frictional drag coefficient is
In dimensionless form, we have
III. VARIATIONAL ITERATION METHOD (VIM)
The basic notion of the technique is illustrated by considering the problem of finding z(η) such that
where L and N are linear and nonlinear operators respectively and g(η) is the source term. For a given z0, approximate solution zp+1 of equation (16) can be found as follows:
where λ is named as the Lagrange multiplier. This Lagrange multiplier is obtained by taking variation δ on both sides of equation (17) with respect to the variable zp.
where is a term being restricted which in turn gives . An unknown parameter λ(s, η) is found by making use of the optimal conditions, see Refs. 24 and 25. This yield an exact solution z(η), where
This method of getting the approximate solution is known as the method of variational iteration abbreviated as (VIM). The VIM was introduced by Inokuti et al.16 However, He14 has established this method for giving solution to wider range of problems, emerging in several fields of pure and applied sciences. In VIM, proper selection initial approximation leads to the fast converging solution, see Refs. 24 and 25 and the references therein.
A. Optimal variational iteration method (OVIM)
Consider a second order nonlinear ordinary differential equation,
subject to the boundary conditions
where is a linear differential operator, N(η) represents the nonlinear operator and g(η) is an inhomogeneous term and is set of real numbers. According to standard VIM, the correction functional used is given as
To make stable correction functional, the lagrange multiplier is defined as λ = s − η, see Refs. 24 and 25, we get the following iterative formula
An unknown auxiliary parameter h is inserted into the iterative formula (20), for n = 0, (20) becomes
For optimal variational iteration method, we will proceed as follows
where a0 and a1 are unknown parameters which can be determined by using the given boundary conditions. Substituting the values of a0 and a1 in (24), we obtain
From the above equation, we develop the following sequence
The iterative algorithm (23)–(25) does not contain undetermined constants except an auxiliary parameter h which is used to control the convergence of approximate solution optimally by minimising the norm 2 of residual function over the domain of the given problem. For further details, see Ref. 22. This algorithm is called optimal variational iteration formula for solving boundary value problem (17).
IV. NUMERICAL EXAMPLE
A. Example 4.1
Using VIM, we construct the correction functional for (26) as follows
Using the optimality conditions,25 one has
We assume
as an initial approximation, where ai, i = 0, 1, 2, 3, are unknowns and will be determined by boundary conditions (27). Using (31) in (30), we obtain the following iterations
Applying the boundary conditions given by Eq. 27 on approximation f3(η), we obtain
We define the following residual error function
Now using OVIM, we proceed as follows
Using the boundary conditions, we have
The OVIM algorithm is given as follows
We obtain the following approximations
We define the following residual function
and the norm 2 of above residual error function is
The above residual function can be used to approximate e3(h), and the optimal value of h can be determined by minimizing the e3(h). The value of h is found to be 0.0027467499661031 when minimum value of e3(h) is 0.0000609419536125150826. We may obtain more accuracy by increasing the number of iterations.
The graphs of residual errors r3(η) and r3(η, h) given by (32) and (33) are plotted in Figure 1. The graph of r3(η) obtained by VIM diverges but the graph of r3(η, h) obtained by OVIM converges to zero. It is concluded that OVIM is more reliable and accurate as compared to VIM.
Residual Error comparison of VIM and OVIM Solutions when α = 0 and S = 0.2.
B. Example 4.2
We construct the correction functional for (34) as follows
Using the optimality conditions,25 one has
Applying the boundary conditions (35) on f2(η), the unknowns are
We define the following residual error function
Now using OVIM, we proceed as follows
We obtain the following approximations
We define the following residual error function
and the norm 2 of residual function (40) is
The optimal value of h can be determined by minimizing e2(h). At h = 0.810519126906262, the above norm 2 of residual error function is minimum.
The graphs of residual errors r2(η) and r2(η, h) given by (39) and (40) are plotted in Figure 2. The graph of r2(η) obtained by VIM diverges but the graph of r2(η, h) obtained by OVIM converges to zero. It is concluded that OVIM is more reliable and accurate as compared to VIM.
Residual Error Comparison of VIM and OVIM Solutions when α = 0.9 and S = 0.2.
V. RESULTS AND DISCUSSIONS
Here we discuss the sway of different parameters like S and α on the components of velocity f, acceleration and skin-friction coefficients . The change in the coefficient of the skin-friction for α and S is shown in Table I.
Values of skin-friction coefficient for different parameters at η = 0.
α . | S . | . |
---|---|---|
0.0 | 0.5 | 2.24676880475314 |
0.2 | 0.5 | 4.09820409923631 |
0.5 | 0.5 | 6.93689203353080 |
0.7 | 0.5 | 9.40436723010559 |
1.0 | 0.5 | 13.8596936496045 |
2.0 | 0.5 | 59.2327101256866 |
0.5 | 0 | 6.15999848254210 |
0.5 | 0.2 | 4.81951263596841 |
0.5 | 0.5 | 7.11968336959140 |
0.5 | 1 | 4.98044449104253 |
0.5 | 1.5 | 2.83681213996240 |
0.5 | 2 | 0.742432181172167 |
α . | S . | . |
---|---|---|
0.0 | 0.5 | 2.24676880475314 |
0.2 | 0.5 | 4.09820409923631 |
0.5 | 0.5 | 6.93689203353080 |
0.7 | 0.5 | 9.40436723010559 |
1.0 | 0.5 | 13.8596936496045 |
2.0 | 0.5 | 59.2327101256866 |
0.5 | 0 | 6.15999848254210 |
0.5 | 0.2 | 4.81951263596841 |
0.5 | 0.5 | 7.11968336959140 |
0.5 | 1 | 4.98044449104253 |
0.5 | 1.5 | 2.83681213996240 |
0.5 | 2 | 0.742432181172167 |
Figure 3 and Figure 4 are displaying the changes made by α on the component of velocity f and . These figures depict that f and the function which are increasing depending upon α. The variation in is considerably small as compared to f. The boundary layer thickness gets increased with increase in α. Fig. 5 and Fig. 6 are depicting the changes made by S on f and . It can be observed from the given figures that these velocities f and are the functions which are increasing depending upon S. The boundary layer thickness gets decreased as we take larger values of S. In Table I, the change of the coefficient of skin-friction for several values of α and S has been shown. It has been noted that the coefficient of skin-friction get increased by increasing α. In this table, it is also indicated that the coefficient of skin-friction get decreased when S < 2, for S = 2 the coefficient of skin-friction is almost zero and for S > 2 it is increased.
VI. CONCLUSION
In this article we analysed the spectacle of boundary layer flow and heat transfer with respect to the motion of second grade fluid over a stretchable surface. Flow we considered was unsteady and fluid was viscous. Technique of OVIM is implemented to solve the governing differential equation. OVIM is the modified version of VIM. The convergence of solution we sought as well as the clout of pertinent parameters on the velocity profile are discussed, tabulated, and graphed. Behaviour of coefficient of skin-friction is also tabularized and elucidated for different parameters. Comparison of residual errors sought via use of VIM and OVIM is also presented. Amazingly OVIM gave much stable and reliable results as compared with VIM.