English

Weighted-norm preconditioners for a multi-layer tide model

Numerical Analysis 2022-07-06 v1 Numerical Analysis

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.

Keywords

Cite

@article{arxiv.2207.02116,
  title  = {Weighted-norm preconditioners for a multi-layer tide model},
  author = {Colin J. Cotter and Robert C. Kirby and Hunter Morris},
  journal= {arXiv preprint arXiv:2207.02116},
  year   = {2022}
}

Comments

23 pages, 7 figures

R2 v1 2026-06-24T12:14:40.175Z