Skip to Main content Skip to Navigation
Journal articles

Compact symmetric Poisson equation discretization for non-hydrostatic sigma coordinates ocean model

Abstract : In anticipation of relaxing the hydrostatic assumption in a sigma coordinates primitive equations ocean model, we show how a projection method can be designed with the use of a compact symmetric 15-point stencil for the Poisson equation. This is achieved by recognizing that, owing the non-orthogonality of the grid, the velocity has a contravariant and covariant set of components. The two sets play a different role in the primitive equations: the contravariant components enter the definition of the model fluxes, whereas the covariant components experience the forces and in particular the pressure gradient. By treating these two sets separately, the discretized gradient and divergence operators are simple finite differences. The two sets of components are related via a linear transformation, the metric tensor, which is entirely determined by the kinetic energy. We show how the spatial discretization of the kinetic energy fully controls the Poisson equation discretization, including its boundary conditions. The discretization of the Poisson equation is shown to converge at second order and to behave as well or better than alternative methods. This approach is a prerequisite to implement later an efficient Poisson solver, such as a multigrid algorithm.
Complete list of metadata

https://hal-insu.archives-ouvertes.fr/insu-03682732
Contributor : Nathalie POTHIER Connect in order to contact the contributor
Submitted on : Tuesday, May 31, 2022 - 11:41:30 AM
Last modification on : Thursday, June 2, 2022 - 3:38:11 AM

Links full text

Identifiers

Citation

Roullet Guillaume, Molemaker M. Jeroen, Ducousso Nicolas, Dubos Thomas. Compact symmetric Poisson equation discretization for non-hydrostatic sigma coordinates ocean model. OCEAN MODELLING, 2017, 118, pp.107-117. ⟨10.1016/j.ocemod.2017.09.001⟩. ⟨insu-03682732⟩

Share

Metrics

Record views

4