A Comparison between Modified Sumudu Decomposition Method and Homotopy Perturbation Method ()
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 [1] and used in [2] [3] [4] [5] [6] among many others. The second method is homotopy perturbation method (HPM) developed by He [7] and used from many researchers [8] [9] [10] . 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:
, (1)
with the initial condition;
. (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:
(3)
where D is the second order linear differential operator
, R is other
linear differential operator of less order than D, N represent the general nonlinear operator and
is the source term.
Taking the Sumudu transform of both sides of Equation (3), we get:
; (4)
Using the differentiation property of the Sumudu transform and given initial conditions, we have:
. (5)
If we apply the inverse operator
to both sides of the Equation (5), we obtain:
. (6)
where
represents the term arising from the source term and the prescribed initial conditions. Now, apply the Adomain decomposition method [11] [12] [13] [14] :
, (7)
The nonlinear term can be decomposed as:
, (8)
For some Adomian polynomials
that are given by:
Substituting Equation (7) and Equation (8) in Equation (6), we get:
. (9)
Accordingly, the formal recursive relation is defined by:
(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 [15] [16] . For specific cases, the function
can be set as the sum of two partial functions, namely
and
. In other words, we can set:
; (11)
Using (11), we introduce a qualitative change in the formation of recursive relation (10). The modified recursive relation can be identified by:
(12)
(13)
(14)
The solution through the modified Sumudu decomposition method is highly depend upon the choice of
and
.
2.2. Basic Idea of HPM
Consider the following general nonlinear differential equation,
, (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
, and
as a series with components
, homotopy polynomials, which can be calculated using the formula:
(16)
To illustrate the homotopy perturbation method (HPM), we consider (15) as;
(17)
with solution
. As a possible remedy, we can define homotopy
as follows:
where
is an integral operator with known solution
, which can be obtained easily. Typically, we may choose a convex homotopy in the form;
(18)
and continuously trace an implicitly defined curve from a starting point
to a solution function
. The embedding parameter p monotonically increase from zero to unit as the trivial problem
is continuously deformed to original problem
.
(19)
when
, Equation (18) corresponds to Equation (17) and Equation (19) becomes the approximate solution of Equation (17), i.e.
(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 [17] :
; (21)
With the initial condition:
. (22)
1) Using HPM
To solve Equations (21)-(22) by homotopy perturbation method, we construct the following homotopy:
(23)
Assume the solution of Equation (23) to be in the form:
(24)
Substituting (24) into (23) and comparing coefficients of terms with identical powers of p, leads to:
(25)
(26)
The given initial value admits the use of:
(27)
The solution reads:
(28)
, (29)
. (30)
Examining the components
and
in Equation (29) and Equation (30), we can easily observe that the last two terms in
and the first two terms in
are the self-canceling (noise terms) [18] . Hence, the non-noise terms in
yields the exact solution of Equations (21)-(22), given by:
. (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:
. (32)
Applying
to both sides of (32) we obtain;
. (33)
Substituting;
; (34)
And the nonlinear terms of;
. (35)
Into (33) gives;
(36)
This gives the modified recursive relation;
(37)
, (38)
(39)
The first few of the components are given by;
(40)
(41)
(42)
The solution in a closed form is given by;
. (43)
Example 2: Consider the following inhomogeneous nonlinear KdV equation [17] :
; (44)
With the initial condition:
. (45)
1) Using HPM
Using homotopy perturbation method like in Example 1, we obtain the following components:
(46)
(47)
(48)
It is obvious that the last three terms in
and the first three terms in
are the self-canceling (noise terms). Keeping the remaining non-noise terms in
leads to the exact solution of Equations (44)-(45), given by:
. (49)
2) Using ADSTM
Proceeding as in Example 1, Equation (44) becomes:
(50)
The modified Sumudu decomposition method admits the of a modified recursive relation given by:
(51)
(52)
(53)
Consequently, we obtain:
(54)
(55)
(56)
The exact solution is given by:
(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.