Abstract

Compatible finite elements provide a framework for preserving important structures in equations of geophysical fluid dynamics, and are becoming important in their use for building atmosphere and ocean models. We survey the application of compatible finite element spaces to geophysical fluid dynamics, including the application to the nonlinear rotating shallow water equations, and the three-dimensional compressible Euler equations. We summarise analytic results about dispersion relations and conservation properties, and present new results on approximation properties in three dimensions on the sphere, and on hydrostatic balance properties.

Highlights

  • Finite element methods have recently become a popular discretization approach for numerical weather prediction, mostly using spectral elements or discontinuous Galerkin methods (Bao et al, 2015; Brdar et al, 2013; Dennis et al, 2011; Fournier et al, 2004; Giraldo et al 2013; Kelly and Giraldo, 2012; Marras et al, 2013; Thomas and Loft, 2005); the use of these methods is reviewed in the work by Marras et al (2015)

  • Compatible finite element methods are built from finite element spaces that have differential operators that map from one space to another

  • They have not been considered much in the numerical weather prediction context, the lowest Raviart–Thomas element has been proposed for ocean modelling (Walters and Casulli, 1998), and the RT0 and lowest order Brezzi–Douglas–Marini element were both analysed in Rostand and Le Roux (2008) as part of the quest for a mixed finite element pair with good dispersion properties when applied to geophysical fluid dynamics

Read more

Summary

Introduction

Finite element methods have recently become a popular discretization approach for numerical weather prediction, mostly using spectral elements or discontinuous Galerkin methods (Bao et al, 2015; Brdar et al, 2013; Dennis et al, 2011; Fournier et al, 2004; Giraldo et al 2013; Kelly and Giraldo, 2012; Marras et al, 2013; Thomas and Loft, 2005); the use of these methods is reviewed in the work by Marras et al (2015). Finite element methods provide the opportunity to use more general grids in numerical weather prediction, either to improve load balancing in massively parallel computation or to facilitate adaptive mesh refinement They allow the development of higher order discretizations. Compatible finite element methods are built from finite element spaces that have differential operators that map from one space to another They have a long history in both numerical analysis (see Boffi et al, 2013, for a summary of contributions) and applications, including aerodynamics, structural mechanics, electromagnetism and porous media flows.

Hilbert complexes
Compatible finite element spaces in two dimensions
Compatible finite element methods for the linear shallow water equations
Stability of pressure gradient operator
Compatible finite element methods for the non-linear shallow water equations
Three-dimensional compatible finite element spaces
Tensor product finite element spaces
Approximation properties and curvilinear elements
Representation of temperature in compatible finite element spaces
Hydrostatic balance
Summary and outlook
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call