Abstract:This paper proposes an $O(N)$ fast direct solver for two-dimensional elastic wave scattering problems. The proxy surface method is extended to elastodynamics to obtain shared coefficients for low-rank approximations from discretized integral operators. The proposed method is a variant of the Martinsson-Rokhlin-type fast direct solver. Our variant avoids the explicit computation of the inverse of the coefficient matrix, thereby reducing the required number of matrix-matrix multiplications. Numerical experiments demonstrate that the proposed solver has a complexity of $O(N)$ in the low-frequency range and has a highly parallel computation efficiency with a strong scaling efficiency of 70\%. Furthermore, multiple right-hand sides can be solved efficiently; specifically, when solving problems with 180 right-hand side vectors, the processing time per vector from the second vector onward was approximately 28,900 times faster than that for the first vector. This is a key advantage of fast direct methods.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to develop a fast direct solver for two - dimensional elastic wave scattering problems. Specifically, the paper proposes a fast direct solver with O(N) complexity based on the proxy surface method to solve the two - dimensional elastic wave scattering problem in cavities. Due to the density of the coefficient matrix, the traditional boundary element method (BEM) has a high time complexity, usually O(N^3), when dealing with large - scale problems, which limits its efficiency in practical applications. The proposed fast direct solver avoids the explicit calculation of the inverse of the coefficient matrix through the low - rank approximation technique, thereby reducing the number of required matrix - matrix multiplications and achieving a solution efficiency with linear time complexity O(N).
### Main contribution points:
1. **Low - rank approximation**: Extended to elastodynamics through the proxy surface method, shared coefficients are extracted from the discretized integral operator to achieve low - rank approximation.
2. **Computational efficiency**: Numerical experiments show that this solver has a time complexity of O(N) in the low - frequency range and has a high parallel computing efficiency, with a strong scaling efficiency reaching 70%.
3. **Multiple right - hand - side problems**: For multiple right - hand - side problems, starting from the second vector, the processing time of each vector is approximately 28,900 times faster than that of the first vector, which is a key advantage of the fast direct method.
4. **Avoidance of spurious eigenfrequencies**: Using the Burton - Miller - type integral equation, the problem of spurious eigenfrequencies existing in the traditional boundary integral equation is solved, improving the reliability of the numerical results.
### Method overview:
- **Proxy surface method**: By introducing a local virtual boundary (proxy surface), the far - field interaction is approximated as a low - rank matrix, thereby reducing the computational complexity.
- **Galerkin discretization**: The Galerkin method is used to discretize the boundary integral equation to handle the hypersingular kernel problem and ensure high - precision solutions.
- **Binary tree decomposition**: The computational domain is decomposed into a binary tree to form a hierarchical structure, facilitating the recursive calculation of low - rank approximation.
- **Fast direct solver**: Based on a variant of the Martinsson - Rokhlin - type fast direct solver, it avoids the explicit calculation of the inverse matrix, reduces the number of matrix - matrix multiplications, and improves the computing speed.
### Application background:
- **Seismic wave analysis**: In seismic wave propagation analysis, the structure is large relative to the wavelength, and linear equations with a large number of degrees of freedom need to be processed, so an efficient solver is very necessary.
- **Nondestructive testing**: When using ultrasonic waves to detect internal defects in materials, an efficient and accurate numerical method is also required.
In conclusion, this paper proposes an efficient fast direct solver, which is suitable for two - dimensional elastic wave scattering problems, especially showing significant advantages when dealing with large - scale problems.