Non-stiff boundary and interface penalties for narrow-stencil finite difference approximations of the Laplacian on curvilinear multiblock grids

The Laplacian appears in several partial differential equations used to model wave propagation. Summation-by-parts-simultaneous approximation term (SBP-SAT) finite difference methods are often used for such equations, as they combine computational efficiency with provable stability on curvilinear multiblock grids. However, the existing SBP-SAT discretization of the Laplacian quickly becomes prohibitively stiff as grid skewness increases. The stiffness stems from the SATs that impose inter-block couplings and Dirichlet boundary conditions. We resolve this issue by deriving stable SATs whose stiffness is almost insensitive to grid skewness. The new discretization thus allows for large time steps in explicit time integrators, even on very skewed grids. It also applies to the variable-coefficient generalization of the Laplacian. We demonstrate the efficacy and versatility of the new SATs by applying them to acoustic wave propagation problems inspired by marine seismic exploration and infrasound monitoring of volcanoes. (C) 2020 Elsevier Inc. All rights reserved.