Hydrostatic Consistency in Sigma Coordinate Ocean Models
Chu, Peter C.
MetadataShow full item record
Truncation error and hydrostatic inconsistency at steep topography are two concerns in sigma coordinate ocean models due to the horizontal pressure gradient being difference of two large terms. A consensus is reached in the ocean modeling community on the first concern (truncation error), but not on the second concern (hydrostatic inconsistency). Since the integration of the pressure gradient over a finite volume equals the integration of the pressure over the surface of that volume (always dynamically consistent), dynamical analysis on finite volumes is used to determine the hydrostatic consistency of a sigma coordinate ocean model. A discrete, hydrostatically consistent scheme is obtained for the sigma coordinate ocean models. Comparison between finite-volume and finite-difference approaches leads to the conclusion that a Boussinesq, hydrostatic, sigma coordinate ocean model with second-order staggered scheme is always hydrostatically consistent. Guidance for improving numerical accuracy is also provided.
RightsThis publication is a work of the U.S. Government as defined in Title 17, United States Code, Section 101. Copyright protection is not available for this work in the United States.
Showing items related by title, author, creator and subject.
Verification of a non-hydrostatic dynamical core using the horizontal spectral element method and vertical finite difference method: 2-D aspects Giraldo, F.X.; Choi, S.-J.; Kim, J.; Shin, S. (Copernicus Publications, 2014);The non-hydrostatic (NH) compressible Euler equations for dry atmosphere were solved in a simplified two-dimensional (2-D) slice framework employing a spectral element method (SEM) for the horizontal discretization and ...
Efficient Construction of Unified Continuous and Discontinuous Galerkin Formulations for the 3D Euler Equations Abdi, Daniel S.; Giraldo, Francis X. (2016-05-28);A unified approach for the numerical solution of the 3D hyperbolic Euler equations using high order methods, namely continuous Galerkin (CG) and discontinuous Galerkin (DG) methods, is presented. First, we examine how ...
Verification of a non-hydrostatic dynamical core using horizontally spectral element vertically finite difference method: 2D Aspects Choi, Suk-Jin; Giraldo, Francis X.; Kim, Junghan; Shin, Seoleun (2014-04);The non-hydrostatic (NH) compressible Euler equations of dry atmosphere are solved in a simplified two dimensional (2D) slice (X-Z) framework employing a spectral element method (SEM) for the horizontal discretization and ...