The present paper is devoted to the convergence control and accelerating the traditional Decomposition Method of Adomian (ADM). By means of perturbing the initial or early terms of the Adomian iterates by adding a parameterized term, containing an embedded parameter, new modified ADM is constructed. The optimal value of this parameter is later determined via squared residual minimizing the error. The failure of the classical ADM is also prevented by a suitable value of the embedded parameter, particularly beneficial for the Duan–Rach modification of the ADM incorporating all the boundaries into the formulation. With the presented squared residual error analysis, there is no need to check out the results against the numerical ones, as usually has to be done in the traditional ADM studies to convince the readers that the results are indeed converged to the realistic solutions. Physical examples selected from the recent application of ADM demonstrate the validity, accuracy and power of the presented novel approach in this paper. Hence, the highly nonlinear equations arising from engineering applications can be safely treated by the outlined method for which the classical ADM may fail or be slow to converge.
Researchers prefer an easily accessible and user friendly method requiring less computational labor while accurately approximating highly nonlinear equations resulting from mathematical modeling of real-life phenomena. The Adomian decomposition method (ADM) is one such popular technique capable of dealing with the prevailing nonlinearities by means of Adomian polynomials [1,2]. A modification of the classical ADM is proposed within the current study based on the recent publications [3,4] successfully generating fast convergent ADM series solutions with as small Adomian polynomials as possible in the solution series.
A quick literature survey exhibits that ADM has been applied to many nonlinear equations [5]. To classify some of the recent bibliography, algebraic equations were contained within the references [6,7]. The ordinary differential equations were dealt within the citations [8–14]. The articles [15–18] covered the efforts to partial differential equations. Mathematical analysis of the convergence of ADM to certain nonlinear equations was fulfilled in the publications [19–21]. It can be successfully used to gain correct physical parameters domain [22]. A traffic model was also very recently treated in [23] via the Adomian method. The publications by [24,25] present investigation of some nonlinear problems via different numerical approaches.
It is now well-known that an inadequate arrangement of the classical ADM series may lead to non-convergent solutions or solutions with a poor convergence rate. To avoid these shortcomings, a parameter is generally inserted at the leading term of the Adomian series and later it is subtracted at the first order term not to break down the equation structure. This procedure was pursued by the recent publications [11,17,18]. However, how a proper value of the inserted parameter will be determined was not mentioned in these references. Instead, a randomly chosen value was assigned to it. A variety of modifications were also offered in the articles [26–31]. A successful formulation of the ADM was made in the recent work of [3] which was named as the optimal ADM. Further applications of the homotopy analytic approximate method may be found in the literature [32–35].
The motivation of the current work is, benefiting from the idea in [11], to devise a method that greatly improves the mathematical property of classical ADM. Within this aim, a reorganization of the ADM series is proposed by altering the early terms so that they incorporate extra controllable terms. The reason of such a treatment is to get a rapidly converging ADM solutions with the least Adomian polynomials. In place of randomly selecting, an optimum value of the introduced parameter is later determined through error on the grounds of total residual. With this value at our disposal, there is no doubt that the ADM method is convergent to the true solution in a most rapid way, not demanding a verification of the ADM solutions by numerical ones. The failure of the classical ADM in the usual form or in the Duan–Rach formulation is also prevented by a suitable value of the embedded parameter. The present approach can also extend the region of convergence of the traditional method. Examples of physical value are provided to justify and validate the given procedure.
Traditional Decomposition Method of Adomian
The usual steps of traditional ADM can be inferred from the aforementioned citations. The methodology in brief is such that under an invertible linear operator L and a forcing function f, it is desired to approximate the function u having the nonlinearity N(u) and satisfying the general nonlinear equation
L(u)=N(u)+f,
with the initial and/or boundary restrictions
B(u)=b.
Having inverted (1) under the restrictions (2) generally leads to
u=g+L-1[N(u)]+L-1[f],
where g is due to the conditions in (2). If u is a single scalar parameter like (1) representing an algebraic equation, then there is no such g in (3), whereas, in the case of a variable u, L−1 denotes an integral operator giving rise to g in (3). Then (3) is a mixed Volterra-Fredholm type equation so-called as the Duan–Rach formulation in the recent literature, see for instance [13,14].
The subsequent series decompositions of u and N(u)
u=∑n=0∞un,N(u)=∑n=0∞An,
in which An’s are the classical Adomian polynomials, are later substituted into (3). The solution u of (1) is finally generated from the recurrence relation
u0=g+L-1[f],un+1=L-1[An],n≥0.
As a result, by means of the relations from (5), an approximate series solution of order M is obtained as
u=∑n=0Mun,
which serves for practical purposes.
In general, the procedure in (5) yields convergent ADM series solutions, see for instance [19–21]. If not, to achieve convergent solutions or for computational conveniences some modifications in the terms ui in (5) are implemented as in the articles [11,12], without a proper mathematical evidence and support.
A Modified Decomposition Method of Adomian
To overcome the divergence of classical ADM or to speed up the convergence rate of the ADM series, the leading order term u0 in (5) (which is in compliance with the previous implementations, in for instance [11]) or some of the early terms, call ue, 0≤e≤Me, in (5) will be modified so that the modification will account for the change in the ADM by incorporating parameterized terms ũe(h) with h an embedded unknown parameter [3]. To exemplify, consider the modified version of (5) in the manner
The following conditions for parameterized terms in the new algorithm (7) should be added
ũ0(0)=ũ1(0)=⋯=ũe(0)=0,
so that it can be reduced to the traditional ADM for h = 0.
It is remarked that there is no a unique way of selecting the ũe(h) terms, but they must be as simple as possible in line with the simplicity inspiration of the ADM method itself. For instance, constants depending on h or simple power functions (or simple integrable functions) in combination with h can be used, refer to the illustrations belove for more details. We should remark that in the absence of tilde terms (or h) (7) dully conforms with (5). Having inserted such tilde terms in the modified ADM (7), the question now arises, how to determine the correct or proper value of h? So far, in ADM applications of this kind, only random or trial values are preassigned, see for instance [11] amongst others. On the other hand, we may either observe the least change in the plot of some physical quantities, the method is so-called as the constant h-level curves [36] giving rise to rough estimates for h, or we outline the following rigorous algorithm to determine the best or optimal h.
Algorithm. Consider the squared residual error corresponding to (1) defined by
Res(h)=∥L(u)-N(u)-f∥2,
where either L2(Ω) or L(Ω) norms are employed. Moreover, u in (8) is owing to the modified ADM from (7). The optimum value of h is such that it minimizes the squared residual error (8) at the approximation level M.
As a consequence, the above Algorithm will generate the best value of h which will ensure the convergence of ADM series solution (7) in a fastest rate of convergence. The minimization task of (8) may be fulfilled by means of contemporary softwares, such as MAPPLE or MATHEMATICA.
Applications
Potential applications of the introduced ADM in (7) are given here. To control the error, we use the norm
err=∥ue-u(t)∥
with the exact ue and ADM solution u.
An Algebraic Equation
As stated by Adomian [1] the classical ADM method (5) fails to result in a convergent solution of
u2-2u-2=0,
for the solution u = −0.73205080757. On the other hand, when the new ADM is built via
u0=-1-h,u1=h+12u02,un+1=12∑k=0nukun-k,n≥1,
Fig. 1 displays h-level curves at selected truncation orders M. The interval h∈(-0.35,-0.25) is observed to yield convergent ADM solutions, which excludes the traditional ADM with h = 0 [1].
Through the residual minimization
Res(h)=u2-2u-2
at the approximation level M = 8, h = −0.2679492 is obtained as the optimum. The history and why this value is the best for the convergence control, as compared to the failure of classical ADM can be visualized in Tab. 1.
Convergence control parameter h regarding (11)
Convergence of modified ADM (12) for Eq. (11) with different h. Paranthesis is for the absolute error (10)
M
h = 0.0
h = −0.20
h = −0.2679492
h = −0.40
2
−1.000000000
−0.7760000000
−0.7320508035 (4.06×10-9)
−0.6880000000
4
−1.250000000
−0.7707200000
−0.7320508054 (2.17×10-9)
−0.7390400000
8
−4.234375000
−0.7794992000
−0.7320508069 (6.24×10-10)
−0.7295606400
12
−33.22070313
−0.8094297223
−0.7320508074 (1.79×10-10)
−0.7332243215
16
−349.4215698
−0.8769006132
−0.7320508075 (5.15×10-11)
−0.7314507918
20
−4071.603386
−1.0260654160
−0.7320508076 (1.48×10-11)
−0.7323668046
The convergent solution of (11) with the new modified ADM (12) at the approximation level M = 8 is found to be
We find h∈(-1500,1500) for the convergence interval in Fig. 2, refer also to [3]. A comparison and convergence accelerating feature of the present ADM can be visualized from the Tab. 2. It is noticed from Tab. 2 that even the 4th-order modified ADM (16) is able to produce the solution accurate of order 10−10, which is adequate in practical purposes.
Convergence control parameter h regarding (15)
Errors in (15) regarding (14). Parenthesis denotes the optimum values of h
M
|ue(1) − u(1)| in [30]
|ue(1) − u(1)| present
2
3.568238353×10-4
2.585085279×10-6
(0.8164225489)
4
8.462701192×10-8
4.525954643×10-11
(0.7993559975)
5
1.233325816×10-9
1.783106924×10-13
(0.7958113139)
10
6.276712343×10-19
1.264455451×10-25
(0.7889813873)
A Fin with Porosity Feature
As taken from [14], a porous fin can be modelled via
u′′=su2+β2u,0<η<1,u′(0)=0,u(1)=1.
In (17), temperature along the fin is u, and s and β are physical parameters.
At the selected values s = 5 and β=1, the history of convergence is tabulated in Tab. 3. It is seen how useful the modified ADM over the classical one by more than doubling the convergence rate at the same number of iterations. The CPU times evaluated by MATHEMATICA through the residual are shown in Tab. 3. The list proves the improved efficiency of the current method over the classical ADM.
Residual errors of classical ADM and modified ADM (18) for Eq. (17) with s = 5 and β=1. Parenthesis denote the optimum values of h
M
Res(h) ADM (h = 0)
CPU
Res(h) modified ADM
CPU
2
9.275340168×10-1
0.031
2.9362768064×10-2
(−0.5306610965)
0.031
4
1.570940630×10-1
0.156
4.4652370216×10-4
(−0.4863076732)
0.178
8
2.113503722×10-3
0.625
5.6679833941×10-8
(−0.4699160217)
0.703
12
1.711030679×10-5
2.026
4.3828152149×10-12
(−0.4655830387)
2.203
The values of u(0) and u′(1) are listed in Tab. 4 with M = 12 from the ADM and also with the optimums hopt. The ADM (18) excellently compares with the full numerical solutions and hence, they can be used without consulting a numerical verification.
Numerical and modified ADM results regarding (17) [3]
s
θ0
u′(1)
θ0 (OADM)
u′(1) (OADM)
hopt
1
0.5227380936
1.1393789158
0.5227380936
1.1393789158
−0.4071445322
2
0.4487571832
1.4182574617
0.4487571832
1.4182574617
−0.4410040002
3
0.3977684724
1.6480019963
0.3977684724
1.6480019963
−0.4567006260
4
0.3596753295
1.8474903623
0.3596753295
1.8474903623
−0.4635604524
5
0.3297469637
2.0261010880
0.3297469637
2.0261010880
−0.4655830387
Instead of the modification of ADM in (18), we may use the Duan–Rach formulation involving no unknown parameters within it except the embedded parameter h
Choosing s=β=1, Fig. 3 shows h −level curves, indicating that the prescription of h must be made within the range (0.39, 0.49) to ensure the convergence of the modified ADM 19. The exact values of u(0) and u′(1) are respectively, 0.522738093570 and 1.13937891581.
Convergence control parameter h regarding (17)
To demonstrate the power of the modified ADM (19), Tab. 5 shows the squared residual error (9) from both the novel and classical ADM. It is unfortunate to observe that the Duan–Rach formulation (19) with h = 0 fails to converge, however, the optimum embedded parameter h insures that the modified ADM is convergent for the present physical problem, even if the convergence is not as fast as the modified formulation in (18).
Values of Res(h) for (17) from the classical and modified ADM methods
M
Res(0)
Res(h)
h
4
0.91×101
9.379592752×10-3
0.5717083134
8
5.90×102
8.485113566×10-5
0.4466738418
12
1.90×105
1.712138085×10-5
0.4205569783
16
9.50×107
8.888688579×10-6
0.4108495377
20
5.68×1010
3.872408211×10-6
0.4061899784
Gelfand Equation
The Gelfand equation [5] involves exponential nonlinearity [8]
u′′=eu,u(η=0)=u(η=1)=0,
with u′(0)=-0.4636325917 [3].
In line with the publications [5,8] when h = 0, the modified ADM is
where An(h,η) are the polynomials of Adomian for eu(η), encompassing all boundary conditions in (20) within the Duan–Rach approach, and we have
u′(0)=-∫ 01∫ 0ηAndηdη.
Fig. 4 shows the predicted convergence control parameters. With M = 12, an optimum value for the embedding parameter h is found to be −0.01274 from the Algorithm in (9). We find that the residual error is 1.6835117258×10-15 with u′(0)=-0.4636325917 from the present approach. On the other hand, and error of 4.813740934×10-9 with u′(0)=-0.4636325899 is calculated from the traditional ADM. Our value is given by
The success of the present modified ADM (21) is thus obvious.
Convergence control parameter h regarding (20)
Electrostatic Cantilever Micro-Electromechanical System
The beam-type electrostatic actuators for the nonlinear cantilever micro-electro mechanical systems are modelled by the fourth-order boundary value problem from [11]
where the Adomian polynomials An(h,η) are due to the negative-power nonlinearities in (23), see [11] for more details.
For the fixed parameters K = 3, αK=0.2, β=0.5 and γ=0.25, Fig. 5 produced from (24) helps us guess the proper values of embedded parameter h. Small values of h are seen to be adequate to get convergent ADM solutions.
Convergence control parameter h regarding (23)
Fig. 6 demonstrates different approximation levels M, and it signifies to h = 0.1045421730 as the optimum h when M = 10. With this optimum value of the embedding parameter, the squared residual error for the current problem is Res=3.2727623195×10-9, whereas the classical residual error with h = 0 is Res=8.7520447833×10-5. This implies that much more Adomian series terms are required for the classical ADM to reach the accuracy of the modified ADM here. Moreover, from [11] the physical values are u′′(0)=-0.56764 and u′′′(0)=1.07585 evaluated with 12 Adomian polynomials. The present values are u′′(0)=-0.5676684138 and u′′′(0)=1.0758898885 evaluated with only 10 Adomian polynomials, correct to 9 decimal places. Hence, the advantage of the present modification is clear for the present highly nonlinear physical problem. We should remark that with the present Algorithm, the validation of the ADM results against the numerical ones is no longer a prerequisite. Also, there is no need to evaluate error remainder functions, nor to consult to Pade-approximates to increase the accuracy, both of which as implemented in [11].
Error regarding (23) for various M
The convergence accelerating feature of the present modified ADM (24) as compared to the classical ADM is better visualized from the Tab. 6. Table also shows the comparable CPU times.
Convergence history of modified ADM vs. classical ADM for Eq. (23)
M
Res(0)
CPU
Res(h) (hopt)
CPU
2
4.2251049880×10-2
0.023
1.4791163690×10-2 (0.0777275068)
0.034
4
7.0270444089×10-3
0.041
3.9824140521×10-4 (0.0961306006)
0.058
6
1.4789254421×10-3
0.070
9.1320572997×10-6 (0.1020430905)
0.089
8
3.4800777077×10-4
0.105
1.8150755783×10-7 (0.1040090218)
0.131
10
8.7520447833×10-5
0.136
3.2727623195×10-9 (0.1045421730)
0.173
Electrostatic Cantilever Nano-Electromechanical System
Nonlinear model for the electrostatic double cantilever nano-electromechanical system in the case of Casimir force (K = 4) is given by [11]
where the Adomian polynomials An(h,η) are due to the negative-power nonlinearities in (25), see [11] for more details.
For the specific parameters αK=1, β=1.5 and γ=0.5, Fig. 7 displays the constant h-level curves drawn at the approximation level M = 8. Similar to the previous example, very small values of h are seen to be adequate to get convergent ADM.
Convergence control parameter h regarding (25)
In order to evaluate the performance of modified ADM over the classical one, Tab. 7 shows the unknown physical quantities u′′(0)=u′′(1) and u′′′(1) at several truncation orders M. The faster convergence rate of the present ADM is apparent.
Values of u′′(0) and u′′′(1) for (25) from the classical and modified ADM methods
M
u′′(0) (h = 0)
u′′′(1) (h = 0)
u′′(0)
u′′′(1)
h
2
−0.2533482143
−1.5156250000
−0.2534302817
−1.5159601707
0.0053058312
4
−0.2534436938
−1.5160577858
−0.2534438205
−1.5160583465
0.0061391414
6
−0.2534438213
−1.5160583545
−0.2534438215
−1.5160583555
0.0063873039
8
−0.2534438215
−1.5160583555
−0.2534438215
−1.5160583555
0.0064510177
To illustrate, the analytical formula computed via the present algorithm (26) at M = 4 for the value of u′′(0) is given by
which is of almost nine degree of accuracy as seen from Tab. 7.
Lane–Emden Equation
We consider the Lane–Emden type boundary value problem from [9]
u′′+kxu′-rus+u=0,0<x<1,u′(0)=0,5u(1)+u′(1)=5,
that models the oxygen diffusion in a spherical cell with Michaelis–Menten uptake kinetics. We take into account the subsequent constants to comply with the literature [9]
k=4,r=0.76129,s=0.03119.
The modified ADM that is offered for the present problem is then
which conforms with the classical ADM of [9] in the limit h→0 and the Adomian polynomials An(h,η) are to account for the nonlinearity rus+u. The interest is to determine the physical value of u(0) = a.
We present Tab. 8 to demonstrate the performance of the modified ADM (29) versus the classical ADM. The expected practical accuracy is met at lower Adomian series approximations via the modified method.
Values of u(0) = a and Res(h) for (28) from the classical and modified ADM methods
M
Res(0)
Res(h)
a (h = 0)
a
h
2
5.05920434×10-5
1.82258949×10-5
0.896870943665
0.896877048532
−0.0289267765
4
2.24898539×10-7
2.18219055×10-8
0.896877026028
0.896877046709
−0.0321993310
6
1.12229262×10-9
2.83931728×10-11
0.896877046624
0.896877046709
−0.0335012655
The Fluid Flow of Jeffery–Hamel
The Jeffery–Hamel fluid flow problem is modelled via [22]
where An(h,η)=-(4-Ha)α2un′(η)-2αRe∑k=0nuk(η)un-k′(η) are the Adomian polynomials.
For the diverging channel, considering the specific parameters α=50, Re = 50 and Ha = 1000 to be in line with [13], the constant h-level curves are depicted in Fig. 8 at the approximation order M = 8. Accordingly, we expect the convergency of the modified ADM (31) in the vicinity of zero.
Convergence control parameter h regarding (30)
The performance of modified ADM (31) is next measured by computing the centerline velocity u(0.5) (numerical value is 0.764064240111) at different approximation levels M as shown in Tab. 9. It is observed that 10 digits of accuracy is quickly reached by the present ADM, whereas the classical ADM falls behind. Hence, even though it was not clearly mentioned in [13] (see Tab. 1 therein), the accuracy of order 10−8 as obtained via the classical ADM demands at least 15–20 Adomian polynomials, whereas only 6 Adomian polynomials are sufficient to gain the same accuracy with the present modification.
Values of u(0.5) and Res(h) for (30) from the classical and modified ADM methods
M
Res(0)
Res(h)
u(0.5) (h = 0)
u(0.5)
h
2
6.92050730×10-2
1.13520536×10-2
0.763699786377
0.764097626541
−0.0180187189
4
1.94819917×10-3
1.87302552×10-5
0.764051346961
0.764064279681
−0.0179816107
6
7.53651632×10-5
3.95453415×10-8
0.764063696627
0.764064239975
−0.0178807933
8
3.34939742×10-6
9.72265613×10-11
0.764064214850
0.764064240111
−0.0178880579
Squeezing Two Parallel Plates
The flow squeezed between two parallel plates are modelled by the nonlinear equations [22]
where
An(h,η)=SA(1-ϕ)2.5(ηun′′′(η)+3un′′(η)+∑k=0n(uk′(η)un-k′′(η)-uk(η)un-k′′′(η)))
are the Adomian polynomials.
To make a comparison with the classical ADM in [14], we set the parameters S = 1, ϕ=0.02 and concentrate on Cu-Water nanofluid with ρs=8933 and ρf=997.1. The corresponding constant h-level curves are plotted in Fig. 9. It appears that the convergence of the modified ADM (32) is guaranteed if h is selected in the overlapping interval (-1,1).
Convergence control parameter h regarding (32)
The effects of iterative number M on the skin friction u′′(1) are next demonstrated in Tab. 10. The advantage of the modified ADM (33) with optimum values of embedding parameter h is clearly observed against the classical ADM with h = 0. It appears that the accuracy of classical ADM of O(10−8) as displayed in Tab. 2 of [18] clearly requires more Adomian series terms (see M = 8 in Tab. 10) as compared to the less terms needed in the modified ADM here.
Values of u′′(1) and Res(h) for (32) from the classical and modified ADM methods
M
Res(0)
Res(h)
u′′(1)(h=0)
u′′(1)
h
2
2.67740258×10-1
6.68835495×10-2
−3.6764267174
−3.6917310584
−0.0991340307
4
1.61057990×10-2
7.92481375×10-4
−3.6952425904
−3.6965216730
−0.1097576248
6
1.09939458×10-3
6.61056361×10-6
−3.6964553942
−3.6965506329
−0.1111585030
8
8.10446170×10-5
7.71479681×10-8
−3.6965431988
−3.6965505168
−0.1112323962
10
6.28790630×10-6
1.09803715×10-9
−3.6965499276
−3.6965505117
−0.1112270087
The following fourth-order modified ADM series solution for the skin friction may serve good to the purpose of engineering applications if not high accuracy is required
The classical Adomian method with h = 0 in (36) is not convergent, whereas with h = 0.68981924, the residual error becomes Res(h)=0.02589439, and the convergency is satisfied as revealed in Fig. 10 for the domain of definition.
Convergence control parameter h regarding (35)
Diffusion Equation
Let us consider the nonlinear diffusion equation, see [18] and [3]
ut+u+u3=(u2ux)x,u(x,0)=ex3,0≤x≤1,0≤t≤1,
for which [18] presents an exact solution.
The form of modified ADM for the partial differential equation (37) is
where
An(h,x,t)
are the Adomian polynomials corresponding to the function uxx − uux.
The traditional ADM with h = 0 in the domain 0≤x≤1,0≤t≤1 is noticed to fail to produce a convergent series as also clear in Fig. 12, because h = 0 is not lying in the overlapping region.
With the definition
Res(h)=||ut+uux-uxx||2=∫ 01∫ 01(ut+uux-uxx)2dxdt,
Tab. 12 justifies the success of the present modified ADM over the classical divergent one, both in terms of accuracy and computational cost.
The residual errors Res(h) for (41) from the classical and modified ADM methods
M
Res(0)
CPU
Res(h)(hopt)
CPU
5
1.75×102
0.012
4.674547981×10-2 (2.2068296567)
0.039
10
1.67×105
0.053
4.082713247×10-3 (2.3147042279)
0.078
15
1.39×108
0.101
6.755704561×10-4 (2.3159080855)
0.145
20
1.24×1011
0.242
2.741368603×10-5 (2.3820449852)
0.289
Concluding Remarks
The aim of the present work is to present superiority over the well-known Adomian decomposition method (ADM) often employed in the recent literature to analytically approximate solutions to highly nonlinear algebraic and differential equations of some real physical motions. Within this aim, a reformulation of the ADM is targeted to prevent first the failure and then convergence acceleration of the classical Adomian polynomials.
To accomplish the objective, the classical ADM is modified by inserting some simple parameterized terms into the early iterates involving an embedded parameter to control and pacing the convergence of the generated ADM series. In order to determine the best suitable value or the optimum value of this parameter, squared residual minimizing of the governing equation is proposed. This enables us to overcome the divergence of the classical ADM, and more importantly, there is no need to check out the results against the numerical ones, as usually has to be done in traditional ADM studies, since the optimum embedded parameter obtained is an insurance for ADM series convergence in a most rapid manner.
Physical examples selected from the recent application of ADM demonstrate the validity, accuracy and power of the present approach in terms of generating the convergent solution within the least number of iterations. In particular, the Duan-Rach modification of the ADM incorporating all the boundaries mostly used in the recent ADM applications takes great benefit of the present proposal, otherwise there is always the inevitable danger that it may lead to non physical solutions. The present approach successfully extends the convergence interval of the studied problem. In conclusion, the present formulation of ADM offers a promising tool to treat more strongly nonlinear equations/systems of real life phenomena.
Funding Statement: The author received no specific funding for this study.
Conflicts of Interest: The author declares that he has no conflicts of interest to report regarding the present study.
ReferencesAdomian, G., Rach, R. (1985). On the solution of algebraic equations by the decomposition method. ,105(1),141–166. DOI 10.1016/0022-247X(85)90102-7.Adomian, G. (1988). A review of the decomposition method in applied mathematics. ,135(2),501–544. DOI 10.1016/0022-247X(88)90170-9.Turkyilmazoglu, M. (2017). Parametrized Adomian decomposition method with optimum convergence. ,27(4),1–22. DOI 10.1145/3106373.Turkyilmazoglu, M. (2019). Accelerating the convergence of decomposition method of Adomian. ,31,54–59. DOI 10.1016/j.jocs.2018.12.014.Duan, J. J., Rach, R., Baleanu, D., Wazwaz, A. M. (2012). A review of the Adomian decomposition method and its applications to fractional differential equations. ,3(2),73–99.Abbasbandy, S. (2003). Improving newton-raphson method for nonlinear equations by modified Adomian decomposition method. ,145(2–3),887–893. DOI 10.1016/S0096-3003(03)00282-0.Chun, C. (2006). A new iterative method for solving nonlinear equations. ,178(2),415–422. DOI 10.1016/j.amc.2005.11.055.Wazwaz, A. M. (2012). A reliable study for extensions of the bratu problem with boundary conditions. ,35(7),845–856. DOI 10.1002/mma.1616.Wazwaz, A. M., Rach, R., Duan, J. S. (2013). Adomian decomposition method for solving the volterra integral form of the lane-emden equations with initial values and boundary conditions. ,219(10),5004–5019. DOI 10.1016/j.amc.2012.11.012.Bhanja, D., Kundu, B., Mandal, P. K. (2013). Thermal analysis of porous pin fin used for electronic cooling. ,64,956–965. DOI 10.1016/j.proeng.2013.09.172.Duan, J. S., Rach, R., Wazwaz, A. M. (2013). Solution of the model of beam-type micro-and nano-scale electrostatic actuators by a new modified Adomian decomposition method for nonlinear boundary value problems. ,49,159–169. DOI 10.1016/j.ijnonlinmec.2012.10.003.Duan, J. S., Chaolu, T., Rach, R., Lu, L. (2013). The Adomian decomposition method with convergence acceleration techniques for nonlinear fractional differential equations. ,66(5),728–736. DOI 10.1016/j.camwa.2013.01.019.Dib, A., Haiahem, A., Said, B. B. (2014). An analytical solution of the MHD jeffery-hamel flow by the modified Adomian decomposition method. ,102,111–115. DOI 10.1016/j.compfluid.2014.06.026.Dib, A., Haiahem, A., Said, B. B. (2015). Approximate analytical solution of squeezing unsteady nanofluid flow. ,269,193–199. DOI 10.1016/j.powtec.2014.08.074.Wazwaz, A. M. (2007). A comparison between the variational iteration method and Adomian decomposition method. ,207(1),129–136. DOI 10.1016/j.cam.2006.07.018.Olivares, A. G. (2003). Analytic solution of partial differential equations with Adomian’s decomposition. ,32(3–4),354–368. DOI 10.1108/03684920310458584.Song, L., Wang, W. (2010). Approximate rational jacobi elliptic function solutions of the fractional differential equations via the enhanced Adomian decomposition method. ,374(31–32),3190–3196. DOI 10.1016/j.physleta.2010.05.057.Song, L., Wang, W. (2013). A new improved Adomian decomposition method and its application to fractional differential equations. ,37(3),1590–1598. DOI 10.1016/j.apm.2012.03.016.Babolian, E., Biazar, J. (2002). On the order of convergence of Adomian method. ,130(2),383–387. DOI 10.1016/S0096-3003(01)00103-5.El-Kalla, I. L. (2008). Convergence of the Adomian method applied to a class of nonlinear integral equations. ,21(4),372–376. DOI 10.1016/j.aml.2007.05.008.Abdelrazec, A., Pelinovsky, D. (2009). Convergence of the Adomian decomposition method for initial-value problems. ,27(4),749–766. DOI 10.1002/num.20549.Turkyilmazoglu, M. (2016). Determination of the correct range of physical parameters in the approximate analytical solutions of nonlinear equations using the Adomian decomposition method. ,13(6),4019–4037. DOI 10.1007/s00009-016-0730-8.Senturk, E., Coskun, S. B., Atay, M. T. (2018). Solution of jamming transition problem using Adomian decomposition method. ,35(5),1950–1964. DOI 10.1108/EC-12-2016-0437.Zare, F., Heydari, M., Loghmani, G. B., Wazwaz, A. M. (2017). Numerical investigation of the beam-type nano-electrostatic actuator model by using the birkhoff interpolation method. ,3(4),129–146. DOI 10.1007/s40819-017-0417-2.Heydari, M., Loghmani, G. B., Wazwaz, A. M. (2017). A numerical approach for a class of astrophysics equations using piecewise spectral-variational iteration method. ,27(2),358–378. DOI 10.1108/HFF-09-2015-0379.Wazwaz, A. M. (1999). A reliable modification of Adomian decomposition method. ,102(1),77–86. DOI 10.1016/S0096-3003(98)10024-3.Wazwaz, A. M., El-Sayed, S. M. (2001). A new modification of the Adomian decomposition method for linear and nonlinear operators. ,122(3),393–405. DOI 10.1016/S0096-3003(00)00060-6.Bakodah, H. O., Banaja, M. A., Alrigi, B. A., Ebaid, A., Rach, R. (2019). An efficient modification of the decomposition method with aconvergence parameter for solving korteweg de vries equations. ,31(4),1424–1430. DOI 10.1016/j.jksus.2018.11.010.Nuruddeen, R. I., Muhammad, L., Nass, A. M., Sulaiman, T. A. (2018). A review of the integral transforms-based decomposition methods and their applications in solving nonlinear pdes. ,7(1),262–280.Hamoud, A. A., Ghadle, K. P. (2017). The combined modified laplace with Adomian decomposition method for solving the nonlinear volterra-fredholm integro-differential equations. ,21,17–28. DOI 10.12941/jksiam.2017.21.017.Hamoud, A. A., Ghadle, K. P. (2018). Modified Adomian decomposition method for solving fuzzy volterra-fredholm integral equations. ,85(1–2),52–69. DOI 10.18311/jims/2018/16260.Rashidi, M. M., Domairry, G., Dinarvand, S. (2009). Approximate solutions for the burger and regularized long wave equations by means of the homotopy analysis method. ,14(3),708–717. DOI 10.1016/j.cnsns.2007.09.015.Rashidi, M. M., Shahmohamadi, H. (2009). Analytical solution of three-dimensional navier stokes equations for the flow near an infinite rotating disk. ,14(7),2999–3006. DOI 10.1016/j.cnsns.2008.10.030.Rashidi, M. M., Mohimanian, P. S. A. (2010). Analytic approximate solutions for unsteady boundary-layer flow and heat transfer due to a stretching sheet by homotopy analysis method. ,15(1),83–95. DOI 10.15388/NA.2010.15.1.14366.Rashidi, M. M., Freidoonimehr, N., Hosseini, A., Beg, O. A., Hung, T. K. (2014). Homotopy simulation of nanofluid dynamicsfrom a non-linearly stretching isothermal permeable sheet with transpiration. ,49(2),469–482. DOI 10.1007/s11012-013-9805-9.Liao, S. J. (2014). . Singapore: World Scientific. DOI 10.1142/8939.