Asynchronous Teams and Tasks in a Message Passing Environment

As the discipline of scientific computing grows, so too does the "skills gap" between the increasingly complex scientific applications and the efficient algorithms required. Increasing demand for computational power on the march towards exascale requires innovative approaches. Closing the skills gap avoids the many pitfalls that lead to poor utilisation of resources and wasted investment. This thesis tackles two challenges: asynchronous algorithms for parallel computing and fault tolerance. First I present a novel asynchronous task invocation methodology for Discontinuous Galerkin codes called enclave tasking. The approach modifies the parallel ordering of tasks that allows for efficient scaling on dynamic meshes up to 756 cores. It ensures high levels of concurrency and intermixes tasks of different computational properties. Critical tasks along domain boundaries are prioritised for an overlap of computation and communication. The second contribution is the teaMPI library, forming teams of MPI processes exchanging consistency data through an asynchronous "heartbeat". In contrast to previous approaches, teaMPI operates fully asynchronously with reduced overhead. It is also capable of detecting individually slow or failing ranks and inconsistent data among replicas. Finally I provide an outlook into how asynchronous teams using enclave tasking can be combined into an advanced team-based diffusive load balancing scheme. Both concepts are integrated into and contribute towards the ExaHyPE project, a next generation code that solves hyperbolic equation systems on dynamically adaptive cartesian grids

Numerical simulations of seismo-acoustic nuisance patterns from an induced M1.8 earthquake in the Helsinki, southern Finland, metropolitan area

Irritating earthquake sounds, reported also at low ground shaking levels, can negatively impact the social acceptance of geo-engineering applications. Concurringly, earthquake sound patterns have been linked to faulting mechanisms, thus opening possibilities for earthquake source characterisation. Inspired by consistent reports of felt and heard disturbances associated with the weeks-long stimulation of a 6 km-deep geothermal system in 2018 below the Otaniemi district of Espoo, Helsinki, we conduct fully-coupled 3D numerical simulations of wave propagation in solid Earth and the atmosphere. We assess the sensitivity of ground shaking and audible noise distributions to the source geometry of small induced earthquakes, using the largest recorded event in 2018 of magnitude ML=1.8. Utilizing recent computational advances, we are able to model seismo-acoustic frequencies up to 25 Hz therefore reaching the lower limit of human sound sensitivity. Our results provide for the first time synthetic spatial nuisance distributions of audible sounds at the 50-100 m scale for a densely populated metropolitan region. In five here presented 3D coupled elastic-acoustic scenario simulations, we include the effects of topography and subsurface structure, and analyse the audible loudness of earthquake generated acoustic waves. We can show that in our region of interest, S-waves are generating the loudest sound disturbance. We compare our sound nuisance distributions to commonly used empirical relationships using statistical analysis. We find that our 3D synthetics are generally smaller than predicted empirically, and that the interplay of source-mechanism specific radiation pattern and topography can add considerable non-linear contributions. Our study highlights the complexity and information content of spatially variable audible effects, even when caused by very small earthquakes.Comment: 29 pages, 9 figures. This paper has been submitted to the Bulletin of the Seismological Society of America for publicatio

Verification of an ADER-DG method for complex dynamic rupture problems

We present results of thorough benchmarking of an arbitrary high-order derivative discontinuous Galerkin (ADER-DG) method on unstructured meshes for advanced earthquake dynamic rupture problems. We verify the method by comparison to well-established numerical methods in a series of verification exercises, including dipping and branching fault geometries, heterogeneous initial conditions, bimaterial interfaces and several rate-and-state friction laws. We show that the combination of meshing flexibility and high-order accuracy of the ADER-DG method makes it a competitive tool to study earthquake dynamics in geometrically complicated setups