Preconditioning mixed finite elements for tide models
Numerical Analysis
2020-03-04 v1 Numerical Analysis
Abstract
We describe a fully discrete mixed finite element method for the linearized rotating shallow water model, possibly with damping. While Crank-Nicolson time-stepping conserves energy in the absence of drag or forcing terms and is not subject to a CFL-like stability condition, it requires the inversion of a linear system at each step. We develop weighted-norm preconditioners for this algebraic system that are nearly robust with respect to the physical and discretization parameters in the system. Numerical experiments using Firedrake support the theoretical results.
Keywords
Cite
@article{arxiv.2003.01632,
title = {Preconditioning mixed finite elements for tide models},
author = {Tate Kernell and Robert C. Kirby},
journal= {arXiv preprint arXiv:2003.01632},
year = {2020}
}