TY - JOUR
T1 - Thetis coastal ocean model
T2 - Discontinuous Galerkin discretization for the three-dimensional hydrostatic equations
AU - Kärnä, Tuomas
AU - Kramer, Stephan C.
AU - Mitchell, Lawrence
AU - Ham, David A.
AU - Piggott, Matthew D.
AU - Baptista, António M.
N1 - Publisher Copyright:
© Author(s) 2018.
PY - 2018/10/30
Y1 - 2018/10/30
N2 - Unstructured grid ocean models are advantageous for simulating the coastal ocean and river-estuary-plume systems. However, unstructured grid models tend to be diffusive and/or computationally expensive, which limits their applicability to real-life problems. In this paper, we describe a novel discontinuous Galerkin (DG) finite element discretization for the hydrostatic equations. The formulation is fully conservative and second-order accurate in space and time. Monotonicity of the advection scheme is ensured by using a strong stability-preserving time integration method and slope limiters. Compared to previous DG models, advantages include a more accurate mode splitting method, revised viscosity formulation, and new second-order time integration scheme. We demonstrate that the model is capable of simulating baroclinic flows in the eddying regime with a suite of test cases. Numerical dissipation is well-controlled, being comparable or lower than in existing state-of-the-art structured grid models.
AB - Unstructured grid ocean models are advantageous for simulating the coastal ocean and river-estuary-plume systems. However, unstructured grid models tend to be diffusive and/or computationally expensive, which limits their applicability to real-life problems. In this paper, we describe a novel discontinuous Galerkin (DG) finite element discretization for the hydrostatic equations. The formulation is fully conservative and second-order accurate in space and time. Monotonicity of the advection scheme is ensured by using a strong stability-preserving time integration method and slope limiters. Compared to previous DG models, advantages include a more accurate mode splitting method, revised viscosity formulation, and new second-order time integration scheme. We demonstrate that the model is capable of simulating baroclinic flows in the eddying regime with a suite of test cases. Numerical dissipation is well-controlled, being comparable or lower than in existing state-of-the-art structured grid models.
UR - http://www.scopus.com/inward/record.url?scp=85055848412&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85055848412&partnerID=8YFLogxK
U2 - 10.5194/gmd-11-4359-2018
DO - 10.5194/gmd-11-4359-2018
M3 - Article
AN - SCOPUS:85055848412
SN - 1991-959X
VL - 11
SP - 4359
EP - 4382
JO - Geoscientific Model Development
JF - Geoscientific Model Development
IS - 11
ER -