Newmark local time stepping on high-performance computing architectures

JOURNAL OF COMPUTATIONAL PHYSICS

卷:334

Article

Rietmann, Max^1,3  Grote, Marcus²  Peter, Daniel^1,3,4  Schenk, Olaf¹ 

[1] Univ Svizzera italiana, Inst Computat Sci, Lugano, Switzerland

[2] Univ Basel, Dept Math & Comp Sci, CH-4003 Basel, Switzerland

[3] Swiss Fed Inst Technol, Inst Geophys, Zurich, Switzerland

[4] King Abdullah Univ Sci & Technol, Div Phys Sci & Engn, Extreme Comp Res Ctr, Thuwal, Saudi Arabia

关键词: Seismic wave propagation;    Local-time stepping;    Multi-rate time stepping;   

DOI  :  10.1016/j.jcp.2016.11.012

来源: Elsevier

【 摘 要 】

In multi-scale complex media, finite element meshes often require areas of local refinement, creating small elements that can dramatically reduce the global time-step for wave propagation problems due to the CFL condition. Local time stepping (LTS) algorithms allow an explicit time-stepping scheme to adapt the time-step to the element size, allowing near-optimal time-steps everywhere in the mesh. We develop an efficient multilevel LTS-Newmark scheme and implement it in a widely used continuous finite element seismic wave-propagation package. In particular, we extend the standard LTS formulation with adaptations to continuous finite element methods that can be implemented very efficiently with very strong element-size contrasts (more than 100x). Capable of running on large CPU and GPU clusters, we present both synthetic validation examples and large scale, realistic application examples to demonstrate the performance and applicability of the method and implementation on thousands of CPU cores and hundreds of GPUs. (C) 2016 Elsevier Inc. All rights reserved.

【 授权许可】

Free