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.
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. Applications4.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.
The approximate values of the skin-friction coefficient f′′(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.331775
Regarding the stream function f(η), 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 ϵ.
The stream function for the class of the Blasius problem at different values of ϵ using 15 terms of the current method.
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 solution f(η)=η/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].
The approximate values of the skin-friction coefficient f′′(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.443088
The stream function for the class of the Falkner-Skan problem at different values of ϵ using 15 terms of the current method.
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.
BelhachmiZ.BrighiB.TaousK.On the concave solutions of the Blasius equation2000692199214MR1819521ZBL0972.34015DattaB. K.Analytic solution for the Blasius equation2003342237240MR1964524ZBL1054.34011KuikenH. K.On boundary layers in fluid mechanics that decay algebraically along stretches of wall that are not vanishingly small198127438740510.1093/imamat/27.4.387MR637502ZBL0472.76045KuikenH. K.A “backward” free-convective boundary layer198134339741310.1093/qjmam/34.3.397MR627276ZBL0483.76046HeJ.Approximate analytical solution of Blasius' equation1998342602632-s2.0-1684437883510.1016/S1007-5704(98)90046-6WazwazA.-M.The variational iteration method for solving two forms of Blasius equation on a half-infinite domain2007188148549110.1016/j.amc.2006.10.009MR2327376ZBL1114.76055YunB. I.Intuitive approach to the approximate analytical solution for the Blasius problem2010215103489349410.1016/j.amc.2009.09.050MR2578930ZBL1183.65098AbbasbandyS.BervillierC.Analytic continuation of Taylor series and the boundary value problems of some nonlinear ordinary differential equations201121852178219910.1016/j.amc.2011.07.035MR2831493ZBL1234.34009KhanM.GondalM. A.Homotopy perturbation Padé transform method for Blasius flow equation using He’s polynomials2011121-817KhanM.GondalM. A.HussainI.VananiS. K.A new comparative study between homotopy analysis transform method and homotopy perturbation transform method on a semi infinite domain2012553-41143115010.1016/j.mcm.2011.09.038MR2887449ZBL1255.65186FazioR.Blasius problem and Falkner-Skan model: Töpfer's algorithm and its extension201373202209AdomianG.1994Boston, Mass, USAKluwer Academic Publishersxiv+352MR1282283WazwazA.-M.RachR.DuanJ.-S.Adomian decomposition method for solving the Volterra integral form of the Lane-Emden equations with initial values and boundary conditions2013219105004501910.1016/j.amc.2012.11.012MR3009461WazwazA. M.RachR.DuanJ. S.A study on the systems of the Volterra integral forms of the Lane-Emden equations by the Adomian decomposition method201310.1002/mma.2776EbaidA. E.-H.Approximate periodic solutions for the non-linear relativistic harmonic oscillator via differential transformation method20101571921192710.1016/j.cnsns.2009.07.003MR2585058ZBL1222.65084EbaidA. E.A reliable aftertreatment for improving the differential transformation method and its application to nonlinear oscillators with fractional nonlinearities201116152853610.1016/j.cnsns.2010.03.012MR2679201ZBL1221.34009VenkataranganS. N.RajalakshmiK.Modification of Adomian's decomposition method to solve equations containing radicals1995296758010.1016/0898-1221(95)00009-NMR1317861ZBL0818.34007WazwazA.-M.The modified decomposition method and Padé approximants for solving the Thomas-Fermi equation19991051111910.1016/S0096-3003(98)10090-5MR1706059ZBL0956.65064WazwazA.-M.The modified decomposition method and Padé approximants for a boundary layer equation in unbounded domain2006177273774410.1016/j.amc.2005.09.102MR2291999ZBL1096.65072KechilS. A.HashimI.Non-perturbative solution of free-convective boundary-layer equation by Adomian decomposition method20073631-21101142-s2.0-3384664078710.1016/j.physleta.2006.11.054Awang KechilS.HashimI.Series solution of flow over nonlinearly stretching sheet with chemical reaction and magnetic field200837213225822632-s2.0-3984909138310.1016/j.physleta.2007.11.027AlizadehE.SedighiK.FarhadiM.Ebrahimi-KebriaH. R.Analytical approximate solution of the cooling problem by Adomian decomposition method20091424624722-s2.0-4774911873810.1016/j.cnsns.2007.09.008HayatT.HussainQ.JavedT.The modified decomposition method and Padé approximants for the MHD flow over a non-linear stretching sheet200910296697310.1016/j.nonrwa.2007.11.020MR2474272ZBL1167.76385RashidiM. M.Mohimanian PourS. A.Explicit solution of axisymmetric stagnation flow towards a shrinking sheet by DTM-Padé2010453–5626172632MR2729560RashidiM. M.KeimaneshM.Using differential transform method and padé approximant for solving mhd flow in a laminar liquid film from a horizontal stretching surface20102010142-s2.0-7795251588010.1155/2010/491319491319EbaidA.A new analytical and numerical treatment for singular two-point boundary value problems via the Adomian decomposition method201123581914192410.1016/j.cam.2010.09.007MR2763114ZBL1209.65077KudenattiR. B.A new exact solution for boundary layer flow over a stretching plate201247727733TajvidiT.RazzaghiM.DehghanM.Modified rational Legendre approach to laminar viscous flow over a semi-infinite flat plate2008351596610.1016/j.chaos.2006.05.031MR2355238ZBL1135.76040BairstowL.Skin friction192529323GoldsteinS.Concerning some solutions of the boundary layer equations in hydrodynamics193026130FalknerV. M.A method of numerical solution of differential equations193621624640HorwarthL.On the solution of the laminar boundary layer equations1938164547579FazioR.The Blasius problem formulated as a free boundary value problem1992951–41710.1007/BF01170800MR1179914ZBL0753.76051BoydJ. P.The Blasius function in the complex plane199984381394MR173723310.1080/10586458.1999.10504626ZBL0980.34053AbbasbandyS.A numerical solution of Blasius equation by Adomian's decomposition method and comparison with homotopy perturbation method20073112572602-s2.0-3374585975510.1016/j.chaos.2005.10.071