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 - Funding Information:
Acknowledgements. The National Science Foundation partially supported this research through cooperative agreement OCE-0424602. The National Oceanic and Atmospheric Administration (NA11NOS0120036 and AB-133F-12-SE-2046), Bonneville Power Administration (00062251), and Corps of Engineers (W9127N-12-2-007 and G13PX01212) provided partial motivation and additional support. This work was supported by the UK’s Engineering and Physical Science Research Council (grant numbers EP/M011054/1, EP/L000407/1) and the Natural Environment Research Council (grant number NE/K008951/1). This work used the Extreme Science and Engineering Discovery Environment (XSEDE), which is supported by National Science Foundation grant number ACI-1053575. The authors acknowledge the Texas Advanced Computing Center (TACC) at the University of Texas at Austin for providing HPC resources that have contributed to the research results reported within this paper.
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
VL - 11
SP - 4359
EP - 4382
JO - Geoscientific Model Development
JF - Geoscientific Model Development
SN - 1991-959X
IS - 11
ER -