2022/07/05 by Colin J. Cotter, Cotter, Colin J., Robert C. Kirby +3
Engineering · Physics and Astronomy · Earth and Planetary Sciences · #Advanced Numerical Methods in Computational Mathematics #Electromagnetic Scattering and Analysis #Geophysical and Geoelectrical Methods
paper · pdf · doi:10.48550/arxiv.2207.02116
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~\citekirby2021preconditioning 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.