Journal
OCEAN MODELLING
Volume 26, Issue 1-2, Pages 86-90Publisher
ELSEVIER SCI LTD
DOI: 10.1016/j.ocemod.2008.09.002
Keywords
Shallow-water equations; Finite element method; Geostrophic balance; Pressure modes; Discontinuous Galerkin
Categories
Funding
- NERC [NE/C52101X/1]
- NERC [NE/F012594/1, ESM010001] Funding Source: UKRI
- Natural Environment Research Council [NE/F012594/1, NER/T/S/2002/00459, NE/C521036/1, NE/C52101X/1, ESM010001] Funding Source: researchfish
Ask authors/readers for more resources
We introduce a mixed discontinuous/continuous finite element pair for ocean modelling, with continuous quadratic layer thickness and discontinuous velocity. We investigate the finite element pair applied to the linear shallow-water equations on an f-plane. The element pair has the property that all geostrophically balanced states which strongly satisfy the boundary conditions have discrete divergence equal to exactly zero and hence are exactly steady states of the discretised equations. This means that the finite element pair has excellent geostrophic balance properties. We also show that the element pair applied to the non-rotating linear shallow-water equations does not have any spurious small eigenvalues. We illustrate these properties using numerical tests and provide convergence calculations which show that the numerical solutions have errors which decay quadratically with element edge length for both velocity and layer thickness. (c) 2008 Elsevier Ltd. All rights reserved.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available