Heinlein, Alexander ORCID: 0000-0003-1578-8104, Perego, Mauro and Rajamanickam, Sivasankaran (2021). FROSch Preconditioners for Land Ice Simulations of Greenland and Antarctica. Technical Report.


Download (4MB) | Preview


Numerical simulations of Greenland and Antarctic ice sheets involve the solution of large-scale highly nonlinear systems of equations on complex shallow geometries. This work is concerned with the construction of Schwarz preconditioners for the solution of the associated tangent problems, which are challenging for solvers mainly because of the strong anisotropy of the meshes and wildly changing boundary conditions that can lead to poorly constrained problems on large portions of the domain. Here, two-level GDSW (Generalized Dryja–Smith–Widlund) type Schwarz preconditioners are applied to different land ice problems, i.e., a velocity problem, a temperature problem, as well as the coupling of the former two problems. We employ the MPI-parallel implementation of multi-level Schwarz preconditioners provided by the package FROSch (Fast and Robust Schwarz) from the Trilinos library. The strength of the proposed preconditioner is that it yields out-of-the-box scalable and robust preconditioners for the single physics problems. To our knowledge, this is the first time two-level Schwarz preconditioners are applied to the ice sheet problem and a scalable preconditioner has been used for the coupled problem. The preconditioner for the coupled problem differs from previous monolithic GDSW preconditioners in the sense that decoupled extension operators are used to compute the values in the interior of the subdomains. Several approaches for improving the performance, such as reuse strategies and shared memory OpenMP parallelization, are explored as well. In our numerical study we target both uniform meshes of varying resolution for the Antarctic ice sheet as well as non uniform meshes for the Greenland ice sheet are considered. We present several weak and strong scaling studies confirming the robustness of the approach and the parallel scalability of the FROSch implementation. Among the highlights of the numerical results are a weak scaling study for up to 32 K processor cores (8 K MPI-ranks and 4 OpenMP threads) and 566 M degrees of freedom for the velocity problem as well as a strong scaling study for up to 4 K processor cores (and MPI-ranks) and 68 M degrees of freedom for the coupled problem.

Item Type: Preprints, Working Papers or Reports (Technical Report)
CreatorsEmailORCIDORCID Put Code
Heinlein, Alexanderalexander.heinlein@uni-koeln.deorcid.org/0000-0003-1578-8104UNSPECIFIED
Perego, Mauromperego@sandia.govUNSPECIFIEDUNSPECIFIED
Rajamanickam, Sivasankaransrajama@sandia.govUNSPECIFIEDUNSPECIFIED
URN: urn:nbn:de:hbz:38-306684
Series Name at the University of Cologne: Technical report series. Center for Data and Simulation Science
Volume: 2021,1
Date: 25 January 2021
Language: English
Faculty: Central Institutions / Interdisciplinary Research Centers
Divisions: Weitere Institute, Arbeits- und Forschungsgruppen > Center for Data and Simulation Science (CDS)
Subjects: Natural sciences and mathematics
Technology (Applied sciences)
Uncontrolled Keywords:
domain decomposition methodsEnglish
monolithic Schwarz preconditionersEnglish
GDSW coarse spacesEnglish
multiphysics simulationsEnglish
parallel computingEnglish
Refereed: No
URI: http://kups.ub.uni-koeln.de/id/eprint/30668


Downloads per month over past year


Actions (login required)

View Item View Item