Energy-conserving Galerkin approximations for quasigeostrophic dynamics

Matthew Watwood, Ian Grooms, Keith A. Julien, Kendall Shafer-Smith

Research output: Contribution to journalArticle

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 et al. (2016) [1] 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) [2] 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.

Original languageEnglish (US)
Pages (from-to)23-40
Number of pages18
JournalJournal of Computational Physics
Volume388
DOIs
StatePublished - Jul 1 2019

Fingerprint

Vorticity
vorticity
collocation
degrees of freedom
approximation
baroclinic instability
energy methods
Legendre functions
energy conservation
Buoyancy
buoyancy
Nonlinear equations
Finite difference method
nonlinear equations
energy
Energy conservation
derivation
Polynomials

Keywords

  • Galerkin
  • Legendre
  • Quasigeostrophic
  • Spectral

ASJC Scopus subject areas

  • Physics and Astronomy (miscellaneous)
  • Computer Science Applications

Cite this

Energy-conserving Galerkin approximations for quasigeostrophic dynamics. / Watwood, Matthew; Grooms, Ian; Julien, Keith A.; Shafer-Smith, Kendall.

In: Journal of Computational Physics, Vol. 388, 01.07.2019, p. 23-40.

Research output: Contribution to journalArticle

Watwood, Matthew ; Grooms, Ian ; Julien, Keith A. ; Shafer-Smith, Kendall. / Energy-conserving Galerkin approximations for quasigeostrophic dynamics. In: Journal of Computational Physics. 2019 ; Vol. 388. pp. 23-40.
@article{5b0af6f0f27a40359170334168217987,
title = "Energy-conserving Galerkin approximations for quasigeostrophic dynamics",
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 et al. (2016) [1] 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) [2] 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.",
keywords = "Galerkin, Legendre, Quasigeostrophic, Spectral",
author = "Matthew Watwood and Ian Grooms and Julien, {Keith A.} and Kendall Shafer-Smith",
year = "2019",
month = "7",
day = "1",
doi = "10.1016/j.jcp.2019.03.029",
language = "English (US)",
volume = "388",
pages = "23--40",
journal = "Journal of Computational Physics",
issn = "0021-9991",
publisher = "Academic Press Inc.",

}

TY - JOUR

T1 - Energy-conserving Galerkin approximations for quasigeostrophic dynamics

AU - Watwood, Matthew

AU - Grooms, Ian

AU - Julien, Keith A.

AU - Shafer-Smith, Kendall

PY - 2019/7/1

Y1 - 2019/7/1

N2 - 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 et al. (2016) [1] 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) [2] 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.

AB - 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 et al. (2016) [1] 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) [2] 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.

KW - Galerkin

KW - Legendre

KW - Quasigeostrophic

KW - Spectral

UR - http://www.scopus.com/inward/record.url?scp=85063405843&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85063405843&partnerID=8YFLogxK

U2 - 10.1016/j.jcp.2019.03.029

DO - 10.1016/j.jcp.2019.03.029

M3 - Article

AN - SCOPUS:85063405843

VL - 388

SP - 23

EP - 40

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

ER -