TY - GEN
T1 - Extreme scale multi-physics simulations of the tsunamigenic 2004 sumatra megathrust earthquake
AU - Uphoff, Carsten
AU - Rettenberger, Sebastian
AU - Bader, Michael
AU - Madden, Elizabeth H.
AU - Ulrich, Thomas
AU - Wollherr, Stephanie
AU - Gabriel, Alice Agnes
N1 - Publisher Copyright:
© 2017 Copyright held by the owner/author(s).
PY - 2017/11/12
Y1 - 2017/11/12
N2 - We presenta high-resolution simulation of the 2004 Sumatra-Andaman earthquake, including non-linear frictional failure on a megathrust-splay fault system. Our method exploits unstructured meshes capturing the complicated geometries in subduction zones that are crucial to understand large earthquakes and tsunami generation. These up-to-date largest and longest dynamic rupture simulations enable analysis of dynamic source effects on the seafloor displacements. To tackle the extreme size of this scenario an end-to-end optimization of the simulation code SeisSol was necessary. We implemented a new cache-aware wave propagation scheme and optimized the dynamic rupture kernels using code generation. We established a novel clustered local-time-stepping scheme for dynamic rupture. In total, we achieved a speed-up of 13.6 compared to the previous implementation. For the Sumatra scenario with 221 million elements this reduced the time-to-solution to 13.9 hours on 86,016 Haswell cores. Furthermore, we used asynchronous output to overlap I/O and compute time.
AB - We presenta high-resolution simulation of the 2004 Sumatra-Andaman earthquake, including non-linear frictional failure on a megathrust-splay fault system. Our method exploits unstructured meshes capturing the complicated geometries in subduction zones that are crucial to understand large earthquakes and tsunami generation. These up-to-date largest and longest dynamic rupture simulations enable analysis of dynamic source effects on the seafloor displacements. To tackle the extreme size of this scenario an end-to-end optimization of the simulation code SeisSol was necessary. We implemented a new cache-aware wave propagation scheme and optimized the dynamic rupture kernels using code generation. We established a novel clustered local-time-stepping scheme for dynamic rupture. In total, we achieved a speed-up of 13.6 compared to the previous implementation. For the Sumatra scenario with 221 million elements this reduced the time-to-solution to 13.9 hours on 86,016 Haswell cores. Furthermore, we used asynchronous output to overlap I/O and compute time.
KW - ADER-DG
KW - Asynchronous output
KW - Dynamic rupture
KW - Earthquake simulation
KW - Hybrid parallelization
KW - Local time stepping
KW - Petascale performance
KW - Tsunami coupling
UR - http://www.scopus.com/inward/record.url?scp=85040180617&partnerID=8YFLogxK
U2 - 10.1145/3126908.3126948
DO - 10.1145/3126908.3126948
M3 - Conference contribution
AN - SCOPUS:85040180617
T3 - Proceedings of the International Conference for High Performance Computing, Networking, Storage and Analysis, SC 2017
BT - Proceedings of the International Conference for High Performance Computing, Networking, Storage and Analysis, SC 2017
PB - Association for Computing Machinery, Inc
T2 - International Conference for High Performance Computing, Networking, Storage and Analysis, SC 2017
Y2 - 12 November 2017 through 17 November 2017
ER -