Summation-by-parts (SBP) finite difference methods have several desirable properties for second-order wave equations. They combine the computational efficiency of narrow-stencil finite difference operators with provable stability on curvilinear multiblock grids. While several techniques for boundary and interface conditions exist, weak imposition via simultaneous approximation terms (SATs) is perhaps the most flexible one. Although SBP methods have been applied to elastic wave equations many times, an SBP-SAT method for general anisotropic elastic wave equations has not yet been presented in the literature. We fill this gap by deriving energy-stable self-adjoint SBP-SAT methods for general anisotropic materials on curvilinear multiblock grids. The methods are based on fully compatible SBP operators. Although this paper focuses on classical SBP finite difference operators, the presented boundary and interface treatments are general and apply to a range of methods that satisfy an SBP property. We demonstrate the stability and accuracy properties of a particular set of fully compatible SBP-SAT schemes using the method of manufactured solutions. We also demonstrate the utility of the new method in elastodynamic cloaking and seismic imaging in mountainous regions.

零件加法（SBP）有限差分方法对于二阶波动方程具有几个理想的属性。他们将窄模板有限差分算子的计算效率与曲线多块网格上的可证明的稳定性相结合。尽管存在几种边界和界面条件的技术，但通过同时逼近项（SAT）进行的弱拼版可能是最灵活的一种。尽管SBP方法已多次应用于弹性波方程，但文献中尚未提出用于一般各向异性弹性波方程的SBP-SAT方法。我们通过为曲线多块网格上的一般各向异性材料导出能量稳定的自伴SBP-SAT方法来填补这一空白。这些方法基于完全兼容的SBP运算符。尽管本文着重于经典SBP有限差分算子，但所提出的边界和界面处理是通用的，适用于满足SBP特性的一系列方法。我们使用制造的解决方案的方法演示了一组特定的完全兼容的SBP-SAT方案的稳定性和准确性。我们还演示了该新方法在山区的弹性动力隐身和地震成像中的实用性。