skip to main content

Title: Discretizing singular point sources in hyperbolic wave propagation problems

Here, we develop high order accurate source discretizations for hyperbolic wave propagation problems in first order formulation that are discretized by finite difference schemes. By studying the Fourier series expansions of the source discretization and the finite difference operator, we derive sufficient conditions for achieving design accuracy in the numerical solution. Only half of the conditions in Fourier space can be satisfied through moment conditions on the source discretization, and we develop smoothness conditions for satisfying the remaining accuracy conditions. The resulting source discretization has compact support in physical space, and is spread over as many grid points as the number of moment and smoothness conditions. In numerical experiments we demonstrate high order of accuracy in the numerical solution of the 1-D advection equation (both in the interior and near a boundary), the 3-D elastic wave equation, and the 3-D linearized Euler equations.
 [1] ;  [2] ;  [1] ;  [2]
  1. Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
  2. Stanford Univ., Stanford, CA (United States)
Publication Date:
OSTI Identifier:
Report Number(s):
Journal ID: ISSN 0021-9991
Grant/Contract Number:
Accepted Manuscript
Journal Name:
Journal of Computational Physics
Additional Journal Information:
Journal Volume: 321; Journal Issue: C; Journal ID: ISSN 0021-9991
Research Org:
Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
Sponsoring Org:
Country of Publication:
United States
58 GEOSCIENCES; 97 MATHEMATICS, COMPUTING, AND INFORMATION SCIENCE; singular sources; hyperbolic wave propagation; moment conditions; smoothness conditions; summation by parts