An energetically balanced, quasi-Newton integrator for non-hydrostatic vertical atmospheric dynamics
arXiv:2004.14090 · doi:10.1016/j.jcp.2020.109988
Abstract
An energetically balanced, implicit integrator for non-hydrostatic vertical atmospheric dynamics on the sphere is presented. The integrator allows for the exact balance of energy exchanges in space and time for vertical atmospheric motions by preserving the skew-symmetry of the non-canonical Hamiltonian formulation of the compressible Euler equations. The performance of the integrator is accelerated by a preconditioning strategy that reduces the dimensionality of the inner linear system. Here we reduce the four component velocity, density, density weighted potential temperature and Exner pressure system into a single equation for the density weighted potential temperature via repeated Schur complement decomposition and the careful selection of coupling terms. As currently formulated, the integrator is based on a horizontal-vertical spatial splitting that does not permit bottom topography. The integrator is validated for standard test cases for baroclinic instability and a non-hydrostatic gravity wave on the sphere and a rising bubble in a high-resolution planar geometry, and shows robust convergence across all of these regimes.
References in corpus (11)
- Preserving energy resp. dissipation in numerical PDEs using the "Average Vector Field" method
- An energy consistent discretization of the nonhydrostatic equations in primitive variables
- Energy-conserving time propagation for a geometric particle-in-cell Vlasov--Maxwell solver
- Mixed Mimetic Spectral Element Method for Stokes Flow: A Pointwise Divergence-Free Solution
- Compatible finite element spaces for geophysical fluid dynamics
- Energy-enstrophy conserving compatible finite element schemes for the rotating shallow water equations with slip boundary conditions
- Discrete conservation properties for shallow water flows using mixed mimetic spectral elements
- Energy conserving upwinded compatible finite element schemes for the rotating shallow water equations
- A Mixed Mimetic Spectral Element Model of the Rotating Shallow Water Equations on the Cubed Sphere
- A Mixed Mimetic Spectral Element Model of the 3D Compressible Euler Equations on the Cubed Sphere
- Energy conserving SUPG methods for compatible finite element schemes in numerical weather prediction