Parallel Finite Element Earthquake Rupture Simulations on Large-scale Multicore Supercomputers
November 3, 2011
2:30 p.m.
Xingfu Wu
Abstract
In this talk, we use the 2008 Ms 8.0 Wenchuan earthquake occurred in Wenchuan county, Sichuan province in China on May 12th, 2008 as an example to address our earthquake rupture simulations. We integrate a 3D mesh generator into the simulation, and use MPI to parallelize the 3D mesh generator, illustrate an element-based partitioning scheme for explicit finite element methods, and based on the partitioning scheme and what we learned from our previous work, we implement our hybrid MPI/OpenMP finite element earthquake simulation code in order to not only achieve multiple levels of parallelism of the code but also to reduce the communication overhead of MPI within a multicore node by taking advantage of the shared address space and on-chip high inter-core bandwidth and low inter-core latency. We evaluate the hybrid MPI/OpenMP finite element earthquake rupture simulations on quad- and hex-core Cray XT 4/5 systems from Oak Ridge National Laboratory using the Southern California Earthquake Center (SCEC) benchmark TPV 210. Our experimental results indicate that the parallel finite element earthquake rupture simulation obtains the accurate output results and has good scalability on these Cray XT systems.