Hybridised multigrid preconditioners for a compatible finite-element dynamical core

Jack Betteridge

doi:10.1002/qj.4515

Hybridised multigrid preconditioners for a compatible finite-element dynamical core

2023-07-12

1

OA

AI

DOI：10.1002/qj.4515

Original request for help

Save

Abstract

En 中文

Compatible finite-element discretisations for the atmospheric equations of motion have recently attracted considerable interest. Semi-implicit timestepping methods require the repeated solution of a large saddle-point system of linear equations. Preconditioning this system is challenging, since the velocity mass matrix is nondiagonal, leading to a dense Schur complement. Hybridisable discretisations overcome this issue: weakly enforcing continuity of the velocity field with Lagrange multipliers leads to a sparse system of equations, which has a similar structure to the pressure Schur complement in traditional approaches. We describe how the hybridised sparse system can be preconditioned with a non-nested two-level preconditioner. To solve the coarse system, we use the multigrid pressure solver that is employed in the approximate Schur complement method previously proposed by the some of the authors. Our approach significantly reduces the number of solver iterations. The method shows excellent performance and scales to large numbers of cores in the Met Office next-generation climate and weather prediction model LFRic.

Keywords:

atmosphere

dynamics

finite-element discretisation

hybridisation

linear solvers

multigrid

numerical methods and NWP

parallel scalability

Hybridised multigrid preconditioners for a compatible finite-element dynamical core

Abstract

Journal

Organization