Abstract:We develop a stable finite difference method for the elastic wave equation in bounded media, where the material properties can be discontinuous at curved interfaces. The governing equation is discretized in second order form by a fourth or sixth order accurate summation-by-parts operator. The mesh size is determined by the velocity structure of the material, resulting in nonconforming grid interfaces with hanging nodes. We use order-preserving interpolation and the ghost point technique to couple adjacent mesh blocks in an energy-conserving manner, which is supported by a fully discrete stability analysis. In our previous work for the wave equation, two pairs of order-preserving interpolation operators are needed when imposing the interface conditions weakly by a penalty technique. Here, we only use one pair in the ghost point method. In numerical experiments, we demonstrate that the convergence rate is optimal, and is the same as when a globally uniform mesh is used in a single domain. In addition, with a predictor-corrector time integration method, we obtain time stepping stability with stepsize almost the same as given by the usual Courant-Friedrichs-Lewy condition.
What problem does this paper attempt to address?
This paper aims to solve the numerical solution problem of the elastic wave equation in bounded media with non - conformal interfaces, especially when the material properties are discontinuous on the curved interface. Specifically, the paper proposes a stable finite - difference method to handle this situation.
### Problem Background
1. **Elastic Wave Equation**: It studies elastic waves propagating in a bounded domain, where the material properties can change abruptly on the curved interface.
2. **Non - conformal Interface**: Due to the change of material properties, the grid has different resolutions in different regions, forming a non - conformal interface (that is, the grid points at the grid interface do not match).
3. **High - order Accuracy Requirement**: In order to improve computational efficiency and accuracy, a high - order accuracy numerical method is required.
### Problems Solved in the Paper
1. **Develop High - order Finite - difference Method**: The paper proposes a method based on fourth - or sixth - order accurate Summation - by - Parts (SBP) operators for discretizing the elastic wave equation.
2. **Handle Non - conformal Interface**: By introducing the ghost - point technique and the order - preserving interpolation method, the energy - conserving coupling between adjacent grid blocks is ensured.
3. **Stability and Convergence**: Through a fully - discrete stability analysis, it is proved that this method has a stability similar to the traditional Courant–Friedrichs–Lewy condition in the time step, and the optimal convergence rate is shown in numerical experiments.
4. **Reduce the Influence of Interpolation Pairs**: Compared with previous work, this method only requires one pair of order - preserving interpolation operators instead of two pairs, thus simplifying the algorithm and improving efficiency.
### Main Contributions
- **High - order Accuracy**: By using fourth - or sixth - order accurate SBP operators, the high - order accuracy of the numerical method is ensured.
- **Non - conformal Interface Handling**: The non - conformal interface problem is effectively handled, ensuring the accuracy and stability of the numerical solution.
- **Ghost - point Technique**: By using the ghost - point technique, the boundary conditions and interface conditions can be imposed, avoiding the introduction of low - order errors.
- **Order - preserving Interpolation**: By constructing order - preserving interpolation operators, the overall convergence rate is ensured to be unaffected.
### Conclusion
This paper successfully solves the numerical solution problem of the elastic wave equation in bounded media with non - conformal interfaces, providing an efficient, stable and high - precision finite - difference method, and providing new tools and ideas for research in related fields.