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.

Implementation of sparse matrix algorithms in an advection–diffusion–chemistry module

SPARSE SQP ALGORITHM FOR CHEMICAL PROCESS OPTIMIZATION

Smoldyn on Graphics Processing Units: Massively Parallel Brownian Dynamics Simulations

An open and parallel multiresolution framework using block-based adaptive grids

On Solving Groundwater Flow and Transport Models with Algebraic Multigrid Preconditioning

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

An Efficient Algorithm for Astrochemical Systems Using Stoichiometry Matrices

Development and acceleration of parallel chemical transport models

Accelerating and parallelizing Lagrangian simulations of mixing-limited reactive transport

Using SIMD and SIMT vectorization to evaluate sparse chemical kinetic Jacobian matrices and thermochemical source terms

Splitting Approach for Solving Multi-Component Transport Models with Maxwell-Stefan-Diffusion

A modular massively parallel computing environment for three-dimensional multiresolution simulations of compressible flows

Sparse grid time-discontinuous Galerkin method with streamline diffusion for transport equations

Parallel SPH modeling using dynamic domain decomposition and load balancing displacement of Voronoi subdomains

A Semi-Lagrangian Adaptive-Rank (SLAR) Method for Linear Advection and Nonlinear Vlasov-Poisson System

Scalable explicit implementation of anisotropic diffusion with Runge-Kutta-Legendre super-time-stepping

A scalable Euler–Lagrange approach for multiphase flow simulation on spectral elements

Parallel Domain Decomposition Methods with Mixed Order Discretization for Fully Implicit Solution of Tracer Transport Problems on the Cubed-Sphere.

Parallel-in-time Multilevel Krylov Methods: A Prototype

Performance Analysis of a Multicore Implementation for Solving a Two-Dimensional Inverse Anomalous Diffusion Problem