# Analytical Solution of Flow and Heat Transfer over a Permeable Stretching Wall in a Porous Medium

Authors: M. Dayyan; S. M. Seyyedi; G. G. Domairry; M. Gorji Bandpy
Journal: Mathematical Problems in Engineering
 (2013)

---

## Abstract

Boundary layer flow through a porous medium over a stretching porous wall  has seen solved with analytical solution. It has been considered two wall boundary conditions which are power-law distribution of either wall temperature or heat flux. These are general enough to cover the isothermal and isoflux cases. In addition to momentum, both first and second laws of thermodynamics analyses of the problem are investigated. The governing equations are transformed into a system of ordinary differential equations. The transformed ordinary equations are solved analytically using homotopy analysis method. A comprehensive parametric study is presented, and it is shown that the rate of heat transfer increases with Reynolds number, Prandtl number, and suction to the surface.


---

## Body

## 1. Introduction

Boundary layer flows over a stretching surface have application in engineering processes such as liquid composite molding, extrusion of plastic sheets, paper production, glass blowing, metal spinning, wire drawing, and hot rolling [1–3]. More importantly, the quality of the products, in the abovementioned processes, depends on the kinematics of stretching and the simultaneous heat and mass transfer rates during the fabrication process. Sakiadis [4, 5] and Crane [6] were the pioneers in the investigations of boundary layer flow over continuously moving surfaces that are quite different from the free stream flow over stationary flat plates. Pop and Na [7] studied free convection heat transfer of non-Newtonian fluids along a vertical wavy surface in a porous medium.The flow field of a stretching surface with a power-law velocity variation was discussed by Banks [8]. Elbashbeshy [9] investigated heat transfer over a stretching surface with variable surface heat flux. Elbashbeshy and Bazid [10] analyzed heat transfer in a porous medium over a stretching surface with internal heat generation and suction or injection. This work was extended by Cortell [11] to include power-law temperature distribution. Steady flow and heat transfer of a viscous incompressible fluid flow through porous medium over a stretching sheet were studied by Sriramalu et al. [12]. Ali [13] investigated thermal boundary layer on a power-law stretched surface with suction or injection for uniform and variable surface temperatures. Recently, Pantokratoras [14] published analytical solution for velocity distribution inside a Darcy-Brinkman porous medium with a stretching boundary. Elbashbeshy [15] included thermal radiation effects in a stretching surface problem. Tamayol and Bahrami [16] understood that porous materials can be used to enhance the heat transfer rate from stretching surfaces to improve processes such as hot rolling and composite fabrication.In view of the above, an analytical solution is developed in the present study to solve fluid flow, heat transfer in a porous medium over a porous plate with linear velocity, subjected to different power-law thermal boundary conditions. The analytical solution is successfully validated in comparison with numerical analysis. There are many effective methods for obtaining the solutions of nonlinear equation such as variational iteration method [17], Adomian method [18], and homotopy perturbation method [19, 20], and we use one of a powerful technique, namely, the homotopy analysis method (HAM) which was expressed by Liao [21–24]. This method has been successfully applied to solve many types of nonlinear problems [25–29] and provides us with great freedom to express solutions of a given nonlinear problem by means of different base functions. Secondly, unlike all previous analytic techniques, the homotopy analysis method always provides us with a family of solution expressions in the auxiliary parameter ℏ, even if a nonlinear problem has a unique solution. Thirdly, unlike perturbation techniques, the homotopy analysis method is independent of any small or large quantities. So, the homotopy analysis method can be applied no matter if governing equations and boundary/initial conditions of a given nonlinear problem contain small or large quantities or not.Above all, there are no rigorous theories to direct us to choose the initial approximations, auxiliary linear operators, auxiliary functions, and auxiliary parameterℏ. From the practical viewpoints, there are some fundamental rules such as the rule of solution expression, the rule of coefficient ergodicity, and the rule of solution existence, which play important roles within the homotopy analysis method. Unfortunately, the rule of solution expression implies such an assumption that we should have, more or less, some knowledge about a given nonlinear problem a prior. So, theoretically, this assumption impairs the homotopy analysis method, although we can always attempt some base functions even if a given nonlinear problem is completely new for us.


## 2. Governing Equation

Consider a steady, constant property, two-dimensional flow through a homogenous porous medium of permeabilityK, over a stretching surface with linear velocity distribution; that is, uw=u0x/L (Figure 1). The transport properties of the medium can be assumed independent of temperature when the difference between wall and ambient temperatures is not significant [4]. The origin is kept fixed, while the wall is stretching, and the y-axis is perpendicular to the surface. The governing equations are [4, 5]
(1)∂u∂x+∂υ∂y=0,ρ(u∂u∂x+υ∂u∂y)=μeff∂2u∂y2-μK′u,u∂T∂x+υ∂T∂y=αeff∂2T∂y2,
where uand υ are velocity components in the x and y directions, respectively. μeff is the effective viscosity which for simplicity in the present study is considered to be identical to the dynamic viscosity, μ. The transport properties of the porous medium such as permeability depend on their microstructure and can be calculated either using existing correlations in the literature or through experimental measurements. αeff is the effective thermal diffusivity of the medium. The hydrodynamic boundary conditions are
(2)u(x*,0)=u0x*,υ(x*,0)=υw,u(x*,∞)=0,
where x*=x/L is the nondimensional x-coordinate and L is the length of the porous plate. The following thermal boundary conditions are considered:
(3)T(x*,0)=T∞+T0(x*)n,T(x*,∞)=T∞,-κ∂T∂y|(x*,0)=q0(x*)n,T(x*,∞)=T∞,
where κ is the effective thermal conductivity of the medium and is a function of thermal conductivities of the fluid and solid phases and the porous medium microstructure.Figure 1
Schematic diagram of problem.Using dimensionless parameters,(4)η=yK,u=u0x*f′(η),υ=-u0LKf(η),ψ=u0x*Kf(η),
where f′ is df/dη.The transformed nonlinear ordinary differential equations are(5)f′′′+Re(ff′′-f′2)-f′=0,θ′′+Re·Pr(fθ′)-n·Re·Pr(f′θ)=0,
where Re=ρu0K/Lμ is the Reynolds number. Equation (5) should be solved subject to the following boundary conditions:
(6)f′(0)=1,f(0)=-υwLu0K=fw,f′(∞)=0,(7a)θ(0)=1,θ(∞)=0Power-law temperature,(7b)θ′(0)=-1,θ′(∞)=0Power-law heat flux,where fw is the injection parameter. Positive/negative values of fw show suction/injection into/from the porous surface, respectively. The wall shear stress term can be calculated as
(8)τL=-μu0Lf′′(0)2K.For power-law fluid wit constant temperature and heat flux boundary conditions, respectively. Employing the definition of convective heat transfer coefficient, the local Nusselt numbers become(9)Nux=hxk={-θ′(0)xK,Power-law temperature,qwxk(Tw-T∞)=xθ(0)K,Power-law heat flux.


## 3. Solution of Problem by Homotopy Analysis Method (HAM)

As mentioned by Liao, a solution may be expressed with different base functions, among which some converge to the exact solution of the problem faster than others. Such base functions are obviously better suited for the final solution to be expressed in terms of. Noting these facts, we have decided to expressg(η) by a set of base functions of the following form:
(10)f(η)=∑m=0∞∑n=0∞∑k=0∞bm,nkηkexp(-nη),θ(η)=∑m=0∞∑n=0∞∑k=0∞bm,nkηkexp(-nη).
The rule of solution expression provides us with a starting point. It is under the rule of solution expression that initial approximations, auxiliary linear operators, and the auxiliary functions are determined. So, according to the rule of solution expression, we choose the initial guess and auxiliary linear operator in the following form:
(11)f0(η)=fw+1-exp(-η),θ0(η)=exp(-η).
We choose linear operator as follows:
(12)ℒ1(f)=f′′′+f′′,ℒ2(θ)=θ′′+θ′,ℒ1(c1+c2η+c3exp(-η))=0,ℒ2(c4+c5exp(-η))=0,
where ci(i=1–5) are constants. Let P∈[0,1] denote the embedding parameter, and let ℏ indicate nonzero auxiliary parameters. Then, we construct the following equations.


### 3.1. Consider Zeroth–Order Deformation Equations

(13)
(
1
-
P
)
ℒ
1
[
f
(
η
;
p
)
-
f
0
(
η
)
]
=
p
ℏ
1
H
(
η
)
N
1
[
f
(
η
;
p
)
]
,
f
(
0
;
p
)
=
f
w
,
f
′
(
0
;
p
)
=
1
,
f
′
(
∞
;
p
)
=
0
,
(
1
-
P
)
ℒ
2
[
θ
(
η
;
p
)
-
θ
0
(
η
)
]
=
p
ℏ
2
H
(
η
)
N
2
[
θ
(
η
;
p
)
]
,
θ
(
0
;
p
)
=
1
,
θ
(
∞
;
p
)
=
0
Power-law  temperature,
θ
′
(
0
;
p
)
=
-
1
,
θ
(
∞
;
p
)
=
0
Power-law  heat  flux,
N
1
[
f
(
η
;
p
)
,
θ
(
η
;
p
)
]
=
∂
3
f
(
η
;
p
)
∂
η
3
+
Re
(
f
(
η
;
p
)
∂
2
f
(
η
;
p
)
∂
η
2
-
(
∂
f
(
η
;
p
)
∂
η
)
2
)
-
∂
f
(
η
;
p
)
∂
η
,
N
2
[
f
(
η
;
p
)
,
θ
(
η
;
p
)
]
=
∂
2
θ
(
η
;
p
)
∂
η
2
+
Re
·
Pr
(
f
(
η
;
p
)
∂
θ
(
η
;
p
)
∂
η
)
-
n
·
Re
·
Pr
(
θ
(
η
;
p
)
∂
f
(
η
;
p
)
∂
η
)
.
For p=0 and p=1, we have
(14)f(η;0)=f0(η),f(η;1)=f(η),θ(η;0)=θ0(η),θ(η;1)=θ(η).
When p increases from 0 to 1, then f(η;p) and θ(η;p) vary from f0(η) and θ0(η) to f(η) and θ(η). By Taylor’s theorem and using (14), we can write the following:
(15)f(η;p)=f0(η)+∑m-1∞fm(η)pm,fm(η)=1m!∂m(f(η;p))∂pm,θ(η;p)=θ0(η)+∑m-1∞θm(η)pm,θm(η)=1m!∂m(θ(η;p))∂pm.
In which ℏ1 and ℏ2 are chosen in such a way that these two series are convergent at p=1; therefore we have the following through (15):
(16)f(η)=f0(η)+∑m=1∞fm(η),θ(η)=θ0(η)+∑m=1∞θm(η).


### 3.2. ConsiderMth-Order Deformation Equations

(17)
ℒ
1
[
f
m
(
η
)
-
χ
m
f
m
-
1
(
η
)
]
=
ℏ
1
H
(
η
)
R
m
f
(
η
)
,
(18)
f
m
(
0
)
=
0
,
f
m
′
(
0
)
=
0
,
f
m
′
(
∞
)
=
0
,
(19)
R
m
f
(
η
)
=
f
m
-
1
′′′
+
∑
k
=
0
m
-
1
f
k
f
m
-
1
-
k
′′
-
∑
k
=
0
m
-
1
f
k
′
f
m
-
1
-
k
′
-
K
f
m
-
1
′
,
(20)
ℒ
2
[
θ
m
(
η
)
-
χ
m
θ
m
-
1
(
η
)
]
=
ℏ
2
H
(
η
)
R
m
θ
(
η
)
,
(21)
θ
m
(
0
)
=
0
,
θ
m
(
∞
)
=
0
,
(22)
θ
m
′
(
0
)
=
0
,
θ
m
(
∞
)
=
0
,
(23)
R
m
θ
(
η
)
=
θ
m
-
1
′′
+
Re
·
Pr
(
∑
k
=
0
m
-
1
f
k
θ
m
-
1
-
k
′
)
-
n
·
Re
·
Pr
(
∑
k
=
0
m
-
1
θ
k
f
m
-
1
-
k
′
)
,
(24)
χ
m
=
{
0
,
m
≤
1
,
1
,
m
>
1
.
The general solutions of (17)–(23) are
(25)fm(η)-χmfm-1(η)=fm*(η)+C1m+C2mη+C3mexp(-η),θm(η)-χmθm-1(η)=θm*(η)+C4m+C5mexp(-η),
where C1m to C5m are constants that can be obtained by applying the boundary condition in (18), (21), and (22).As discussed by Liao, the rule of coefficient ergodicity and the rule of solution existence play important roles in determining the auxiliary function and ensuring that the high-order deformation equations are closed and have solutions. In many cases, by means of the rule of solution expression and the rule of coefficient ergodicity, auxiliary functions can be uniquely determined. So we define the auxiliary functionH(η) in the following form:
(26)H(η)=exp(-η).


## 3.1. Consider Zeroth–Order Deformation Equations

(13)
(
1
-
P
)
ℒ
1
[
f
(
η
;
p
)
-
f
0
(
η
)
]
=
p
ℏ
1
H
(
η
)
N
1
[
f
(
η
;
p
)
]
,
f
(
0
;
p
)
=
f
w
,
f
′
(
0
;
p
)
=
1
,
f
′
(
∞
;
p
)
=
0
,
(
1
-
P
)
ℒ
2
[
θ
(
η
;
p
)
-
θ
0
(
η
)
]
=
p
ℏ
2
H
(
η
)
N
2
[
θ
(
η
;
p
)
]
,
θ
(
0
;
p
)
=
1
,
θ
(
∞
;
p
)
=
0
Power-law  temperature,
θ
′
(
0
;
p
)
=
-
1
,
θ
(
∞
;
p
)
=
0
Power-law  heat  flux,
N
1
[
f
(
η
;
p
)
,
θ
(
η
;
p
)
]
=
∂
3
f
(
η
;
p
)
∂
η
3
+
Re
(
f
(
η
;
p
)
∂
2
f
(
η
;
p
)
∂
η
2
-
(
∂
f
(
η
;
p
)
∂
η
)
2
)
-
∂
f
(
η
;
p
)
∂
η
,
N
2
[
f
(
η
;
p
)
,
θ
(
η
;
p
)
]
=
∂
2
θ
(
η
;
p
)
∂
η
2
+
Re
·
Pr
(
f
(
η
;
p
)
∂
θ
(
η
;
p
)
∂
η
)
-
n
·
Re
·
Pr
(
θ
(
η
;
p
)
∂
f
(
η
;
p
)
∂
η
)
.
For p=0 and p=1, we have
(14)f(η;0)=f0(η),f(η;1)=f(η),θ(η;0)=θ0(η),θ(η;1)=θ(η).
When p increases from 0 to 1, then f(η;p) and θ(η;p) vary from f0(η) and θ0(η) to f(η) and θ(η). By Taylor’s theorem and using (14), we can write the following:
(15)f(η;p)=f0(η)+∑m-1∞fm(η)pm,fm(η)=1m!∂m(f(η;p))∂pm,θ(η;p)=θ0(η)+∑m-1∞θm(η)pm,θm(η)=1m!∂m(θ(η;p))∂pm.
In which ℏ1 and ℏ2 are chosen in such a way that these two series are convergent at p=1; therefore we have the following through (15):
(16)f(η)=f0(η)+∑m=1∞fm(η),θ(η)=θ0(η)+∑m=1∞θm(η).


## 3.2. ConsiderMth-Order Deformation Equations

(17)
ℒ
1
[
f
m
(
η
)
-
χ
m
f
m
-
1
(
η
)
]
=
ℏ
1
H
(
η
)
R
m
f
(
η
)
,
(18)
f
m
(
0
)
=
0
,
f
m
′
(
0
)
=
0
,
f
m
′
(
∞
)
=
0
,
(19)
R
m
f
(
η
)
=
f
m
-
1
′′′
+
∑
k
=
0
m
-
1
f
k
f
m
-
1
-
k
′′
-
∑
k
=
0
m
-
1
f
k
′
f
m
-
1
-
k
′
-
K
f
m
-
1
′
,
(20)
ℒ
2
[
θ
m
(
η
)
-
χ
m
θ
m
-
1
(
η
)
]
=
ℏ
2
H
(
η
)
R
m
θ
(
η
)
,
(21)
θ
m
(
0
)
=
0
,
θ
m
(
∞
)
=
0
,
(22)
θ
m
′
(
0
)
=
0
,
θ
m
(
∞
)
=
0
,
(23)
R
m
θ
(
η
)
=
θ
m
-
1
′′
+
Re
·
Pr
(
∑
k
=
0
m
-
1
f
k
θ
m
-
1
-
k
′
)
-
n
·
Re
·
Pr
(
∑
k
=
0
m
-
1
θ
k
f
m
-
1
-
k
′
)
,
(24)
χ
m
=
{
0
,
m
≤
1
,
1
,
m
>
1
.
The general solutions of (17)–(23) are
(25)fm(η)-χmfm-1(η)=fm*(η)+C1m+C2mη+C3mexp(-η),θm(η)-χmθm-1(η)=θm*(η)+C4m+C5mexp(-η),
where C1m to C5m are constants that can be obtained by applying the boundary condition in (18), (21), and (22).As discussed by Liao, the rule of coefficient ergodicity and the rule of solution existence play important roles in determining the auxiliary function and ensuring that the high-order deformation equations are closed and have solutions. In many cases, by means of the rule of solution expression and the rule of coefficient ergodicity, auxiliary functions can be uniquely determined. So we define the auxiliary functionH(η) in the following form:
(26)H(η)=exp(-η).


## 4. Convergence of the HAM Solution

As was mentioned in introduction, the convergence and the rate of approximation of the HAM solution strongly depend on the values of auxiliary parameterℏ. By means of the so-called ℏ curves, it is easy to find out the so-called valid regions of ℏ to gain a convergent solution series. According to Figures 2 and 3, the acceptable range of auxiliary parameter for Pr=1,n=0,Re=1,andfw=0 is -0.3<ℏ1<-0.1 and -1.8<ℏ2<-0.3.Figure 2
Theℏ1 validity for Pr=1, n=0, Re=1, and fw=0.Figure 3
Theℏ2 validity for Pr=1, n=0, Re=1, and fw=0.Figures4 and 5 show how auxiliary parameters ℏ1 and ℏ2 vary with changing fw. If fw increases, the range of convergency of solution is restricted and then decreased.Figure 4
Theℏ1 validity for various fw, Pr=1, n=0, and Re=1.Figure 5
Theℏ2 validity for various fw, Pr=1, n=0, Re=1, and fw=0.


## 5. Results and Discussion

In the present study, the Homotopy analysis method has been used to solve the nonlinear equations of heat transfer and fluid flow over a permeable stretching wall in a porous medium. The nondimensional numbers introduced in the present analysis are Reynolds number (Re), Prandtl number (Pr), and the injection number (fw). Another important parameter is the power of the surface temperature/heat flux distribution (n), which is considered in the present study.In order to ensure the convergence of the solution series, the variation off has been plotted at different orders of approximation in Figure 6. The comparison between results of HAM and numerical solution (Runge-Kutta) has been displayed in Table 1. It can be observed that there is a good agreement between HAM method and numerical solution. Figures 7, 8, and 9 illustrate effect of variation of wall injection/suction parameter (fw) on velocity and temperature distribution. It can be observed that all boundary layer thickness decreased by increasing fw from negative to positive (from injection to suction). Figures 10, 11, and 12 show the effect of Reynolds number on the velocity profile when Pr=1, n=0, and fw=0. This figure shows that the boundary layer thickness and thermal boundary layer thickness are quite opposite to that of Reynolds number. The influence of n on temperature field for both types of the thermal boundary conditions considered has been studied in Figures 13 and 14. We notice that increasing n reduces the thermal boundary layer thickness regardless of the boundary condition type leading to a heat transfer augmentation. Effect of Prandtl number (isothermal) on the temperature field is plotted in Figure 15. This figure shows that thermal boundary layer thickness directly depends on Prandtl number.Table 1
The results of HAM and NS forf(η), f′(η), and θ(η) when fw=0, Pr=1, n=0, and Re=1.
η
f
(
η
)
f
′
(
η
)
θ
(
η
) (temperature)
θ
(
η
) (heat flux)
HAM
NS
HAM
NS
HAM
NS
HAM
NS
0.0
0.00000
0.00000
0.99999
1.00000
0.999999
1.00000
1.9565
1.9615
0.2
0.17042
0.17420
0.75376
0.75363
0.89633
0.89996
1.7677
1.7627
0.4
0.30564
0.30549
0.56813
0.75363
0.80205
0.80327
1.5658
1.5706
0.6
0.40571
0.40443
0.42820
0.42804
0.71301
0.71225
1.3856
1.3898
0.8
0.47505
0.47900
0.32272
0.32259
0.62792
0.62826
1.2210
1.2229
1.0
0.53802
0.53519
0.24322
0.24311
0.55226
0.55191
1.0690
1.0712
1.2
0.58010
0.57754
0.18330
0.18322
0.48134
0.48329
0.9313
0.9340
1.4
0.61161
0.60946
0.13814
0.13808
0.42056
0.42214
0.8094
0.8134
1.6
0.64096
0.63352
0.10410
0.10406
0.36595
0.36800
0.6943
0.7058
1.8
0.65374
0.65165
0.07845
0.07842
0.31817
0.32029
0.6095
0.6110
2.0
0.66952
0.66531
0.05912
0.05910
0.27657
0.27841
0.5193
0.5278
2.2
0.67812
0.67560
0.04455
0.04454
0.23934
0.24175
0.4480
0.4549
2.4
0.68567
0.68336
0.03357
0.03357
0.20813
0.20971
0.3861
0.3913
2.6
0.69054
0.68921
0.02530
0.02529
0.17972
0.18177
0.3323
0.3358
2.8
0.69675
0.69362
0.01907
0.01906
0.15519
0.15743
0.2858
0.2874
3.0
0.69842
0.69694
0.01437
0.01436
0.13407
0.13623
0.2456
0.2453Figure 6
The variation off at different orders of approximations, Pr=1, n=0, Re=1, and fw=0.Figure 7
Velocity profilef′ for various fw when Pr=1, n=0, and Re=1.Figure 8
Temperature profile (isothermal)θ for various fw when Pr=1, n=0, and Re=1.Figure 9
Temperature profile (isoflux)θ for various fw when Pr=1, n=0, and Re=1.Figure 10
Velocity profilef′ for various Re when Pr=1, n=0, and fw=0.Figure 11
Temperature profile (isothermal)θ for various Re when Pr=1, n=0, and fw=0.Figure 12
Temperature profile (isoflux)θ for various Re when Pr=1, n=0, and fw=0.Figure 13
Temperature profile (isothermal)θ for various n when Pr=1, n=0, and fw=0.Figure 14
Temperature profile (isoflux)θ for various n when Pr=1, n=0, and fw=0.Figure 15
Temperature profileθ (isothermal) for various Pr when Pr=1, n=0, and fw=0.Tables2 and 3 compare the results of HAM and numerical solution when the Reynolds number varies for fw=0, Pr=1, and n=0 for isothermal and heat flux boundary conditions, respectively. It can be seen that for isothermal boundary condition with increasing Reynolds number, the wall shear stress for both boundary conditions considered here (isothermal and heat flux) and Nusselt number increase with Reynolds number for isothermal state and independent of Reynolds number for heat flux boundary condition.Table 2
The results of HAM and NS forf′′(0), θ′(0) when fw=0, Pr=1, and n=0 for temperature.
Re
ℏ
1
ℏ
2
-
f
′
′
(
0
)
-
θ
′
(
0
)
HAM
NS
HAM
NS
1
-
0.2
-
0.4
1.4198
1.4242
0.5030
0.5033
1.5
-
0.25
-
0.4
1.5799
1.5811
0.6456
0.6422
2
-
0.2
-
0.2
1.7234
1.7320
0.7518
0.7592
5
-
0.1
-
0.1
2.4394
2.4494
1.2636
1.2576Table 3
The results of HAM and NS forf′′(0), θ′(0) when fw=0, Pr=1, and n=0 for heat flux.
Re
ℏ
1
ℏ
2
-
f
′
′
(
0
)
-
θ
′
(
0
)
HAM
NS
HAM
NS
1
-
0.25
-
0.2
1.4191
1.4242
0.9997
1
1.5
-
0.2
-
0.4
1.5791
1.5811
1.0000
1
2
-
0.15
-
0.3
1.7298
1.7320
0.9985
1
2.5
-
0.1
-
0.1
2.4445
2.4494
0.9939
1Comparison of the results of HAM and numerical solution has been shown in Tables4 and 5 for various Prandtl numbers when Re=1, n=0 and fw=-0.3,fw=0.3, respectively. It can be observed that the Prandtl number has no effect on the wall shear stress for both boundary conditions and Nusselt number for heat flux boundary condition.Table 4
The results of HAM and NS forf′′(0), θ′(0) for various Pr when, Re=1, n=0, and  fw=-0.3.
Pr
-
f
′
′
(
0
)
-
θ
′
(
0
)
f
w
=
-
0.3
f
w
=
-
0.3
Temperature
Heat flux
Temperature
Heat flux
HAM
NS
HAM
NS
HAM
NS
HAM
NS
0.5
1.2699
1.2721
1.2711
1.2721
0.216
0.219
1.000
0.999
1
1.2699
1.2721
1.2711
1.2721
0.3128
0.313
1.000
1.000
5
1.2699
1.2721
1.2711
1.2721
0.636
0.623
1.000
0.999Table 5
The results of HAM and NS forf′′(0), θ′(0) for various Pr when, Re=1, n=0, and  fw=0.3.
Pr
-
f
′
′
(
0
)
-
θ
′
(
0
)
f
w
=
0.3
f
w
=
0.3
Temperature
Heat flux
Temperature
Heat flux
HAM
NS
HAM
NS
HAM
NS
HAM
NS
0.5
1.5699
1.5721
1.5695
1.5721
0.4145
0.4100
1.0002
1.0000
1
1.5699
1.5721
1.5695
1.5721
0.7258
0.7216
1.0002
1.0000
5
1.5699
1.5721
1.5695
1.5721
2.5821
2.5826
1.0002
1.0000


## 6. Summary and Conclusions

Homotopy analysis method (HAM) is applied to compute wall driven flow through a porous medium over a stretching permeable surface subjected to power-law temperature and heat flux boundary conditions. The proper range of the auxiliary parameterℏ to ensure the convergency of the solution series was obtained through the so-called ℏ curves. When compared with other analytical methods, it is clear that HAM provides highly accurate analytical solutions for nonlinear problems. Moreover, second-law (of thermodynamics) aspects of the problem are investigated. The highlights of this study are the following. (i)
The nondimensional viscous boundary layer thickness has a direct relationship with Reynolds number; thus Nusselt number rate increases withRe.(ii)
Nusselt number, wall shear stress have a reverse relationship with and mass transfer from the wall,fw.(iii)
Increasing the Prandtl number results in reduction of thermal boundary layer thickness.Consequently, Nusselt number increase withPr.



---

*Source: 682795-2013-09-16.xml*