Abstract:In this paper, we present and analyze fully discrete finite difference schemes designed for solving the initial value problem associated with the fractional Korteweg-de Vries (KdV) equation involving the fractional Laplacian. We design the scheme by introducing the discrete fractional Laplacian operator which is consistent with the continuous operator, and posses certain properties which are instrumental for the convergence analysis. Assuming the initial data (u_0 \in H^{1+\alpha}(\mathbb{R})), where (\alpha \in [1,2)), our study establishes the convergence of the approximate solutions obtained by the fully discrete finite difference schemes to a classical solution of the fractional KdV equation. Theoretical results are validated through several numerical illustrations for various values of fractional exponent $\alpha$. Furthermore, we demonstrate that the Crank-Nicolson finite difference scheme preserves the inherent conserved quantities along with the improved convergence rates.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to develop a fully discrete finite - difference scheme for solving the initial - value problems related to the fractional - order Korteweg - de Vries (KdV) equation. Specifically, the form of the equation considered in the paper is:
\[
\begin{cases}
u_t+\left(\frac{u^2}{2}\right)_x-(-\Delta)^{\alpha / 2}u_x = 0, &(x, t)\in Q_T:=\mathbb{R}\times(0, T],\\
u(x, 0)=u_0(x), &x\in\mathbb{R},
\end{cases}
\]
where \(T>0\) is fixed, \(u_0\) is the given initial condition, and \(u:Q_T\rightarrow\mathbb{R}\) represents the unknown solution. The non - local operator \(-(-\Delta)^{\alpha / 2}\) corresponds to the fractional - order Laplacian operator, introducing non - locality in the dynamics of the equation. The parameter \(\alpha\in[1, 2)\) determines the order of the fractional - order and plays a crucial role in the behavior of the solution.
The main objective of the paper is to design a discrete fractional - order Laplacian operator that is consistent with the continuous operator and has certain properties that are essential for the convergence analysis. Assuming that the initial data \(u_0\in H^{1 +\alpha}(\mathbb{R})\), the study establishes the convergence of the approximate solution obtained by the fully discrete finite - difference scheme to the classical solution of the fractional - order KdV equation. The theoretical results are verified by numerical experiments, which are carried out for different fractional - order parameters \(\alpha\).
Furthermore, the paper also shows that the Crank - Nicolson finite - difference scheme not only preserves the conserved quantities (such as mass, momentum, and energy) in the classical solution but also has an improved convergence rate. These conserved quantities are defined as:
\[
C_1(u):=\int_{\mathbb{R}}u(x, t)dx,\quad C_2(u):=\int_{\mathbb{R}}u^2(x, t)dx,
\]
\[
C_3(u):=\int_{\mathbb{R}}\left((-(-\Delta)^{\alpha / 4}u)^2-\frac{u^3}{3}\right)(x, t)dx,\quad\alpha\in[1, 2).
\]
Through numerical experiments, the paper proves that the Crank - Nicolson scheme is superior in performance to the implicit Euler scheme, which is as expected.