Optimizing CAD and Mesh Generation Workflow for SeisSol

SeisSol is a simulation software for seismic wave propagation and earthquake scenarios. It solves the fully elastic wave equations in heterogeneous media. Incorporating dynamic rupture simulation it performs complex multiphysics earthquake simulations. To account for complicated geometries SeisSol uses a fully unstructured tetrahedral mesh. Recent publications [1], [2] have shown that SeisSol is able run at petascale performance on modern supercomputers. In this poster we present SeisSol workflow optimizations for three scalability-critical stages: CAD model generation, parallel mesh generation, and initialization of the mesh database in SeisSol. CAD model generation is automated through a combination of point cloud surface meshing techniques and discrete surface boolean opertions. On these CAD models the PUMGen tool couples Simmetrix1 unstructured parallel mesh generation to a novel parallel mesh file format. Using this format we were able to initialize a 3.6×10 element tetrahedral mesh within seconds on 9216 processes in SeisSol.