Skip to Main content Skip to Navigation
Journal articles

Assessing the robustness and scalability of the accelerated pseudo-transient method

Abstract : The development of highly efficient, robust and scalable numerical algorithms lags behind the rapid increase in massive parallelism of modern hardware. We address this challenge with the accelerated pseudo-transient (PT) iterative method and present a physically motivated derivation. We analytically determine optimal iteration parameters for a variety of basic physical processes and confirm the validity of theoretical predictions with numerical experiments. We provide an efficient numerical implementation of PT solvers on graphical processing units (GPUs) using the Julia language. We achieve a parallel efficiency of more than 96 % on 2197 GPUs in distributed-memory parallelisation weak-scaling benchmarks. The 2197 GPUs allow for unprecedented tera-scale solutions of 3D variable viscosity Stokes flow on 49953 grid cells involving over 1.2 trillion degrees of freedom (DoFs). We verify the robustness of the method by handling contrasts up to 9 orders of magnitude in material parameters such as viscosity and arbitrary distribution of viscous inclusions for different flow configurations. Moreover, we show that this method is well suited to tackle strongly nonlinear problems such as shear-banding in a visco-elasto-plastic medium. A GPU-based implementation can outperform direct-iterative solvers based on central processing units (CPUs) in terms of wall time, even at relatively low spatial resolution. We additionally motivate the accessibility of the method by its conciseness, flexibility, physically motivated derivation and ease of implementation. This solution strategy thus has a great potential for future high-performance computing (HPC) applications, and for paving the road to exascale in the geosciences and beyond.
Document type :
Journal articles
Complete list of metadata

https://hal-insu.archives-ouvertes.fr/insu-03752043
Contributor : Isabelle Dubigeon Connect in order to contact the contributor
Submitted on : Tuesday, August 16, 2022 - 11:10:47 AM
Last modification on : Thursday, September 1, 2022 - 8:44:06 AM

File

gmd-15-5757-2022.pdf
Publisher files allowed on an open archive

Licence


Distributed under a Creative Commons Attribution 4.0 International License

Identifiers

Citation

Ludovic Räss, Ivan Utkin, Thibault Duretz, Samuel Omlin, Yuri Y Podladchikov. Assessing the robustness and scalability of the accelerated pseudo-transient method. Geoscientific Model Development, European Geosciences Union, 2022, 15 (14), pp.5757 - 5786. ⟨10.5194/gmd-15-5757-2022⟩. ⟨insu-03752043⟩

Share

Metrics

Record views

21

Files downloads

1