Energy-conserving Galerkin approximations for quasigeostrophic dynamics

DOI10.1016/J.JCP.2019.03.029zbMATH Open1452.86008arXiv1810.10628OpenAlexW2898242084MaRDI QIDQ2220596FDOQ2220596

Keith Julien, Ian Grooms, Matthew Watwood, K. Shafer Smith

Publication date: 25 January 2021

Published in: Journal of Computational Physics (Search for Journal in Brave)

Abstract: A method is presented for constructing energy-conserving Galerkin approximations in the vertical coordinate of the full quasigeostrophic model with active surface buoyancy. The derivation generalizes the approach of Rocha emph{et al.} (2016) to allow for general bases. Details are then presented for a specific set of bases: Legendre polynomials for potential vorticity and a recombined Legendre basis from Shen (1994) for the streamfunction. The method is tested in the context of linear baroclinic instability calculations, where it is compared to the standard second-order finite-difference method and to a Chebyshev collocation method. The Galerkin scheme is quite accurate even for a small number of degrees of freedom $N$, and growth rates converge much more quickly with increasing $N$ for the Galerkin scheme than for the finite-difference scheme. The Galerkin scheme is at least as accurate as finite differences and can in some cases achieve the same accuracy as the finite difference scheme with ten times fewer degrees of freedom. The energy-conserving Galerkin scheme is of comparable accuracy to the Chebyshev collocation scheme in most linear stability calculations, but not in the Eady problem where the Chebyshev scheme is significantly more accurate. Finally the three methods are compared in the context of a simplified version of the nonlinear equations: the two-surface model with zero potential vorticity. The Chebyshev scheme is the most accurate, followed by the Galerkin scheme and then the finite difference scheme. All three methods conserve energy with similar accuracy, despite not having any a priori guarantee of energy conservation for the Chebyshev scheme. Further nonlinear tests with non-zero potential vorticity to assess the merits of the methods will be performed in a future work.

Full work available at URL: https://arxiv.org/abs/1810.10628

zbMATH Keywords

spectralGalerkin schemeLegendrequasigeostrophic

Mathematics Subject Classification ID

Finite element, Rayleigh-Ritz and Galerkin methods for initial value and initial-boundary value problems involving PDEs (65M60) Applications to the sciences (65Z05) Hydrology, hydrography, oceanography (86A05)

Cites Work

• Title not available (Why is that?)
• Title not available (Why is that?)
• Spectral Methods in MATLAB
• Atmospheric and Oceanic Fluid Dynamics
• A Practical Guide to Pseudospectral Methods
• Spectral Methods
• Numerical algorithms for stationary statistical properties of dissipative dynamical systems
• Applied analysis
• Long time stability of a classical efficient scheme for two-dimensional Navier-Stokes equations
• Efficient Spectral-Galerkin Method I. Direct Solvers of Second- and Fourth-Order Equations Using Legendre Polynomials
• Formation of strong fronts in the 2-D quasigeostrophic thermal active scalar
• Efficient multi-dimensional solution of PDEs using Chebyshev spectral methods
• Computational design for long-term numerical integration of the equations of fluid motion: two-dimensional incompressible flow. Part I
• Stochastic superparameterization in quasigeostrophic turbulence
• Global in time classical solutions to the 3D quasi-geostrophic system for large initial data
• The periodic quasigeostrophic equations: existence and uniqueness of strong solutions
• Validity of the Quasigeostrophic Model for Large-Scale Flow in the Atmosphere and Ocean
• The Cauchy Problem and the Continuous Limit for the Multilayer Model in Geophysical Fluid Dynamics
• Balanced and unbalanced routes to dissipation in an equilibrated Eady flow
• Ocean turbulence at meso and submesoscales: connection between surface and interior dynamics

Cited In (2)

• On energy exchanges between eddies and the mean flow in quasigeostrophic turbulence
• Energy conserving SUPG methods for compatible finite element schemes in numerical weather prediction

Uses Software

• Matlab
• Q-GCM

Recommendations

• Title not available (Why is that?) 👍 👎
• Title not available (Why is that?) 👍 👎
• Title not available (Why is that?) 👍 👎
• An energy-conserving Galerkin scheme for a class of nonlinear dispersive equations 👍 👎
• On the energy and helicity conservations for the 2-D quasi-geostrophic equation 👍 👎
• Energy conserving Galerkin finite element schemes for the primitive equations of numerical weather prediction 👍 👎
• Enslaved finite difference approximations for quasigeostrophic shallow flows 👍 👎
• Quasi-geostrophic approximation of anelastic convection 👍 👎
• Energy conservation for the generalized surface quasi-geostrophic equation 👍 👎
• Energy and potential enstrophy flux constraints in quasi-geostrophic models 👍 👎

This page was built for publication: Energy-conserving Galerkin approximations for quasigeostrophic dynamics

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2220596)