A Comparison between Modified Sumudu Decomposition Method and Homotopy Perturbation Method

Abstract

In this paper, we present a comparative study between the modified Sumudu decomposition method (MSDM) and homotopy perturbation method (HPM). The study outlines the important features of the two methods. The analysis will be explained by discussing the nonhomogeneous Kortewege-de Vries (KdV) problems.

Share and Cite:

Ahmed, S. (2018) A Comparison between Modified Sumudu Decomposition Method and Homotopy Perturbation Method. Applied Mathematics, 9, 199-206. doi: 10.4236/am.2018.93014.

1. Introduction

This paper outlines effective comparison between two powerful methods that were newly advanced. The first method is modified form of the Sumudu decomposition method (MSDM) developed by Devedra Kumar  and used in      among many others. The second method is homotopy perturbation method (HPM) developed by He  and used from many researchers    . The two methods give quickly convergent series and are successfully used to solve a large number of linear and nonlinear problems.

In the present study, we consider the nonhomogeneous KdV equation of the following type:

${U}_{t}+aU{U}_{x}+b{U}_{xxx}=f\left(x,t\right)$ , (1)

with the initial condition;

$U\left(x,0\right)=g\left(x\right)$ . (2)

where a and b are constants.

The nonlinear KdV Equation (1) is an important mathematical model with wide applications in quantum mechanics and nonlinear optics. The KdV equation has several applications to physical problems. It approximately describes the evolution of long water waves. In addition, it is used in various fields such as, shallow water waves, acoustic waves in plasma, and long internal waves in a density.

The main objective of this contribution is to introduce a comparative study to solve KdV equation by using two of the most recently developed methods, the (MSDM), and the (HPM). The two methods accurately compute the solution in a series form or in an exact form.

2. Analysis of Numerical Methods

2.1. Basic Idea of the MSDM

To illustrate the basic idea of this method, we consider a general non-homogeneous partial differential equation with the initial conditions of the form:

$\begin{array}{l}DU\left(x,t\right)+RU\left(x,t\right)+NU\left(x,t\right)=g\left(x,t\right)\\ U\left(x,0\right)=h\left(x\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}{U}_{t}\left(x,0\right)=f\left(x\right).\end{array}$ (3)

where D is the second order linear differential operator $D=\frac{{\partial }^{2}}{\partial {t}^{2}}$ , R is other

linear differential operator of less order than D, N represent the general nonlinear operator and $g\left(x,t\right)$ is the source term.

Taking the Sumudu transform of both sides of Equation (3), we get:

$S\left[LU\left(x,t\right)\right]+S\left[RU\left(x,t\right)\right]=S\left[g\left(x,t\right)\right]$ ; (4)

Using the differentiation property of the Sumudu transform and given initial conditions, we have:

$S\left[U\left(x,t\right)\right]={u}^{2}S\left[g\left(x,t\right)\right]+h\left(x\right)+uf\left(x\right)-{u}^{2}S\left[RU\left(x,t\right)+NU\left(x,t\right)\right]$ . (5)

If we apply the inverse operator ${S}^{-1}$ to both sides of the Equation (5), we obtain:

$U\left(x,t\right)=G\left(x,t\right)-{S}^{-1}\left[{u}^{2}S\left[RU\left(x,t\right)+NU\left(x,t\right)\right]\right]$ . (6)

where $G\left(x,t\right)$ represents the term arising from the source term and the prescribed initial conditions. Now, apply the Adomain decomposition method     :

$U\left(x,t\right)=\underset{n=0}{\overset{\infty }{\sum }}{U}_{n}\left(x,t\right)$ , (7)

The nonlinear term can be decomposed as:

$NU\left(x,t\right)=\underset{n=0}{\overset{\infty }{\sum }}{A}_{n}\left(U\right)$ , (8)

For some Adomian polynomials ${A}_{n}\left(U\right)$ that are given by:

${A}_{n}\left({U}_{0},{U}_{1},{U}_{2},\cdots ,{U}_{n}\right)=\frac{1}{n!}\frac{{\text{d}}^{n}}{\text{d}{\lambda }^{n}}{\left[N\left(\underset{n=0}{\overset{\infty }{\sum }}{\lambda }^{n}{U}_{n}\right)\right]}_{\lambda =0},\text{\hspace{0.17em}}\text{\hspace{0.17em}}n=0,1,2,\cdots .$

Substituting Equation (7) and Equation (8) in Equation (6), we get:

$\underset{n=0}{\overset{\infty }{\sum }}{U}_{n}\left(x,t\right)=G\left(x,t\right)-{S}^{-1}\left[{u}^{2}S\left[R\underset{n=0}{\overset{\infty }{\sum }}{U}_{n}\left(x,t\right)+\underset{n=0}{\overset{\infty }{\sum }}{A}_{n}\left(U\right)\right]\right]$ . (9)

Accordingly, the formal recursive relation is defined by:

$\begin{array}{l}{U}_{0}\left(x,t\right)=G\left(x,t\right),\\ {U}_{k+1}\left(x,t\right)=-{S}^{-1}\left[{u}^{2}S\left[R{U}_{k}+{A}_{k}\right]\right],\text{\hspace{0.17em}}k\ge 0.\end{array}$ (10)

The modified Sumudu decomposition method introduces a slight variation to the recursive relation (10) that will lead to the determination of the components of U in a faster and easier way   . For specific cases, the function $G\left(x,t\right)$ can be set as the sum of two partial functions, namely ${f}_{1}\left(x,t\right)$ and ${f}_{2}\left(x,t\right)$ . In other words, we can set:

$G\left(x,t\right)={f}_{1}+{f}_{2}$ ; (11)

Using (11), we introduce a qualitative change in the formation of recursive relation (10). The modified recursive relation can be identified by:

${U}_{0}={f}_{1},$ (12)

${U}_{1}={f}_{2}-{S}^{-1}\left[{u}^{2}S\left[R{U}_{0}+{A}_{0}\right]\right],$ (13)

${U}_{k+1}=-{S}^{-1}\left[{u}^{2}S\left[R{U}_{k}+{A}_{k}\right]\right];\text{\hspace{0.17em}}\text{\hspace{0.17em}}k\ge 1$ (14)

The solution through the modified Sumudu decomposition method is highly depend upon the choice of ${f}_{1}\left(x,t\right)$ and ${f}_{2}\left(x,t\right)$ .

2.2. Basic Idea of HPM

Consider the following general nonlinear differential equation,

$u-N\left(u\right)=f$ , (15)

where N is a nonlinear operator from Hilbert space H to H, u is an unknown function, and f is a known function in H.

The homotopy perturbation method u as a series with components ${u}_{n}$ , and $N\left(u\right)$ as a series with components ${H}_{n}$ , homotopy polynomials, which can be calculated using the formula:

${H}_{n}=\frac{1}{n!}\frac{{\text{d}}^{n}}{\text{d}{\lambda }^{n}}N{\left(\underset{i=0}{\overset{\infty }{\sum }}{u}_{i}{\lambda }^{i}\right)|}_{\lambda =0}$ (16)

To illustrate the homotopy perturbation method (HPM), we consider (15) as;

$L\left(v\right)=v\left(x\right)-f\left(x\right)-N\left(v\right)=0$ (17)

with solution $u\left(x\right)$ . As a possible remedy, we can define homotopy $H\left(v,p\right)$ as follows:

$H\left(v,0\right)=F\left(v\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}H\left(v,1\right)=L\left(v\right)$

where $F\left(v\right)$ is an integral operator with known solution ${v}_{0}$ , which can be obtained easily. Typically, we may choose a convex homotopy in the form;

$H\left(v,p\right)=\left(1-p\right)F\left(v\right)+pL\left(v\right)=0$ (18)

and continuously trace an implicitly defined curve from a starting point $H\left({v}_{0},0\right)$ to a solution function $H\left(u,0\right)$ . The embedding parameter p monotonically increase from zero to unit as the trivial problem $F\left(v\right)=0$ is continuously deformed to original problem $L\left(v\right)=0$ .

$v={v}_{0}+p{v}_{1}+{p}^{2}{v}_{2}+{p}^{3}{v}_{3}+\cdots$ (19)

when $p\to 1$ , Equation (18) corresponds to Equation (17) and Equation (19) becomes the approximate solution of Equation (17), i.e.

$v=\underset{p\to 1}{\mathrm{lim}}{v}_{0}+{v}_{1}+{v}_{2}+{v}_{3}+\cdots$ (20)

3. Application

In this section, we demonstrate the analysis of two methods by applying two methods to the following two Kortewege-de Vries (KdV) partial differential equations.

Example 1: Consider the following inhomogeneous nonlinear KdV equation  :

${u}_{t}+u{u}_{x}+{u}_{xxx}=\mathrm{sin}x+t\mathrm{cos}x\left(t\mathrm{sin}x-1\right)$ ; (21)

With the initial condition:

$u\left(x,0\right)=0$ . (22)

1) Using HPM

To solve Equations (21)-(22) by homotopy perturbation method, we construct the following homotopy:

$\frac{\partial v}{\partial t}-\frac{\partial {u}_{0}}{\partial t}=p\left(-u\frac{\partial u}{\partial x}-\frac{{\partial }^{3}u}{\partial {x}^{3}}+\mathrm{sin}x+t\mathrm{cos}x\left(t\mathrm{sin}x-1\right)-\frac{\partial {u}_{0}}{\partial t}\right)$ (23)

Assume the solution of Equation (23) to be in the form:

$v={v}_{0}+p{v}_{1}+{p}^{2}{v}_{2}+{p}^{3}{v}_{3}+\cdots$ (24)

Substituting (24) into (23) and comparing coefficients of terms with identical powers of p, leads to:

${p}^{0}:\frac{\partial {v}_{0}}{\partial y}-\frac{\partial {u}_{0}}{\partial y}=0.$ (25)

${p}^{1}:\frac{\partial {v}_{1}}{\partial t}=-{u}_{0}\frac{\partial {u}_{0}}{\partial x}-\frac{{\partial }^{3}{u}_{0}}{\partial {x}^{3}}+\mathrm{sin}x+{t}^{2}\mathrm{cos}x\mathrm{sin}x-t\mathrm{cos}x-\frac{\partial {u}_{0}}{\partial t}$ (26)

$⋮$

The given initial value admits the use of:

${u}_{0}\left(x,0\right)=0$ (27)

${u}_{0}\left(x,t\right)=0$ (28)

${u}_{1}\left(x,t\right)=t\mathrm{sin}x+\frac{{t}^{3}}{3}\mathrm{sin}x\mathrm{cos}x-\frac{{t}^{2}}{2}\mathrm{cos}x$ , (29)

${u}_{2}\left(x,y\right)=-\frac{{t}^{3}}{3}\mathrm{sin}x\mathrm{cos}x+\frac{{t}^{2}}{2}\mathrm{cos}x-\frac{{t}^{4}}{3}{\mathrm{sin}}^{2}x-\frac{{t}^{4}}{3}{\mathrm{sin}}^{3}x+\cdots$ . (30)

Examining the components ${u}_{1}$ and ${u}_{2}$ in Equation (29) and Equation (30), we can easily observe that the last two terms in ${u}_{1}$ and the first two terms in ${u}_{2}$ are the self-canceling (noise terms)  . Hence, the non-noise terms in ${u}_{1}$ yields the exact solution of Equations (21)-(22), given by:

$u\left(x,t\right)=t\mathrm{sin}x$ . (31)

Notes on HPM:

From the previous analysis, we can observe that:

・ HPM can be applied it to various nonlinear problems. The main disadvantage is that we should suitably choose an initial guess.

・ HPM needs some modification to the rapid convergence of the series solution.

To overcome these disadvantages of HPM, the following ADSTM method is suggested.

2) Using MSDM

By taking Sumudu transform for (21) and using (22) we obtain:

$S\left[U\left(x,t\right)\right]=u\mathrm{sin}x+2{u}^{3}\mathrm{cos}x\mathrm{sin}x-{u}^{2}\mathrm{cos}x-uS\left[U{U}_{x}+{U}_{xxx}\right]$ . (32)

Applying ${S}^{-1}$ to both sides of (32) we obtain;

$U\left(x,t\right)=t\mathrm{sin}x+\frac{{t}^{3}}{3}\mathrm{cos}x\mathrm{sin}x-\frac{{t}^{2}}{2}\mathrm{cos}x-{S}^{-1}\left[uS\left[U{U}_{x}+{U}_{xxx}\right]\right]$ . (33)

Substituting;

$U\left(x,t\right)=\underset{n=0}{\overset{\infty }{\sum }}{U}_{n}\left(x,t\right)$ ; (34)

And the nonlinear terms of;

$U{U}_{x}=\underset{n=0}{\overset{\infty }{\sum }}{A}_{n}$ . (35)

Into (33) gives;

$\begin{array}{c}\underset{n=0}{\overset{\infty }{\sum }}{U}_{n}\left(x,t\right)=t\mathrm{sin}x+\frac{{t}^{3}}{3}\mathrm{cos}x\mathrm{sin}x-\frac{{t}^{2}}{2}\mathrm{cos}x\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-{S}^{-1}\left[uS\left[\underset{n=0}{\overset{\infty }{\sum }}{A}_{n}+{\left(\underset{n=0}{\overset{\infty }{\sum }}{U}_{n}\left(x,t\right)\right)}_{xxx}\right]\right]\end{array}$ (36)

This gives the modified recursive relation;

${U}_{0}\left(x,t\right)=t\mathrm{sin}x,$ (37)

${U}_{1}\left(x,t\right)=\frac{{t}^{3}}{3}\mathrm{cos}x\mathrm{sin}x-\frac{{t}^{2}}{2}\mathrm{cos}x-{S}^{-1}\left[uS\left[{A}_{0}+{\left({U}_{0}\right)}_{xxx}\right]\right]$ , (38)

${U}_{k+1}\left(x,t\right)=-{S}^{-1}\left({A}_{k}+{U}_{k}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}k\ge 1.$ (39)

The first few of the components are given by;

${U}_{0}\left(x,t\right)=t\mathrm{sin}x,$ (40)

${U}_{1}\left(x,t\right)=\frac{{t}^{3}}{3}\mathrm{cos}x\mathrm{sin}x-\frac{{t}^{2}}{2}\mathrm{cos}x-{S}^{-1}\left[uS\left[{A}_{0}+{\left({U}_{0}\right)}_{xxx}\right]\right]=0,$ (41)

${U}_{k+1}\left(x,t\right)=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}k\ge 1.$ (42)

The solution in a closed form is given by;

$U\left(x,t\right)=t\mathrm{sin}x$ . (43)

Example 2: Consider the following inhomogeneous nonlinear KdV equation  :

${u}_{t}+u{u}_{x}+{u}_{xxx}=\left(xt+\mathrm{sin}x\right)\left(t+\mathrm{cos}x\right)+x+\mathrm{cos}x$ ; (44)

With the initial condition:

$u\left(x,0\right)=\mathrm{sin}x$ . (45)

1) Using HPM

Using homotopy perturbation method like in Example 1, we obtain the following components:

${u}_{0}\left(x,t\right)=\mathrm{sin}x$ (46)

${u}_{1}\left(x,t\right)=xt+\frac{1}{3}x{t}^{3}+\frac{1}{2}x{t}^{2}\mathrm{cos}x+t\mathrm{cos}x+t\mathrm{sin}x\mathrm{cos}x+\frac{1}{2}{t}^{2}\mathrm{sin}x,$ (47)

${u}_{2}\left(x,t\right)=-t\mathrm{cos}x-t\mathrm{sin}x\mathrm{cos}x-\frac{1}{2}{t}^{2}\mathrm{sin}x-\cdots .$ (48)

It is obvious that the last three terms in ${u}_{1}$ and the first three terms in ${u}_{2}$ are the self-canceling (noise terms). Keeping the remaining non-noise terms in ${u}_{1}$ leads to the exact solution of Equations (44)-(45), given by:

$u\left(x,t\right)=xt+\mathrm{sin}x$ . (49)

Proceeding as in Example 1, Equation (44) becomes:

$\begin{array}{c}\underset{n\text{\hspace{0.17em}}=\text{\hspace{0.17em}}0}{\overset{\infty }{\sum }}{U}_{n}\left(x,t\right)=xt+\mathrm{sin}x+\frac{1}{3}x{t}^{3}+\frac{1}{2}x{t}^{2}\mathrm{cos}x+t\mathrm{cos}x+t\mathrm{sin}x\mathrm{cos}x\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\frac{1}{2}{t}^{2}\mathrm{sin}x-{S}^{-1}\left[uS\left[\underset{n=0}{\overset{\infty }{\sum }}{A}_{n}+{\left(\underset{n=0}{\overset{\infty }{\sum }}{U}_{n}\left(x,t\right)\right)}_{xxx}\right]\right]\end{array}$ (50)

The modified Sumudu decomposition method admits the of a modified recursive relation given by:

${U}_{0}\left(x,t\right)=xt+\mathrm{sin}x,$ (51)

$\begin{array}{c}{U}_{1}\left(x,t\right)=\frac{1}{3}xt+\frac{1}{2}x{t}^{2}\mathrm{cos}x+t\mathrm{cos}x+t\mathrm{sin}x\mathrm{cos}x+\frac{1}{2}{t}^{2}\mathrm{sin}x\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-{S}^{-1}\left[uS\left[{A}_{0}+{\left({U}_{0}\left(x,t\right)\right)}_{xxx}\right]\right],\end{array}$ (52)

${U}_{k+1}\left(x,t\right)=-{S}^{-1}\left[uS\left[{A}_{k}+{\left({U}_{k}\left(x,t\right)\right)}_{xxx}\right]\right],\text{\hspace{0.17em}}\text{\hspace{0.17em}}k\ge 1.$ (53)

Consequently, we obtain:

${U}_{0}\left(x,t\right)=xt+\mathrm{sin}x,$ (54)

${U}_{1}\left(x,t\right)=0,$ (55)

${U}_{k+1}\left(x,t\right)=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}k\ge 1.$ (56)

The exact solution is given by:

$U\left(x,t\right)=xt+\mathrm{sin}x.$ (57)

4. Conclusions

The main objective of this paper is to introduce a comparative study between modified Sumudu decomposition method and Homotopy perturbation method. The two methods are strong and effective methods that both hand approximations of higher reliability of series solution. A remarkable conclusion can made here. Modified Sumudu decomposition method gives rapid convergence of the series solution without showing noise terms. However, homotopy perturbation method provides the components of the series solution, where these components should have noise terms.

More importantly, the modified Sumudu decomposition method (MSDM) gives better efficiency in many cases, and this implies that the modified Sumudu decomposition method has an advantage over the homotopy perturbation method.

Conflicts of Interest

The authors declare no conflicts of interest.     customer@scirp.org +86 18163351462(WhatsApp) 1655362766  Paper Publishing WeChat 