A family of mixed finite element pairs with optimal geostrophic balance

C. J. Cotter

DOI: https://doi.org/10.48550/arXiv.0905.0126

2009-05-02

Abstract: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.

Numerical Analysis

What problem does this paper attempt to address?

### Problems the paper attempts to solve

This paper aims to introduce a new type of mixed finite element pairs for geographical grids and adaptive grid refinement in numerical weather prediction and ocean models. Specifically, the paper attempts to solve the following problems:

1. **Optimizing geostrophic balance**: The authors prove that when these finite element pairs are applied to the linear rotating shallow - water equations, the geostrophic equilibrium state is completely stationary, which means that the numerical scheme does not introduce any spurious inertial - gravity waves. This makes these finite element pairs optimal in a sense in numerical weather prediction and ocean modeling applications.

2. **Avoiding spurious pressure modes**: By proving an inf - sup condition, the authors show that these finite element pairs do not have spurious pressure modes that would contaminate the numerical solutions over the long - time scales required for large - scale geophysical applications.

3. **Extension to three - dimensional rotating stratified incompressible flow**: The authors further extend these properties to the three - dimensional rotating stratified incompressible flow in non - hydrostatic ocean modeling and show that for the linearized equations, the equilibrium state remains completely stationary on any unstructured grid again.

4. **Optimal stencil properties of the discrete pressure Poisson equation**: The authors also show that the discrete pressure Poisson equation obtained by these discretization methods satisfies the optimal stencil properties.

5. **Application of unstructured grids**: All of these characteristics make these discretization methods excellent candidates for numerical weather prediction and large - scale ocean modeling applications that require unstructured grids.

### Key formulas

- Linearized shallow - water equations:
\[
u_t+fu^{\perp}+g\nabla\eta = 0,
\]
\[
\eta_t+\overline{D}\nabla\cdot u = 0.
\]

- Geostrophic balance condition:
\[
fu^{\perp}+g\nabla\eta = 0.
\]

- Dispersion relation:
\[
\omega(\omega^2 - f^2 - g\overline{D}(k^2 + l^2)) = 0.
\]

- Inf - sup condition:
\[
\sup_{w_{\delta}\in V}\frac{\int_{\Omega}w_{\delta}\cdot\nabla\phi_{\delta}\,dV}{\sqrt{\int_{\Omega}|w_{\delta}|^2\,dV}}\geq\beta\sqrt{\int_{\Omega}(\phi_{\delta})^2\,dV},\quad\forall\phi_{\delta}\in H,
\]
where \(\beta\) does not depend on the side lengths in the grid \(M\).

### Summary

This paper solves the common problems of spurious inertial - gravity waves and spurious pressure modes in numerical weather prediction and ocean modeling by introducing new mixed finite element pairs, ensures the stability of the geostrophic equilibrium state, and is applicable to any unstructured grid. These improvements are of great significance for improving the accuracy and stability of numerical simulations.