A Numerical Technique for Coupling the Momentum and the Continuity Equations for Semi-Implicit 3D Ocean Models

TL;DR

This paper introduces an efficient semi-implicit finite difference method for 3D ocean modeling that couples momentum and continuity equations with linear complexity, improving computational efficiency in simulating coastal and oceanic processes.

Contribution

It presents a novel O(N) coupling technique for momentum and continuity equations in semi-implicit 3D ocean models, enhancing computational efficiency.

Findings

Method achieves linear complexity O(N).

Couples momentum and continuity equations efficiently.

Applicable to 3D hydrostatic primitive equations.

Abstract

Semi-implicit methods are powerful and efficient tools for the three-dimensional modeling of coastal and oceanic processes. A semi-implicit finite difference method for 3D hydrostatic primitive equations is presented in this paper. The governing equations are time-discretized with an implicit treatment of barotropic pressure gradient and vertical viscosity. The discretized momentum equations along a water column are coupled with the depth-integrated continuity equation of the column to construct a linear system for free-surface elevations. The novelty of this work lies in formulating an efficient method with the order of complexity O(N) for coupling the momentum and the continuity equations. In this method, the horizontal velocity components are expressed in terms of neighbor free-surface elevations by some simple recursive formulas and then are substituted in the integrated continuity equation.