Abstract:The Navier equation is the governing equation of elastic waves, and computing its solution accurately and rapidly has a wide range of applications in geophysical exploration, materials science, etc. In this paper, we focus on the efficient and high-precision numerical algorithm for the time harmonic elastic wave scattering problems from cornered domains via the boundary integral equations in two dimensions. The approach is based on the combination of Nyström discretization, analytical singular integrals and kernel-splitting method, which results in a high-order solver for smooth boundaries. It is then combined with the recursively compressed inverse preconditioning (RCIP) method to solve elastic scattering problems from cornered domains. Numerical experiments demonstrate that the proposed approach achieves high accuracy, with stabilized errors close to machine precision in various geometric configurations. The algorithm is further applied to investigate the asymptotic behavior of density functions associated with boundary integral operators near corners, and the numerical results are highly consistent with the theoretical formulas.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to develop a numerical algorithm with high efficiency and high precision for solving the time - harmonic elastic wave scattering problem in two - dimensional angular domains. Specifically, the author focuses on how to accurately and rapidly calculate the solution of the Navier equation in angular domains (i.e., regions whose boundaries contain sharp corners or edges).
### Background and Problem Description of the Paper
The Navier equation is the basic governing equation for elastic waves and has wide applications in geophysical exploration, materials science, etc. However, for the elastic wave scattering problem in angular domains, due to the interaction between the geometric shape and the properties of the wave equation, complex phenomena such as corner singularities and diffraction effects will occur, which makes it very difficult to design effective numerical methods.
### Main Research Contents
1. **Problem Definition**:
- Consider a two - dimensional angular domain filled with a homogeneous isotropic elastic medium.
- This domain is constrained by Dirichlet or Neumann boundary conditions.
- The goal is to develop a robust and high - precision numerical method to solve these scattering problems.
2. **Method Overview**:
- Use a combination of Nyström discretization, analytic singular integral, and kernel splitting methods to handle high - order solvers on smooth boundaries.
- For non - smooth boundaries (i.e., when there are corner points), combine the Recursive Compressed Inverse Preconditioning (RCIP) method to handle corner singularities.
3. **Innovation Points**:
- Propose a new singular kernel decomposition method, which is compact and easy to calculate.
- Demonstrate the effectiveness of singular integral discretization through detailed numerical implementation.
- Verify the asymptotic behavior of the boundary integral operator density function near corner points, and the results are highly consistent with theoretical predictions.
### Mathematical Formulas
- Navier Equation:
\[
\mu \Delta \mathbf{u}+(\lambda + \mu)\nabla(\nabla\cdot\mathbf{u})+\rho\omega^{2}\mathbf{u} = 0,\quad x\in\Omega
\]
where \(\mu\) and \(\lambda\) are Lamé coefficients, \(\rho\) is the mass density, and \(\omega\) is the angular frequency.
- Wave Numbers of Compressional and Shear Waves:
\[
k_{p}=\frac{\omega}{\sqrt{\frac{\lambda + 2\mu}{\rho}}},\quad k_{s}=\frac{\omega}{\sqrt{\frac{\mu}{\rho}}}
\]
- Free - Space Green Function:
\[
G(x,y)=\frac{i}{4\mu}\left\{H^{(1)}_{0}(k_{s}|x - y|)I+\frac{1}{k_{s}^{2}}\nabla_{x}\nabla_{x}^{T}\left(H^{(1)}_{0}(k_{s}|x - y|)-H^{(1)}_{0}(k_{p}|x - y|)\right)\right\}
\]
where \(H^{(1)}_{0}\) is the zero - order Hankel function of the first kind, and \(I\) is a 2×2 identity matrix.
### Summary
This paper proposes a novel and efficient algorithm by combining multiple numerical techniques, which solves the computational challenges of elastic wave scattering problems in angular domains. This method not only improves the computational precision but also significantly reduces the demand for computational resources, providing strong support for applications in related fields.