A hybrid finite element-spectral boundary integral approach: Applications to dynamic rupture modeling in unbounded domains

The finite element method (FEM) and the spectral boundary integral method (SBI) have both been widely used in the study of dynamic rupture simulations along a weak interface. Here, we present a hybrid method that combines FEM and SBI through the consistent exchange of displacement and traction boundary conditions, thereby benefiting from the flexibility of FEM in handling problems with nonlinearities or small-scale heterogeneities and from the superior performance and accuracy of SBI. We validate the hybrid method using a benchmark problem from the Southern California Earthquake Center's dynamic rupture simulation validation exercises.We further demonstrate the capability and computational efficiency of the hybrid scheme for resolving off-fault heterogeneities by studying a 2D in-plane shear crack in two different settings: one where the crack is embedded in a high-velocity zone and another where it is embedded in a low-velocity zone. Finally, we discuss the potential of the hybrid method for addressing a wide range of problems in geophysics and engineering.