Approximate Analytical Solution for the Forced Korteweg-de Vries Equation

1 Department of Mathematical Sciences, Faculty of Science, Universiti Teknologi Malaysia, 81310 Johor Bahru, Johor, Malaysia 2 Faculty of Computer and Mathematical Sciences, Universiti Teknologi MARA, 40450 Shah Alam, Selangor, Malaysia 3 UTM Centre for Industrial and Applied Mathematics, Universiti Teknologi Malaysia, 81310 Johor Bahru, Johor, Malaysia 4Department of Mathematics, Faculty of Science, University of Kordofan, North Kordofan State, El-Obied 51111, Sudan


Introduction
An analytical model of Tsunami propagation was proposed by Pelinovsky et al. [1] as follows: with where  = (, ) refers to the elevation of free water surface,  = (, ) represents the solid bottom, ℎ is assumed to be the constant mean water depth, and  ≈ √ℎ is the long wave speed with  being gravity acceleration.If the right-hand side of (1), that is, / = 0, then (1) will be the Kortewegde Vries equation which is completely integrable [2].But if / ̸ = 0, then the equation is difficult to be integrable where the term / is called forcing term.Therefore, (1) is known as "KdV equation with forcing term" or "forced KdV (fKdV)" equation [3].
The fKdV equation to be considered in this paper is as follows: where  and  are constant.This form (3) is an alternative approach in describing the governing equations for the basic hydrodynamic model of Tsunami generation, for example, by atmospheric disturbances (Akylas [4] and Nosov and Skachko [5]).The appropriate evolution equation, which asymptotically approximates the Boussinesq equation, leads to this fKdV equation (3).The forcing term in the fKdV can be assumed to be derivable from atmospheric disturbances.Various forms of this equation have been extensively studied (see Grimshaw et al. [6], Pelinovsky et al. [7]) and numerical results show that the solution contains the set of solitary waves.Evidently, the behaviour of Tsunami waves on the open ocean is considerably more complicated than the solitary wave model; nonetheless this fKdV equation has been used in Yaacob et al. [8] as a simple mathematical model that could describe the Tsunami waves modelling.In the present paper, the above-mentioned considerations motivate us to investigate the approximate analytical solutions of the fKdV equation by using HAM.
Up to now, to the best of our knowledge, the solutions of fKdV equation of (3) can only be obtained by numerical or perturbation techniques [3,9,10].Recently, the analytical solution of fKdV equation with a certain form of forcing term has been solved by Zhao and Guo [3], using Hirota direct method.Thus the main reason of this work is to solve fKdV equation by using the homotopy analysis method (HAM) for various forcing terms, including in [3].This approximate analytical HAM solution will be used to be compared with the analytical soliton-type solution of Zhao and Guo [3].This analysis would in a way demonstrate the effectiveness and the potential of HAM to solve difficult nonlinear equations.
The homotopy analysis method (HAM) is an analytic method to solve nonlinear partial differential problems, which was first introduced by Liao in 1992 [11].It is known that HAM has greater flexibility in the selection of a proper set of base functions for the solution and a much simpler way in the control of the convergence rate and region compared to perturbation approach [12][13][14][15][16][17][18][19].Furthermore, this analytic technique does not have restriction of nonperturbation methods, such as Lyapunov's artificial small parameter method, the -expansion method, and Adomian's decomposition method [20].The analytic technique had been applied successfully in many nonlinear problems in engineering and science [21,22], for example, the magnetohydrodynamic flows of non-Newtonian fluids over a stretching sheet [23], nonlinear progressive waves [24], free oscillations of positively damped systems with algebraically decaying amplitude [25], free oscillations of self-excited systems [26], and similarity boundary layer equations [27].
The paper is organized as follows.The HAM of the fKdV equation ( 3) is presented in Section 2. Section 3 briefly illustrates the HAM solution of KdV equation.Section 4 is divided into 2 parts.The first (Section 4.1) discusses the concept of HAM and is followed by some examples of external forces employed on the fKdV equation.Meanwhile, Section 4.2 compares the results of the analytical solution of [3] with the generated HAM solution.Conclusion is concisely laid out in Section 5.

Basic Ideas of HAM
We consider differential equation in the form of (, ), where  is the nonlinear operator and (, ) is the unknown function.
By Taylor expansion of (, ; ), in a power series with respect to , we obtain where With appropriate property of  0 (, ), ℎ, , and ℓ, the series of ( 6) is convergence at  = 1, and one has Let us define the following vectors: Differentiating the zeroth-order deformation equation ( 5) -times with respect to the embedding parameter  and dividing it by ! and finally letting  = 0, we obtain the following th-order deformation equation: where It is noted that HAM contains control-convergence parameter ℎ where it provides an easy way to adjust and control the convergence region.This can be done by plotting the ℎcurves.

HAM Solution on KdV
When the forcing term of (3) becomes zero, then the equation becomes the Korteweg-de Vries (KdV) equation.Wazwaz [2] gave a form of the exact solution of KdV equation.Nazari et al. [28] used HAM to obtain a comparable approximate analytical solution of the KdV equation, which is in good agreement with Wazwaz's exact solution of KdV equation.Figure 1 is given the graphical representation of the HAM solution of KdV equation, where the forcing term of (3) is zero; / = 0.
Further to this research output, we will now apply HAM to solve different fKdV equations based on various forcing terms.

Approximate Analytical Solution for fKdV Using HAM
For HAM, the zeroth-order deformation equation is given by and we use the th-order deformation problem now becomes In this paper, MATHEMATICA © is used to solve the nonlinear differential equation ( 16) for the case / =  2 .
It is found that the solution for 10th-order approximation is Liao [11] pointed out that the valid region of ℎ is a horizontal line segment.Looking at ℎ-curves in Figure 2, the valid region of convergence interval is −1.5 < ℎ < −0.4.We choose three different values of ℎ in the valid interval.The outcomes for the first case are shown in Figure 3.
The HAM solution at 10th-order approximation is presented with various values of ℎ, where convergence interval is shown in Figure 3. Figures 3(a   It is noted that the wave is normal with no turbulence in Figure 3(e) where then ℎ value is chosen to be 1.4.Comparing Figures 3(e) and 3(f) with the same value of ℎ, it is observed that the wave exhibits anomalous behaviour when time is increased.For all cases above, we found that the wave displays certain chaotic manner when we increase the time.Figure 3 shows clearly that the wave changes in abnormal mode corresponding to the quadratic force in the right-hand side of (18).

Forcing Term, 𝜕𝑓/𝜕𝑥 = sin(x). Consider the nonlinear
KdV equation with the quadratic term sin() as the forcing term.The fKdV equation is as follows: Similarly MATHEMATICA © is used to solve the nonlinear differential equation ( 16) for the case / = sin().The solution at 9th-order approximation is Liao [11] pointed out that the horizontal line segment will be the valid region of ℎ.Based on Figure 4, the admissible interval will be −1.7 < ℎ < −0.3, so we choose ℎ = −0.3,ℎ = −1, and ℎ = −1.7.
The outcomes for the sin() force term are shown in Figure 5.
Figure 5 shows the HAM solution at 9th-order approximation with three different convergent values of ℎ.Figures 5(a), 5(b), 5(c), and 5(d) are the outcomes for different intervals of time less than ℎ = −0.3.The interval for Figure 5(a) is 0 ≤  ≤ 0.1 and the wave portrays a normal wave without any disturbance.The wave seems to be disturbed and the transition seen clearly in Figures 5(b) and 5(c) after the time is increased.In Figure 5(d), the time is increased to 0.7 s and it is observed that the wave starts to fluctuate.Upon considering ℎ = −1, the output is shown in Figures 5(e) and 5(f) where the suddenness and large energy of waves are created in short period of time.Figures 5(g) and 5(h) are the outcomes of ℎ = −1.7 for 0 ≤  ≤ 0.06 and 0 ≤  ≤ 0.1 where the waves are compressed and propagate in short time.As a whole it can be concluded that the wave changes tremendously over time with the sine force incorporated in fKdV equation.The fKdV equation is as follows: Similarly MATHEMATICA © is used to solve the nonlinear differential equation ( 16) for the case / =   .The solution at 17th-order approximation is Based on Figure 6, the accepted convergence interval will be −1.7 ≤ ℎ ≤ −0.5.For further deliberation, the convergence values of ℎ are chosen as, ℎ = −0.5, ℎ = −1, and ℎ = −1.7.
The outcomes for the   force term are shown in Figure 7. Figure 7 presents the solution of HAM at 17th-order approximationwith different convergence values of ℎ.Figures 7(a) and 7(b) are observed for ℎ = −0.5 under the time interval of 0 ≤  ≤ 0.001 and 0 ≤  ≤ 0.05, respectively.The wave is undisturbed when  is 0.001 s but the wave changes over time.This similar scenario occurs in Figure 7(c), upon using ℎ = −1.0.When ℎ = −1.7 is used the outcome is depicted in Figure 7(d) for the interval of 0 ≤  ≤ 0.01.It is found that the wave steadily rises over at a point in Figure 7(d).

Comparison of fKdV HAM Solution with Zhao and
Guo's Analytical Solution.In 2009, Zhao and Guo [3] gave a number of analytical solutions corresponding to the fKdV equation by using Hirota's direct method.In this paper, we will compare the analytical solution of [3] with the corresponding HAM solution.This comparison would again reveal the effectiveness and potential of HAM to solve difficult nonlinear equations.
The fKdV equation used in [3] is as follows: with () =  arctan () +  1  +  0 .The parameters chosen by Zhao and Guo [3] and the same input are used in our HAM procedure and are shown in Figure 8, that are Figure 8 shows the comparison between the graphs of analytical solution of [3] with the HAM solution.Due to very high computational and some limitations on technologies being used, we have simplified (25), by making the term  as constant in the forcing term.For purpose of comparison with the result of [3], we have chosen to show four cases at  = 0.25 s, 0.5 s, 0.75 s, and 1 s.It should be noted that, since the time has been fixed as constant in the right-hand side term, then for each case of different ℎ, the ensuing different orders of approximation are obtained.HAM solution shows perfect agreement at (, 1) with the analytical solution of [3].As a whole, based on the graphs obtained, HAM solution shows good approximation with the analytical solution of [3].

Conclusion
In this paper, HAM is used to get the approximate analytical solutions of forced KdV equations for different applied forcing terms.With the existence of auxiliary parameter ℎ, HAM provides us with a straightforward way to control the convergence of approximation series whereby it makes a noticeable difference between HAM and the existing perturbation methods.Figures 2-7 show distinct fluctuations in the wave profiles of the HAM solutions, with various forcing terms being incorporated in the fKdV equations.Nonetheless, Figure 8 shows very close approximation of HAM solution with the analytical solution of Zhao and Guo [3].These outcomes bring to light that HAM has flexibility and potential to solve complicated nonlinear problems.Since fKdV equation is associated with modelling of Tsunami propagation, as proposed by Pelinovsky et al. [1], we are hopeful that this work would have some added insight into more complicated Tsunami propagation models.In order to achieve more practical oriented results, there are suggestions that further work could be considered under data-driven (measurements) framework, for example, a recent overview can be seen from [29].

Figure 1 :
Figure 1: The HAM solution of KdV whereby forcing term is zero.
) and 3(b) obtained for ℎ = −0.6,where the interval for Figure3(a) is 0 ≤  ≤ 0.2 and for Figure3(b) is 0 ≤  ≤ 0.3, respectively.It is observed that there is a slight formation on the waves when time had been increased from  is 0.2 s to 0.3 s.The phenomenon is similar to the result for ℎ = −1 of Figures3(c) and 3(d).