Optimal Energy Conserving Local Discontinuous Galerkin Methods for Elastodynamics: Semi and Fully Discrete Error Analysis

Abstract

We present an arbitrary high-order local discontinuous Galerkin (LDG) method with alternating fluxes for solving linear elastodynamics problems in isotropic media. Both the semi-discrete analysis and fully discrete analysis for a leap-frog LDG method are given to show that the proposed method simultaneously enjoys the energy conserving property and optimal convergence rates in both the displacement and stress, when the tensor product polynomials of the degree k are used on Cartesian meshes. Numerical experiments demonstrate that the proposed method has several advantages including the exact energy conservation, slow-growing errors in long time simulation, and subtle dependence on the first Lame parameter $$\lambda $$ .