A three-dimensional hybrid finite element -- spectral boundary integral method for modeling earthquakes in complex unbounded domains

DOI10.1002/NME.6816arXiv2102.08756OpenAlexW3194109659MaRDI QIDQ6360862FDOQ6360862

Authors: Gabriele Albertini, Ahmed Elbanna, David S. Kammer 

Publication date: 15 February 2021

Abstract: We present a 3D hybrid method which combines the Finite Element Method (FEM) and the Spectral Boundary Integral method (SBIM) to model nonlinear problems in unbounded domains. The flexibility of FEM is used to model the complex, heterogeneous, and nonlinear part -- such as the dynamic rupture along a fault with near fault plasticity -- and the high accuracy and computational efficiency of SBIM is used to simulate the exterior half spaces perfectly truncating all incident waves. The exact truncation allows us to greatly reduce the domain of spatial discretization compared to a traditional FEM approach, leading to considerable savings in computational cost and memory requirements. The coupling of FEM and SBIM is achieved by the exchange of traction and displacement boundary conditions at the computationally defined boundary. The method is suited to implementation on massively parallel computers. We validate the developed method by means of a benchmark problem. Three more complex examples with a low velocity fault zone, low velocity off-fault inclusion, and interaction of multiple faults, respectively, demonstrate the capability of the hybrid scheme in solving problems of very large sizes. Finally, we discuss potential applications of the hybrid method for problems in geophysics and engineering.

Full work available at URL: https://doi.org/10.1002/nme.6816

Recommendations

• Simulation of earthquake rupture dynamics in complex geometries using coupled finite difference and finite volume methods
• Dynamic earthquake rupture simulations on nonplanar faults embedded in 3D geometrically complex, heterogeneous elastic solids
• Low and high order finite element method: experience in seismic modeling
• High-performance finite-element simulations of seismic wave propagation in three-dimensional nonlinear inelastic geological media
• A nitsche-extended finite element method for earthquake rupture on complex fault systems

Mathematics Subject Classification ID

Fracture and damage (74Rxx) Numerical and other methods in solid mechanics (74Sxx) Special subfields of solid mechanics (74Lxx)

This page was built for publication: A three-dimensional hybrid finite element -- spectral boundary integral method for modeling earthquakes in complex unbounded domains

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6360862)