Abstract
High order finite difference approximations are derived for the second order wave equation with discontinuous coefficients, on rectangular geometries. The discontinuity is treated by splitting the domain at the discontinuities in a multi block fashion. Each sub-domain is discretized with compact second derivative summation by parts operators and the blocks are patched together to a global domain using the projection method. This guarantees a conservative, strictly stable and high order accurate scheme. The analysis is verified by numerical simulations in one and two spatial dimensions.
| Original language | English |
|---|---|
| Pages (from-to) | 249-269 |
| Number of pages | 21 |
| Journal | Journal of Computational Physics |
| Volume | 220 |
| Issue number | 1 |
| DOIs | |
| Publication status | Published - 20 Dec 2006 |
| Externally published | Yes |
Keywords
- Discontinuous media
- High order finite difference methods
- Numerical stability
- Second derivatives
- Wave equation
ASJC Scopus subject areas
- Numerical Analysis
- Modeling and Simulation
- Physics and Astronomy (miscellaneous)
- General Physics and Astronomy
- Computer Science Applications
- Computational Mathematics
- Applied Mathematics