Abstract:The demand for substantial increases in the spatial resolution of global weather and climate prediction models makes it necessary to use numerically efficient and highly scalable algorithms to solve the equations of large‐scale atmospheric fluid dynamics. For stability and efficiency reasons, several of the operational forecasting centres, in particular the Met Office and the European Centre for Medium‐Range Weather Forecasts (ECMWF) in the UK, use semi‐implicit semi‐Lagrangian time‐stepping in the dynamical core of the model. The additional burden with this approach is that a three‐dimensional elliptic partial differential equation (PDE) for the pressure correction has to be solved at every model time step and this often constitutes a significant proportion of the time spent in the dynamical core. In global models, this PDE must be solved in a thin spherical shell. To run within tight operational time‐scales, the solver has to be parallelized and there seems to be a (perceived) misconception that elliptic solvers do not scale to large processor counts and hence implicit time‐stepping cannot be used in very high‐resolution global models. After reviewing several methods for solving the elliptic PDE for the pressure correction and their application in atmospheric models, we demonstrate the performance and very good scalability of Krylov subspace solvers and multigrid algorithms for a representative model equation with more than 10 10 unknowns on 65 536 cores on the High‐End Computing Terascale Resource (HECToR), the UK's national supercomputer. For this, we tested and optimized solvers from two existing numerical libraries (the Distributed and Unified Numerics Environment (DUNE) and Parallel High Performance Preconditioners (hypre)) and implemented both a conjugate gradient solver and a geometric multigrid algorithm based on a tensor‐product approach, which exploits the strong vertical anisotropy of the discretized equation. We study both weak and strong scalability and compare the absolute solution times for all methods; in contrast to one‐level methods, the multigrid solver is robust with respect to parameter variations.

Massively parallel solvers for elliptic partial differential equations in numerical weather and climate prediction

High level implementation of geometric multigrid solvers for finite element problems: Applications in atmospheric modelling

Multigrid preconditioners for the mixed finite element dynamical core of the LFRic atmospheric model

Improving Met Office Weather and Climate Forecasts with Bespoke Multigrid Solvers

Parallel Semi-Implicit Spectral Element Methods for Atmospheric General Circulation Models

10M-core scalable fully-implicit solver for nonhydrostatic atmospheric dynamics.

A Scalable Fully Implicit Compressible Euler Solver for Mesoscale Nonhydrostatic Simulation of Atmospheric Flows.

Parallel Preconditioned GMRES Solvers for 3-D Helmholtz Equations in Regional Non-hydrostatic Atmosphere Model

Accelerating the 3D Euler Atmospheric Solver Through Heterogeneous CPU-GPU Platforms

Scalability of high-performance PDE solvers

Hybridised multigrid preconditioners for a compatible finite element dynamical core

Semi-Lagrangian shallow water equations solver on the cubed-sphere grid as a prototype of new-generation global atmospheric model

A GPU-accelerated Cartesian grid method is proposed for solving the heat, wave, and Schrodinger equations on irregular domains

A Scalable Helmholtz Solver in GRAPES over Large‐scale Multicore Cluster

Helmholtz Solving and Performance Optimization in Global/Regional Assimilation and Prediction System

Parallel-in-time Multilevel Krylov Methods: A Prototype

Scalable robust solvers for unstructured FE geodynamic modeling applications: Solving the Stokes equation for models with large localized viscosity contrasts

Matrix-free multigrid block-preconditioners for higher order Discontinuous Galerkin discretisations

Strong Scaling for Numerical Weather Prediction at Petascale with the Atmospheric Model NUMA

Parallel Multilevel Methods for Implicit Solution of Shallow Water Equations with Nonsmooth Topography on the Cubed-Sphere

Towards a Scalable Hierarchical High-order CFD Solver