Numerical model of nonhydrostatic ocean dynamics based on methods of artificial compressibility and multicomponent splitting

An algorithm is proposed for solving three-dimensional ocean hydrodynamics equations without hydrostatic approximation and traditional simplification of Coriolis acceleration. It is based on multicomponent splitting of the modified model with artificial compressibility. The original system of equations is split into two subsystems describing the transport of three velocity components and adjustment of the density and velocity fields. At the adjustment stage, the horizontal velocity components are represented as a sum of the depth means and deviations; the two corresponding subsystems are derived. For barotropic dynamics, the compressibility effect is represented as the boundary condition at the free surface, while for the baroclinic subsystem, it is introduced as ε-regularization of the continuity equation. Then, the baroclinic equations are split into two subsystems describing the hydrostatic and nonhydrostatic dynamics. The nonhydrostatic dynamics is computed at a separate splitting stage. The algorithm is included into the Institute of Numerical Mathematics of the Russian Academy of Sciences model based on “primitive” equations and verified by solving the hydrodynamics problem for the Sea of Marmara.