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.

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, 79, 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.

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

(1)

The constitutive equation for a second grade incompressible homogenous fluid27 is

(2)

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

(3)

where ddt 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

(4)
(5)

where ν and ρ denote the kinematic viscosity and density of the fluid respectively. We assume the corresponding boundary conditions are

(6)

where c1 and c2 are positive real numbers and 1 − c2t > 0. Consider the dimensionless transformations19 

(7)

where η is the dimensionless local similarity variable and ψ(x, y) is the stream function defined by (v1,v2)=(ψy,ψx) and ν is the kinematic viscosity of the second grade fluid. The continuity equation (4) is automatically satisfied and from (2)–(6), one can write

(8)
(9)

Here S=c2c1 denotes the unsteadiness parameter, α=c1α1μ(1c2t) 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 

(10)

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 τw on the surface of the thin liquid film sheet is

(11)

and the local skin-friction coefficient or frictional drag coefficient is

(12)

In dimensionless form, we have

(13)

where Rex1/2=c1x2/ν(1c2t) is the local Reynold’s number. We will solve (10) and (11) analytically using VIM and OVIM in the next section.

The basic notion of the technique is illustrated by considering the problem of finding z(η) such that

(14)

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:

(15)

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 zp(s)̃ is a term being restricted which in turn gives δzp(s)̃=0. 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

(16)

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.

Consider a second order nonlinear ordinary differential equation,

(17)

subject to the boundary conditions

(18)

where L(η)=d2dη2 is a linear differential operator, N(η) represents the nonlinear operator and g(η) is an inhomogeneous term and R is set of real numbers. According to standard VIM, the correction functional used is given as

(19)

To make stable correction functional, the lagrange multiplier is defined as λ = sη, see Refs. 24 and 25, we get the following iterative formula

(20)

An unknown auxiliary parameter h is inserted into the iterative formula (20), for n = 0, (20) becomes

(21)

For optimal variational iteration method, we will proceed as follows

(22)

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

(23)
(24)
(25)

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).

In this section, we use VIM and OVIM to find an approximate solution of boundary value problem (7) and (8)for different values of S and α, while keeping β fixed.

For S = 0.2, α = 0, β = 1, Equations (8) and (9) become

(26)
(27)

Using VIM, we construct the correction functional for (26) as follows

(28)

Using the optimality conditions,25 one has

(29)

Using (28) and (29), we have

(30)

We assume

(31)

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

(32)

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

(33)

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.

FIG. 1.

Residual Error comparison of VIM and OVIM Solutions when α = 0 and S = 0.2.

FIG. 1.

Residual Error comparison of VIM and OVIM Solutions when α = 0 and S = 0.2.

Close modal

For S = 0.2, α = 0.9, β = 1, (8) and (9) become

(34)
(35)

We construct the correction functional for (34) as follows

(36)

Using the optimality conditions,25 one has

(37)

Using (37) in (36), We have the iterative scheme as

(38)

Using (31) and (38), we obtain the following iterations

Applying the boundary conditions (35) on f2(η), the unknowns are

We define the following residual error function

(39)

Now using OVIM, we proceed as follows

We obtain the following approximations

We define the following residual error function

(40)

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.

FIG. 2.

Residual Error Comparison of VIM and OVIM Solutions when α = 0.9 and S = 0.2.

FIG. 2.

Residual Error Comparison of VIM and OVIM Solutions when α = 0.9 and S = 0.2.

Close modal

Here we discuss the sway of different parameters like S and α on the components of velocity f, acceleration dfdη and skin-friction coefficients Rex1/2Cf. The change in the coefficient of the skin-friction Rex1/2Cf for α and S is shown in Table I.

TABLE I.

Values of skin-friction coefficient Rex1/2Cf for different parameters at η = 0.

αSRex1/2Cf
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 6.15999848254210 
0.5 0.2 4.81951263596841 
0.5 0.5 7.11968336959140 
0.5 4.98044449104253 
0.5 1.5 2.83681213996240 
0.5 0.742432181172167 
αSRex1/2Cf
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 6.15999848254210 
0.5 0.2 4.81951263596841 
0.5 0.5 7.11968336959140 
0.5 4.98044449104253 
0.5 1.5 2.83681213996240 
0.5 0.742432181172167 

Figure 3 and Figure 4 are displaying the changes made by α on the component of velocity f and dfdη. These figures depict that f and dfdη the function which are increasing depending upon α. The variation in dfdη 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 dfdη. It can be observed from the given figures that these velocities f and dfdη 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 Rex1/2Cf 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 Rex1/2Cf get decreased when S < 2, for S = 2 the coefficient of skin-friction Rex1/2Cf is almost zero and for S > 2 it is increased.

FIG. 3.

Influence of second grade parameter α on f.

FIG. 3.

Influence of second grade parameter α on f.

Close modal
FIG. 4.

Influence of second grade parameter α on dfdη.

FIG. 4.

Influence of second grade parameter α on dfdη.

Close modal
FIG. 5.

Influence of unsteadiness parameter S on dfdη.

FIG. 5.

Influence of unsteadiness parameter S on dfdη.

Close modal
FIG. 6.

Influence of unsteadiness parameter S on dfdη.

FIG. 6.

Influence of unsteadiness parameter S on dfdη.

Close modal

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.

1.
T.
Akbar
,
R.
Nawaz
,
M.
Kamran
, and
A.
Rasheed
, “
Magnetohydrodynamic (MHD) flow analysis of second grade fluids in a porous medium with prescribed vorticity
,”
AIP Advances
5
(
11
),
117133
(
2015
).
2.
H.
Andersson
,
J.
Aarseth
,
N.
Braud
, and
B.
Dandapat
, “
Flow of a power-law fluid film on an unsteady stretching surface
,”
Journal of Non-Newtonian Fluid Mechanics
62
(
1
),
1
8
(
1996
).
3.
H. I.
Andersson
,
J. B.
Aarseth
, and
B. S.
Dandapat
, “
Heat transfer in a liquid film on an unsteady stretching surface
,”
International Journal of Heat and Mass Transfer
43
(
1
),
69
74
(
2000
).
4.
A.
Aziz
, “
A similarity solution for laminar thermal boundary layer over a flat plate with a convective surface boundary condition
,”
Communications in Nonlinear Science and Numerical Simulation
14
(
4
),
1064
1068
(
2009
).
5.
A.
Aziz
and
T.
Aziz
, “
MHD flow of a third grade fluid in a porous half space with plate suction or injection: An analytical approach
,”
Applied Mathematics and Computation
218
(
21
),
10443
10453
(
2012
).
6.
C.-H.
Chen
, “
Heat transfer in a power-law fluid film over a unsteady stretching sheet
,”
Heat and Mass transfer
39
(
1
),
791
796
(
2003
).
7.
R.
Cortell
, “
Flow and heat transfer of a fluid through a porous medium over a stretching surface with internal heat generation/absorption and suction/blowing
,”
Fluid Dynamics Research
37
(
4
),
231
(
2005
).
8.
R.
Cortell
, “
A note on magnetohydrodynamic flow of a power-law fluid over a stretching sheet
,”
Applied Mathematics and Computation
168
(
1
),
557
566
(
2005
).
9.
R.
Cortell
, “
Effects of viscous dissipation and work done by deformation on the MHD flow and heat transfer of a viscoelastic fluid over a stretching sheet
,”
Physics Letters A
357
(
4
),
298
305
(
2006
).
10.
V.
Coscia
and
G. P.
Galdi
, “
Existence, uniqueness and stability of regular steady motions of a second-grade fluid
,”
International Journal of Non-Linear Mechanics
29
(
4
),
493
506
(
1994
).
11.
J. E.
Dunn
and
R. L.
Fosdick
, “
Thermodynamics, stability, and boundedness of fluids of complexity 2 and fluids of second grade
,”
Archive for Rational Mechanics and Analysis
56
(
3
),
191
252
(
1974
).
12.
H.
Ghaneai
and
M.
Hosseini
, “
Variational iteration method with an auxiliary parameter for solving wave-like and heat-like equations in large domains
,”
Computers and Mathematics with Applications
69
(
5
),
363
373
(
2015
).
13.
T.
Hayat
,
M. R.
Mohyuddin
, and
S.
Asghar
, “
Some inverse solutions for unsteanian fluid
,”
Tamsui Oxford Journal of Mathematical Sciences
21
(
1
),
1
20
(
2005
).
14.
J.-H.
He
, “
Variational iteration method a kind of non-linear analytical technique: some examples
,”
International Journal of Non-Linear Mechanics
34
(
4
),
699
708
(
1999
).
15.
J.-H.
He
, “
Notes on the optimal variational iteration method
,”
Applied Mathematics Letters
25
(
10
),
1579
1581
(
2012
).
16.
M.
Inokuti
,
H.
Sekine
, and
T.
Mura
, “
General use of the Lagrange multiplier in nonlinear mathematical physics1
,” in
S.
Nemat-Nasser
, editor
Variational Methods in the Mechanics of Solids
, pages
156
162
,
Pergamon
,
1980
.
17.
A.
Ishak
,
R.
Nazar
, and
I.
Pop
, “
Heat transfer over an unsteady stretching permeable surface with prescribed wall temperature
,”
Nonlinear Analysis: Real World Applications
10
(
5
),
2909
2913
(
2009
).
18.
N. S.
Khan
,
S.
Zuhra
,
Z.
Shah
,
E.
Bonyah
,
W.
Khan
, and
S.
Islam
, “
Slip flow of eyring-powell nanoliquid film containing graphene nanoparticles
,”
AIP Advances
8
(
11
),
115302
(
2018
).
19.
I.-C.
Liu
and
H. I.
Andersson
, “
Heat transfer in a liquid film on an unsteady stretching sheet
,”
International Journal of Thermal Sciences
47
(
6
),
766
772
(
2008
).
20.
O.
Makinde
and
A.
Aziz
, “
MHD mixed convection from a vertical plate embedded in a porous medium with a convective boundary condition
,”
International Journal of Thermal Sciences
49
(
9
),
1813
1820
(
2010
).
21.
S.
Mohammad Mehdi Hosseini
,
S.
Tauseef Mohyud-Din
, and
H.
Ghaneai
, “
Variational iteration method for Hirota-Satsuma coupled KdV equation using auxiliary parameter
,”
International Journal of Numerical Methods for Heat & Fluid Flow
,
22
(
3
),
277
286
(
2012
).
22.
S. T.
Mohyud-Din
,
W.
Sikander
,
U.
Khan
, and
N.
Ahmed
, “
Optimal variational iteration method for nonlinear problems
,”
Journal of the Association of Arab Universities for Basic and Applied Sciences
24
,
191
197
(
2017
).
23.
S. T.
Mohyud-Din
,
W.
Sikander
,
U.
Khan
, and
N.
Ahmed
, “
Optimal variational iteration method using Adomian’s polynomials for physical problems on finite and semi-infinite intervals
,”
The European Physical Journal Plus
132
(
5
),
236
(
2017
).
24.
M. A.
Noor
,
K. I.
Noor
,
M.
Rafiq
, and
E. A.
Al-said
, “
Variational iteration method for solving a system of second order boundary value problems
,”
International Journal of Nonlinear Sciences and Numerical Simulation
11
(
12
),
1109
1120
(
2010
).
25.
M. A.
Noor
,
M.
Rafiq
,
K. I.
Noor
,
E.
Al-Said
, and
A.
Khalifa
, “
Solving system of fourth-order boundary value problems by variational iteration method
,”
World Applied Sciences Journal
13
(
11
),
2312
2322
(
2011
).
26.
M.
Rafiq
,
H.
Ahmad
, and
S. T.
Mohyud-Din
, “
Variational iteration method with an auxiliary parameter for solving Volterras population model
,”
Nonlinear Science Letter A
8
(
4
),
389
396
(
2017
).
27.
K. R.
Rajagopal
and
A. S.
Gupta
, “
An exact solution for the flow of a non-Newtonian fluid past an infinite porous plate
,”
Meccanica
19
(
2
),
158
160
(
1984
).
28.
B.
Sakiadis
, “
Boundary-layer behavior on continuous solid surfaces: I. Boundary-layer equations for two-dimensional and axisymmetric flow
,”
7
,
26
28
(
1961
).
29.
B. C.
Sakiadis
, “
Boundary-layer behavior on continuous solid surfaces: II. The boundary layer on a continuous flat surface
,”
AIChE Journal
7
(
2
),
221
225
(
1961
).
30.
A.
Shahzad
and
R.
Ali
, “
Approximate analytic solution for magneto-hydrodynamic flow of a non-Newtonian fluid over a vertical stretching sheet
,”
Can J Appl Sci
2
(
1
),
202
215
(
2012
).
31.
A.
Shahzad
and
R.
Ali
, “
MHD flow of a non-Newtonian power law fluid over a vertical stretching sheet with the convective boundary condition
,”
Walailak Journal of Science and Technology (WJST)
10
(
1
),
43
56
(
2012
).
32.
A.
Shahzad
,
R.
Ali
,
M.
Hussain
, and
M.
Kamran
, “
Unsteady axisymmetric flow and heat transfer over time-dependent radially stretching sheet
,”
Alexandria Engineering Journal
56
(
1
),
35
41
(
2017
).
33.
M.
Sheikholeslami
,
H.
Ashorynejad
,
D.
Ganji
, and
A.
Y
ld
r
m
ııı, “
Homotopy perturbation method for three-dimensional problem of condensation film on inclined rotating disk
,”
Scientia Iranica
19
(
3
),
437
442
(
2012
).
34.
M.
Sheikholeslami
and
D.
Ganji
, “
Heat transfer of cu-water nanofluid flow between parallel plates
,”
Powder Technology
235
,
873
879
(
2013
).
35.
M.
Sheikholeslami
,
D.
Ganji
, and
H.
Ashorynejad
, “
Investigation of squeezing unsteady nanofluid flow using adm
,”
Powder Technology
239
,
259
265
(
2013
).
36.
M.
Sheikholeslami
and
D. D.
Ganji
, “
Nanofluid flow and heat transfer between parallel plates considering brownian motion using DTM
,”
Computer Methods in Applied Mechanics and Engineering
283
,
651
663
(
2015
).
37.
M.
Sheikholeslami
,
D. D.
Ganji
,
H. R.
Ashorynejad
, and
H. B.
Rokni
, “
Analytical investigation of Jeffery-Hamel flow with high magnetic field and nanoparticle by Adomian decomposition method
,”
Applied Mathematics and Mechanics
33
(
1
),
25
36
(
2012
).
38.
M.
Sheikholeslami
,
M.
Rashidi
,
D. M.
Al Saad
,
F.
Firouzi
,
H. B.
Rokni
, and
G.
Domairry
, “
Steady nanofluid flow between parallel plates considering thermophoresis and Brownian effects
,”
Journal of King Saud University-Science
28
(
4
),
380
389
(
2016
).
39.
M.
Shekholeslami
,
H.
Ashorynejad
,
D.
Domairry
, and
I.
Hashim
, “
Investigation of the laminar viscous flow in a semi-porous channel in the presence of uniform magnetic field using optimal homotopy asymptotic method
,”
Sains Malaysiana
41
(
10
),
1281
1285
(
2012
).
40.
A.
Siddiqui
and
P.
Kaloni
, “
Certain inverse solutions of a non-Newtonian fluid
,”
International journal of non-linear mechanics
21
(
6
),
459
473
(
1986
).
41.
T. W.
Ting
, “
Certain non-steady flows of second-order fluids
,”
Archive for Rational Mechanics and Analysis
14
(
1
),
1
26
(
1963
).
42.
M.
Turkyilmazoglu
, “
An optimal variational iteration method
,”
Applied Mathematics Letters
24
(
5
),
762
765
(
2011
).
43.
C.
Wang
, “
Analytic solutions for a liquid film on an unsteady stretching surface
,”
Heat and Mass transfer
42
(
8
),
759
766
(
2006
).
44.
C.
Wang
and
I.
Pop
, “
Analysis of the flow of a power-law fluid film on an unsteady stretching surface by means of homotopy analysis method
,”
Journal of Non-Newtonian Fluid Mechanics
138
(
2
),
161
172
(
2006
).
45.
C. Y.
Wang
, “
Liquid film on an unsteady stretching surface
,”
Quarterly of Applied Mathematics
48
(
4
),
601
610
(
1990
).