Journal Article FZJ-2022-02207

http://join2-wiki.gsi.de/foswiki/pub/Main/Artwork/join2_logo100x88.png
A suite of Richardson preconditioners for semi-implicit all-scale atmospheric models

 ;

2022
Elsevier Amsterdam

Journal of computational physics 463, 111296 () [10.1016/j.jcp.2022.111296]

This record in other databases:  

Please use a persistent id in citations:   doi:

Abstract: The paper documents a suite of preconditioners for Krylov-subspace solvers of elliptic boundary-value problems (BVPs) that underlie semi-implicit integrations of the nonhydrostatic equations governing the dynamics of all-scale atmospheric flows. Effective preconditioning of the linear operators inherent in the semi-implicit models lies at the heart of the state-of-the-art multiscale-flow simulation. This is especially evident in simulations of global weather and climate—posed on a thin spherical shell—where some form of direct tridiagonal inversion of the operator in the vertical is crucial to relax the often enormous stiffness of the problem. The documented preconditioners stem from the Richardson's (1910) idea of augmenting an elliptic BVP with a transient diffusion equation. Exploiting this idea for mixed explicit-implicit pseudo-time-stepping schemes leads to a broad suite of stationary-iteration solvers, including the many classical algorithms. Here, the high-performance all-scale EULAG model (Smolarkiewicz et al. (2014) [58]), with a flexible three-dimensional decomposition of MPI tasks, is furnished with the preconditioners akin to the classical alternating-direction-implicit (ADI) algorithms, generalized to optional permutations of parallel tridiagonal inversions. The utility of various options is found to be problem dependent, in terms of computational accuracy as well as efficiency. The main thrust of the work is on the long-range forecasts using large anisotropic grids. The relative efficiency and/or accuracy gains attainable with the developed preconditioners are illustrated for idealized scenarios representative of atmospheric flows from planetary to a single-cloud and laboratory scales. The key insight that best encapsulates the significance and novelty of the present work is that there is no single “super-preconditioner” that will perform best in all cases, yet the suite as a whole offers substantial gains in the model performance.

Classification:

Contributing Institute(s):
  1. Jülich Supercomputing Center (JSC)
Research Program(s):
  1. 5111 - Domain-Specific Simulation & Data Life Cycle Labs (SDLs) and Research Groups (POF4-511) (POF4-511)

Appears in the scientific report 2022
Database coverage:
Medline ; OpenAccess ; Clarivate Analytics Master Journal List ; Current Contents - Physical, Chemical and Earth Sciences ; Ebsco Academic Search ; Essential Science Indicators ; IF < 5 ; JCR ; NationallizenzNationallizenz ; SCOPUS ; Science Citation Index Expanded ; Web of Science Core Collection
Click to display QR Code for this record

The record appears in these collections:
Document types > Articles > Journal Article
Workflow collections > Public records
Institute Collections > JSC
Publications database
Open Access

 Record created 2022-05-19, last modified 2023-01-23


OpenAccess:
Download fulltext PDF
External link:
Download fulltextFulltext by OpenAccess repository
Rate this document:

Rate this document:
1
2
3
 
(Not yet reviewed)