Energy stable and high-order-accurate finite difference methods on staggered grids

JOURNAL OF COMPUTATIONAL PHYSICS

卷:346

Article

O'Reilly, Ossian^1,2  Lundquist, Tomas²  Dunham, Eric M.^1,3  Nordstrom, Jan² 

[1] Stanford Univ, Dept Geophys, Stanford, CA 94305 USA

[2] Linkoping Univ, Dept Math, Div Computat Math, SE-58183 Linkoping, Sweden

[3] Stanford Univ, Inst Computat & Math Engn, Stanford, CA 94305 USA

关键词: Staggered grids;    High order finite difference methods;    Summation-by-parts;    Weakly enforced boundary conditions;    Energy stability;    Wave propagation;   

DOI  :  10.1016/j.jcp.2017.06.030

来源: Elsevier

【 摘 要 】

For wave propagation over distances of many wavelengths, high-order finite difference methods on staggered grids are widely used due to their excellent dispersion properties. However, the enforcement of boundary conditions in a stable manner and treatment of interface problems with discontinuous coefficients usually pose many challenges. In this work, we construct a provably stable and high-order-accurate finite difference method on staggered grids that can be applied to a broad class of boundary and interface problems. The staggered grid difference operators are in summation-by-parts form and when combined with a weak enforcement of the boundary conditions, lead to an energy stable method on multiblock grids. The general applicability of the method is demonstrated by simulating an explosive acoustic source, generating waves reflecting against a free surface and material discontinuity. (C) 2017 Elsevier Inc. All rights reserved.

【 授权许可】

Free