A fast direct solver for two-dimensional transmission problems of elastic waves

This paper describes a fast direct boundary element method for elastodynamic transmission problems in two dimensions,which can be used for analyzing elastic wave scattering by an inclusion.This paper reports the development of an efficient solver based on a discretization method that isbroadly applicable regardless of the inclusion shape.From the smoothness of the solutions of the Navier--Cauchy equation,it is reasonable thatthe displacement is approximated by the piecewise linear basesand the traction is approximated by the piecewise constant bases.However, in this mixed bases strategy, Calder\'on preconditioning, thatis, an analytical preconditioning with excellent performance, cannot be applied,which means we cannot adopt iterative solvers.To address this challenge, we developed a fast direct solver formulated using both Burton--Miller and Poggio--Miller--Chang--Harrington--Wu--Tsai (PMCHWT) boundary integral equations.Our method uses a technique based on the proxy method for low-rank approximation of the coefficient matrix's off-diagonal blocks. To handle transmission problems, the proposed fast direct solver uses separate binary tree partitions for nodes and elements.Numerical examples demonstrate that our solver achieves linear computational complexityand can efficiently handle problems with multiple right-sides.Notably, the solver based on the Burton--Miller formulation is approximately 20% faster than the one using the PMCHWT formulation.Our new method provides a versatile, fast solver, whose performance is independent of the shape of inclusions and computational parameters, such as frequency and density, for elastodynamic transmission problems.