Premium
Simulation of shallow‐water jets with a unified element‐based continuous/discontinuous Galerkin model with grid flexibility on the sphere
Author(s) -
Marras S.,
Kopera M. A.,
Giraldo F. X.
Publication year - 2014
Publication title -
quarterly journal of the royal meteorological society
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.744
H-Index - 143
eISSN - 1477-870X
pISSN - 0035-9009
DOI - 10.1002/qj.2474
Subject(s) - grid , inviscid flow , shallow water equations , galerkin method , discontinuous galerkin method , mathematics , finite element method , geometry , mathematical analysis , computer science , physics , mechanics , thermodynamics
We test the behaviour of a unified continuous/discontinuous Galerkin (CG/DG) shallow‐water model in spherical geometry with curved elements on three different grids of ubiquitous use in atmospheric modelling: (i) the cubed‐sphere, (ii) the reduced latitude–longitude, and (iii) the icosahedral grid. Both conforming and non‐conforming grids are adopted including static and dynamically adaptive grids for a total of twelve mesh configurations. The behaviour of CG and DG on the different grids are compared for a nonlinear midlatitude perturbed jet and for a linear case that admits an analytic solution. Because the inviscid solution on certain grids shows a high sensitivity to the resolution, the viscous counterpart of the governing equations is also solved and the results compared. The logically unstructured element‐based CG/DG model described in this article is flexible with respect to arbitrary grids. However, we were unable to define a best grid configuration that could possibly minimize the error regardless of the characteristic geometry of the flow. This is especially true if the governing equations are not regularized by the addition of a sufficiently large, fully artificial, diffusion mechanism, as will be shown. The main novelty of this study lies in the unified implementation of two element‐based Galerkin methods that share the same numerical machinery and do not rely on any specific grid configuration to solve the shallow‐water equation on the sphere.