Abstract:We develop a novel and efficient discontinuous Galerkin spectral element method (DG-SEM) for the spherical rotating shallow water equations in vector invariant form. We prove that the DG-SEM is energy stable, and discretely conserves mass, vorticity, and linear geostrophic balance on general curvlinear meshes. These theoretical results are possible due to our novel entropy stable numerical DG fluxes for the shallow water equations in vector invariant form. We experimentally verify these results on a cubed sphere mesh. Additionally, we show that our method is robust, that is can be run stably without any dissipation. The entropy stable fluxes are sufficient to control the grid scale noise generated by geostrophic turbulence without the need for artificial stabilisation.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to develop a new discontinuous Galerkin spectral element method (DG - SEM) for solving the vector - invariant form of the rotating shallow - water equations (RSWE) on the sphere. Specifically, the author aims to prove that the proposed method can discretely preserve mass, vorticity, and linear geostrophic balance on general curvilinear grids and is energy - stable. In addition, the author also hopes to verify the robustness and stability of this method in practical applications, especially its ability to control grid - scale noise generated by geostrophic turbulence without introducing artificial dissipation.
### Decomposition of the core problems in the paper:
1. **Mass conservation**:
- The author hopes to ensure the local conservation of total mass through numerical methods.
- The formula is expressed as:
\[
\frac{d}{dt} \int_{\Omega} D \, d\Omega = 0
\]
2. **Conservation of absolute vorticity**:
- The conservation of absolute vorticity is crucial for accurately simulating mid - latitude weather dynamics.
- The formula is expressed as:
\[
\frac{d}{dt} \int_{\Omega} \omega \, d\Omega = 0
\]
- Here, \(\omega\) is the absolute vorticity and satisfies the continuity equation:
\[
\omega_t+\nabla\cdot(\omega u) = 0
\]
3. **Energy conservation and stability**:
- The method should minimize or avoid energy dissipation as much as possible to ensure nonlinear stability.
- The expression for energy \(E\) is:
\[
E=\frac{1}{2}D u\cdot u+\frac{1}{2}g D^{2}
\]
- The energy conservation formula is:
\[
\frac{d}{dt} \int_{\Omega} E \, d\Omega = 0
\]
4. **Linear geostrophic balance**:
- Under a constant Coriolis parameter \(f\), the linearized rotating shallow - water equations should be able to maintain a steady - state geostrophic pattern.
- The steady - state solution is obtained by choosing an arbitrary stream function \(\psi\) and setting \(D = -\frac{f}{g}\psi\) and \(u=\nabla\times\psi k\).
5. **Entropy stability**:
- By using entropy - stable numerical fluxes, ensure the entropy conservation within the volume term, while the boundary fluxes dissipate entropy, thereby achieving local entropy stability.
### Summary
The main objective of this paper is to develop an efficient and stable DG - SEM method for solving the rotating shallow - water equations on the sphere. This method not only theoretically proves the discrete conservation properties of mass, vorticity, energy, and linear geostrophic balance, but also demonstrates its robustness and stability in experiments, especially when dealing with geostrophic turbulence without introducing artificial dissipation.