_{1}

Finite Element Method (FEM), based on p and h versions approach, and the Adomians decomposition algorithm (ADM) are introduced for solving the Emden-Fowler Equation. A number of special cases of p and h versions of FEM are introduced. Several iterated forms of the ADM are considered also. To demonstrate the efficiency of both methods, the numerical solutions of different examples are compared for both methods with the analytical solutions. It is observed that the results obtained by FEM are quite satisfactory and more accurate than ADM. Moreover, the FEM method is applicable for a wide range of classes including the singularity cases with the given special treatments by the FEM. Comparing the results with the existing true solutions shows that the FEM approach is highly accurate and converges rapidly.

In this work, we present an alternate algorithm to solve the Emden-Fowler Equation [

Adomians decomposition algorithm has been recently employed to solve a wide range of problems (see [

d d t ( q ( t ) d u d t ) + q ( t ) u n = 0 , t ≥ 0 (1)

with specified initial condition

u ( 0 ) = α (2)

d u d t = β

The balance of this paper is as follows. In Section 2.1, we give a brief description of Adomians method and then introduce a modified version of this algorithm. We apply the modified scheme to Equations (1), (2). While in Section 2.2, we seek a finite element solution for solving Emden-Fowler Equation. We consider the mesh point u_{i} over the interval [0, 1], with x_{0} = 0 and x_{n} = 1, noting that the mesh points distance = h.

In Section 3, several interesting examples that arise in applications are used to illustrate the algorithm with error estimates.

And tables are not prescribed, although the various table text styles are provided. The formatter will need to create these components, incorporating the applicable criteria that follow.

In this section, we first describe algorithm of Adomians decomposition method as it applies to general nonlinear equation of the form

u − N ( u ) = f (3)

where N is a nonlinear operator on a Hilbert space H and f is a known element of H. We assume that for a given f a unique solution u of Equation (3) exists. We then introduce a modified version of this algorithm to handle equations of the form Equation (1).

The Adomians algorithm assumes a series solution for u given by

u = ∑ n = 0 ∞ u n (4)

And the nonlinear operator N to be decomposed into

N ( u ) = ∑ n = 0 ∞ A n (5)

where A_{n}’s are Adomians polynomials of u 0 , u 1 , ⋯ , u n given by

A n = 1 n ! d n d λ n [ N ∑ i = 0 ∞ λ i u i ] λ = 0 , n = 0 , 1 , ⋯ (6)

Substituting Equations (4) and (5) into the functional Equation (3) yields.

∑ n = 0 ∞ u n − ∑ n = 0 ∞ A n = f (7)

The convergence of the series in Equation (7) will yield

u 0 = f u 1 = A 0 u 2 = A 1 ⋮ u n = A n − 1 (8)

Thus, one can, recurrently determine every term of the series ∑ n = 0 ∞ u n the convergence of this series has been established (see [

To illustrate the scheme, let the nonlinear operator N(u) be a nonlinear function of u, say g(u) then the first four Adomians Polynomials are given by

A 0 = g ( u 0 ) A 1 = u 1 g ′ ( u 0 ) A 2 = u 2 g ′ ( u 0 ) + 1 2 ! u 1 2 g ″ ( u 0 ) ⋮ A 3 = u 3 g ′ ( u 0 ) + u 1 u 2 g ″ ( u 0 ) + 1 3 ! u 1 3 + g ‴ ( u 0 ) (9)

How do we interpret and solve Emden-Fowler Equation in this setting?

Following the Adomians decomposition analysis [

L t = ∂ ∂ t (10)

Equation (1) can be rewritten in terms of the linear operator

L t [ p ( t ) L t ( u ) ] + q ( t ) u n = 0 (11)

It was shown in [

Thus, the inverse operator of L_{t}, namely L t − 1 exists and is the twofold indefinite integral; i.e.

[ L t − 1 f ] ( t ) = ∫ 0 t d u ∫ 0 t d v f ( v ) (12)

Operating on both sides of (11) with L t − 1 yields

L t − 1 L t [ p ( t ) L t ( u ) ] + L t − 1 [ q ( t ) u n ] = 0 (13)

From which it follows, upon using the initial conditions given in Equation (2).

u ( t ) = u ( 0 ) + L t − 1 [ p ( 0 ) β p ( t ) ] − L t − 1 [ L t − 1 ( q ( t ) u n ) P ( t ) ] (14)

The Adomians decomposition method yields the solution in the form given in (4) i.e. u ( t ) = u 0 + u 1 + u 2 + ⋯ + u n ,then from Equation (14)

u 0 + u 1 + u 2 + ⋯ = 1 + L t − 1 [ p ( 0 ) β β ( t ) ] − L t − 1 [ L t − 1 ( q ( t ) u n ) p ( t ) ] (15)

For later numerical computation, let the expression

ϕ n = ∑ i = 0 n − 1 u i ( t ) (16)

can serve as a practical solution. We will show through several examples, that this Adomians decomposition method which converts the given equation to recurrences relation whose terms are computed using maple 15.

We seek a finite element solution for solving Emden-Fowler Equation, we consider the mesh point u_{i} over the interval [0, 1], with x_{0} = 0 and x_{n} = 1, note that the mesh points distance = h, with x n = x 0 + i h . Given

u ( x ) = ∑ n = 0 n − 1 α i φ i (17)

The finite element method minimizes the integral (17) over a class of Piecewise polynomials. The idea is to choose a finite number of trial functions φ i ,and among all their linear combinations ∑ n = 0 n − 1 α i φ i to find the one, which is the minimum, the unknown α i determined by a system of N discrete algebraic equations, which the computer can handle. Therefore, the goal is to choose trial functions φ i ,which are convenient enough for the given integral (17) to be compute and minimized, and at the same time general enough to approximate closely the unknown u. The software TWEPEP starts by a subdivision of the given region into smaller pieces which are triangles with standard six-node with a quadratic basis function, and with one edge curved when adjacent to a curved boundary according to the isoperimetric method. It is also optional to use 10-piont cubic (3rd degree) or the 15-points quartic (4th degree) isoperimetric triangular elements for greater accuracy. Each time a triangle is partitioned, it is divided by a line from the midpoint of its longest side to the opposite vertex. If this side is not on the boundary, the triangle which shares that side must also be divided to avoid non-conforming elements with discontinuous basis functions. An initial triangulation with sufficient triangles to define the region is supplied by the user, then the refinement and grading of this triangulation is guided by a user supplied function D3EST which should be largest where the final triangulation is to be densest. The Cuthill-McKee algorithm ( [_{ }

In all cases, the algebraic system solved by Newton’s method. One iteration per time step is done for parabolic problems and one iteration is sufficient for linear elliptic problems. The linear system is solved directly by block Gaussian elimination, without row interchanges since pivoting is unnecessary when the matrix is positive definite. Symmetry also taken advantage of in the elimination process if it is present then the storage and computational work halved. If the Jacobian matrix is too large to keep in core, the frontal method is used efficiently organize its storage out of core.

In this section, we shall consider three examples the first is the linearized Emden-Fowler Equation, and the other two are nonlinear. The out-come of modified Adomians decomposition method compared with and the finite element method the p, the h versions, and any known solution to the underlying Emden-Fowler Equation. Then solution obtained generalized by maple [

Example 1. Consider Emden-Fowler Equation of the form

1 t 2 d d t ( t 2 d d t ( u ) ) + t 2 u = 0

u ( 0 ) = 1 and d u ( 0 ) d t = 0 (18)

Substituting these values into the general formula Equation (13), we obtain

u 0 = 1 u 1 = − L t − 1 ( t 2 u 0 t 2 ) = − 1 6 t 2 u 2 = − L t − 1 ( t 2 u 1 t 2 ) = 1 120 t 4 u 3 = − L t − 1 ( t 2 u 2 t 2 ) = − 1 5040 t 6 ⋮ (19)

The exact solution to Equation (18) given by u ( t ) = 1 t sin t .

Example 2. Consider the Emden Fowler of the form

1 t 2 ( d d t ( t 2 d d t u ) ) + u 5 = 0

u ( 0 ) = 1 , d u ( 0 ) d t = 0 (20)

Using Equation (13), the various u_{i} are given by

u 0 = 1 u 1 = − L t − 1 ( t 2 A 0 t 2 ) = − 1 6 t 2 u 2 = − L t − 1 ( t 2 A 1 t 2 ) = − 1 24 t 4 u 3 = − L t − 1 ( t 2 A 2 t 2 ) = − 5 432 t 6 (21)

where A 0 = f ( u 0 ) = 1

A 1 = u 1 f ′ ( u 0 ) = 5 u 0 5 = − 5 6 t 2

A 2 = u 2 f ′ ( u 0 ) + 1 2 u 1 2 f ″ ( u 0 ) = 5 u 0 4 u 2 + 10 u 0 3 u 1 2 = 35 71 t 4 (22)

A 3 = u 3 f ′ ( u 0 ) + 1 2 ( 2 u 1 u 2 ) f ″ ( u 0 ) + 1 3 ! ( u 0 u 1 3 ) f ‴ ( u 0 ) = − 35 144 t 6

The exact solution of Equation (20) is given by

u ( t ) = 1 1 + 1 3 t

The exact solution compared with the numerical solution using Adomians method.

Example 3. Our final Example deals with Emden-Fowler Equation of the form

d d t ( 1 1 + t ( d d t ( u ) ) ) + 1 2 ( 1 + t ) 3 u 2 = 0 (23)

u ( 0 ) = 1 , d u ( 0 ) d t = 1 2 (24)

Using Equation (13) the various iterates are given by

u 0 = 1 + 1 2 L t − 1 [ L t − 1 [ 1 + t ] 1 / 2 ] = 2 3 + 1 3 ( 1 + t ) 3 / 2 u 1 = − 1 2 L t − 1 [ L t − 1 ( A 0 ( 1 + t ) 3 ) ( 1 + t ) 1 / 2 ] u 2 = − 1 2 L t − 1 [ L t − 1 ( A 1 ( 1 + t ) 3 ) ( 1 + t ) 1 / 2 ] u 3 = − 1 2 L t − 1 [ L t − 1 ( A 2 ( 1 + t ) 3 ) ( 1 + t ) 1 / 2 ] ⋮

where

A 0 = f ( u 0 ) = u 0 2 = [ 2 3 + 1 3 ( 1 + t ) 3 / 2 ] 2 (25)

A 1 = 2 u 0 u 1 A 2 = 2 u 0 u 1 + u 1 2 A 3 = 2 u 0 u 3 + 2 u 0 u 1 u 2 ⋮ _{ }

The exact solution of Equation (23) and (24) is given by

u ( t ) = 1 + t (26)

The exact solution compared with the numerical solution using Adomians method, shown in

In the next section, we will give the numerical results arising from the implementation of this adaptive collocation approach over the Emden-Fowler problem.

In this section, the ADM collocation method used to solve the Emden-Fowler problem for different values of the u_{i} using the computer algebra system Maple 15. In Tables 1-3, the numerical solutions obtained by the above method at the mesh points t = 0 , 0.01 , 0.02 , ⋯ , 0.1 . Respectively, are compared with the given analytic solutions the resulting relative error is also compared with h version FEM numerical techniques, namely, 50, 100 and 150 elements.

In Tables 4-6, the numerical solution obtained by the FEM collocation at the mesh points t = 0 , 0.01 , 0.02 , ⋯ , 0.1 ,is compared with the exact solutions given. A Fortran code called TWOPEP [

t | ADM Relative Error | FEM Relative Error | ||||
---|---|---|---|---|---|---|

Four terms | Five terms | Six terms | 50 | 100 | 150 | |

0.0 | 0.0 | 0.0 | 0.0 | 0.0 | 0.0 | 0.0 |

0.02 | 0.453 × 10^{−3} | 1.123 × 10^{−4} | 1.013 × 10^{−4} | 0.234 × 10^{−3} | 3.325 × 10^{−4} | 1.120 × 10^{−4} |

0.04 | 1.987 × 10^{−4} | 1.347 × 10^{−4} | 0.047 × 10^{−4} | 1.126 × 10^{−3} | 1.034 × 10^{−4} | 2.613 × 10^{−5} |

0.06 | 2.798 × 10^{−4} | 0.340 × 10^{−4} | 0.434 × 10^{−5} | 2.023 × 10^{−2} | 1.122 × 10^{−2} | 1.106 × 10^{−5} |

0.08 | 1.001 × 10^{−5} | 0.150 × 10^{−5} | 0.504 × 10^{−6} | 4.750 × 10^{−3} | 2.322 × 10^{−4} | 1.457 × 10^{−6} |

0.10 | 1.338 × 10^{−5} | 2.243 × 10^{−5} | 1.243 × 10^{−5} | 6.445 × 10^{−3} | 7.805 × 10^{−5} | 3.345 × 10^{−6} |

t | ADM Relative Error | FEM Relative Error | ||||
---|---|---|---|---|---|---|

Four terms | Five terms | Six terms | 50 | 100 | 150 | |

0.0 | 0.0 | 0.0 | 0.0 | 0.0 | 0.0 | 0.0 |

0.02 | 1.331 × 10^{−3} | 3.408 × 10^{−5} | 1.408 × 10^{−4} | 0.435 × 10^{−3} | 3.025 × 10^{−4} | 1.420 × 10^{−4} |

0.04 | 9.034 × 10^{−3} | 1.000 × 10^{−4} | 9.00 × 10^{−4} | 1.512 × 10^{−3} | 1.034 × 10^{−4} | 2.152 × 10^{−5} |

0.06 | 1.141 × 10^{−3} | 1.092 × 10^{−4} | 9.092 × 10^{−5} | 3.023 × 10^{−2} | 0.122 × 10^{−3} | 1.102 × 10^{−5} |

0.08 | 3.114 × 10^{−3} | 2.843 × 10^{−4} | 2.843 × 10^{−4} | 4.750 × 10^{−3} | 3.421 × 10^{−5} | 1.456 × 10^{−6} |

0.10 | 7.081 × 10^{3} | 0.774 × 10^{−4} | 0.075 × 10^{−4} | 6.445 × 10^{−3} | 6.505 × 10^{−5} | 3.345 × 10^{−6} |

t | Relative Error | Number of elements | ||||
---|---|---|---|---|---|---|

Four terms | Five terms | Six terms | 50 | 100 | 150 | |

0.0 | 0.0 | 0.0 | 0.0 | 0.0 | 0.0 | 0.0 |

0.02 | 4.335 × 10^{−2} | 3.425 × 10^{−2} | 1.425 × 10^{−2} | 0.335 × 10^{−3} | 3.325 × 10^{−4} | 1.420 × 10^{−4} |

0.04 | 1.512 × 10^{−2} | 1.034 × 10^{−2} | 1.634 × 10^{−2} | 1.512 × 10^{−3} | 1.034 × 10^{−4} | 2.653 × 10^{−5} |

0.06 | 3.023 × 10^{−2} | 1.122 × 10^{−2} | 1.002 × 10^{−2} | 3.023 × 10^{−2} | 1.122 × 10^{−2} | 1.102 × 10^{−5} |

0.08 | 4.750 × 10^{−3} | 2.322 × 10^{−4} | 1.022 × 10^{−4} | 4.750 × 10^{−3} | 2.322 × 10^{−4} | 1.456 × 10^{−6} |

0.10 | 6.445 × 10^{−3} | 7.805 × 10^{−4} | 2.805 × 10^{−4} | 6.445 × 10^{−3} | 7.805 × 10^{−5} | 3.345 × 10^{−6} |

F.E.M p version | F.E.M h version | ||||
---|---|---|---|---|---|

50 elements | 75 elements | 100 elements | 125 elements | 150 elements | |

Quadratic | 1.071 × 10^{−3} | 2.138 × 10^{−3} | 6.306 × 10^{−4} | 4.006 × 10^{−4} | 1.303 × 10^{−5} |

Cubic | 2.402 × 10^{−4} | 1.485 × 10^{−4} | 2.759 × 10^{−5} | 2.051 × 10^{−5} | 2.759 × 10^{−5} |

Quartic | 1.541 × 10^{−5} | 7.683 × 10^{−6} | 1.481 × 10^{−6} | 1.6211 × 10^{−6} | 0.381 × 10^{−6} |

F.E.M p version | F.E.M h version | ||||
---|---|---|---|---|---|

100 elements | 125 elements | 150 elements | 175 elements | 200 elements | |

Quadratic | 3.971 × 10^{−3} | 2.338 × 10^{−3} | 7.378 × 10^{−4} | 4.006 × 10^{−4} | 1.303 × 10^{−5} |

Cubic | 2.502 × 10^{−4} | 1.485 × 10^{−4} | 2.759 × 10^{−5} | 2.051 × 10^{−5} | 2.759 × 10^{−5} |

Quartic | 1.541 × 10^{−5} | 8.683 × 10^{−6} | 0.683 × 10^{−6} | 1.621 × 10^{−6} | 0.381 × 10^{−6} |

F.E.M p version | F.E.M h version | ||||
---|---|---|---|---|---|

100 elements | 125 elements | 150 elements | 175 elements | 200 elements | |

Quadratic | 3.271 × 10^{−3} | 2.138 × 10^{−3} | 6.306 × 10^{−4} | 4.006 × 10^{−4} | 1.303 × 10^{−5} |

Cubic | 2.402 × 10^{−4} | 1.485 × 10^{−4} | 2.759 × 10^{−5} | 2.051 × 10^{−5} | 2.059 × 10^{−5} |

Quartic | 1.541 × 10^{−5} | 7.683 × 10^{−6} | 1.681 × 10^{−6} | 1.621 × 10^{−6} | 0.331 × 10^{−7} |

From Tables 4-6, we observe that the h and the p collocation FEM method yields a almost exact solution for elements in order = 50, 100, and 150, respectively [

We observe generally notable change in the numerical results when the number of elements and the degree increased, we can say almost exact.

From Tables 1-3, it is evident that the relative error for all given examples decreases as the number of terms of the ADM increase, convergence is rapid and the relative error is small.

Likewise, result in Tables 4-6 shows the F.E.M error as the number of elements is subdivided (the h version) also as the degree of the polynomial is increased (p version). With more dense elements near t = 0 has the upper hand results.

In this work, the Adomians decomposition method, and the finite element method p and the h versions have been successfully applied to find the solution of nonlinear Emden-Fowler Equation. For general cases, Adomians decomposition method usually fails to solve singular initial value problems of Emden-Fowler type, while an effective modification of FEM that only requires denser elements around the singularity point is introduced. We conclude that FEM is a powerful method for the solution of nonlinear Emden-Fowler Equation. However, the computations by the ADM are simpler and faster than FEM techniques. The FME offers great advantages of straightforward applicability, computational efficiency and high accuracy.

The author declares no conflicts of interest regarding the publication of this paper.

Hossen, K.A.A. (2020) On Solutions of Emden-Fowler Equation. American Journal of Computational Mathematics, 10, 90-99. https://doi.org/10.4236/ajcm.2020.101006