Numerical methods for the generalized Zakharov system

  • Authors:

  • Weizhu Bao;Fangfang Sun;G. W. Wei

  • Affiliations:

  • Department of Computational Science, National University of Singapore, Singapore 117543, Singapore;Department of Computational Science, National University of Singapore, Singapore 117543, Singapore;Department of Mathematics, Michigan State University, East Lansing, MI

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2003

Quantified Score

Hi-index 31.47

Visualization

Abstract

We present two numerical methods for the approximation of the generalized Zakharov system (ZS). The first one is the time-splitting spectral (TSSP) method, which is explicit, time reversible, and time transverse in variant if the generalized ZS is, keeps the same decay rate of the wave energy as that in the generalized ZS, gives exact results for the plane-wave solution, and is of spectral-order accuracy in space and second-order accuracy in time. The second one is to use a local spectral method, the discrete singular convolution (DSC) for spatial derivatives and the fourth-order Runge-Kutta (RK4) for time integration, which is of high (the same as spectral)-order accuracy in space and can be applied to deal with general boundary conditions. In order to test accuracy and stability, we compare these two methods with other existing methods: Fourier pseudospectral method (FPS) and wavelet-Galerkin method (WG) for spatial derivatives combining with the RK4 for time integration, as well as the standard finite difference method (FD) for solving the ZS with a solitary-wave solution. Furthermore, extensive numerical tests are presented for plane waves, solitary-wave collisions in 1d, as well as a 2d problem of the generalized ZS. Numerical results show that TSSP and DSC are spectralorder accuracy in space and much more accurate than FD, and for stability, TSSP requires k = O(h), DSC-RK4 requires k = O(h2) for fixed acoustic speed, where k is the time step and h is the spatial mesh size.