Semi-analytical Exponential RKN Integrators for Efficiently Solving High-Dimensional Nonlinear Wave Equations Based on FFT Techniques.
Lijie Mei,Li Huang,Xinyuan Wu,Shixiang Huang
DOI: https://doi.org/10.1016/j.cpc.2019.04.017
IF: 4.717
2019-01-01
Computer Physics Communications
Abstract:In the light of the analytical expression of solutions to nonlinear wave equations (Wu et al. (2015)), using FFT techniques in this paper, we present a class of semi-analytical exponential/extended Runge–Kutta–Nyström (ERKN) integrators, which can nearly preserve the spatial continuity as well as the oscillation of the nonlinear waves equations. The standard ERKN methods require, in every time step, the computation of the matrix–vector product whose computational complexity in terms of basic multiplication is O(N2), once a direct calculation procedure is implemented. To achieve the near preservation of spatial continuity in the implementation of ERKN integrators, the number of spatial mesh grids N is required to be large enough, which may result in the massive calculation in the magnitude of O(N2) for each time step and the possible overflow of memory once a direct calculation procedure is implemented. To avoid this potential obstacle in applications, we design and analyse two novel algorithms which are incorporated with the Fast Fourier Transform (FFT) in the implementation of ERKN integrators. It is shown that the two algorithms reduce the magnitude of calculation amount from O(N2) to O(NlogN) in terms of basic multiplication. Meanwhile, they can greatly decrease the occupied memory, which allow the sufficient mesh grids to nearly preserve the spatial continuity especially for high-dimensional wave equations. Finally, numerical experiments including a two-dimensional sine–Gordon equation are carried out, and the numerical results soundly support the superiority and potentiality of the new algorithms.