Abstract:To improve the computational accuracy and efficiency of long-time wavefield simulations, we have developed a so-called symplectic interior penalty discontinuous Galerkin (IPDG) method for 2D acoustic equation. For the symplectic IPDG method, the scalar wave equation is first transformed into a Hamiltonian system. Then, the high-order IPDG formulations are introduced for spatial discretization because of their high accuracy and ease of dealing with computational domains with complex boundaries. The time integration is performed using an explicit third-order symplectic partitioned Runge-Kutta scheme so that it preserves the Hamiltonian structure of the wave equation in long-term simulations. Consequently, the symplectic IPDG method combines the advantages of discontinuous Galerkin method and the symplectic time integration. We investigate the properties of the method in detail for high-order spatial basis functions, including the stability criteria, numerical dispersion and dissipation relationships, and numerical errors. The analyses indicate that the symplectic SIPG method is nondissipative and retains low numerical dispersion. We also find that different symplectic IPDG methods have different convergence behaviors. It is indicated that using coarse meshes with a high-order method produces smaller errors and retains high accuracy. We have applied our method to simulate the scalar wavefields for different models, including layered models, a rough topography model, and the Marmousi model. The numerical results show that the symplectic IPDG method can suppress numerical dispersion effectively and provide accurate information on the wavefields. We also conduct a long-term experiment that verifies the capability of symplectic IPDG method for long-time simulations.

A low dispersion precise integration time domain method based on wavelet galerkin scheme

An Improved MOD Discontinuous Galerkin Method for Solving Time-Domain Electric Field Integral Equation

Discontinuous Galerkin Method Using Laguerre Polynomials for Solving a Time-Domain Electric Field Integral Equation

An Efficient Method to Reduce the Numerical Dispersion in the LOD-FDTD Method Based on the (2, 4) Stencil

A Fast Time Domain PMCHW Equations for Predicting Transient Responses of Dielectric Objects

An Improved Time Integration Method Based on Galerkin Weak Form with Controllable Numerical Dissipation

Efficient and Parallel Solution of High-order Continuous Time Galerkin for Dissipative and Wave Propagation Problems

Improved Estimation of Time Step Bound for Discontinuous Galerkin Time-Domain Method

Numerical investigation of the element-free Galerkin method (EFGM) with appropriate temporal discretization techniques for transient wave propagation problems

A Time-Domain Discontinuous Galerkin Trefftz Method

A TIME INTEGRATION METHOD BASED ON THE WEAK FORM GALERKIN METHOD

Penalty Factor Threshold and Time Step Bound Estimations for Discontinuous Galerkin Time-Domain Method Based on Helmholtz Equation

A Hybrid Explicit-Implicit Scheme for the Time-Dependent Wigner Equation

Phase-Informed Discontinuous Galerkin Method for Extremely High-Frequency Wave Modeling

An novel finite difference dispersion error elimination mechanism in the Lax–Wendroff high‐order time discretization

A wavelet multi-resolution enabled interpolation Galerkin method for two-dimensional solids

A Global-Local Error-Driven p -Adaptive Scheme For Discontinuous Galerkin Time Domain Method With Local Time Stepping Strategy

Wavelet-Galerkin Method for the Reaction-Diffusion Equation

Symplectic Interior Penalty Discontinuous Galerkin Method for Solving the Seismic Scalar Wave Equation

Highly accurate wavelet solution for the two-dimensional Bratu's problem

A mixed discontinuous Galerkin method for the wave equation