skip to main content

SciTech ConnectSciTech Connect

This content will become publicly available on June 1, 2017

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.
Authors:
 [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:
1343000
Report Number(s):
LLNL-JRNL--679293
Journal ID: ISSN 0021-9991
Grant/Contract Number:
AC52-07NA27344
Type:
Accepted Manuscript
Journal Name:
Journal of Computational Physics
Additional Journal Information:
Journal Volume: 321; Journal Issue: C; Journal ID: ISSN 0021-9991
Publisher:
Elsevier
Research Org:
Lawrence Livermore National Laboratory (LLNL), Livermore, CA (United States)
Sponsoring Org:
USDOE
Country of Publication:
United States
Language:
English
Subject:
58 GEOSCIENCES; 97 MATHEMATICS, COMPUTING, AND INFORMATION SCIENCE singular sources; hyperbolic wave propagation; moment conditions; smoothness conditions; summation by parts