Discretizing singular point sources in hyperbolic wave propagation problems
2016 (English)In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 321, 532-555 p.Article in journal (Refereed) Published
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.
Place, publisher, year, edition, pages
2016. Vol. 321, 532-555 p.
Singular sources, Hyperbolic wave propagation, Moment conditions, Smoothness conditions, Summation by parts
Computational Mathematics Mathematical Analysis
IdentifiersURN: urn:nbn:se:liu:diva-132630DOI: 10.1016/j.jcp.2016.05.060OAI: oai:DiVA.org:liu-132630DiVA: diva2:1047470