Parallel Octree-Based Finite Element Method for Large-Scale Earthquake Ground Motion Simulation

Jacobo Bielak, & Omar Ghattas

Published 2005, SCEC Contribution #947

We present a parallel octree-based finite element method for large-scale earthquake ground motion simulation in realistic basins. The octree representation combines the low memory per node and good cache performance of finite difference methods with the spatial adaptivity to local seismic wavelengths characteristic of unstructured finite element methods. Several tests are provided to verify the numerical performance of the method against Green's function solutions for homogeneous and piecewise homogeneous media, both with and without anelastic attenuation. A comparison is also provided against a finite difference code and an unstructured tetrahedral finite element code for a simulation of the 1994 Northridge Earthquake. The numerical tests all show very good agreement with analytical solutions and other codes. Finally, performance evaluation indicates excellent single-processor performance and parallel scalability over a range of 1 to 2048 processors for Northridge simulations with up to 300 million degrees of freedom.

Citation
Bielak, J., & Ghattas, O. (2005). Parallel Octree-Based Finite Element Method for Large-Scale Earthquake Ground Motion Simulation. Computer Modeling in Engineering and Sciences, 10(2), 99-112.