Loading…

A finite element method for solving the shallow water equations on the sphere

Within the framework of ocean general circulation modeling, the present paper describes an efficient way to discretize partial differential equations on curved surfaces by means of the finite element method on triangular meshes. Our approach benefits from the inherent flexibility of the finite eleme...

Full description

Saved in:
Bibliographic Details
Published in:Ocean modelling (Oxford) 2009, Vol.28 (1), p.12-23
Main Authors: Comblen, Richard, Legrand, Sébastien, Deleersnijder, Eric, Legat, Vincent
Format: Article
Language:English
Subjects:
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:Within the framework of ocean general circulation modeling, the present paper describes an efficient way to discretize partial differential equations on curved surfaces by means of the finite element method on triangular meshes. Our approach benefits from the inherent flexibility of the finite element method. The key idea consists in a dialog between a local coordinate system defined for each element in which integration takes place, and a nodal coordinate system in which all local contributions related to a vectorial degree of freedom are assembled. Since each element of the mesh and each degree of freedom are treated in the same way, the so-called pole singularity issue is fully circumvented. Applied to the shallow water equations expressed in primitive variables, this new approach has been validated against the standard test set defined by [Williamson, D.L., Drake, J.B., Hack, J.J., Jakob, R., Swarztrauber, P.N., 1992. A standard test set for numerical approximations to the shallow water equations in spherical geometry. Journal of Computational Physics 102, 211–224]. Optimal rates of convergence for the P 1 NC - P 1 finite element pair are obtained, for both global and local quantities of interest. Finally, the approach can be extended to three-dimensional thin-layer flows in a straightforward manner.
ISSN:1463-5003
1463-5011
DOI:10.1016/j.ocemod.2008.05.004