Collocation Method for Solving the Generalized KdV Equation

Show more

1. Introduction

Several physical processes for example dispersion of long waves in shallow water waves under gravity, bubble-liquid mixtures, ion acoustic plasma waves, fluid mechanics, nonlinear optics and wave phenomena in enharmonic crystals can be expressed by the KdV equation which was first introduced by Korteweg and de Vries [1]. The equation was solved analytically by Zabusky, Fornberg and Whitham, [2] [3]. Zabusky and Kruskal [4] were first obtained numerical solutions of the equation with finite difference method. Gardner et al. [5] demonstrated existence and uniqueness of solutions of the KdV equation. Several scientists have used various numerical methods including pseudospectral method [4], finite difference method [6] [7], finite element method [8] - [15] and heat balance integral method [15] to solve the equation. Numerical solutions of the KdV equation were obtained using differential quadrature method based on cosine expansion by B. Saka [16]. Like the KdV equation, in recent years, various numerical methods have been improved for the solution of the MKdV equation. Kaya [17] calculated the explicit solutions of the higher order modified Korteweg de-Vries equation by Adomian decomposition method. MKdV equation has been solved by using Galerkins’ method with quadratic B-spline finite elements by Biswas et al. [18]. Raslan and Baghdady [19] [20] indicated the accuracy and stability of the difference solution of the MKdV equation and they obtained the numerical aspects of the dynamics of shallow water waves along lakes’ shores and beaches modeled by the MKdV equation. A new variety of (3 + 1)-dimensional MKdV equations and multiple soliton solutions for each new equation were established by Wazwaz [21] [22]. Lumped Petrov-Galerkin and Galerkin methods were practiced to the MKdV equation by Ak et al. [23] [24]. GKdV equation has received much less attention, presumably because of its higher nonlinearity for $p>2$. The symmetry group was calculated for the equation and several classes of solutions were obtained in [25]. Liu and Yi [26] developed and analyzed a Hamiltonian preserving DG method for solving the generalized KdV equation. The initial value problem of a kind of GKdV equations is considered by using Sobolev space theory and finite element method by Lai et al. [27]. Alvarado and Omel’yanov [28] create a finite differences scheme to simulate the solution of the Cauchy problem and present some numerical results for the problem of the solitary waves interaction. A class of fully discrete scheme for the generalized Korteweg-de Vries equation in a bounded domain $\left(\mathrm{0,}L\right)$ has studied by Sepulveda and Villagran [29]. Collocation finite element method based on quintic B-spline functions is applied to the generalized KdV equation by Ak et al. [30]. Solitary wave solution for the GKdV equation by using ADM has been obtained by Ismail et al. [31].

In this article, we will take in consideration for the following GKdV equation

${U}_{t}+\epsilon {U}^{2}{U}_{x}+\mu {U}_{xxx}=0,$ (1)

with the homogeneous boundary conditions

$\begin{array}{l}U\left(a,t\right)=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}U\left(b,t\right)=0,\\ {U}_{x}\left(a,t\right)=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}{U}_{x}\left(b,t\right)=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}t>0\end{array}$ (2)

and an initial condition

$U\left(x,0\right)={U}_{0}\left(x\right),\mathrm{}a\le x\le b,$ (3)

where t is time, x is the space coordinate, $\epsilon $ and $\mu $ are positive parameters. One of the primary mathematical models for describing the theory of water waves in shallow channels is the following Korteweg de-Vries (KdV) equation:

${U}_{t}+\epsilon U{U}_{x}+\mu {U}_{xxx}=0.$ (4)

The terms $U{U}_{x}$ and ${U}_{xxx}$ in the Equation (4) stand for the nonlinear convection and dispersion, respectively. In this paper, we have numerically solved the GKdV equation using collocation method with septic B-spline finite elements. We have investigated the motion of a single soliton wave to show the performance and profiency of the proposed method. Also we have showed the suggested method is unconditionally stable applying the von-Neumann stability analysis.

2. Septic B-Spline Collocation Method

We think of a mesh $a={x}_{0}<{x}_{1}<\cdots <{x}_{N}=b$ as a uniform divide of the solution area $a\le x\le b$ by the points ${x}_{m}$ with $h=\frac{b-a}{N}={x}_{m+1}-{x}_{m}$. The septic B-splines ${\varphi}_{m}\left(x\right)$, ( $m=-\mathrm{3,}-\mathrm{2,}\cdots \mathrm{,}N+3$ ) at the knots ${x}_{m}$ are given by [32]

${\varphi}_{m}\left(x\right)=\frac{1}{{h}^{7}}\{\begin{array}{ll}{\left(x-{x}_{m-4}\right)}^{7}\hfill & \left[{x}_{m-4},{x}_{m-3}\right]\hfill \\ {\left(x-{x}_{m-4}\right)}^{7}-8{\left(x-{x}_{m-3}\right)}^{7}\hfill & \left[{x}_{m-3},{x}_{m-2}\right]\hfill \\ {\left(x-{x}_{m-4}\right)}^{7}-8{\left(x-{x}_{m-3}\right)}^{7}+28{\left(x-{x}_{m-2}\right)}^{7}\hfill & \left[{x}_{m-2},{x}_{m-1}\right]\hfill \\ {\left(x-{x}_{m-4}\right)}^{7}-8{\left(x-{x}_{m-3}\right)}^{7}+28{\left(x-{x}_{m-2}\right)}^{7}-56{\left(x-{x}_{m-1}\right)}^{7}\hfill & \left[{x}_{m-1},{x}_{m}\right]\hfill \\ {\left({x}_{m+4}-x\right)}^{7}-8{\left({x}_{m+3}-x\right)}^{7}+28{\left({x}_{m+2}-x\right)}^{7}-56{\left({x}_{m+1}-x\right)}^{7}\hfill & \left[{x}_{m},{x}_{m+1}\right]\hfill \\ {\left({x}_{m+4}-x\right)}^{7}-8{\left({x}_{m+3}-x\right)}^{7}+28{\left({x}_{m+2}-x\right)}^{7}\hfill & \left[{x}_{m+1},{x}_{m+2}\right]\hfill \\ {\left({x}_{m+4}-x\right)}^{7}-8{\left({x}_{m+3}-x\right)}^{7}\hfill & \left[{x}_{m+2},{x}_{m+3}\right]\hfill \\ {\left({x}_{m+4}-x\right)}^{7}\hfill & \left[{x}_{m+3},{x}_{m+4}\right]\hfill \\ 0\hfill & \text{otherwise}\hfill \end{array}$ (5)

The set of septic B-spline functions $\left\{{\varphi}_{-3}\left(x\right)\mathrm{,}{\varphi}_{-2}\left(x\right)\mathrm{,}\cdots \mathrm{,}{\varphi}_{N+2}\left(x\right)\mathrm{,}{\varphi}_{N+3}\left(x\right)\right\}$ forms a basis for the problem region of solution $\left[a\mathrm{,}b\right]$. The approximate solution ${U}_{N}\left(x\mathrm{,}t\right)$ to the exact solution $U\left(x\mathrm{,}t\right)$ in the form:

${U}_{N}\left(x,t\right)={\displaystyle \underset{m=-3}{\overset{N+3}{\sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{\varphi}_{m}\left(x\right){\delta}_{m}\left(t\right)$ (6)

where ${\varphi}_{m}\left(x\right)$ are septic B-splines and ${\delta}_{m}\left(t\right)$ are time dependent parameters to be identified from the boundary and collocation conditions. A characteristic finite interval $\left[{x}_{m}\mathrm{,}{x}_{m+1}\right]$ is turn into the interval $\left[\mathrm{0,1}\right]$ by a domestic coordinate conversion described by $h\xi =x-{x}_{m}$, $0\le \xi \le 1$. So septic B-splines (5) in terms of $\xi $ over $\left[\mathrm{0,1}\right]$ can be written as

$\begin{array}{l}{\varphi}_{m-3}=1-7\xi +21{\xi}^{2}-35{\xi}^{3}+35{\xi}^{4}-21{\xi}^{5}+7{\xi}^{6}-{\xi}^{7},\\ {\varphi}_{m-2}=120-392\xi +504{\xi}^{2}-280{\xi}^{3}+84{\xi}^{5}-42{\xi}^{6}+7{\xi}^{7},\\ {\varphi}_{m-1}=1191-1715\xi +315{\xi}^{2}+665{\xi}^{3}-315{\xi}^{4}-105{\xi}^{5}+105{\xi}^{6}-21{\xi}^{7},\\ {\varphi}_{m}=2416-1680\xi +560{\xi}^{4}-140{\xi}^{6}+35{\xi}^{7},\end{array}$

$\begin{array}{l}{\varphi}_{m+1}=1191+1715\xi +315{\xi}^{2}-665{\xi}^{3}-315{\xi}^{4}+105{\xi}^{5}+105{\xi}^{6}-35{\xi}^{7},\\ {\varphi}_{m+2}=120+392\xi +504{\xi}^{2}+280{\xi}^{3}-84{\xi}^{5}-42{\xi}^{6}+21{\xi}^{7},\\ {\varphi}_{m+3}=1+7\xi +21{\xi}^{2}+35{\xi}^{3}+35{\xi}^{4}+21{\xi}^{5}+7{\xi}^{6}-{\xi}^{7},\\ {\varphi}_{m+4}={\xi}^{7}.\end{array}$ (7)

Using Equation (5) and Equation (6), the nodal values of ${U}_{m}\mathrm{,}{{U}^{\prime}}_{m}\mathrm{,}{{U}^{\u2033}}_{m}\mathrm{,}{{U}^{\u2034}}_{m}$ and ${U}_{m}^{iv}$ at the knots ${x}_{m}$ are obtained as the following:

$\begin{array}{l}{U}_{N}\left({x}_{m},t\right)={U}_{m}={\delta}_{m-3}+120{\delta}_{m-2}+1191{\delta}_{m-1}+2416{\delta}_{m}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.05em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}+1191{\delta}_{m+1}+120{\delta}_{m+2}+{\delta}_{m+3},\\ {{U}^{\prime}}_{m}=\frac{7}{h}\left(-{\delta}_{m-3}-56{\delta}_{m-2}-245{\delta}_{m-1}+245{\delta}_{m+1}+56{\delta}_{m+2}+{\delta}_{m+3}\right),\\ {{U}^{\u2033}}_{m}=\frac{42}{{h}^{2}}\left({\delta}_{m-3}+24{\delta}_{m-2}+15{\delta}_{m-1}-80{\delta}_{m}+15{\delta}_{m+1}+24{\delta}_{m+2}+{\delta}_{m+3}\right),\\ {{U}^{\u2034}}_{m}=\frac{210}{{h}^{3}}\left(-{\delta}_{m-3}-8{\delta}_{m-2}+19{\delta}_{m-1}-19{\delta}_{m+1}+8{\delta}_{m+2}+{\delta}_{m+3}\right),\\ {U}_{m}^{iv}=\frac{840}{{h}^{4}}\left({\delta}_{m-3}-9{\delta}_{m-1}+16{\delta}_{m}-9{\delta}_{m+1}+{\delta}_{m+3}\right)\end{array}$ (8)

where the symbols $\text{'},\text{'}\text{'},\text{'}\text{'}\text{'}$ symbolize differentiation according to x, respectively. Using (5) and (8) in the Equation (1) this guides to a set of ordinary differential equations of the form

$\begin{array}{l}{\stackrel{\dot{}}{\delta}}_{m-3}+120{\stackrel{\dot{}}{\delta}}_{m-2}+1191{\stackrel{\dot{}}{\delta}}_{m-1}+2416{\stackrel{\dot{}}{\delta}}_{m}+1191{\stackrel{\dot{}}{\delta}}_{m+1}+120{\stackrel{\dot{}}{\delta}}_{m+2}+{\stackrel{\dot{}}{\delta}}_{m+3}\\ +\frac{7\epsilon {Z}_{m}}{h}\left(-{\delta}_{m-3}-56{\delta}_{m-2}-245{\delta}_{m-1}+245{\delta}_{m+1}+56{\delta}_{m+2}+{\delta}_{m+3}\right)\\ +\frac{210\mu}{{h}^{3}}\left(-{\delta}_{m-3}-8{\delta}_{m-2}+19{\delta}_{m-1}-19{\delta}_{m+1}+8{\delta}_{m+2}+{\delta}_{m+3}\right)=0,\end{array}$ (9)

where

${Z}_{m}={\left({\delta}_{m-3}+120{\delta}_{m-2}+1191{\delta}_{m-1}+2416{\delta}_{m}+1191{\delta}_{m+1}+120{\delta}_{m+2}+{\delta}_{m+3}\right)}^{2}\mathrm{.}$

If time parameters ${\delta}_{i}$ and its time derivatives ${\stackrel{\dot{}}{\delta}}_{i}$ in Equation (9) are separated by the Crank-Nicolson form and finite difference approach, respectively:

${\delta}_{i}=\frac{{\delta}_{i}^{n+1}+{\delta}_{i}^{n}}{2}\mathrm{,}$ (10)

and usual finite difference approximation

${\stackrel{\dot{}}{\delta}}_{i}=\frac{{\delta}_{i}^{n+1}-{\delta}_{i}^{n}}{\Delta t}$ (11)

we acquired a repetition relationship between two time levels n and $n+1$ relating two unknown parameters ${\delta}_{i}^{n+1}\mathrm{,}{\delta}_{i}^{n}$ for $i=m-3,m-2,\cdots ,m+2,m+3$

$\begin{array}{l}{\gamma}_{1}{\delta}_{m-3}^{n+1}+{\gamma}_{2}{\delta}_{m-2}^{n+1}+{\gamma}_{3}{\delta}_{m-1}^{n+1}+{\gamma}_{4}{\delta}_{m}^{n+1}+{\gamma}_{5}{\delta}_{m+1}^{n+1}+{\gamma}_{6}{\delta}_{m+2}^{n+1}+{\gamma}_{7}{\delta}_{m+3}^{n+1}\\ ={\gamma}_{7}{\delta}_{m-3}^{n}+{\gamma}_{6}{\delta}_{m-2}^{n}+{\gamma}_{5}{\delta}_{m-1}^{n}+{\gamma}_{4}{\delta}_{m}^{n}+{\gamma}_{3}{\delta}_{m+1}^{n}+{\gamma}_{2}{\delta}_{m+2}^{n}+{\gamma}_{1}{\delta}_{m+3}^{n},\end{array}$ (12)

where

$\begin{array}{l}{\gamma}_{1}=\left[1-E{Z}_{m}-M\right],\\ {\gamma}_{2}=\left[120-56E{Z}_{m}-8M\right],\\ {\gamma}_{3}=\left[1191-245E{Z}_{m}+19M\right],\\ {\gamma}_{4}=\left[2416\right],\end{array}$

$\begin{array}{l}{\gamma}_{5}=\left[1191+245E{Z}_{m}-19M\right],\\ {\gamma}_{6}=\left[120+56E{Z}_{m}+8M\right],\\ {\gamma}_{7}=\left[1+E{Z}_{m}+M\right],\\ m=0,1,\cdots ,N,E=\frac{7\epsilon}{2h}\Delta t,M=\frac{105\mu}{{h}^{3}}\Delta t.\end{array}$ (13)

The system (12) contains of $\left(N+1\right)$ linear equations containing $\left(N+7\right)$ unknown coefficients ${\left({\delta}_{-3}\mathrm{,}{\delta}_{-2}\mathrm{,}{\delta}_{-1}\mathrm{,}\cdots \mathrm{,}{\delta}_{N+1}\mathrm{,}{\delta}_{N+2}\mathrm{,}{\delta}_{N+3}\right)}^{\text{T}}$. To acquire a solution of this system, we require six additional restrictions. These are obtained from the boundary conditions (2) and can be used to remove ${\delta}_{-3}\mathrm{,}{\delta}_{-2}\mathrm{,}{\delta}_{-1}$ and ${\delta}_{N+1}\mathrm{,}{\delta}_{N+2}\mathrm{,}{\delta}_{N+3}$ from the systems (12) which occurs a matrix equation for the $N+1$ unknowns ${d}^{n}={\left({\delta}_{0}\mathrm{,}{\delta}_{1}\mathrm{,}\cdots \mathrm{,}{\delta}_{N}\right)}^{\text{T}}$ of the form

$A{d}^{n+1}=B{d}^{n}.$ (14)

The resulting system is effectively solved with a version of the Thomas algorithm and we implement an inner iteration ${\delta}^{n\ast}={\delta}^{n}+\frac{1}{2}\left({\delta}^{n}-{\delta}^{n-1}\right)$ at each time step to overcome the non-linearity caused by ${Z}_{m}$. Before the beginning of the solution procedure, initial parameters ${d}^{0}$ are established by using the initial condition and following derivatives at the boundaries;

${U}_{N}\left(x,0\right)=U\left({x}_{m},0\right);\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}m=0,1,2,\cdots ,N$ (15)

${\left({U}_{N}\right)}_{x}\left(a,0\right)=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\left({U}_{N}\right)}_{x}\left(b,0\right)=0,$ (16)

${\left({U}_{N}\right)}_{xx}\left(a,0\right)=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\left({U}_{N}\right)}_{xx}\left(b,0\right)=0,$ (17)

${\left({U}_{N}\right)}_{xxx}\left(a,0\right)=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\left({U}_{N}\right)}_{xxx}\left(b,0\right)=0.$ (18)

So, by taking account (18), we obtain the following matrix form for the initial vector ${d}^{0}$ ;

$W{d}^{0}=b,$

where

$W=\left[\begin{array}{cccccccccc}1536& 2712& 768& 24& & & & & & \\ \frac{82731}{81}& \frac{210568.5}{81}& \frac{104796}{81}& \frac{10063.5}{81}& 1& & & & & \\ \frac{9600}{81}& \frac{96597}{81}& \frac{195768}{81}& \frac{96474}{81}& 120& 1& & & & \\ & & \ddots & & & & & & & \\ & & & & & & & & & \\ & & & 1& 120& 1191& 2416& 1191& 120& 1\\ & & & & 1& 120& \frac{96474}{81}& \frac{195768}{81}& \frac{96597}{81}& \frac{9600}{81}\\ & & & & & 1& \frac{10063.5}{81}& \frac{104796}{81}& \frac{210568.5}{81}& \frac{82731}{81}\\ & & & & & & 24& 768& 2712& 1536\end{array}\right]$

${d}^{0}={\left({\delta}_{0}\mathrm{,}{\delta}_{1}\mathrm{,}{\delta}_{2}\mathrm{,}\cdots \mathrm{,}{\delta}_{N-2}\mathrm{,}{\delta}_{N-1}\mathrm{,}{\delta}_{N}\right)}^{\text{T}}$

and

$b={\left(U\left({x}_{0},0\right),U\left({x}_{1},0\right),\cdots ,U\left({x}_{N-1},0\right),U\left({x}_{N},0\right)\right)}^{\text{T}}$.

3. Stability Analysis

To implement the von Neumann stability analysis, GKdV equation is linearized by thinking about the quantity ${U}^{p}$ in the nonlinear term ${U}^{p}{U}_{x}$ is locally invariable. Substituting the Fourier mode ${\delta}_{m}^{n}={\xi}^{n}{\text{e}}^{imkh}$, ( $i=\sqrt{-1}$ ) in which k is a mode number and h is the element size, into the Equation (12) gives the growth factor $\xi $ of the form

$\xi =\frac{A-iB}{A+iB}$ (19)

where

$\begin{array}{l}A=2\mathrm{cos}\left(3kh\right)+240\mathrm{cos}\left(2kh\right)+2382\mathrm{cos}\left(kh\right)+2416\\ B=2\left(E{Z}_{m}+M\right)\mathrm{sin}\left(3kh\right)+2\left(56E{Z}_{m}+8M\right)\mathrm{sin}\left(2kh\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}+2\left(245E{Z}_{m}-19M\right)\mathrm{sin}\left(kh\right)\end{array}$ (20)

The modulus of the (19) is found 1, hence the linearized algorithm is unconditionally stable.

4. Test Problems

In this section, we introduce some numerical examples including: motion of single soliton wave whose exact solution is known to test validity of our algorithm for solving GKdV equation. The initial boundary value problem (1)-(2) possesses following conservative quantities;

$\begin{array}{l}{I}_{1}={\displaystyle {\int}_{-\infty}^{\infty}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}U\left(x,t\right)\text{d}x,\\ {I}_{2}={\displaystyle {\int}_{-\infty}^{\infty}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{U}^{2}\left(x,t\right)\text{d}x,\\ {I}_{3}={\displaystyle {\int}_{-\infty}^{\infty}}\left[{U}^{p+2}\left(x,t\right)-\frac{\mu \left(p+1\right)\left(p+2\right)}{2\epsilon}{U}_{x}^{2}\left(x,t\right)\right]\text{d}x\end{array}$ (21)

which correspond to the mass, momentum and energy of the shallow water waves, respectively [33] [34]. To compare the numerical solution with the exact solution we use the following error norms:

${L}_{2}={\Vert {U}^{exact}-{U}_{N}\Vert}_{2}\simeq \sqrt{h{\displaystyle \underset{j=0}{\overset{N}{\sum}}}{\left|{U}_{j}^{exact}-{\left({U}_{N}\right)}_{j}\right|}^{2}}\mathrm{,}$

and

${L}_{\infty}={\Vert {U}^{exact}-{U}_{N}\Vert}_{\infty}\simeq \underset{j}{\mathrm{max}}\left|{U}_{j}^{exact}-{\left({U}_{N}\right)}_{j}\right|\mathrm{.}$

The Motion of Single Solitary Wave

For this test problem, Equation (1) is examined with the boundary conditions $U\to 0$ as $x\to \pm \infty $ and the initial condition $U\left(x,0\right)=A{\text{sech}}^{\frac{2}{p}}\left[k\left(x-{x}_{0}\right)\right]$ where $A={\left[\frac{c\left(p+1\right)\left(p+2\right)}{2\epsilon}\right]}^{\frac{1}{p}}$ is amplitude and $k=\frac{p}{2}\sqrt{\frac{c}{\mu}}$ is width of the single soliton. The exact soliton solution of the GKdV equation is

$U\left(x,t\right)=A{\text{sech}}^{\frac{2}{p}}\left[k\left(x-{x}_{0}-ct\right)\right],$

where c and ${x}_{0}$ are arbitrary constants. In order to exemplify the validity of our numerical algorithm, we conceive the first case of a single soliton solution for the parameters $p=1$, $\epsilon =1$, $\mu =4.84\times {10}^{-4}$, $h=0.01$, $\Delta t=0.005$, $c=0.3$ and $\epsilon =3$, $\mu =1$, $h=0.1$, $\Delta t=0.01$ through the interval $\left[\mathrm{0,80}\right]$ to compare with that of previous papers [8] [10] [14] [16] [30]. All parameters are given in all refarans. For these parameters, the single solitary wave has the amplitude 0.9 and 0.3, respectively. The three invariants ${I}_{1}\mathrm{,}{I}_{2}$ and ${I}_{3}$ together with the ${L}_{2}$, ${L}_{\infty}$ error norms for the problem are documented and compared in Table 1 for times up to $t=1$. As seen from the table that ${L}_{2}$ and ${L}_{\infty}$ error norms are found small enough and the conservation of the invariants can be seen to be almost constant.

Solitary wave profiles are demonstrated at $t=0,0.1,0.2,\cdots ,1$ in Figure 1 in which the soliton moves to the right at a nearly unchanged speed and amplitude as time increases, as expected.

Table 1. Comparisons of results for invariants and error norms with $p=1$, $\epsilon =1$, $\mu =4.84\times {10}^{-4}$, $c=0.3$, $h=0.01$, $\Delta t=0.005$, $x\in \left[\mathrm{0,2}\right]$ and $\epsilon =3$, $\mu =1$, $c=0.3$, $h=0.1$, $\Delta t=0.01$.

(a) (b)

Figure 1. Motion of single solitary wave for (a) $p=1$, $\epsilon =1$, $\mu =4.84\times {10}^{-4}$, $c=0.3$, $h=0.01$, $\Delta t=0.005$ and (b) $\epsilon =3$, $\mu =1$, $c=0.3$, $h=0.1$, $\Delta t=0.01$.

For the second set, we choose the parameters $p=2$, $\epsilon =3$, $\mu =1$, $h=0.1$, $\Delta t=0.01$, $c=0.845$ and 0.3 throughout the interval $\left[\mathrm{0,80}\right]$ just to be able to compare them with earlier papers [19] [23] [24]. These values yield the amplitude 1.3416 and 0.7746 and the computations are done until time $t=20$ and $t=1$. We calculate the values of the error norms ${L}_{2}$ and ${L}_{\infty}$ and invariants ${I}_{1},{I}_{2}$ and ${I}_{3}$ for different time levels and compare them with earlier papers in Table 2. This table indicates that the error norms obtained by our method are found much better than the others and the calculated values of invariants are in good conformity with the others. The motion of solitary wave using our scheme are plotted at times $t=0,5,10,15,20$ and $t=0,0.1,0.2,\cdots ,1$ in Figure 2 in which the soliton moves to the right at a nearly unchanged speed and amplitude as time increases, as expected.

Finally, we have taken the parameters $p=3$, $\epsilon =1$, $\mu =4.84\times {10}^{-4}$, $h=0.01$, $\Delta t=0.005$, $c=0.3$ and $\epsilon =6$, $\mu =1$, $h=0.1$, $\Delta t=0.01$, $c=0.6$ over the region $\left[\mathrm{0,80}\right]$. Thereby, solitary wave has amplitude 1.44 and 1.0, respectively. Simulations are executed to time $t=1$ to invent the error norms ${L}_{2}$ and ${L}_{\infty}$ and the numerical invariants ${I}_{1}\mathrm{,}{I}_{2}$ and ${I}_{3}$. The calculated values are presented in Table 3. As can be seen in Table 3, the error norms ${L}_{2}$ and ${L}_{\infty}$ are sensibly small and the quantities of the invariants remain almost constant during the computer run. The behaviors of solutions for values of $h=t=0.125$ at times $t=0,20$ and 40 are depicted in Figure 3.

5. Conclusion

In this paper, a septic B-spline collocation method has been successfully applied to the GKdV equation to examine the motion of a single solitary wave whose analytical solution is known. To show how good and accurate the numerical solutions of the test problems, we have computed the error norms ${L}_{2}$ and ${L}_{\infty}$ and conserved quantities ${I}_{1}\mathrm{,}{I}_{2}$ and ${I}_{3}$. According to the tables in the paper, one can have easily seen that our error norms are enough small and the invariants are well conserved. Stability analysis has been done and the linearized numerical scheme has been obtained unconditionally stable. Thus, we can say

(a) (b)

Figure 2. Motion of single solitary wave for (a) $p=2$, $\epsilon =3$, $\mu =1$, $h=0.1$, $\Delta t=0.01$, $c=0.845$ and (b) $c=0.3$, $h=0.1$, $\Delta t=0.01$.

(a) (b)

Figure 3. Motion of single solitary wave for (a) $p=3$, $\epsilon =3$, $\mu =1$, $h=0.01$, $\Delta t=0.005$, $c=0.845$ and (b) $c=0.3$, $h=0.1$, $\Delta t=0.01$.

Table 2. Comparisons of results for invariants and error norms with $p=2$, $\epsilon =3$, $\mu =1$, $h=0.1$, $\Delta t=0.01$, $c=0.845$ and $c=0.3$, $h=0.1$, $\Delta t=0.01$.

Table 3. Values of the invariants and error norms for $p=3$, $\epsilon =3$, $\mu =1$, $h=0.01$, $\Delta t=0.005$, $c=0.845$ and $c=0.3$, $h=0.1$, $\Delta t=0.01$.

that our numerical scheme is useful to obtain the numerical solutions of other important nonlinear problems in various fields.

References

[1] Korteweg, D.J. and de Vries, G. (1895) On the Change of Form of Long Waves Advancing in a Rectangular Canal and on a New Type of Long Stationary Waves. The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, 39, 422-443.

https://doi.org/10.1080/14786449508620739

[2] Zabusky, N.J. (1967) A Synergetic Approach to Problem of Nonlinear Dispersive Wave Propagation and Interaction. In: Ames, W., Ed., Proceedings of the Symposium Nonlinear Partial Differential Equation, Academic Press, 223-258.

[3] Fornberg, B. and Whitham, G.B. (1978) A Numerical and Theoretical Study of Certain Nonlinear Wave Phenomena. Philosophical Transactions of the Royal Society A, 289, 373-404.

https://doi.org/10.1098/rsta.1978.0064

[4] Zabusky, N.J. and Kruskal, M.D. (1965) Interaction of Solitons in a Collisionless plasma and the Recurrence of Initial States. Physical Review Letters, 15, 240-243.

https://doi.org/10.1103/PhysRevLett.15.240

[5] Gardner, C.S., Green, J.M., Kruskal, M.D. and Miura, R.M. (1967) Method for Solving Korteweg-de Vries Equation. Physical Review Letters, 19, 1095-1097.

https://doi.org/10.1103/PhysRevLett.19.1095

[6] Goda, K. (1975) On Instability of Some Finite Difference Schemes for Korteweg-de Vries Equation. Journal of the Physical Society of Japan, 39, 229-236.

https://doi.org/10.1143/JPSJ.39.229

[7] Vliengenthart, A.C. (1971) On Finite Difference Methods for the Korteweg-de Vries Equation. Journal of Engineering Mathematics, 5, 137-155.

https://doi.org/10.1007/BF01535405

[8] Soliman, A.A. (2004) Collocation Solution of the Korteweg-de Vries Equation Using Septic Splines. International Journal of Computer Mathematics, 81, 325-331.

https://doi.org/10.1080/00207160410001660817

[9] Irk, D., Dağ, I. and Saka, B. (2006) A Small Time Solutions for the Korteweg--de Vries Equation Using Spline Approximation. Applied Mathematics and Computation, 173, 834-846.

https://doi.org/10.1016/j.amc.2005.04.018

[10] Canıvar, A., Sari, M. and Dağ, I. (2010) A Taylor-Galerkin Finite Element Method for the KdV Equation Using Cubic B-Splines. Physica B: Condensed Matter, 405, 3376-3383.

https://doi.org/10.1016/j.physb.2010.05.008

[11] Aksan, E.N. and Özdeş, A. (2006) Numerical Solution of Korteweg-de Vries Equation by Galerkin B-Spline Finite Element Method. Applied Mathematics and Computation, 175, 1256-1265.

https://doi.org/10.1016/j.amc.2005.08.038

[12] Irk, D. (2017) Quintic B-Spline Galerkin Method for the KdV Equation. Anadolu University Journal of Science and Technology B-Theoritical Sciences, 5, 111-119.

[13] Ersoy, O. and Dağ, I. (2015) The Exponential Cubic B-Spline Algorithm for Korteweg-de Vries Equation. Advances in Numerical Analysis, 2015, Article ID: 367056.

https://doi.org/10.1155/2015/367056

[14] Dağ, I. and Dereli, Y. (2008) Numerical Solutions of KdV Equation Using Radial Basis Function. Applied Mathematical Modelling, 32, 535-546.

https://doi.org/10.1016/j.apm.2007.02.001

[15] Kutluay, S., Bahadir, A.R. and Özdeş, A. (2000) A Small Time Solutions for the Korteweg-de Vries Equation. Applied Mathematics and Computation, 107, 203-210.

https://doi.org/10.1016/S0096-3003(98)10119-4

[16] Saka, B. (2009) Cosine Expansion-Based Differential Quadrature Method for Numerical Solution of the KdV Equation. Chaos, Solitons and Fractals, 40, 2181-2190.

https://doi.org/10.1016/j.chaos.2007.10.004

[17] Kaya, D. (2005) An Application for the Higher Order Modified KdV Equation by Decomposition Method. Communications in Nonlinear Science and Numerical Simulation, 10, 693-702.

https://doi.org/10.1016/j.cnsns.2003.12.009

[18] Biswas, A. (2010) 1-Soliton Solution of the K(m, n) Equation with Generalized Evolution and Time-Dependent Damping and Dispersion. Computers and Mathematics with Applications, 59, 2536-2540.

[19] Raslan, K.R. and Baghdady, H.A. (2015) A Finite Difference Scheme for the Modified Korteweg-de Vries Equation. General Mathematics Notes, 27, 101-113.

[20] Raslan, K.R. and Baghdady, H.A. (2014) New Algorithm for Solving the Modified Korteweg-de Vries (mKdV) Equation. International Journal of Research and Reviews in Applied Sciences, 18, 59-64.

[21] Wazwaz, A.M. (2014) A Variety of (3+1)-Dimensional mKdV Equations Derived by Using the mKdV Recursion Operator. Computers and Fluids, 93, 41-45.

https://doi.org/10.1016/j.compfluid.2014.01.010

[22] Wazwaz, A.M. (2015) New (3+1)-Dimensional Nonlinear Evolution Equations with mKdV Equation Constituting Its Main Part: Multiple Soliton Solutions. Chaos, Solitons and Fractals, 76, 93-97.

https://doi.org/10.1016/j.chaos.2015.03.018

[23] Ak, T., Karakoc, S.B.G. and Biswas, A. (2017) Application of Petrov-Galerkin Finite Element Method to Shallow Water Waves Model: Modified Korteweg-de Vries Equation. Scientia Iranica B, 24, 1148-1159.

https://doi.org/10.24200/sci.2017.4096

[24] Ak, T., Karakoc, S.B.G. and Biswas, A. (2017) A New Approach for Numerical Solution of Modified Korteweg-de Vries Equation. Iranian Journal of Science and Technology, Transactions A: Science, 41, 1109-1121.

https://doi.org/10.1007/s40995-017-0238-5

[25] Bracken, P. (2005) Specific Solutions of the Generalized Korteweg-de Vries Equation with Possible Physical Applications. Central European Journal of Physics, 3, 127-138.

https://doi.org/10.2478/BF02476511

[26] H. Liu and Yi, N. (2016) A Hamiltonian Preserving Discontinuous Galerkin Method for the Generalized Korteweg-de Vries Equation. Journal of Computational Physics, 321, 776-796.

https://doi.org/10.1016/j.jcp.2016.06.010

[27] Lai, X., Cao, Q. and Twizell, E.H. (2005) The Global Domain of Attraction and the Initial Value Problems of a Kind of GKdV Equations. Chaos, Solitons and Fractals, 23, 1613-1628.

https://doi.org/10.1016/S0960-0779(04)00413-8

[28] Garcia Alvarado, M.G. and Omel’yanov, G.A. (2012) Interaction of Solitary Waves for the Generalized KdV Equation. Communications in Nonlinear Science and Numerical Simulation, 17, 3204-3218.

https://doi.org/10.1016/j.cnsns.2011.12.001

[29] Sepulveda, M. and Villagran, O.V. (2008) Numerical Methods for Generalized KdV Equations. In: Anais do XXXI Congresso Nacional de Matematica Aplicada e Computacional, 894-900.

[30] Ak, T., Triki, H., Dhawan, S., Bhowmik, S.K., Moshokoae, S.P., Ullah, M.Z. and Biswas, A. (2017) Computational Analysis of Shallow Water Waves with Korteweg-de Vries Equation. Scientia Iranica B, 25, 2582-2597.

[31] Ismail, H.N.A., Raslan, K.R. and Salem, G.S.E. (2004) Solitary Wave Solutions for the General KDV Equation by Adomian Decomposition Method. Applied Mathematics and Computation, 154, 17-29.

https://doi.org/10.1016/S0096-3003(03)00686-6

[32] Prenter, P.M. (1975) Splines and Variational Methods. John Wiley & Sons, New York.

[33] Miura, R.M., Gardner, C.S. and Kruskal, M.D. (1968) Korteweg-de Vries Equation and Generalizations. II. Existence of Conservation Laws and Constants of Motion. Journal of Mathematical Physics, 9, 1204.

https://doi.org/10.1063/1.1664701

[34] Miura, R.M. (1976) The Korteweg-de Vries Equation: A Survey of Results. SIAM Review, 18, 412-459.

https://doi.org/10.1137/1018076