C. J. Cotter

C. J. Cotter, J. Thuburn

We describe discretisations of the shallow water equations on the sphere using the framework of finite element exterior calculus, which are extensions of the mimetic finite difference framework presented in Ringler, Thuburn, Klemp, and Skamarock (Journal of Computational Physics, 2010). The exterior calculus notation provides a guide to which finite element spaces should be used for which physical variables, and unifies a number of desirable properties. We present two formulations: a ``primal'' formulation in which the finite element spaces are defined on a single mesh, and a ``primal-dual'' formulation in which finite element spaces on a dual mesh are also used. Both formulations have velocity and layer depth as prognostic variables, but the exterior calculus framework leads to a conserved diagnostic potential vorticity. In both formulations we show how to construct discretisations that have mass-consistent (constant potential vorticity stays constant), stable and oscillation-free potential vorticity advection.

C. J. Cotter, J. Shipton

We show how two-dimensional mixed finite element methods that satisfy the conditions of finite element exterior calculus can be used for the horizontal discretisation of dynamical cores for numerical weather prediction on pseudo-uniform grids. This family of mixed finite element methods can be thought of in the numerical weather prediction context as a generalisation of the popular polygonal C-grid finite difference methods. There are a few major advantages: the mixed finite element methods do not require an orthogonal grid, and they allow a degree of flexibility that can be exploited to ensure an appropriate ratio between the velocity and pressure degrees of freedom so as to avoid spurious mode branches in the numerical dispersion relation. These methods preserve several properties of the C-grid method when applied to linear barotropic wave propagation, namely: a) energy conservation, b) mass conservation, c) no spurious pressure modes, and d) steady geostrophic modes on the $f$-plane. We explain how these properties are preserved, and describe two examples that can be used on pseudo-uniform grids: the recently-developed modified RT0-Q0 element pair on quadrilaterals and the BDFM1-\pdg element pair on triangles. All of these mixed finite element methods have an exact 2:1 ratio of velocity degrees of freedom to pressure degrees of freedom. Finally we illustrate the properties with some numerical examples.

C. J. Cotter, D. A. Ham, C. C. Pain et al.

We introduce a new mixed discontinuous/continuous Galerkin finite element for solving the 2- and 3-dimensional wave equations and equations of incompressible flow. The element, which we refer to as P1dg-P2, uses discontinuous piecewise linear functions for velocity and continuous piecewise quadratic functions for pressure. The aim of introducing the mixed formulation is to produce a new flexible element choice for triangular and tetrahedral meshes which satisfies the LBB stability condition and hence has no spurious zero-energy modes. We illustrate this property with numerical integrations of the wave equation in two dimensions, an analysis of the resultant discrete Laplace operator in two and three dimensions, and a normal mode analysis of the semi-discrete wave equation in one dimension.

J. Shipton, T. H. Gibson, C. J. Cotter

We describe a compatible finite element discretisation for the shallow water equations on the rotating sphere, concentrating on integrating consistent upwind stabilisation into the framework. Although the prognostic variables are velocity and layer depth, the discretisation has a diagnostic potential vorticity that satisfies a stable upwinded advection equation through a Taylor-Galerkin scheme; this provides a mechanism for dissipating enstrophy at the gridscale whilst retaining optimal order consistency. We also use upwind discontinuous Galerkin schemes for the transport of layer depth. These transport schemes are incorporated into a semi-implicit formulation that is facilitated by a hybridisation method for solving the resulting mixed Helmholtz equation. We illustrate our discretisation with some standard rotating sphere test problems.

C. J. Cotter

Diffeomorphic matching (only one of several names for this technique) is a technique for non-rigid registration of curves and surfaces in which the curve or surface is embedded in the flow of a time-series of vector fields. One seeks the flow between two topologically-equivalent curves or surfaces which minimises some metric defined on the vector fields, \emph{i.e.} the flow closest to the identity in some sense. In this paper, we describe a new particle-mesh discretisation for the evolution of the geodesic flow and the embedded shape. Particle-mesh algorithms are very natural for this problem because Lagrangian particles (particles moving with the flow) can represent the movement of the shape whereas the vector field is Eulerian and hence best represented on a static mesh. We explain the derivation of the method, and prove conservation properties: the discrete method has a set of conserved momenta corresponding to the particle-relabelling symmetry which converge to conserved quantities in the continuous problem. We also introduce a new discretisation for the geometric current matching condition of (Vaillant and Glaunes, 2005). We illustrate the method and the derived properties with numerical examples.

C. J. Cotter, D. Kuzmin

Motivated by finite element spaces used for representation of temperature in the compatible finite element approach for numerical weather prediction, we introduce locally bounded transport schemes for (partially-)continuous finite element spaces. The underlying high-order transport scheme is constructed by injecting the partially-continuous field into an embedding discontinuous finite element space, applying a stable upwind discontinuous Galerkin (DG) scheme, and projecting back into the partially-continuous space; we call this an embedded DG scheme. We prove that this scheme is stable in L2 provided that the underlying upwind DG scheme is. We then provide a framework for applying limiters for embedded DG transport schemes. Standard DG limiters are applied during the underlying DG scheme. We introduce a new localised form of element-based flux-correction which we apply to limiting the projection back into the partially-continuous space, so that the whole transport scheme is bounded. We provide details in the specific case of tensor-product finite element spaces on wedge elements that are discontinuous P1/Q1 in the horizontal and continuous P2 in the vertical. The framework is illustrated with numerical tests.

C. J. Cotter, D. A. Ham, C. C. Pain

We introduce a mixed discontinuous/continuous finite element pair for ocean modelling, with continuous quadratic pressure/layer depth 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 illustrate these properties using numerical tests and provide convergence calculations which show that the discretisation has quadratic errors, indicating that the element pair is stable.

C. J. Cotter, D. A. Ham, A. T. T. McRae et al.

We provide an approach to implementing the shallow atmosphere approximation in three dimensional finite element discretisations for dynamical cores. The approach makes use of the fact that the shallow atmosphere approximation metric can be obtained by writing equations on a three-dimensional manifold embedded in $\mathbb{R}^4$ with a restriction of the Euclidean metric. We show that finite element discretisations constructed this way are equivalent to the use of a modified three dimensional mesh for the construction of metric terms. We demonstrate our approach via a convergence test for a prototypical elliptic problem.

C. J. Cotter

Inertia-gravity mode and Rossby mode dispersion properties are examined for discretisations of the linearized rotating shallow-water equations using the $P1_{DG}$-$P2$ finite element pair on arbitrary triangulations in planar geometry. A discrete Helmholtz decomposition of the functions in the velocity space based on potentials taken from the pressure space is used to provide a complete description of the numerical wave propagation for the discretised equations. In the $f$-plane case, this decomposition is used to obtain decoupled equations for the geostrophic modes, the inertia-gravity modes, and the inertial oscillations. As has been noticed previously, the geostrophic modes are steady. The Helmholtz decomposition is used to show that the resulting inertia-gravity wave equation is third-order accurate in space. In general the \pdgp finite element pair is second-order accurate, so this leads to very accurate wave propagation. It is further shown that the only spurious modes supported by this discretisation are spurious inertial oscillations which have frequency $f$, and which do not propagate. The Helmholtz decomposition also allows a simple derivation of the quasi-geostrophic limit of the discretised $P1_{DG}$-$P2$ equations in the $β$-plane case, resulting in a Rossby wave equation which is also third-order accurate.

S. C. Kramer, C. J. Cotter, C. C. Pain

We discuss the ill conditioning of the matrix for the discretised Poisson equation in the small aspect ratio limit, and motivate this problem in the context of nonhydrostatic ocean modelling. Efficient iterative solvers for the Poisson equation in small aspect ratio domains are crucial for the successful development of nonhydrostatic ocean models on unstructured meshes. We introduce a new multigrid preconditioner for the Poisson problem which can be used with finite element discretisations on general unstructured meshes; this preconditioner is motivated by the fact that the Poisson problem has a condition number which is independent of aspect ratio when Dirichlet boundary conditions are imposed on the top surface of the domain. This leads to the first level in an algebraic multigrid solver (which can be extended by further conventional algebraic multigrid stages), and an additive smoother. We illustrate the method with numerical tests on unstructured meshes, which show that the preconditioner makes a dramatic improvement on a more standard multigrid preconditioner approach, and also show that the additive smoother produces better results than standard SOR smoothing. This new solver method makes it feasible to run nonhydrostatic unstructured mesh ocean models in small aspect ratio domains.

C. J. Cotter

We introduce a family of mixed finite element pairs for use on geodesic grids and with adaptive mesh refinement for numerical weather prediction and ocean modelling. We prove that when these finite element pairs are applied to the linear rotating shallow water equations, the geostrophically balanced states are exactly steady, which means that the numerical schemes do not introduce any spurious inertia-gravity waves; this makes these finite element pairs in some sense optimal for numerical weather prediction and ocean modelling applications. We further prove that these finite element pairs satisfy an inf-sup condition which means that they are free of spurious pressure modes which would pollute the numerical solution over the timescales required for large-scale geophysical applications. We then discuss the extension to incompressible Euler-Boussinesq equations with rotation, and show that for the linearised equations the balanced states are again exactly steady on arbitrary unstructured meshes. We also show that the discrete pressure Poisson equation resulting from these discretisations satisfies an optimal stencil property. All these properties make the discretisations in this family excellent candidates for numerical weather prediction and large-scale ocean modelling applications when unstructured grids are required.