# A New Approach for a Class of the Blasius Problem via a Transformation and Adomian’s Method

Authors: Abdelhalim Ebaid; Nwaf Al-Armani
Journal: Abstract and Applied Analysis
 (2013)

---

## Abstract

The main feature of the boundary layer flow problems is the inclusion of the boundary conditions at infinity. Such boundary conditions cause difficulties for any of the series methods when applied to solve such problems. To the best of the authors’ knowledge, two procedures were used extensively in the past two decades to deal with the boundary conditions at infinity, either the Padé approximation or the direct numerical codes. However, an intensive work is needed to perform the calculations using the Padé technique. Regarding this point, a new idea is proposed in this paper. The idea is based on transforming the unbounded domain into a bounded one by the help of a transformation. Accordingly, the original differential equation is transformed into a singular differential equation with classical boundary conditions. The current approach is applied to solve a class of the Blasius problem and a special class of the Falkner-Skan problem via an improved version of Adomian’s method (Ebaid, 2011). In addition, the numerical results obtained by using the proposed technique are compared with the other published solutions, where good agreement has been achieved. The main characteristic of the present approach is the avoidance of the Padé approximation to deal with the infinity boundary conditions.


---

## Body

## 1. Introduction

During the past two decades much effort has been spent in the numerical treatment of boundary value problems over an unbounded domain. In fact, these problems arise very frequently in many fields such as in fluid dynamics, aerodynamics, and quantum mechanics. A few notable examples are the Blasius and Falkner-Skan equations. The Blasius equation is one of the basic equations in fluid dynamics. It describes the velocity profile of the fluid in the boundary layer theory [1, 2] on a half-infinite interval. Several analytical and numerical methods have been proposed in [1–11] to handle this problem. The two forms of the Blasius problem are represented by the same differential equation with different sets of boundary conditions, as will be indicated later. The main feature of the Blasius problem is the existence of the boundary conditions at infinity. Such conditions at infinity cause difficulties for any of the series methods, such as the Adomian decomposition method [12–14] and the differential transformation method (or the Taylor series method) [15, 16]. This is because the infinity boundary condition cannot be imposed directly in the series, where the Padé approximation should be established before applying the boundary condition at infinity. It was observed in the past two decades that many authors [17–25] have been resorted to either the Padé technique or some numerical methods to treat the boundary conditions at infinity. Although the results obtained by using the Padé technique were accurate in many cases, a massive computational work was needed to obtain accurate approximate solutions. A possible way to avoid the Padé technique is to change the boundary conditions at infinity into classical conditions. Therefore, a suggestion is proposed in this paper to transform the domain of the problem from an unbounded domain into a bounded one with the help of a simple transformation.According to the suggested transformation, the original Blasius equation is transformed into a system of two singular differential equations. Hence, the two mentioned forms are described by this system with two different sets of boundary conditions at classical point. The transformed singular system will be solved by a recent version of the ADM [26]. The first form of the original Blasius problem is given by [6]
(1)f′′′(η)+12f(η)f′′(η)=0,
subject to the following boundary conditions:
(2)f(0)=0,f′(0)=1,f′(∞)=0,
while the second form is given by (1) with the following boundary conditions:
(3)f(0)=0,f′(0)=0,f′(∞)=1.
A class of Blasius problem is given by
(4)f′′′(η)+γf(η)f′′(η)=0,
subject to the following boundary conditions:
(5)f(0)=0,f′(0)=1-ϵ,f′(∞)=ϵ,
where γ and ϵ are finite constants. This class will be studied for ϵ∈[0,1]. Here, it is noted that (2) and (3) are special examples of (5) for ϵ=0 and ϵ=1, respectively. In addition, the class (4)-(5) reduces to the two forms of the Blasius problem when (γ=1/2, ϵ=0) and (γ=1/2, ϵ=1), respectively. At the same time, when γ=1, the suggested class reduces to a special class of the Falkner-Skan problem, at δ=0, which is well known as [27]
(6)f′′′(η)+f(η)f′′(η)+δ[ϵ2-(f′(η))2]=0,
with the class of boundary conditions (5), where δ refers to the pressure gradient parameter, while ϵ refers to the velocity ratio parameter, ϵ=U∞/(U∞+Uw). Equation (6) with the boundary conditions (5) is a new version of the Falkner-Skan equation relating free stream velocity U∞ to composite reference velocity, that is, sum of the velocities of stretching boundary Uw and free stream U∞. In order to use the improved Adomian’s method [26] to solve the class (4)-(5), we first transform the governing equation (4) into the following system of differential equations:
(7)f′(η)=u(η),u′′(η)+γf(η)u′(η)=0.
Here, we may indicate that in the theory of the boundary layer, it is usually important to get information about three quantities: the skin-friction coefficient f′′(0), the fluid velocity f′(η), and the stream function f(η). Also, it is well known that at ϵ=1 the problem reduces to one of the two forms of the Blasius problem which has been studied extensively during the past decades.


## 2. A Transformation and a New System

The unbounded domain of the independent variableη∈[0,∞) can be changed into a bounded one by using a new independent variable t (say) ∈[0,1) using the transformation t=1-e-η. Accordingly, the governing system should be expressed in terms of the new variable t. In order to do that, we introduce the following relations between the derivatives with respect to η and the derivatives with respect to t:
(8)ddη(□)=(1-t)ddt(□),d2dη2(□)=(1-t)2d2dt2(□)-(1-t)ddt(□).
The relations given by (8) are obtained by using the chain rule in the differential calculus. Therefore, the system (7) becomes
(9)f′(t)=(11-t)u(t),(10)u′′(t)=(11-t)u′(t)-γ(11-t)f(t)u′(t),
subject to the following set of boundary conditions:
(11)f(0)=0,u(0)=1-ϵ,u(1)=ϵ.Equation (9) with the initial condition f(0)=0 can be easily integrated as an initial value problem, while (10) with the boundary conditions given in (11) should be solved as a two-point boundary value problem. In this regard, the improved Adomian decomposition method is suggested to deal with such a singular two-point boundary value problem. Before launching into the the main idea of this paper, we give an analysis for the improved Adomian decomposition method in the next section to solving (10) with general two-point boundary conditions u(a)=α and u(b)=β.


## 3. The Improved Adomian Decomposition Method

Consider the second order differential equation:(12)u′′(t)+p(t)u′(t)+q(t)f(t)u′(t)=0,
subject to the boundary conditions
(13)u(a)=α,u(b)=β,
where at least one of the functions p(t) and q(t) has a singular point and f(t) is an unspecified function. In order to apply the approach suggested in [26], we first rewrite (12) as
(14)u′′(t)=-p(t)u′(t)-q(t)f(t)u′(t).
Now, suppose that p(t) and q(t) have the same singular point (t=t0, say), Ebaid [26] proposed the following inverse operator to solve (14) with the boundary conditions (13):
(15)L-1[·]=∫atdt′∫ct′[·]dt′′-t-ab-a∫abdt′∫ct′[·]dt′′,a≠b,c(arbitrary)≠t0.
Operating both sides of (14) with this inverse operator, we have
(16)u(t)-u(a)-t-ab-a[u(b)-u(a)]=-L-1[p(t)u′(t)+q(t)f(t)u′(t)],
which can be rewritten as
(17)u(t)=α+t-ab-a(β-α)-L-1[p(t)u′(t)+q(t)f(t)u′(t)].
Based on Adomian’s method, the solutions u(t) and f(t) of system (9)-(10) are assumed in the following form:
(18)u(t)=∑n=0∞un(t),f(t)=∑n=0∞fn(t).
Inserting these series into (17), we obtain
(19)∑n=0∞un(t)=α+t-ab-a(β-α)-L-1[p(t)∑n=0∞un′(t)+q(t)∑n=0∞∑i=0nfi(t)un-i′(t)].
Substituting p(t)=-(1/(1-t)) and q(t)=γ(1/(1-t)) into the last equation yields
(20)∑n=0∞un(t)=α+t-ab-a(β-α)+L-1[-γ(11-t)∑n=0∞∑i=0nfi(t)un-i′(t)(11-t)∑n=0∞un′(t)-γ(11-t)∑n=0∞∑i=0nfi(t)un-i′(t)].
To overcome the difficulty of the singular point, we may replace the function 1/(1-t) with the series form ∑n=0∞tn, where t∈[0,1). Thus, we have
(21)∑n=0∞un(t)=α+t-ab-a(β-α)+L-1[∑n=0∞(∑j=0n∑i=0jtn-jfi(t)uj-i′(t))∑n=0∞∑i=0ntn-iui′(t)-γ∑n=0∞(∑j=0n∑i=0jtn-jfi(t)uj-i′(t))].
According to the modified decomposition method [18], the solution u(t) can be evaluated by using the recurrence scheme:
(22)u0(t)=α,u1(t)=t-ab-a(β-α)+L-1[u0′(t)-γf0(t)u0′(t)],un+1(t)=L-1[∑i=0ntn-iui′(t)-γ(∑j=0n∑i=0jtn-jfi(t)uj-i′(t))],n≥1.
On integrating (9) with respect to t form 0 to t, it then follows that
(23)f(t)=f(0)+∫0t(11-z)u(z)dz.
Hence, f(t) is given by the recurrence scheme:
(24)f0(t)=0,fn+1(t)=∫0tzn-iui(z)dz,n≥0.
The algorithms given by (22) and (24) are applied in the next section to construct the approximate solutions.


## 4. Applications

### 4.1. A Class of the Blasius Problem

Here, we show how to implement (22) and (24) to solve the class of the Blasius problem. On substituting γ=1/2, a=0, b=1, α=1-ϵ, and β=ϵ into (22), and using (24) we obtain
(25)u0(t)=1-ϵ,f0(t)=0,u1(t)=(2ϵ-1)t,f1(t)=(1-ϵ)t,un+1(t)=L-1[∑i=0ntn-iui′(t)-12(∑j=0n∑i=0jtn-jfi(t)uj-i′(t))],fn+1(t)=∫0t(11-z)un(z)dz,n≥1,L-1[·]=∫0tdt′∫ct′[·]dt′′-t∫01dt′∫ct′[·]dt′′,c≠1(c=0,forsimplicity).
The first few terms of the stream function f(t) are evaluated by implementing the previous algorithm and are listed in the following:
(26)f0(t)=0,f1(t)=(1-ϵ)t,f2(t)=(ϵ2)t2,f3(t)=(14-ϵ2)t2+(-16+2ϵ3)t3,f4(t)=(ϵ24-ϵ212)t2+(14-ϵ2)t3+(-316+29ϵ28+ϵ224)t4,f5(t)=(196-ϵ48)t2+(ϵ24-ϵ212)t3+(732-41ϵ96-ϵ248)t4+(-43240+41ϵ40+ϵ215)t5.
The desired mth order approximate solution ϕm(η) obtained by Adomian's method is expressed as
(27)ϕm(η)=∑n=0m-1fn(η).
Hence, the approximate solutions ϕ3(η), ϕ5(η), and ϕ7(η) are, respectively, given in terms of the original variable η as
(28)ϕ3(η)=(1-ϵ)(1-e-η)+14(1-e-η)2+(-16+2ϵ3)(1-e-η)3,ϕ5(η)=(1-ϵ)(1-e-η)+(2596+ϵ48-ϵ212)(1-e-η)2+(112+5ϵ24-ϵ212)(1-e-η)3+(132+17ϵ96+ϵ248)(1-e-η)4+(-43240+41ϵ40+ϵ215)(1-e-η)5,ϕ7(η)=(1-ϵ)(1-e-η)+(171640+ϵ72-553ϵ25760-ϵ3960)(1-e-η)2+(47480+11ϵ60-53ϵ2576+ϵ31440)(1-e-η)3+(31768+451ϵ2304-59ϵ21152-ϵ3288)(1-e-η)4+(11960+137ϵ720-7ϵ2288-ϵ3240)(1-e-η)5+(315760+367ϵ2880+65ϵ21152+ϵ3960)(1-e-η)6+(-138960+5921ϵ40320-83ϵ224192+89ϵ334560-19ϵ4120960)×(1-e-η)7.
Here, we refer to that the series solution obtained previous leads to an exact solution at ϵ=0.5. In this case, the approximate solutions become
(29)ϕ3(η)=12(1-e-η)+14(1-e-η)2+16(1-e-η)3=12∑r=131r(1-e-η)r,ϕ5(η)=12(1-e-η)+14(1-e-η)2+16(1-e-η)3+18(1-e-η)4+110(1-e-η)5=12∑r=151r(1-e-η)r,ϕ7(η)=12(1-e-η)+14(1-e-η)2+16(1-e-η)3+18(1-e-η)4+110(1-e-η)5+112(1-e-η)6+114(1-e-η)7=12∑r=171r(1-e-η)r.
Therefore, the m-term series solution is given by
(30)ϕm(η)=12∑r=1m1r(1-e-η)r,
and thus, the following exact solution is obtained as m→∞:
(31)f(η)=limm→∞ϕm(η)=12∑r=1∞1r(1-e-η)r=-12ln[1-(1-e-η)]=η2.
This exact solution satisfies the boundary conditions and can be easily verified by direct substitution. For more validation, the results obtained by the present technique are checked here via a comparison with those published in the literature. It is well known that at ϵ=1, the problem reduces to one of the two forms of the Blasius problem. In that case, the skin-friction coefficient is computed by many authors as discussed in Section 5.


### 4.2. Special Class of the Falkner-Skan Problem

Here, the proposed approach is applied to a special class of the Falkner-Skan problem. As mentioned before, this special class is given by (4)-(5) at γ=1. Proceeding as in the previous example, the approximate solution can be obtained by using the recurrence scheme:
(32)u0(t)=1-ϵ,f0(t)=0,u1(t)=(2ϵ-1)t,f1(t)=(1-ϵ)t,un+1(t)=L-1[∑i=0ntn-iui′(t)-∑j=0n∑i=0jtn-jfi(t)uj-i′(t)],fn+1(t)=∫0t(11-z)un(z)dz,n≥1,L-1[·]=∫0tdt′∫ct′[·]dt′′-t∫01dt′∫ct′[·]dt′′,c≠1.
The first few terms of the stream function f(t) are evaluated by implementing the algorithm given in (32) and are listed in the following:
(33)f0(t)=0,f1(t)=(1-ϵ)t,f2(t)=(ϵ2)t2,f3(t)=(14-ϵ2)t2+(-16+2ϵ3)t3,f4(t)=(-124+ϵ6-ϵ26)t2+(14-ϵ2)t3+(-16+13ϵ24+ϵ212)t4,f5(t)=(-124+ϵ6-ϵ26)t3+(524-19ϵ48-ϵ224)t4+(-17120+5ϵ12+2ϵ215)t5.
Hence, the approximate solutions ϕ3(η), ϕ5(η), and ϕ7(η) are, respectively, given in terms of the original variable η as
(34)ϕ3(η)=(1-ϵ)(1-e-η)+14(1-e-η)2+(-16+2ϵ3)(1-e-η)3,ϕ5(η)=(1-ϵ)(1-e-η)+(524+ϵ6-ϵ26)(1-e-η)2+(124+ϵ3-ϵ26)(1-e-η)3+(124+7ϵ48+ϵ224)(1-e-η)4+(-17120+5ϵ12+2ϵ215)(1-e-η)5,ϕ7(η)=(1-ϵ)(1-e-η)+(191960+19ϵ96-91ϵ2480-ϵ3240)(1-e-η)2+(13360+17ϵ48-3ϵ216+ϵ3360)(1-e-η)3+(1144+9ϵ32-ϵ212-ϵ372)(1-e-η)4+(1480+17ϵ80-ϵ240-ϵ360)(1-e-η)5+(130+13ϵ288+31ϵ2288+ϵ3240)(1-e-η)6+(-67630+25ϵ84+181ϵ21680+41ϵ32520)(1-e-η)7.
The effectiveness of the present technique is used here not only to obtain the exact solution of the Falkner-Skan equation at δ=0 and ϵ=0.5 but also to get numerical solutions with good accuracy. On inserting ϵ=0.5 into the approximate solutions given by (34), we have
(35)ϕ3(η)=12(1-e-η)+14(1-e-η)2+16(1-e-η)3=12∑r=131r(1-e-η)r,ϕ5(η)=12(1-e-η)+14(1-e-η)2+16(1-e-η)3+18(1-e-η)4+110(1-e-η)5=12∑r=151r(1-e-η)r,ϕ7(η)=12(1-e-η)+14(1-e-η)2+16(1-e-η)3+18(1-e-η)4+110(1-e-η)5+112(1-e-η)6+114(1-e-η)7=12∑r=171r(1-e-η)r.
As indicated in Section 4.1, these approximate solutions lead to the same exact solution given by (31): f(η)=η/2 in the limit.


## 4.1. A Class of the Blasius Problem

Here, we show how to implement (22) and (24) to solve the class of the Blasius problem. On substituting γ=1/2, a=0, b=1, α=1-ϵ, and β=ϵ into (22), and using (24) we obtain
(25)u0(t)=1-ϵ,f0(t)=0,u1(t)=(2ϵ-1)t,f1(t)=(1-ϵ)t,un+1(t)=L-1[∑i=0ntn-iui′(t)-12(∑j=0n∑i=0jtn-jfi(t)uj-i′(t))],fn+1(t)=∫0t(11-z)un(z)dz,n≥1,L-1[·]=∫0tdt′∫ct′[·]dt′′-t∫01dt′∫ct′[·]dt′′,c≠1(c=0,forsimplicity).
The first few terms of the stream function f(t) are evaluated by implementing the previous algorithm and are listed in the following:
(26)f0(t)=0,f1(t)=(1-ϵ)t,f2(t)=(ϵ2)t2,f3(t)=(14-ϵ2)t2+(-16+2ϵ3)t3,f4(t)=(ϵ24-ϵ212)t2+(14-ϵ2)t3+(-316+29ϵ28+ϵ224)t4,f5(t)=(196-ϵ48)t2+(ϵ24-ϵ212)t3+(732-41ϵ96-ϵ248)t4+(-43240+41ϵ40+ϵ215)t5.
The desired mth order approximate solution ϕm(η) obtained by Adomian's method is expressed as
(27)ϕm(η)=∑n=0m-1fn(η).
Hence, the approximate solutions ϕ3(η), ϕ5(η), and ϕ7(η) are, respectively, given in terms of the original variable η as
(28)ϕ3(η)=(1-ϵ)(1-e-η)+14(1-e-η)2+(-16+2ϵ3)(1-e-η)3,ϕ5(η)=(1-ϵ)(1-e-η)+(2596+ϵ48-ϵ212)(1-e-η)2+(112+5ϵ24-ϵ212)(1-e-η)3+(132+17ϵ96+ϵ248)(1-e-η)4+(-43240+41ϵ40+ϵ215)(1-e-η)5,ϕ7(η)=(1-ϵ)(1-e-η)+(171640+ϵ72-553ϵ25760-ϵ3960)(1-e-η)2+(47480+11ϵ60-53ϵ2576+ϵ31440)(1-e-η)3+(31768+451ϵ2304-59ϵ21152-ϵ3288)(1-e-η)4+(11960+137ϵ720-7ϵ2288-ϵ3240)(1-e-η)5+(315760+367ϵ2880+65ϵ21152+ϵ3960)(1-e-η)6+(-138960+5921ϵ40320-83ϵ224192+89ϵ334560-19ϵ4120960)×(1-e-η)7.
Here, we refer to that the series solution obtained previous leads to an exact solution at ϵ=0.5. In this case, the approximate solutions become
(29)ϕ3(η)=12(1-e-η)+14(1-e-η)2+16(1-e-η)3=12∑r=131r(1-e-η)r,ϕ5(η)=12(1-e-η)+14(1-e-η)2+16(1-e-η)3+18(1-e-η)4+110(1-e-η)5=12∑r=151r(1-e-η)r,ϕ7(η)=12(1-e-η)+14(1-e-η)2+16(1-e-η)3+18(1-e-η)4+110(1-e-η)5+112(1-e-η)6+114(1-e-η)7=12∑r=171r(1-e-η)r.
Therefore, the m-term series solution is given by
(30)ϕm(η)=12∑r=1m1r(1-e-η)r,
and thus, the following exact solution is obtained as m→∞:
(31)f(η)=limm→∞ϕm(η)=12∑r=1∞1r(1-e-η)r=-12ln[1-(1-e-η)]=η2.
This exact solution satisfies the boundary conditions and can be easily verified by direct substitution. For more validation, the results obtained by the present technique are checked here via a comparison with those published in the literature. It is well known that at ϵ=1, the problem reduces to one of the two forms of the Blasius problem. In that case, the skin-friction coefficient is computed by many authors as discussed in Section 5.


## 4.2. Special Class of the Falkner-Skan Problem

Here, the proposed approach is applied to a special class of the Falkner-Skan problem. As mentioned before, this special class is given by (4)-(5) at γ=1. Proceeding as in the previous example, the approximate solution can be obtained by using the recurrence scheme:
(32)u0(t)=1-ϵ,f0(t)=0,u1(t)=(2ϵ-1)t,f1(t)=(1-ϵ)t,un+1(t)=L-1[∑i=0ntn-iui′(t)-∑j=0n∑i=0jtn-jfi(t)uj-i′(t)],fn+1(t)=∫0t(11-z)un(z)dz,n≥1,L-1[·]=∫0tdt′∫ct′[·]dt′′-t∫01dt′∫ct′[·]dt′′,c≠1.
The first few terms of the stream function f(t) are evaluated by implementing the algorithm given in (32) and are listed in the following:
(33)f0(t)=0,f1(t)=(1-ϵ)t,f2(t)=(ϵ2)t2,f3(t)=(14-ϵ2)t2+(-16+2ϵ3)t3,f4(t)=(-124+ϵ6-ϵ26)t2+(14-ϵ2)t3+(-16+13ϵ24+ϵ212)t4,f5(t)=(-124+ϵ6-ϵ26)t3+(524-19ϵ48-ϵ224)t4+(-17120+5ϵ12+2ϵ215)t5.
Hence, the approximate solutions ϕ3(η), ϕ5(η), and ϕ7(η) are, respectively, given in terms of the original variable η as
(34)ϕ3(η)=(1-ϵ)(1-e-η)+14(1-e-η)2+(-16+2ϵ3)(1-e-η)3,ϕ5(η)=(1-ϵ)(1-e-η)+(524+ϵ6-ϵ26)(1-e-η)2+(124+ϵ3-ϵ26)(1-e-η)3+(124+7ϵ48+ϵ224)(1-e-η)4+(-17120+5ϵ12+2ϵ215)(1-e-η)5,ϕ7(η)=(1-ϵ)(1-e-η)+(191960+19ϵ96-91ϵ2480-ϵ3240)(1-e-η)2+(13360+17ϵ48-3ϵ216+ϵ3360)(1-e-η)3+(1144+9ϵ32-ϵ212-ϵ372)(1-e-η)4+(1480+17ϵ80-ϵ240-ϵ360)(1-e-η)5+(130+13ϵ288+31ϵ2288+ϵ3240)(1-e-η)6+(-67630+25ϵ84+181ϵ21680+41ϵ32520)(1-e-η)7.
The effectiveness of the present technique is used here not only to obtain the exact solution of the Falkner-Skan equation at δ=0 and ϵ=0.5 but also to get numerical solutions with good accuracy. On inserting ϵ=0.5 into the approximate solutions given by (34), we have
(35)ϕ3(η)=12(1-e-η)+14(1-e-η)2+16(1-e-η)3=12∑r=131r(1-e-η)r,ϕ5(η)=12(1-e-η)+14(1-e-η)2+16(1-e-η)3+18(1-e-η)4+110(1-e-η)5=12∑r=151r(1-e-η)r,ϕ7(η)=12(1-e-η)+14(1-e-η)2+16(1-e-η)3+18(1-e-η)4+110(1-e-η)5+112(1-e-η)6+114(1-e-η)7=12∑r=171r(1-e-η)r.
As indicated in Section 4.1, these approximate solutions lead to the same exact solution given by (31): f(η)=η/2 in the limit.


## 5. Results and Discussion

Atϵ=1, Bairstow [29] found that f′′(0)=0.335 using a power series, whereas Goldstein [30] obtained f′′(0)=0.332. Besides, using a finite difference method, Falkner [31] computed that f′′(0)=0.3325765, and Horwarth [32] yields that f′′(0)=0.332057. In [33], Fazio computed that f′′(0)=0.332057336215. Also, in [34] Boyd used Töpfer’s algorithm to obtain the accurate value f′′(0)=0.33205733621519630. Adomain’s method was implemented in [35] by Abbasbandy, and it was found that f′′(0)=0.333329, whereas a variational iteration method with the Padé approximants allows Wazwaz [6] to calculate the value f′′(0)=0.3732905625. Tajvidi et al. [28] apply the modified rational Legendre functions to get a value of f′′(0)=0.33209. The values of the skin-friction coefficient are given in Table 1 at ϵ=0,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9, and 1 using 11, 13, and 15 terms of the series (27). The current method finds that the skin-friction at ϵ=1 approximately equals ϕ15′′(0)=0.331775, which is very close to those values discussed previously.Table 1
The approximate values of the skin-friction coefficientf′′(0) for the class of the Blasius problem using 11, 13, and 15 terms of Adomian's series.
ϵ
ϕ
11
′
′
(
0
)
ϕ
13
′
′
(
0
)
ϕ
15
′
′
(
0
)
0.0
-
0.456523
-
0.454506
-
0.453122
0.1
-
0.356731
-
0.354995
-
0.353838
0.2
-
0.261082
-
0.259650
-
0.258715
0.3
-
0.169677
-
0.168611
-
0.167920
0.4
-
0.082616
-
0.082014
-
0.081624
0.5
0.000000
0.000000
0.000000
0.6
0.078072
0.077292
0.076775
0.7
0.151503
0.149723
0.148522
0.8
0.220194
0.217151
0.215060
0.9
0.284048
0.279437
0.276206
1.0
0.342969
0.336441
0.331775Regarding the stream functionf(η), it is plotted in Figure 1 using 15 terms, and the fluid velocity f′(η) is depicted in Figure 2 using the approximate solutions ϕ7(η), ϕ9(η), ϕ11(η), ϕ13(η), and ϕ15(η) at ϵ=0,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9, and 1. It is observed from Figure 2 that the approximate solutions using a few terms of Adomian's series converge rapidly to a certain curve at some values of the parameter ϵ.Figure 1
The stream function for the class of the Blasius problem at different values ofϵ using 15 terms of the current method.Figure 2
The fluid velocity for the class of the Blasius problem at different values ofϵ using 7, 9, 11, 13, and 15 terms of the current method.The exact solutionf(η)=η/2 obtained in Section 4.2 for the Falkner-Skan equation at ϵ=0.5 has been reported very recently by Kudenatti [27]. He has derived the exact solution to the Falkner-Skan equation for general values of the pressure gradient parameter δ. In order to check the accuracy of our approach, the values of the skin-friction coefficient are compared in Table 2 with those exactly obtained by Kudenatti [27] in the range 0<ϵ<0.5. The results reveal that a good agreement has been achieved via the present approach. In addition, the stream function f(η) is graphed in Figure 3 at several values of the parameter ϵ by using 15 terms of the decomposition series. At the same values and in Figure 4, the fluid velocity is depicted by using the approximate solutions ϕ11(η), ϕ13(η), and ϕ15(η). It can be concluded from Figure 4 that our results are a coincidence with those exactly obtained in [27] at the values ϵ=0.1,0.2,0.3,0.4,0.5,0.7, and 1, while the fluid velocity at the other values ϵ=0,0.6,0.8, and 0.9 was not discussed by Kudenatti [27].Table 2
The approximate values of the skin-friction coefficientf′′(0) for the class of the Falkner-Skan problem using 11, 15, and 33 terms of Adomian's series.
ϵ
ϕ
11
′
′
(
0
)
ϕ
15
′
′
(
0
)
ϕ
33
′
′
(
0
)
Exact valuesReference  [28]
0.0
-
0.617661
-
0.622494
-
0.625945
-
0.627504
0.1
-
0.479542
-
0.485036
-
0.490729
-
0.492625
0.2
-
0.348276
-
0.353437
-
0.360126
-
0.363901
0.3
-
0.224304
-
0.228296
-
0.234485
-
0.237219
0.4
-
0.108067
-
0.110256
-
0.114249
-
0.115811
0.5
0.000000
0.000000
0.000000
0.000000
0.6
0.099467
0.101753
0.107483
0.7
0.189909
0.194253
0.207152
0.8
0.270907
0.276719
0.297613
0.9
0.342046
0.34835
0.377038
1.0
0.402921
0.408321
0.443088Figure 3
The stream function for the class of the Falkner-Skan problem at different values ofϵ using 15 terms of the current method.Figure 4
The fluid velocity for the class of the Falkner-Skan problem at different values ofϵ using 11, 13, and 15 terms of the current method.


## 6. Conclusion

An approach is presented in this paper to treat the boundary condition at infinity which is the main feature of the boundary layer equations. The suggested approach is based on changing the boundary condition at infinity to a classical one by the help of a transformation. The current approach is applied to solve a class of the Blasius problem and a special class of the Falkner-Skan problem via an improved version of Adomian's method. Moreover, exact solutions are deduced at a certain value of the velocity ratio parameterϵ. In addition, the current numerical results are compared with the other published solutions, where good agreement is found. One of the main advantages of the present approach is the avoidance of the Padé approximation to deal with the infinity boundary condition.



---

*Source: 753049-2013-07-31.xml*