Abstract:We develop a high order accurate numerical method for solving the elastic wave equation in second-order form. We hybridize the computationally efficient Cartesian grid formulation of finite differences with geometrically flexible discontinuous Galerkin methods on unstructured grids by a penalty based technique. At the interface between the two methods, we construct projection operators for the pointwise finite difference solutions and discontinuous Galerkin solutions based on piecewise polynomials. In addition, we optimize the projection operators for both accuracy and spectrum. We prove that the overall discretization conserves a discrete energy, and verify optimal convergence in numerical experiments.
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is to develop an efficient numerical method for solving the elastic wave equation with complex geometries in discontinuous media. Specifically, the authors propose a hybrid numerical method that combines the computationally efficient Cartesian grid finite - difference method (FD) and the geometrically flexible unstructured - grid discontinuous Galerkin method (DG). A projection operator is constructed at the interface between the two methods through a penalty - based technique to achieve high - order accuracy and stability.
### Background and Problem Description of the Paper
The elastic wave equation has wide applications in fields such as seismic wave modeling, especially when dealing with problems of discontinuous material properties due to strata or complex geometric structures. Analytic solutions are usually not obtainable for such problems, so it is particularly important to develop effective numerical schemes to obtain approximate solutions. Time - dependence requires that numerical schemes possess stability and sufficient accuracy.
### Main Contributions
1. **Hybrid FD - DG Discretization**: The authors propose a hybrid discretization scheme based on the SBP (Summation - by - Parts) finite - difference method and the IPDG (Interior Penalty Discontinuous Galerkin) method for solving the second - order form of the elastic wave equation. By using norm - compatible projections at the FD and DG interface, this method is proven to be optimally convergent and satisfies the discrete energy estimate, ensuring stability.
2. **Optimization of the Projection Operator**: The authors not only consider the accuracy of the projection operator but also optimize its spectral characteristics. By minimizing the second - order error and two different spectral metrics (the Frobenius norm and the eigenvalue - flattening norm inspired by [16]), the authors show how the selection of these parameters can significantly improve the stability and accuracy of the system.
### Mathematical Model and Method
- **Elastic Wave Equation**: In the two sub - domains Ωu and Ωv divided by the horizontal interface Γ, the elastic wave equation can be expressed as:
\[
\rho_u\frac{\partial^2u_1}{\partial t^2}=\mu_u\Delta u_1+(\mu_u + \lambda_u)\frac{\partial}{\partial x}\left(\frac{\partial u_1}{\partial x}+\frac{\partial u_2}{\partial y}\right),
\]
\[
\rho_u\frac{\partial^2u_2}{\partial t^2}=\mu_u\Delta u_2+(\mu_u + \lambda_u)\frac{\partial}{\partial y}\left(\frac{\partial u_1}{\partial x}+\frac{\partial u_2}{\partial y}\right),
\]
Similar equations also apply to the displacement fields v1 and v2 in the sub - domain Ωv.
- **Boundary Conditions and Interface Conditions**: To ensure the well - posedness of the problem, initial conditions, boundary conditions, and appropriate interface conditions need to be imposed. For example, on the interface Γ, the continuity of displacement and traction is required:
\[
\mathbf{u}=\mathbf{v}\quad\text{along}\quad\Gamma,
\]
\[
\sigma(\mathbf{v})\cdot\mathbf{n}^++\sigma(\mathbf{u})\cdot\mathbf{n}^- = 0\quad\text{along}\quad\Gamma,
\]
where \(\mathbf{n}^+\) and \(\mathbf{n}^-\) are the unit outer normal vectors on both sides of the interface respectively.
### Discretization Method
- **SBP Finite - Difference Method**: By introducing a difference operator with the SBP property, high - order accuracy can be achieved while maintaining stability.
- **IPDG Method**: By introducing an interior