This paper presents a fully discrete numerical scheme for one-dimensional nonlocal wave equations and provides a rigorous theoretical analysis. To facilitate the spatial discretization, we introduce an auxiliary variable analogous to the gradient field in local discontinuous Galerkin (DG) methods for classical partial differential equations (PDEs) and reformulate the equation into a system of equations. The proposed scheme then uses a DG method for spatial discretization and the Crank-Nicolson method for time integration. We prove optimal L2 error convergence for both the solution and the auxiliary variable under a special class of radial kernels at the semi-discrete level. In addition, for general kernels, we demonstrate the asymptotic compatibility of the scheme, ensuring that it recovers the classical DG approximation of the local wave equation in the zero-horizon limit. Furthermore, we prove that the fully discrete scheme preserves the energy of the nonlocal wave equation. A series of numerical experiments are presented to validate the theoretical findings.
Du et al. (Thu,) studied this question.