# Fourier Splitting Method for Kawahara Type Equations

Authors: Pablo U. Suárez; J. Héctor Morales
Journal: Journal of Computational Methods in Physics
 (2014)

---

## Abstract

In this work, we integrate numerically the Kawahara and generalized Kawahara equation by using an algorithm based on Strang’s
splitting method. The linear part is solved using the Fourier transform and the nonlinear part is solved with the aid of the exponential
operator method. To assess the accuracy of the solution, we compare
known analytical solutions with the numerical solution. Further, we
show that ast increases the conserved quantities remain constant.


---

## Body

## 1. Introduction

The Kawahara equation is a generalized nonlinear dispersive equation which has a form of the Korteweg-de Vries (KdV) equation with an additional fifth-order derivative term [1]. This equation is an evolution partial differential equation (PDE) that describes in one spatial dimension the propagation of shallow water waves with surface tension and magnetoacoustic waves in plasma [1]. Analytical solutions have been investigated [2]; however, no general solutions have been found yet. The inverse scattering transform (IST), a general method used to integrate nonlinear evolution PDEs [3], cannot be applied to this equation because the Painlevé test of integrability fails [4–6]. In this work, we will employ the following generalized version of the Kawahara equation:
(1)
u
t
+
μ
u
x
+
γ
u
n
u
x
+
p
u
3
x
-
q
u
5
x
=
0
,
u
(
x
,
0
)
=
f
(
x
)
,
x
∈
R
,
t
∈
[
0
,
T
]
,
where the rescaled nondimensional quantities μ, γ, p, and q depend on the physical parameters, n is a integer, and the case for n
=
1 corresponds to the magneto-acoustic wave and n
=
2 to the shallow water approximations. When μ
=
0, γ
=
6, n
=
1, p
=
1, and q
=
0, we reduce to a classical representation of the KdV equation.In order to study the equation by means of numerical methods and produce test of such methods, we will make use of the traveling wave ansatz to get a solitary wave (SW) solution to the Kawahara equation. We will adopt the SW solution to be our exact solution to compare with our numerical approximations. Finite differences, finite elements, and radial basis function have been used to get numerical integration of the Kawahara equation [7–9]. The present method has high accuracy; it is fast and easy to implement as seen in other works that deal with nonlinearity—see, for instance, [10]. In this work we employ a split-step (Fourier) method to integrate the equation. The splitting algorithm is a pseudospectral numerical method which consists in approximate iteratively the solution of the nonlinear evolution equation using exponential operators. We advance in time in small steps and treat separately the linear and nonlinear parts. The nonlinear term is handled by using the matrix exponential method, this method was used in [10–12] to solve similar PDEs. To handle the linear part, we take the Fourier and inverse Fourier transform at each time step. This method is explained in Sections 2 and 3, respectively. We demonstrate the implementation of the numerical scheme and we make some comparisons with the exact solutions given in [4]. Finally, in the last section, we make some concluding remarks about the method.


## 2. Numerical Scheme

We consider the generalized Kawahara equation(2)
u
t
+
μ
u
x
+
p
u
3
x
-
q
u
5
x
+
γ
u
n
u
x
=
0
.
In order to implement a splitting method, we rewrite (2) as follows:
(3)
u
t
=
-
(
L
+
N
)
u
,
where the differential operators L and N are the following linear and nonlinear operators, respectively,
(4)
L
u
:
=
μ
u
x
+
p
u
3
x
-
q
u
5
x
,
N
u
:
=
γ
u
n
u
x
.The Kawahara equation can be split into a linear part,(5)
(
u
t
)
L
=
-
(
μ
∂
u
∂
x
+
p
∂
3
u
∂
x
3
-
q
∂
5
u
∂
x
5
)
,
and a nonlinear part,
(6)
(
u
t
)
N
=
-
γ
u
n
-
1
∂
u
∂
x
u
.
When n
=
1, both (5) and (6) have analytical solutions [6]. For instance, a plane wave solution to (5) is given by u
(
x
,
t
)
=
exp
⁡
(
i
(
k
x
-
ω
(
k
)
t
)
), where ω
(
k
)
=
μ
k
+
p
k
3
+
q
k
5 is the dispersion relation. A solution to the nonlinear equation (6) is given by the similarity solution u
(
x
,
t
)
=
x
/
γ
t. However, the Kawahara equation contains both parts and in general does not have an analytical solution. Therefore, we use the traveling wave ansatz [4] to obtain a particular solution to (2).


### 2.1. Nonlinear Part

We first consider the nonlinear operator(7)
N
u
=
u
n
-
1
∂
u
∂
x
this nonlinear operator can also be represented in Fourier space as
(8)
N
u
=
u
n
-
1
F
-
1
(
i
k
F
(
u
)
)
.
We choose this form of the operator because it avoids the harsher numerical error introduced by finite differences. Note that nonlinear operator is applied by simple multiplication.Now, if we consider only a small stepΔ
t along t for the Kawahara equation, one can treat the operators as sequential steps with only a “small” numerical error. To advance in a such “small” nonlinear step, we use the analytical solution of (6):
(9)
u
(
t
+
Δ
t
,
x
)
=
exp
⁡
[
-
γ
u
n
-
1
Δ
t
F
-
1
(
i
k
F
(
u
^
(
t
,
ω
)
)
)
]
u
(
t
,
x
)
.


### 2.2. Linear Part

Next, we advance a linear step by using the analytical solution for (5). The linear equation can be solved by taking the Fourier transform along the x; thus, the solution is
(10)
u
(
t
+
Δ
t
,
x
)
=
F
-
1
(
exp
⁡
[
-
Δ
t
(
i
k
μ
+
p
i
k
3
+
q
i
k
5
)
]
u
^
(
t
,
ω
)
)
.
The Fourier transform is approximated accurately and quickly by the use of the well-known fast Fourier transform (FFT).


### 2.3. Second Order Scheme

In general, an evolution problem can be approximated as [1]
(11)
u
(
t
+
Δ
t
,
x
)
=
exp
⁡
(
Δ
t
(
L
+
N
)
)
u
(
t
,
x
)
iiiiiiiiiiiiiiiiii
≈
exp
⁡
(
Δ
t
L
)
exp
⁡
(
Δ
t
N
)
u
(
t
,
x
)
.
The left hand of the equation is only an approximation because in general the operators L and N do not commute.The above method is only first order accurate. To improve the accuracy, we apply Strang’s splitting scheme [11]:
(12)
u
(
t
+
Δ
t
,
x
)
≈
exp
⁡
(
Δ
t
2
N
)
exp
⁡
(
Δ
t
L
)
exp
⁡
(
Δ
t
2
N
)
u
(
t
,
x
)
.
Using this method, we obtained a second accurate scheme.


## 2.1. Nonlinear Part

We first consider the nonlinear operator(7)
N
u
=
u
n
-
1
∂
u
∂
x
this nonlinear operator can also be represented in Fourier space as
(8)
N
u
=
u
n
-
1
F
-
1
(
i
k
F
(
u
)
)
.
We choose this form of the operator because it avoids the harsher numerical error introduced by finite differences. Note that nonlinear operator is applied by simple multiplication.Now, if we consider only a small stepΔ
t along t for the Kawahara equation, one can treat the operators as sequential steps with only a “small” numerical error. To advance in a such “small” nonlinear step, we use the analytical solution of (6):
(9)
u
(
t
+
Δ
t
,
x
)
=
exp
⁡
[
-
γ
u
n
-
1
Δ
t
F
-
1
(
i
k
F
(
u
^
(
t
,
ω
)
)
)
]
u
(
t
,
x
)
.


## 2.2. Linear Part

Next, we advance a linear step by using the analytical solution for (5). The linear equation can be solved by taking the Fourier transform along the x; thus, the solution is
(10)
u
(
t
+
Δ
t
,
x
)
=
F
-
1
(
exp
⁡
[
-
Δ
t
(
i
k
μ
+
p
i
k
3
+
q
i
k
5
)
]
u
^
(
t
,
ω
)
)
.
The Fourier transform is approximated accurately and quickly by the use of the well-known fast Fourier transform (FFT).


## 2.3. Second Order Scheme

In general, an evolution problem can be approximated as [1]
(11)
u
(
t
+
Δ
t
,
x
)
=
exp
⁡
(
Δ
t
(
L
+
N
)
)
u
(
t
,
x
)
iiiiiiiiiiiiiiiiii
≈
exp
⁡
(
Δ
t
L
)
exp
⁡
(
Δ
t
N
)
u
(
t
,
x
)
.
The left hand of the equation is only an approximation because in general the operators L and N do not commute.The above method is only first order accurate. To improve the accuracy, we apply Strang’s splitting scheme [11]:
(12)
u
(
t
+
Δ
t
,
x
)
≈
exp
⁡
(
Δ
t
2
N
)
exp
⁡
(
Δ
t
L
)
exp
⁡
(
Δ
t
2
N
)
u
(
t
,
x
)
.
Using this method, we obtained a second accurate scheme.


## 3. Numerical Examples and Comparisons

In this section, we present some numerical examples to test the accuracy of the proposed method.We can benchmark the numerical efficiency of the method by computing theL
2 and L
∞ error norms
(13)
L
2
=
Δ
t
∑
j
=
0
N
‍
|
u
j
e
-
u
j
a
|
2
,
L
∞
=
max
⁡
j
|
u
j
e
-
u
j
a
|
,
where u
e is the exact solution and u
a the approximate solution. To complete the analysis, we also calculated numerically the conserved quantities linear momentum (M) and energy (E).


### 3.1. Traveling Wave Solution of the Generalized Kawahara Equation

We show the accuracy of proposed method for generalized Kawahara equation with nonlinearityn
=
1,2 and we let γ
=
p
=
q
=
1 and μ
=
0. For computational purposes, the infinite interval is truncated to x
∈
[
-
100,100
] and t
∈
[
0,50
]. We use periodic boundary conditions to simulate the behaviour at infinity. For all these experiments, we have Δ
t
=
.
005 and the number of space points is 512.


#### 3.1.1. Kawahara Equationn
=
1

For the casen
=
1, the equation reduces to the well known Kawahara equation. An exact solution is given by
(14)
u
(
x
,
t
)
=
105
169
sec
h
4
(
1
2
13
(
x
-
36
169
t
)
)
.In Table1, we see the numerical results obtained by the Fourier splitting technique on the Kawahara equation. As t increases, the error grows. This growth can be attributed to periodic boundary conditions simulating the behaviour at infinity. This is corroborated by the fact that L
∞ is obtained at points near end of the computational domain. However, our numerical solution is reliable because the conserved quantities remain within numerical error throughout the computational domain.Table 1
Absolute errors and conserved quantities for the Kawahara equation (casen
=
1).
t
L
2
L
∞
I
1
I
2
20.0
5.1011
E
-
05
2.1776
E
-
05
5.9737
1.2726
30.0
7.7990
E
-
05
3.2469
E
-
05
5.9738
1.2726
40.0
1.1144
E
-
04
4.5564
E
-
05
5.9738
1.2726
50.0
1.5310
E
-
04
6.2665
E
-
05
5.9739
1.2727


#### 3.1.2. Kawahara Equationn
=
2

The second case we examine is the generalized Kawaharan
=
2. An exact solution is given by
(15)
u
(
x
,
t
)
=
3
10
sec
h
2
(
1
2
5
(
x
-
4
25
t
)
)
.In Table2, we see again an excellent agreement between the numerical results obtained by the Fourier splitting technique and the analytical solution. The errors increase with time but these errors are not numerically significant and the maximum value of the error is always near the boundary.Table 2
Absolute errors and conserved quantities for the Kawahara equation (casen
=
2).
t
L
2
L
∞
I
1
I
2
20.0
1.3882
E
-
04
5.9450
E
-
05
8.4854
2.6835
30.0
2.5212
E
-
04
1.0612
E
-
04
8.4855
2.6836
40.0
4.0950
E
-
04
1.6662
E
-
04
8.4855
2.6837
50.0
6.1195
E
-
04
2.4806
E
-
04
8.4856
2.6837Our method can accurately approximate the analytical solutions for the integer valuesn
>
2. We omit these results as they are not particularly illuminating.


#### 3.1.3. KdV-Kawahara Equation

For our final example, we setμ
=
1. With this additional term, we have the so-called KdV-Kawahara which has applications in magnetohydrodynamics; see [9] for details. An exact traveling wave solution to this equation is given by
(16)
u
(
x
,
t
)
=
105
169
sec
h
4
(
1
2
13
[
x
-
205
169
t
-
x
0
]
)
,
and the initial condition is taken from this solution with x
0
=
0:
(17)
u
(
x
,
0
)
=
105
169
sec
h
4
(
1
2
13
x
)
.
Figure 1 shows the numerical solution. A soliton solution can be seen to propagate throughout the interval. There are no “rough parts” and the solution appears to be smooth. To corroborate what is seen in the figure, we computed the error and the conserved quantities in Table 3. We can see that the analytical and numerical solution agree. Further, the conserved quantities remain constant.Table 3
Table error for the KdV-Kawahara equation.
t
L
2
L
∞
I
1
I
2
40.0
1.1155
E
-
03
4.5900
E
-
04
5.9749
1.2737
50.0
1.5330
E
-
03
6.2752
E
-
04
5.9753
1.2741
60.0
2.0416
E
-
03
8.3314
E
-
04
5.9757
1.2744
70.0
3.6657
E
-
03
2.1694
E
-
03
5.9760
1.2747Figure 1
Pulse propagation of the KdV-Kawahara equation(
μ
≠
0
).


## 3.1. Traveling Wave Solution of the Generalized Kawahara Equation

We show the accuracy of proposed method for generalized Kawahara equation with nonlinearityn
=
1,2 and we let γ
=
p
=
q
=
1 and μ
=
0. For computational purposes, the infinite interval is truncated to x
∈
[
-
100,100
] and t
∈
[
0,50
]. We use periodic boundary conditions to simulate the behaviour at infinity. For all these experiments, we have Δ
t
=
.
005 and the number of space points is 512.


### 3.1.1. Kawahara Equationn
=
1

For the casen
=
1, the equation reduces to the well known Kawahara equation. An exact solution is given by
(14)
u
(
x
,
t
)
=
105
169
sec
h
4
(
1
2
13
(
x
-
36
169
t
)
)
.In Table1, we see the numerical results obtained by the Fourier splitting technique on the Kawahara equation. As t increases, the error grows. This growth can be attributed to periodic boundary conditions simulating the behaviour at infinity. This is corroborated by the fact that L
∞ is obtained at points near end of the computational domain. However, our numerical solution is reliable because the conserved quantities remain within numerical error throughout the computational domain.Table 1
Absolute errors and conserved quantities for the Kawahara equation (casen
=
1).
t
L
2
L
∞
I
1
I
2
20.0
5.1011
E
-
05
2.1776
E
-
05
5.9737
1.2726
30.0
7.7990
E
-
05
3.2469
E
-
05
5.9738
1.2726
40.0
1.1144
E
-
04
4.5564
E
-
05
5.9738
1.2726
50.0
1.5310
E
-
04
6.2665
E
-
05
5.9739
1.2727


### 3.1.2. Kawahara Equationn
=
2

The second case we examine is the generalized Kawaharan
=
2. An exact solution is given by
(15)
u
(
x
,
t
)
=
3
10
sec
h
2
(
1
2
5
(
x
-
4
25
t
)
)
.In Table2, we see again an excellent agreement between the numerical results obtained by the Fourier splitting technique and the analytical solution. The errors increase with time but these errors are not numerically significant and the maximum value of the error is always near the boundary.Table 2
Absolute errors and conserved quantities for the Kawahara equation (casen
=
2).
t
L
2
L
∞
I
1
I
2
20.0
1.3882
E
-
04
5.9450
E
-
05
8.4854
2.6835
30.0
2.5212
E
-
04
1.0612
E
-
04
8.4855
2.6836
40.0
4.0950
E
-
04
1.6662
E
-
04
8.4855
2.6837
50.0
6.1195
E
-
04
2.4806
E
-
04
8.4856
2.6837Our method can accurately approximate the analytical solutions for the integer valuesn
>
2. We omit these results as they are not particularly illuminating.


### 3.1.3. KdV-Kawahara Equation

For our final example, we setμ
=
1. With this additional term, we have the so-called KdV-Kawahara which has applications in magnetohydrodynamics; see [9] for details. An exact traveling wave solution to this equation is given by
(16)
u
(
x
,
t
)
=
105
169
sec
h
4
(
1
2
13
[
x
-
205
169
t
-
x
0
]
)
,
and the initial condition is taken from this solution with x
0
=
0:
(17)
u
(
x
,
0
)
=
105
169
sec
h
4
(
1
2
13
x
)
.
Figure 1 shows the numerical solution. A soliton solution can be seen to propagate throughout the interval. There are no “rough parts” and the solution appears to be smooth. To corroborate what is seen in the figure, we computed the error and the conserved quantities in Table 3. We can see that the analytical and numerical solution agree. Further, the conserved quantities remain constant.Table 3
Table error for the KdV-Kawahara equation.
t
L
2
L
∞
I
1
I
2
40.0
1.1155
E
-
03
4.5900
E
-
04
5.9749
1.2737
50.0
1.5330
E
-
03
6.2752
E
-
04
5.9753
1.2741
60.0
2.0416
E
-
03
8.3314
E
-
04
5.9757
1.2744
70.0
3.6657
E
-
03
2.1694
E
-
03
5.9760
1.2747Figure 1
Pulse propagation of the KdV-Kawahara equation(
μ
≠
0
).


## 3.1.1. Kawahara Equationn
=
1

For the casen
=
1, the equation reduces to the well known Kawahara equation. An exact solution is given by
(14)
u
(
x
,
t
)
=
105
169
sec
h
4
(
1
2
13
(
x
-
36
169
t
)
)
.In Table1, we see the numerical results obtained by the Fourier splitting technique on the Kawahara equation. As t increases, the error grows. This growth can be attributed to periodic boundary conditions simulating the behaviour at infinity. This is corroborated by the fact that L
∞ is obtained at points near end of the computational domain. However, our numerical solution is reliable because the conserved quantities remain within numerical error throughout the computational domain.Table 1
Absolute errors and conserved quantities for the Kawahara equation (casen
=
1).
t
L
2
L
∞
I
1
I
2
20.0
5.1011
E
-
05
2.1776
E
-
05
5.9737
1.2726
30.0
7.7990
E
-
05
3.2469
E
-
05
5.9738
1.2726
40.0
1.1144
E
-
04
4.5564
E
-
05
5.9738
1.2726
50.0
1.5310
E
-
04
6.2665
E
-
05
5.9739
1.2727


## 3.1.2. Kawahara Equationn
=
2

The second case we examine is the generalized Kawaharan
=
2. An exact solution is given by
(15)
u
(
x
,
t
)
=
3
10
sec
h
2
(
1
2
5
(
x
-
4
25
t
)
)
.In Table2, we see again an excellent agreement between the numerical results obtained by the Fourier splitting technique and the analytical solution. The errors increase with time but these errors are not numerically significant and the maximum value of the error is always near the boundary.Table 2
Absolute errors and conserved quantities for the Kawahara equation (casen
=
2).
t
L
2
L
∞
I
1
I
2
20.0
1.3882
E
-
04
5.9450
E
-
05
8.4854
2.6835
30.0
2.5212
E
-
04
1.0612
E
-
04
8.4855
2.6836
40.0
4.0950
E
-
04
1.6662
E
-
04
8.4855
2.6837
50.0
6.1195
E
-
04
2.4806
E
-
04
8.4856
2.6837Our method can accurately approximate the analytical solutions for the integer valuesn
>
2. We omit these results as they are not particularly illuminating.


## 3.1.3. KdV-Kawahara Equation

For our final example, we setμ
=
1. With this additional term, we have the so-called KdV-Kawahara which has applications in magnetohydrodynamics; see [9] for details. An exact traveling wave solution to this equation is given by
(16)
u
(
x
,
t
)
=
105
169
sec
h
4
(
1
2
13
[
x
-
205
169
t
-
x
0
]
)
,
and the initial condition is taken from this solution with x
0
=
0:
(17)
u
(
x
,
0
)
=
105
169
sec
h
4
(
1
2
13
x
)
.
Figure 1 shows the numerical solution. A soliton solution can be seen to propagate throughout the interval. There are no “rough parts” and the solution appears to be smooth. To corroborate what is seen in the figure, we computed the error and the conserved quantities in Table 3. We can see that the analytical and numerical solution agree. Further, the conserved quantities remain constant.Table 3
Table error for the KdV-Kawahara equation.
t
L
2
L
∞
I
1
I
2
40.0
1.1155
E
-
03
4.5900
E
-
04
5.9749
1.2737
50.0
1.5330
E
-
03
6.2752
E
-
04
5.9753
1.2741
60.0
2.0416
E
-
03
8.3314
E
-
04
5.9757
1.2744
70.0
3.6657
E
-
03
2.1694
E
-
03
5.9760
1.2747Figure 1
Pulse propagation of the KdV-Kawahara equation(
μ
≠
0
).


## 4. Conclusions

In this work, we have shown that the split step Fourier method can solve the generalized Kawahara equation very efficiently. The method is fast because of the use of FFTs and the exponential operator method. The numerical solution obtained was shown to agree with known analytical results. We have observed that at the boundaries the error increases. This is attributed to the fact that we use periodic boundary conditions to simulate the behavior at infinity. The value ofn does not make the problem any more complicated for this method. In the future, we will study absorbing boundary conditions to alleviate the error near the boundaries. Finally, the authors would like to comment that this method should be applied to other nonlinear PDEs with high order derivatives.



---

*Source: 894956-2014-08-26.xml*