Detailed Information

Arbitrary-order symplectic time integrator for the acoustic wave equation using the pseudo-spectral method

Cited 0 time in Web of Science Cited 0 time in Scopus


자연과학대학 협동과정 계산과학전공
Issue Date
서울대학교 대학원
Symplectic time integratoracoustic wave equationpseudo-spectral methodextremely accurate simulation
학위논문 (박사)-- 서울대학교 대학원 : 협동과정 계산과학전공, 2017. 2. 신창수.
A Hamiltonian system is symplectic. To simulate a Hamiltonian system, symplectic time integrators are generally applied
otherwise, the energy or the generalized energy is not conserved in the volume of interest. In this study, the symplectic nature of the acoustic wave system is proven. Then, a symplectic scheme that can be extended arbitrarily in temporal dimensions is suggested. The method is based on the Lax-Wendroff expansion of the time differentiation of acoustic wave variables, such as pressure and velocity, existing on the staggered time axis, i.e., one is on the integer grid, and the other is defined on the half integer of the time step. The series can be reduced to the pseudo-differential operator, which enables the application of other approximation techniques, such as the Jacobi-Anger expansion. By virtue of considering the property of the nature of the acoustic wave phenomena, the scheme is more stable and accurate than methods that do not consider symplecticity. Moreover, the phase error per time step can be kept sufficiently small to conduct simulation over long periods of time. According to the analysis of the scheme, the larger the time strides are, the more efficient the simulation is in terms of computing power when a sufficient number of multiplications of the map are accumulated. The effectiveness and accuracy are verified through simulation results using a homogeneous model in which the computed wavefield is equivalent to the analytic solution. The numerical results of the wavefield in the heterogeneous model also yield equivalent results irrespective of the time step lengths. The scheme can be applied to the source problems
however, the time step is confined to describing the entire frequency component of the wavelet.
Files in This Item:
Appears in Collections:


Item View & Download Count

  • mendeley

Items in S-Space are protected by copyright, with all rights reserved, unless otherwise indicated.