### 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 language | English (US) |
---|---|

Pages (from-to) | 23-40 |

Number of pages | 18 |

Journal | Journal of Computational Physics |

Volume | 388 |

DOIs | |

State | Published - Jul 1 2019 |

### Fingerprint

### Keywords

- Galerkin
- Legendre
- Quasigeostrophic
- Spectral

### ASJC Scopus subject areas

- Physics and Astronomy (miscellaneous)
- Computer Science Applications

### Cite this

*Journal of Computational Physics*,

*388*, 23-40. https://doi.org/10.1016/j.jcp.2019.03.029

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

Research output: Contribution to journal › Article

*Journal of Computational Physics*, vol. 388, pp. 23-40. https://doi.org/10.1016/j.jcp.2019.03.029

}

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 -