Improving vertical advection in oceanic modeling

Florian Lemarié 1 Laurent Debreu 1
1 AIRSEA - Mathematics and computing applied to oceanic and atmospheric flows
Inria Grenoble - Rhône-Alpes, LJK - Laboratoire Jean Kuntzmann, UJF - Université Joseph Fourier - Grenoble 1, INPG - Institut National Polytechnique de Grenoble
Abstract : Recent papers by Shchepetkin (2015) and Lemarié et al. (2015) have emphasized that the time-step of an oceanic model with an Eulerian vertical coordinate and an explicit time-stepping scheme is very often restricted by vertical advection in a few hot spots (i.e. most of the grid points are integrated with small Courant numbers, compared to the Courant-Friedrichs-Lewy (CFL) condition, except just few spots where numerical instability of the explicit scheme occurs first). The consequence is that the numerics for vertical advection must have good stability properties while being robust to changes in Courant number in terms of accuracy. An other constraint for oceanic models is the strict control of numerical mixing imposed by the highly adiabatic nature of the oceanic interior (i.e. mixing must be very small in the vertical direction below the boundary layer). We examine in this talk the possibility of mitigating vertical Courant-Friedrichs-Lewy (CFL) restriction, while avoiding numerical inaccuracies associated with standard implicit advection schemes (i.e. large sensitivity of the solution on Courant number, large phase delay, and possibly excess of numerical damping with unphysical orientation), see Baldauf & Skamarock (2009) or Shchepetkin (2015). Most regional oceanic models have been successfully using fourth order compact schemes (e.g. Lele, 1992) for vertical advection. In this talk we present a general framework to derive generic expressions for (one-step) coupled time and space high order compact schemes (see Daru & Tenaud (2004) for a thorough description of coupled time and space schemes). Among other properties, we show that those schemes are unconditionally stable and have very good accuracy properties even for large Courant numbers while having a very reasonable computational cost. We illustrate the properties of those schemes and compare them to existing schemes using linear and nonlinear numerical experiments.
Type de document :
Communication dans un congrès
Partial Differential Equations on the Sphere Workshop, Oct 2015, Seoul, South Korea. 〈〉
Liste complète des métadonnées
Contributeur : Florian Lemarié <>
Soumis le : mardi 5 janvier 2016 - 10:18:09
Dernière modification le : jeudi 11 janvier 2018 - 06:27:25


  • HAL Id : hal-01250645, version 1



Florian Lemarié, Laurent Debreu. Improving vertical advection in oceanic modeling. Partial Differential Equations on the Sphere Workshop, Oct 2015, Seoul, South Korea. 〈〉. 〈hal-01250645〉



Consultations de la notice