|
|
Roullet, Guillaume; Molemaker, Maarten Jeroen; Ducousso, Nicolas; Dubos, Thomas. |
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... |
Tipo: Text |
Palavras-chave: Non-hydrostatic; Sigma coordinates; Poisson equation. |
Ano: 2017 |
URL: http://archimer.ifremer.fr/doc/00399/51045/51818.pdf |
| |