An investigation is presented for the two-dimensional and axisymmetric stagnation flows of a couple stress fluids intrude on a moving plate under partial slip conditions. The governing partial differential equations are converted into ordinary differential equations by a similarity transformation. The important physical parameters of skin friction coefficients of the fluid are also obtained. The homotopy analysis method (HAM) is employed to obtain the analytical solution of the problem. Also, the convergence of the solutions is established by plotting graphs of convergence control parameter. The impacts of couple stresses and slip conditions on the flow and temperature of the fluid have been observed. The numerical comparison for the considered fluid is compared with previous solutions as special case.
1. Introduction
The fluids exhibiting a boundary slip are important in industrial applications, for example, the polishing of artificial heart valves, rarefied fluid problems, and flow on multiple interfaces. There are many cases where no slip condition is replaced with Navier’s partial slip condition. Partial slip condition on solid boundary occurs in many problems such as oscillatory flow channel, transient flow, some coated surfaces, some rough or porous surfaces, and heat transfer on moving plate. The flow on a moving plate is termed as a basic content for convection processes. The partial slip condition on a moving plate was considered by Wang [1]; the steady, laminar, axis-symmetric flow of a Newtonian fluid due to a stretching sheet with partial slip was studied by Ariel [2], Nadeem et al. [3] investigated steady state rotating and MHD flow of a third grade fluid past a rigid plate with slip; flow and heat transfer of a non-Newtonian fluid past a stretching sheet with partial slip are considered by Sahoo [4], and Jamil and Khan [5] considered the slip effects on fractional viscoelastic fluids; the steady boundary layer flow past a moving horizontal flat plate with a slip effect is studied by Kumaran and Pop [6].
The theory of couple stresses, introduced by Stokes [7], explain the rheological behavior of various complex non-Newtonian fluids with body stresses and body couples which cannot be illustrated by the classical theory of continuum mechanics. Due to the rotational interaction of particles, the force-stress tensor is not symmetric and flow behaviors of such fluids are not similar to the Newtonian ones. It draws the researcher’s attention with the growing applications of such fluids in engineering, biomedical, and chemical industries. The peristaltic transport of a couple stress fluid in an asymmetric channel with the effect of the induced magnetic field has been considered by Nadeem and Akram [8], Khan et al. [9] investigate the approximate solution of couple stress fluid with expanding or contracting porous channel, Ramana Murthy and Nagaraju considered the flow of a couple stress fluid generated by a circular cylinder subjected to longitudinal and torsional oscillations [10], Hayat and colleagues observed that the velocity and the boundary layer thickness are decreasing functions of the couple stress fluid parameter in his study of melting heat transfer in the boundary layer flow of a couple stress fluid over a stretching surface [11], and the time dependence of the run up flow of a couple stress fluid between rigid parallel plates is examined by Devakar and Iyengar [12] in which the flow is induced by a constant pressure gradient which is suddenly withdrawn and the parallel plates are set to move instantaneously with different velocities in the direction of the applied pressure gradient.
In present paper, the effect of velocity slip factor λ and thermal slip factor β in two- dimensional and axisymmetric stagnation flows of a couple stress fluid are considered. Heat transfer analysis has also been taken into account. The problem is solved and analyzed with the application of HAM. The HAM provides a simple way to ensure convergence of the solution. It does not only optimize the solution but also gives freedom to choose the base function. The homotopy analysis method contains the auxiliary parameter ħ, which provides a straightforward way to adjust and control the convergence region of the series solution. The HAM can overcome the restrictions and limitations of perturbation methods [13, 14]. It has been successfully applied to many types of nonlinear problems such as the nonlinear equations arising in heat transfer [15], the nonhomogeneous Blasius problem [16], the generalized three-dimensional MHD flow over a porous stretching sheet [17], the axisymmetric flow, and heat transfer of a second grade fluid past a stretching sheet [18], and Kawahara equation [19], and an analytical study for the effects of slip factors on unsteady stagnation point flow and heat transfer towards a stretching sheet has also been carried out through HAM [20]. The solution can be expressed with the help of different linear operators, out of which some converge to the exact solution faster than others.
2. Formulation of the Problem
Let us consider a two-dimensional stagnation flow with heat transfer in the x-z plane moving with velocity U in the x-direction and velocity V in y-direction (as shown in Figure 1). The partial slip condition on the surface of moving plate for both motion of the fluid and heat transfer is also considered. Under the considered assumptions, in the absence of body forces and body couples, the momentum and energy equations of an incompressible couple stress fluid are presented by the following equations:
(1)∇·V=0,ρ(∂V∂t+(V·∇)V)=-∇p+μ∇2V-γ∇4V,ρCp(∂T∂t+(V·∇)T)=k∇2T.
The potential flow far from the plate is given by
(2)u=ax,v=0,w=-az,p=p0-ρa2(x2+z2)2,
where, u,v,w are velocity components in the Cartesian x,y,z directions, a is the strength of the stagnation flow, ρ is the density, p is the pressure, and p0 is the stagnation pressure.
For couple stress flow, we set
(3)u=axf′(η)+Ug(η),v=Vh(η),w=-aνf(η),θ(η)=T-T∞T0-T∞,η≡aνz,(4)p=p0-ρ(a2x22+w22-νwz+γwz,z,z),
where γ is the couple stress fluid parameter and ν is the kinematic viscosity. The subscript z represents differentiation with respect to z. Using (3)-(4), the momentum and energy equations are reduced to ordinary differential equations
(5)f′′′+ff′′-(f′)2-ξf(5)+1=0,g′′+fg′-f′g-ξg(4)=0,h′′+fh′-ξh(4)=0,θ′′+Prfθ′=0,
where Pr=ν/α is Prandtl number and ξ=aγ/(ν2ρ) is the couple stress parameter.
On the plate, velocity slip and thermal slip similar to Navier’s condition is
(6)u-U=N(ρνuz-γuzzz),v-V=N(ρνvz-γvzzz),T-T0=STz,
where N and S are slip constants. Substitution of (3) yields
(7)f′(0)=λ(f′′(0)-ξf′′′′(0)),g(0)=1+λ(g′(0)-ξg′′′(0)),h(0)=1+λ(h′(0)-ξh′′′(0)),θ(0)=1+βθ′(0),
where λ=Nρνa is the nondimensional slip factor and β=Sa/ν is the thermal slip factor. For an impermeable plate and a potential flow which is recovered at infinity, we have
(8)f(0)=0,f′(∞)=1,g(∞)=0,h(∞)=0,θ(∞)=0.
2.2. The Axisymmetric Stagnation Point Flow
The stagnation flow in the x-z plane is moving two-dimensionally with velocity U in the x direction and velocity V in the y direction. The potential flow far from the plate is given by
(9)u=ax,v=ay,w=-2az,p=p0-ρ[(x2+z2)2+2a2z2].
We introduce the following transformations:(10)u=axf′(η)+Ug(η),v=ayf′(η),w=-2aνf(η),p=p0-ρ[a2(x2+y2)2+w22+νwz+γwz,z,z].
The equations of motion and energy are reduced to ordinary differential equations:
(11)f′′′+2ff′′-(f′)2-ξf(5)+1=0,g′′+2fg′-f′g-ξg(4)=0,θ′′+2Prfθ′=0
with the same boundary conditions as in (7)-(8).
The shear stress of the fluid for two-dimensional and axisymmetric flows can be found as
(12)τxz=(μ(∂u∂z)-ξ(∂3u∂z3))|z=0,τyz=(μ(∂v∂z)-ξ(∂3v∂z3))|z=0.
Hence, the local skin friction coefficients are given by
(13)Rex1/2Cf=f′′(0)-ξf′′′′(0),Cg=g′(0)-ξg′′′(0)U2/aν,Cg=h′(0)-ξh′′′(0)V2/aν,
where Rex1/2=ax2/ν is the local Reynolds numbers.
3. Solution of the Problem
For the two-dimensional problem presented in (5), the course of action for the HAM solution, we select
(14)f0=-11+λ+η+e-η1+λ,g0=e-η1+λ,h0=e-η1+λ,θ0=e-η1+β
as the initial approximation of f0,g0,h0, and θ0, respectively, which satisfy the following linear operator and corresponding boundary conditions:
(15)ℒ1=d3dη3+d2dη2,ℒ2=d2dη2-1
such that
(16)ℒ1[c1+c2eη+c3e-η]=0,ℒ2[c4eη+c5e-η]=0,
where c1,c2,…,c5 are arbitrary constants.
The following nonlinear operators are defined as
(17)N1[f~(η,p)]=∂3f~(η,p)∂η3+f~(η,p)∂2f~(η,p)∂η2-(∂f~(η,p)∂η)2-ξ∂5f~(η,p)∂η5+1,N2[f~(η,p),g~(η,p)]=∂2g~(η,p)∂η2+f~(η,p)∂g~(η,p)∂η-∂f~(η,p)∂ηg~(η,p)-ξ∂4g~(η,p)∂η4,N3[f~(η,p),h~(η,p)]=∂2h~(η,p)∂η2+f~(η,p)∂h~(η,p)∂η-ξ∂4h~(η,p)∂η4,N4[f~(η,p),θ~(η,p)]=∂2θ~(η,p)∂η2+Prf~(η,p)∂θ~(η,p)∂η
and then construct the zeroth-order deformation equations:
(18)(1-p)ℒ1[f~(η,p)-f0(η)]-pħH1(η)N1[f~(η,p)]=0,(1-p)ℒ2[g~(η,p)-g0]-pħH2(η)N2[f~(η,p),g~(η,p)]=0,(1-p)ℒ2[h~(η,p)-h0]-pħH3(η)N3[f~(η,p),h~(η,p)]=0,(1-p)ℒ2[θ~(η,p)-θ0(η)]-pħH4(η)N4[f~(η,p),θ~(η,p)]=0,f~′(0,p)=λ(f~′′(0,p)-ξf~′′′′(0,p)),f~(0,p)=0,g~(0,p)=1+λ(g~′(0,p)-ξg~′′′(0)),h~(0,p)=1+λ(h~′(0,p)-ξh~′′′(0)),θ~(0,p)=1+βθ~′(0,p),f~′(∞,p)=1,g~(∞,p)=0,h~(∞,p)=0,θ~(∞,p)=0
in which p∈[0,1] is the embedding parameter and ħ is the auxiliary nonzero parameter and Hi(η),i=1,2,3,4 are the auxiliary functions. Here, we selected the auxiliary functions as Hk=e-η,k=1,2,3,4.
By Taylor’s theorem,
(19)f(η,p)=f0(η)+∑m=1∞fm(η)pm,g(η,p)=g0(η)+∑m=1∞gm(η)pm,h(η,p)=h0(η)+∑m=1∞hm(η)pm,θ(η,p)=θ0(η)+∑m=1∞θm(η)pm.
The mth-order deformation problems with the corresponding boundary conditions are given by
(20)ℒ1[fm-χmfm-1]=ħe-ηR1[fm-1],ℒ2[gm-χmgm-1]ħe-ηR2[fm-1,gm-1],ℒ2[hm-χmhm-1]=ħe-ηR3[fm-1,hm-1],ℒ2[θm-χmθm-1]=ħe-ηR4[fm-1,θm-1],(21)fm′(0)=λ(fm′′(0)-ξf′′′′(0)),fm(0)=0,gm(0)=λ(gm′(0)-ξg′′′(0)),hm(0)=λ(hm′(0)-ξh′′′(0)),θm(0)=βθm′(0),fm′(∞)=0,gm(∞)=0,hm(∞)=0,θm(∞)=0,
where
(22)R1[fm-1]=fm-1′′′+∑i=0m-1fifm-1-i′′-∑i=0m-1fi′fm-1-i′-ξfm-1(5)+(1-χm),R2[fm-1,gm-1]=gm-1′′+∑i=0m-1figm-1-i′-∑i=0m-1fi′gm-1-i-ξgm-1(4),R3[fm-1,hm-1]=hm-1′′+∑i=0m-1fihm-1-i′-ξhm-1(4),R4[fm-1,θm-1]=θm-1′′+Pr∑i=0m-1fiθm-1-i′,χm={0,m≤1,1,m≥2.
According to the procedure defined previously, it is easy to solve the linear equations (20), one after the other in the order m, especially by means of any computation software, such as Mathematica.
In a similar manner the solution of axisymmetric problem defined in (11) can also be determined by taking the same linear operators and initial guesses and proceeding in the same way as for two-dimensional problem.
4. Results and Discussions
The paper presents an analysis of the flow and heat transfer of the stagnation slip flow of a couple stress fluid over a moving plate. Approximate analytical solutions are derived for momentum and energy equations. We found analytic solutions corresponding to the heat transfer of stagnation slip flow over a moving plate by using symbolic software Mathematica 9. The HAM is adopted to achieve the analytic solution of the problem. In HAM, a deformation equation defines a continuous variation from an initial approximation to the solution through an embedding parameter. The convergence region of the obtained solution can be determined with the help of embedding parameter ħ which provides the straight line parallel to x-axis in which the resultant solution is convergent. For the two-dimensional flow, Figures 2(a)–2(c), the ħ curves are plotted for f′′(0), g′(0) and h′(0), and for axisymmetric flow, Figures 3(a) and 3(b) are plotted for f′′(0) and g′(0) that presents the admissible values of ħ under which the solution is convergent.
The ħ curves for 6th order of approximation of HAM when λ=0.2, ξ=1. (a) for f′′(0), (b) for g′(0), (c) for h′(0).
The ħ curves for 6th order of approximation of HAM when λ=0.2, ξ=1. (a) for f′′(0), (b) for g′(0).
For both two-dimensional and axisymmetric stagnation flows, the effect of couple stress parameter ξ on temperature and velocity profiles, the effects of velocity slip factor λ on velocity components, and the effect of thermal slip factor β on temperature, which in turn depends on velocity, are presented in the form of graphical and tabulated results. In Tables 1, 2, and 3 (for two-dimensional) and Tables 4 and 5 (for axisymmetric), the initial values of ReCf, (U2/av)Cg, and (V2/av)Cn are presented for the various values of ξ and λ. At ξ=0 the problem recovers the Newtonian flow problem of Wang [1]. The results in the tables are in agreement with the results of Wang [1] (for ξ=0) and follow the same trend for increasing slip factor λ when there are couple stresses in the fluid. Figures 4 and 5 represent the influence of velocity slip factor λ on f′ and f. An increase in the slip factor λ decelerates the velocity f′ (see Figure 4). The increase in the function f is observed with the increasing in slip factor as shown in Figure 5. It can be noticed from Figure 6 that positive increase in the couple stress parameter ξ has the negative impact on the velocity f′ of the couple stress fluid; it influences more on axisymmetric flow as compared to two-dimensional flow, whereas couple stress parameter ξ decelerates the velocity function f for both two-dimensional and axisymmetric flows (see Figure 7).
Initial values of ReCf for the two-dimensional stagnation flow when ℏ=-0.61 (bracket values from [1]).
λ∖ξ
0
0.1
0.2
0.3
0.4
0
1.23260 (1.23259)
1.43251
2.58260
5.25453
9.72257
0.2
1.04385 (1.04259)
1.17052
1.94111
3.82149
7.15871
0.5
0.82483 (0.82148)
0.90095
1.38118
2.60405
4.90074
1
0.59637 (0.59346)
0.63885
0.90465
1.60141
2.97625
2
0.37692 (0.37589)
0.39660
0.51389
0.82142
1.44177
10
0.094047 (0.094036)
0.09559
0.10423
0.12349
0.14926
20
0.048476 (0.048472)
0.04889
0.05126
0.05608
0.05987
50
0.019753 (0.019752)
0.01982
0.02022
0.02097
0.021125
Initial values of U2/aνCg for the two dimensional stagnation flow when ℏ=-0.35 (bracket values from [1]).
λ∖ξ
0
0.1
0.2
0.3
0.4
0
−0.83406 (−0.81130)
−0.98987
−1.37103
−2.07901
−3.19268
0.2
−0.75213 (−0.77521)
−0.84363
−1.09389
−1.59464
−2.43874
0.5
−0.63674 (−0.67196)
−0.68416
−0.83462
−1.16169
−1.75519
1
−0.49511 (−0.52189)
−0.51574
−0.59259
−0.78238
−1.15175
2
−0.33602 (−0.34911)
−0.34277
−0.37504
−0.45978
−0.64203
10
−0.091548 (−0.09248)
−0.09184
−0.09389
−0.10003
−0.11507
20
−0.047812 (−0.04806)
−0.04788
−0.048423
−0.05008
−0.05418
50
−0.019642 (0.019685)
−0.01965
−0.019744
−0.02002
−0.02070
Initial values of V2/aνCh for the two dimensional stagnation flow when ℏ=-0.5 (bracket values from [1]).
λ∖ξ
0
0.1
0.2
0.3
0.4
0
−0.68635 (−0.57047)
−0.68046
−0.51044
0.02198
1.18096
0.2
−0.62273 (−0.55440)
−0.63793
−0.59093
−0.34999
0.31703
0.5
−0.53707 (−0.50283)
−0.56007
−0.58301
−0.54172
−0.27375
1
−0.43014 (−0.41618)
−0.45097
−0.49846
−0.561178
−0.56149
2
−0.30360 (−0.29951)
−0.31649
−0.35525
−0.43436
−0.54462
10
−0.08882 (−0.08871)
−0.09018
−0.09528
−0.10908
−0.13932
20
−0.047054 (−0.04703)
−0.04744
−0.04896
−0.05324
−0.06318
50
−0.01951 (−0.01951)
−0.019580
−0.01985
−0.020635
−0.02252
Initial values of ReCf for the axisymmetric flow when ℏ=-0.59 (bracket values from [1]).
λ∖ξ
0
0.1
0.2
0.3
0.4
0
1.31166 (1.31193)
1.35335
1.50955
1.61316
1.22802
0.2
1.11399 (1.11097)
1.17608
1.44234
1.85177
2.07303
0.5
0.873726 (0.86688)
0.94076
1.24538
1.82827
2.50949
1
0.622290 (0.6173)
0.67807
0.94791
1.54311
2.44917
2
0.386430 (0.38526)
0.41918
0.59476
1.03573
1.84473
10
0.094590 (0.09459)
0.09752
0.11831
0.18272
0.33557
20
0.048630 (0.04862)
0.04943
0.05557
0.07543
0.12409
50
0.019780 (0.01977)
0.01992
0.02100
0.02463
0.03398
Initial values of U2/aνCg for the axisymmetric flow when ℏ=-0.3 (bracket values from [1]).
λ∖ξ
0
0.1
0.2
0.3
0.4
0
−0.87723 (−0.93873)
−0.98615
−1.17030
−1.43528
−1.77243
0.2
−0.78757 (−0.88766)
−0.86814
−1.02163
−1.26119
−1.58740
0.5
−0.66257 (−0.74987)
−0.72114
−0.84428
−1.05073
−1.35156
1
−0.51106 (−0.56453)
−0.55219
−0.64428
−0.80782
−1.06257
2
−0.343538 (−0.36643)
−0.36980
−0.42972
−0.54027
−0.72341
10
−0.092122 (−0.09356)
−0.09888
−0.11326
−0.13908
−0.18304
20
−0.047969 (−0.04835)
−0.05146
−0.05867
−0.07123
−0.09205
50
−0.019669 (−0.01973)
−0.02109
−0.02396
−0.02883
−0.03662
Effect of slip factor λ on f′(η) (slip and no-slip cases) when ħ=-0.3, ξ=0.5. (a) for two-dimensional and (b) for axisymmetric case.
Effect of slip factor λ on f(η) (slip and no-slip cases) when ħ=-0.2, λ=0.2 (a) for two-dimensional, and (b) for axisymmetric case.
Effect of couple stress parameter ξ on f′(η) (Newtonian and couple stress cases) when ħ=-0.2, λ=0.2. (a) for two-dimensional and (b) for axisymmetric cases.
Effect of couple stress parameter ξ on f(η) (Newtonian and couple stress cases) when ħ=-0.2, λ=0.2. (a) for two-dimensional and (b) for axisymmetric.
Figure 8 illustrates the influence of velocity slip factor λ (for both two-dimensional and axisymmetric flows) on the velocity component g. An increase in the slip factor λ results in an exponential decay of g; the effect of couple stress parameter ξ on the velocity component g reveals that it supports the two-dimensional flow and the axisymmetric flow in axial direction (see Figure 9). From Figure 10, it is noted that the influence of velocity slip factor λ on the translational velocity h is a reducing function of slip factor λ. The effect of couple stress parameter ξ on the translational velocity h is plotted in Figure 11; it is observed that it is also a decreasing function of couple stress parameter.
Effect of slip factor λ on g(η) (slip and no-slip cases) when ħ=-0.2, ξ=0.5. (a) for two-dimensional flow and (b) for axisymmetric flow.
Effects of couple stress parameter ξ on g(η) (Newtonian and couple stress cases) when λ=0.2. (a) for two-dimensional when ħ=0.3 and (b) for axisymmetric flow when ħ=0.3.
Effect of slip factor λ on h(η) (slip and no-slip cases) when ħ=-0.4, ξ=0.3.
Effect of couple stress parameter ξ on h(η) (Newtonian and couple stress cases) when ħ=0.2, λ=0.2.
Figures 12(a), 12(b), and 12(c) for two-dimensional flow and Figures 13(a), 13(b), and 13(c) for axisymmetric flow presented the temperature profiles for various values of thermal slip β. It has been analyzed that the increment in the thermal slip β tends to decrease the temperature of the couple stress fluid. From the definition of Prandtl number, it is quite obvious that a large Prandtl number has a lower thermal diffusivity; therefore, an increase in Pr tends to decrease the temperature.
Temperature profile θ(η) of two-dimensional flow for various values of thermal slip β (slip and no-slip cases) when ħ=-0.5, ξ=0.5. (a) for Pr=7. (b) for Pr=0.7. (c) for Pr=0.07.
Temperature profile θ(η) of axisymmetric flow for various values of thermal slip β (slip and no-slip cases) when ħ=-0.5, ξ=0.5. (a) for Pr=7. (b) for Pr=0.7. (c) for Pr=0.07.
The consequences of an increase in Prandtl number Pr and thermal slip factor β on the heat transfer are visualized and compared with the viscous fluid considered by Wang [1] in Figures 14, 15, and 16. It is seen that the heat transfer is increasing function of Pr and λ and decreasing function of β. It has also been observed that the presence of couple stresses decreases the heat transfer as compared to Newtonian fluid. The heat transfer coefficient for axisymmetric flow is smaller as compared with large value of -θ′(0) for two-dimensional flow.
Profile of temperature gradient at the plate -θ′(0) against λ when Pr=7, ξ=0.3, ħ=0.5. (a) for two-dimensional. (b) for axisymmetric.
Profile of temperature gradient at the plate -θ′(0) against λ when Pr=0.7, ξ=0.3, ħ=0.5. (a) for two-dimensional. (b) for axisymmetric.
Profile of temperature gradient at the plate -θ′(0) against β when Pr=0.07, ξ=0.3, ħ=0.5. (a) for two-dimensional (b) for axisymmetric.
5. Conclusions
The study has been done to explore the influence of couple stresses with partial slip condition on the flow pattern of couple stress fluid with heat transfer.
The trends for the physical parameters of skin frictions, Cf,Cg, and Ch for two-dimensional and axisymmetric flows are decreasing function of slip factor.
The velocity functions f′ and f are increasing functions of slip factor λ(0≤λ<∞).
The axial velocity g and translational velocity h are decreasing functions of slip factor λ.
Couple stress parameter ξ has the negative impact on velocity functions f′, f on both two-dimensional and axisymmetric flows.
The axial velocity g increases in two-dimensional flow and decreases in axisymmetric flow for the increasing values of couple stress parameter ξ.
The velocity function h found no support from the couple stresses in the fluid for the flow in translational direction in both type of flows.
The temperature of the couple stress fluid is the decreasing function of thermal slip β.
Prandtl number has the positive effect on the heat transfer of the fluid.
Notations(x,y,z):
Space coordinates
(u,v,w):
Velocity components
ν:
Kinematic viscosity
ρ:
Density
k:
Thermal conductivity
cp:
Specific heat
T:
Temperature of fluid
γ:
Couple stress parameter
N,S:
Slip constants
Pr:
Prandtl number
ξ:
Dimensionless couple stress parameter
a,U,V:
Dimensional rate constant
λ:
Nondimensional slip factor
β:
Thermal slip factor
Re:
Reynolds number
Cf, Cg, Ch:
Skin friction coefficients.
Conflict of Interests
The authors declare that they have no conflict of interests.
Authors’ Contribution
The work including proofreading was done by all the authors.
Acknowledgment
The author Najeeb Alam Khan is thankful and grateful to the Dean of Faculty of Sciences, University of Karachi, Karachi, Pakistan, for supporting and facilitating this research work.
WangC. Y.Stagnation slip flow and heat transfer on a moving plate20066123766876722-s2.0-3375049046410.1016/j.ces.2006.09.003ArielP. D.Axisymmetric flow due to a stretching sheet with partial slip2007547-8116911832-s2.0-3474888548810.1016/j.camwa.2006.12.063NadeemS.HayatT.AbbasbandyS.AliM.Effects of partial slip on a fourth-grade fluid with variable viscosity: an analytic solution20101128568682-s2.0-7044964353010.1016/j.nonrwa.2009.01.030SahooB.Flow and heat transfer of a non-Newtonian fluid past a stretching sheet with partial slip20101536026152-s2.0-7034980490610.1016/j.cnsns.2009.04.032JamilM.KhanN. A.Slip effects on fractional viscoelastic fluids201120111919381310.1155/2011/193813KumaranV.PopI.Nearly parallel blasius flow with slip20111612461946242-s2.0-7996027661810.1016/j.cnsns.2011.03.045StokesV. K.Couple stresses in fluids196699170917152-s2.0-36849106425NadeemS.AkramS.Peristaltic flow of a couple stress fluid under the effect of induced magnetic field in an asymmetric channel2011811971092-s2.0-7865059239010.1007/s00419-009-0397-8KhanN. A.MahmoodA.AraA.Approximate solution of couple stress fluid with expanding or contracting porous channel201330339940810.1108/02644401311314358Ramana MurthyJ. V.NagarajuG.Flow of a couple stress fluid generated by a circular cylinder subjected to longitudinal and torsional oscillations2009210451461HayatT.MustafaM.IqbalZ.AlsaediA.Stagnation-point flow of couple stress fluid with melting heat transfer2013342167176DevakarM.IyengarT. K. V.Run up flow of a couple stress fluid between parallel plates201015129372-s2.0-77949401655LiangS.JeffreyD. J.Comparison of homotopy analysis method and homotopy perturbation method through an evolution equation20091412405740642-s2.0-6734916253810.1016/j.cnsns.2009.02.016SajidM.HayatT.AsgharS.Comparison between the HAM and HPM solutions of thin film flows of non-Newtonian fluids on a moving belt2007501-227352-s2.0-3454823055410.1007/s11071-006-9140-yAbbasbandyS.The application of homotopy analysis method to nonlinear equations arising in heat transfer200636011091132-s2.0-3375033927010.1016/j.physleta.2006.07.065AllanF. M.SyamM. I.On the analytic solutions of the nonhomogeneous Blasius problem200518223623712-s2.0-1964439486510.1016/j.cam.2004.12.017HayatT.JavedT.On analytic solution for generalized three-dimensional MHD flow over a porous stretching sheet20073703-42432502-s2.0-3514888537410.1016/j.physleta.2007.05.108HayatT.SajidM.Analytic solution for axisymmetric flow and heat transfer of a second grade fluid past a stretching sheet2007501-275842-s2.0-3375071442910.1016/j.ijheatmasstransfer.2006.06.045AbbasbandyS.Homotopy analysis method for the Kawahara equation20101113073122-s2.0-7035073448710.1016/j.nonrwa.2008.11.005KhanN. A.JamilM.NadeemK. A.Effects of slip factors on unsteady stagnation point flow and heat transfer towards a stretching sheet: an analytical study201243877979410.1615/HeatTransRes.2012004962