The ‘recovered space’ advection scheme for lowest-order compatible finite element methods

This is the author accepted manuscript. The final version is available from the publisher via the DOI in this record.We present a new compatible finite element advection scheme for the compressible Euler equations. Unlike the discretisations described in Cotter and Kuzmin (2016) and Shipton et al. (2018), the discretisation uses the lowest-order family of compatible finite element spaces, but still retains second-order numerical accuracy. This scheme obtains this second-order accuracy by first ‘recovering’ the function in higher-order spaces, before using the discontinuous Galerkin advection schemes of Cotter and Kuzmin (2016). As well as describing the scheme, we also present its stability properties and a strategy for ensuring boundedness. We then demonstrate its properties through some numerical tests, before presenting its use within a model solving the compressible Euler equations.Engineering and Physical Science Research Council (EPSRC

On coupling resolved and unresolved physical processes in finite element discretisations of geophysical fluids

At the heart of modern numerical weather forecasting and climate modelling lie simulations of two geophysical fluids: the atmosphere and the ocean. These endeavours rely on numerically solving the equations that describe these fluids. A key challenge is that the fluids contain motions spanning a range of scales. As the small-scale processes (unresolved by the numerical model) affect the resolved motions, they need to be described in the model, which is known as parametrisation. One major class of methods for numerically solving such partial differential equations is the finite element method. This thesis focuses on the coupling of such parametrised processes to the resolved flow within finite element discretisations. Four sets of research are presented, falling under two main categories. The first is the development of a compatible finite element discretisation for use in numerical weather prediction models, so as to avoid the bottleneck in computational scalability associated with the convergence at the poles of latitude-longitude grids. We present a transport scheme for use with the lowest-order function spaces in such a compatible finite element method, which is motivated by the coupling of the resolved and unresolved processes within the model. This then facilitates the use of the lower-order spaces within Gusto, a toolkit for studying such compatible finite element discretisations. Then, we present a compatible finite element discretisation of the moist compressible Euler equations, parametrising the unresolved moist processes. This is a major step in the development of Gusto, extending it to describe its first unresolved processes. The second category with which this thesis is concerned is the stochastic variational framework presented by Holm [Variational principles for stochastic fluid dynamics, P. Roy. Soc. A-Math. Phy. 471 (2176), (2015)]. In this framework, the effect of the unresolved processes and their uncertainty is expressed through a stochastic component to the advecting velocity. This framework ensures the circulation theorem is preserved by the stochastic equations. We consider the application of this formulation to two simple geophysical fluid models. First, we discuss the statistical properties of an enstrophy-preserving finite element discretisation of the stochastic quasi-geostrophic equation. We find that the choice of discretisation and the properties that it preserves affects the statistics of the solution. The final research presented is a finite element discretisation of the stochastic Camassa-Holm equation, which is used to numerically investigate the formation of ‘peakons’ within this set-up, finding that they do still always form despite the noise’s presence.Open Acces

Energy conserving SUPG methods for compatible finite element schemes in numerical weather prediction

We present an energy conserving space discretisation based on a Poisson bracket that can be used to derive the dry compressible Euler as well as thermal shallow water equations. It is formulated using the compatible finite element method, and extends the incorporation of upwinding for the shallow water equations as described in Wimmer, Cotter, and Bauer (2019). While the former is restricted to DG upwinding, an energy conserving SUPG scheme for the (partially) continuous Galerkin thermal field space is newly introduced here. The energy conserving property is validated by coupling the Poisson bracket based spatial discretisation to an energy conserving time discretisation. Further, the discretisation is demonstrated to lead to an improved temperature field development with respect to stability when upwinding is included. An approximately energy conserving full discretisation with a smaller computational cost is also presented.Comment: 27 pages, 9 figures, first version: all comments welcom

Compatible finite element methods for geophysical fluid dynamics

This article surveys research on the application of compatible finite element methods to large scale atmosphere and ocean simulation. Compatible finite element methods extend Arakawa's C-grid finite difference scheme to the finite element world. They are constructed from a discrete de Rham complex, which is a sequence of finite element spaces which are linked by the operators of differential calculus. The use of discrete de Rham complexes to solve partial differential equations is well established, but in this article we focus on the specifics of dynamical cores for simulating weather, oceans and climate. The most important consequence of the discrete de Rham complex is the Hodge-Helmholtz decomposition, which has been used to exclude the possibility of several types of spurious oscillations from linear equations of geophysical flow. This means that compatible finite element spaces provide a useful framework for building dynamical cores. In this article we introduce the main concepts of compatible finite element spaces, and discuss their wave propagation properties. We survey some methods for discretising the transport terms that arise in dynamical core equation systems, and provide some example discretisations, briefly discussing their iterative solution. Then we focus on the recent use of compatible finite element spaces in designing structure preserving methods, surveying variational discretisations, Poisson bracket discretisations, and consistent vorticity transport.Comment: correction of some typo

Postprocessing of Non-Conservative Flux for Compatibility with Transport in Heterogeneous Media

A conservative flux postprocessing algorithm is presented for both steady-state and dynamic flow models. The postprocessed flux is shown to have the same convergence order as the original flux. An arbitrary flux approximation is projected into a conservative subspace by adding a piecewise constant correction that is minimized in a weighted $L^2$ norm. The application of a weighted norm appears to yield better results for heterogeneous media than the standard $L^2$ norm which has been considered in earlier works. We also study the effect of different flux calculations on the domain boundary. In particular we consider the continuous Galerkin finite element method for solving Darcy flow and couple it with a discontinuous Galerkin finite element method for an advective transport problem.Comment: 34 pages, 17 figures, 11 table

Aspects of Unstructured Grids and Finite-Volume Solvers for the Euler and Navier-Stokes Equations

One of the major achievements in engineering science has been the development of computer algorithms for solving nonlinear differential equations such as the Navier-Stokes equations. In the past, limited computer resources have motivated the development of efficient numerical schemes in computational fluid dynamics (CFD) utilizing structured meshes. The use of structured meshes greatly simplifies the implementation of CFD algorithms on conventional computers. Unstructured grids on the other hand offer an alternative to modeling complex geometries. Unstructured meshes have irregular connectivity and usually contain combinations of triangles, quadrilaterals, tetrahedra, and hexahedra. The generation and use of unstructured grids poses new challenges in CFD. The purpose of this note is to present recent developments in the unstructured grid generation and flow solution technology

Hybridizable compatible finite element discretizations for numerical weather prediction: implementation and analysis

There is a current explosion of interest in new numerical methods for atmospheric modeling. A driving force behind this is the need to be able to simulate, with high efficiency, large-scale geophysical flows on increasingly more parallel computer systems. Many current operational models, including that of the UK Met Office, depend on orthogonal meshes, such as the latitude-longitude grid. This facilitates the development of finite difference discretizations with favorable numerical properties. However, such methods suffer from the ``pole problem," which prohibits the model to make efficient use of a large number of computing processors due to excessive concentration of grid-points at the poles. Recently developed finite element discretizations, known as ``compatible" finite elements, avoid this issue while maintaining the key numerical properties essential for accurate geophysical simulations. Moreover, these properties can be obtained on arbitrary, non-orthogonal meshes. However, the efficient solution of the resulting discrete systems depend on transforming the mixed velocity-pressure (or velocity-pressure-buoyancy) system into an elliptic problem for the pressure. This is not so straightforward within the compatible finite element framework due to inter-element coupling. This thesis supports the proposition that systems arising from compatible finite element discretizations can be solved efficiently using a technique known as ``hybridization." Hybridization removes inter-element coupling while maintaining the desired numerical properties. This permits the construction of sparse, elliptic problems, for which fast solver algorithms are known, using localized algebra. We first introduce the technique for compatible finite element discretizations of simplified atmospheric models. We then develop a general software abstraction for the rapid implementation and composition of hybridization methods, with an emphasis on preconditioning. Finally, we extend the technique for a new compatible method for the full, compressible atmospheric equations used in operational models.Open Acces