Weighted-Norm Preconditioners for a Multilayer Tide Model
Abstract
We derive a linearized rotating shallow water system modeling tides, which can be discretized by mixed finite elements. Unlike previous models, this model allows for multiple layers stratified by density. Like the single-layer case~\cite{kirby2021preconditioning} a weighted-norm preconditioner gives a (nearly) parameter-robust method for solving the resulting linear system at each time step, but the all-to-all coupling between the layers in the model poses a significant challenge to efficiency. Neglecting the inter-layer coupling gives a preconditioner that degrades rapidly as the number of layers increases. By a careful analysis of the matrix that couples the layers, we derive a robust method that requires solving a reformulated system that only involves coupling between adjacent layers. Numerical results obtained using Firedrake confirm the theory.
- Publication:
-
SIAM Journal on Scientific Computing
- Pub Date:
- August 2023
- DOI:
- 10.1137/22M150753X
- arXiv:
- arXiv:2207.02116
- Bibcode:
- 2023SJSC...45A1789C
- Keywords:
-
- Mathematics - Numerical Analysis;
- 65F08;
- 65N30;
- G.1.8;
- G.1.3
- E-Print:
- 23 pages, 7 figures