Please use this identifier to cite or link to this item: https://doi.org/10.1137/110855004
Title: An exponential wave integrator sine pseudospectral method for the klein-gordon-zakharov system
Authors: Bao, W. 
Dong, X.
Zhao, X.
Keywords: E-scalability
Error bounds
Exponential wave integrator
High-plasma-frequency limit
Klein-Gordon-Zakharov system
Sine pseudospectral method
Subsonic limit
Issue Date: 2013
Citation: Bao, W., Dong, X., Zhao, X. (2013). An exponential wave integrator sine pseudospectral method for the klein-gordon-zakharov system. SIAM Journal on Scientific Computing 35 (6) : A2903-A2927. ScholarBank@NUS Repository. https://doi.org/10.1137/110855004
Abstract: An exponential wave integrator sine pseudospectral method is presented and analyzed for discretizing the Klein-Gordon-Zakharov (KGZ) system with two dimensionless parameters 0 > e ≤ 1 and 0 > γ = 1 which are inversely proportional to the plasma frequency and the speed of sound, respectively. The main idea in the numerical method is to apply the sine pseudospectral discretization for spatial derivatives followed by using an exponential wave integrator for temporal derivatives in phase space. The method is explicit, symmetric in time, and it is of spectral accuracy in space and second-order accuracy in time for any fixed e = e0 and γ = γ0. In the O(1)-plasma frequency and speed of sound regime, i.e., e = O(1) and γ = O(1), we establish rigorously error estimates for the numerical method in the energy space H1 × L2. We also study numerically the resolution of the method in the simultaneous high-plasma-frequency and subsonic limit regime, i.e., (e, γ) - 0 under e ≥ γ. In fact, in this singular limit regime, the solution of the KGZ system is highly oscillating in time, i.e., there are propagating waves with wavelength of O(e2) and O(1) in time and space, respectively. Our extensive numerical results suggest that, in order to compute "correct" solutions in the simultaneous high-plasma-frequency and subsonic limit regime, the meshing strategy (or e-scalability) is time step t = O(e2) and mesh size h = O(1) independent of e. Finally, we also observe numerically that the method has the property of near conservation of the energy over long time in practical computations.Copyright © by SIAM.
Source Title: SIAM Journal on Scientific Computing
URI: http://scholarbank.nus.edu.sg/handle/10635/102835
ISSN: 10648275
DOI: 10.1137/110855004
Appears in Collections:Staff Publications

Show full item record
Files in This Item:
There are no files associated with this item.

SCOPUSTM   
Citations

14
checked on May 25, 2018

WEB OF SCIENCETM
Citations

11
checked on May 9, 2018

Page view(s)

29
checked on May 4, 2018

Google ScholarTM

Check

Altmetric


Items in DSpace are protected by copyright, with all rights reserved, unless otherwise indicated.