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.

The multi-symplectic Fourier pseudospectral method for solving two-dimensional Hamiltonian PDEs

SYMPLECTIC PSEUDOSPECTRAL TIME-DOMAIN SCHEME FOR SOLVING TIME-DEPENDENT SCHRODINGER EQUATION

High-order Symplectic FDTD Scheme for Solving a Time-Dependent Schrödinger Equation

Energy conserving methods for Hamiltonian PDEs based on spectral space decomposition

General Local Energy-Preserving Integrators for Solving Multi-Symplectic Hamiltonian PDEs.

Direct/split invariant-preserving Fourier pseudo-spectral methods for the rotation-two-component Camassa–Holm system with peakon solitons

An exact in time Fourier pseudospectral method with multiple conservation laws for three-dimensional Maxwell’s equations

Multi-symplectic exponentially fitted RKN methods for wave equations

Conservative Multiquadric Quasi-Interpolation Method for Hamiltonian Wave Equations

Applications of the Multi-Symplectic Euler-Box Scheme

Symplectic Hamiltonian Finite Element Methods for Semilinear Wave Propagation

Three kinds of novel multi-symplectic methods for stochastic Hamiltonian partial differential equations

Pseudospectral method for solving PDEs using Matrix Product States

A Fourier Spectral-Discontinuous Galerkin Method for Time-Dependent 3-D Schrödinger–Poisson Equations with Discontinuous Potentials

Multi-Symplectic Method for the Two-Component Camassa–Holm (2CH) System

A symplectic method in two dimensional steady state heat conduction problems

Numerical implementation of the multisymplectic Preissman scheme and its equivalent schemes

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

Fourier-Gegenbauer Pseudospectral Method for Solving Time-Dependent One-Dimensional Fractional Partial Differential Equations with Variable Coefficients and Periodic Solutions

Combined finite difference and spectral methods for the numerical solution of hyperbolic equation with an integral condition

Explicit Multi-Symplectic Extended Leap-Frog Methods for Hamiltonian Wave Equations