Fourier Methods with Extended Stability Intervals for the Korteweg–de Vries Equation

Abstract

A full leap frog Fourier method for integrating the Korteweg-de Vries (KdV) equation ut + uux-eUxxx 0 results in an O(N-3) stability constraint on the time step, where N is the number of Fourier modes used. This stability limit is much more restrictive than the accuracy limit for many applications.In this paper we propose a method for which the stability limit is extended by treating the linear dispersive ux term implicitly. Thus timesteps can be taken up to an accuracy limit larger than the explicit stability limit. The implicit method is implemented without solving linear systems by integrating in time in the Fourier space and discretizing the nonlinear uu term by leap frog. A second method we propose uses basis functions which solve the linear part of the KdV equation and leap frog for time integration. A linearized stability analysis of the proposed schemes proves that a version of the first scheme possesses a certain kind of unconditional stability and that the second scheme has an O(N-1) stability limit. The accuracy of the schemes for soliton propagation is analyzed by examining the truncation error. In addition, we analyze a linearization of a nonlinear finite element scheme proposed by Winther that treats the ux term by Crank-Nicolson. Numerical experiments on soliton solutions show that the linearized stability analysis gives accurate predictions for all the nonlinear schemes, that the truncation errors give an estimate for the accuracies and that the Fourier methods are more accurate than the finite element method.

A Fast Spectral Algorithm for Nonlinear Wave Equations with Linear Dispersion

“…Since theL-operator is diagonal, the need to invert theL term does not incur any extra computational cost. This idea has appeared many times in the literature for both multistep and Runge-Kutta methods, often for diffusive problems and bearing a name such as "implicit-explicit" or "linearly implicit" [1][2][3][4][5]. Graphically, we illustrate a scheme of the form (4) in Fig.…”

“…2. Employing a slightly different u-stencil than the one for AB/AM, this type of argument led Chan and Kerkhoven [5] …”

Numerical study of acoustic modes in ducted shear flow

“…The first of these two approximations is the classical Crank-Nicolson (trapezoidal) scheme, the second scheme is its more stable modification suggested in Chan & Kerkhoven [24].…”

Numerical Solutions of the Forced Perturbed Korteweg-De Vries Equation With Variable Coefficients

“…It has been used to solve the KdV [13], Burgers [14] and KdV-Burgers equations [15]. The Chan and Kerkhoven [13] scheme is extended for the FpKdV equation (1).…”