Discretizing singular point sources in hyperbolic wave propagation problems

Citations of this article
Mendeley users who have this article in their library.


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.




Petersson, N. A., O’Reilly, O., Sjögreen, B., & Bydlon, S. (2016). Discretizing singular point sources in hyperbolic wave propagation problems. Journal of Computational Physics, 321, 532–555. https://doi.org/10.1016/j.jcp.2016.05.060

Register to see more suggestions

Mendeley helps you to discover research relevant for your work.

Already have an account?

Save time finding and organizing research with Mendeley

Sign up for free