A HIGH-ORDER WEIGHTED RUNGE-KUTTA DISCONTINUOUS GALERKIN METHOD FOR SOLVING 2D ACOUSTIC AND ELASTIC WAVE EQUATIONS IN ISOTROPIC AND ANISOTROPIC MEDIA

A high-order weighted Runge-Kutta Discontinuous Galerkin Method for solving 2D acoustic and elastic wave equations in isotropic and anisotropic media is proposed in this paper, which is an extension of the existing first-order and second-order methods to higher-order cases. For this method, second-order seismic wave equations are first transformed into a first-order hyperbolic system, then local Lax-Friedrichs (LLF) numerical flux discontinuous Galerkin formulations for spatial discretization are employed, directly leading to a semi discrete ordinary differential equation (ODE) system. For time discretization, an implicit diagonal Runge-Kutta method is introduced. To avoid solving a large-scale system of linear equations, a two-step explicit iterative process is implemented. In addition, a weighting factor is introduced for the iteration to enrich the method. The basis functions we use are 1st similar to 5th order polynomials, leading to 2nd - and 6th order of spatial accuracy. Numerical properties of the high-order weighted Runge-Kutta Discontinuous Galerkin Method are investigated in detail, including numerical error, stability criteria and numerical dispersion, which validate the superiority of the high order method. The proposed method is then applied to several 2D wave propagation problems in isotropic and anisotropic media, including acoustic-elastic interface problems. Results illustrate that this method can effectively suppress numerical dispersion and provide accurate information on the wave field on coarse mesh. We also compare the proposed method with the finite difference method to investigate the computational efficiency.