Stability constraints for oceanic numerical models: implications for the formulation of time and space discretizations

Florian Lemarié 1 Laurent Debreu 1 Gurvan Madec 2 Jérémie Demange 3 Jean-Marc Molines 4 Marc Honnorat 3
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
2 NEMO R&D - Nucleus for European Modeling of the Ocean
LOCEAN - Laboratoire d'Océanographie et du Climat : Expérimentations et Approches Numériques
3 MOISE - Modelling, Observations, Identification for Environmental Sciences
Inria Grenoble - Rhône-Alpes, LJK - Laboratoire Jean Kuntzmann, INPG - Institut National Polytechnique de Grenoble
Abstract : Except for vertical diffusion (and possibly the external mode and bottom drag), oceanic models usually rely on explicit time-stepping algorithms subject to Courant-Friedrichs-Lewy (CFL) stability criteria. Implicit methods could be unconditionally stable, but an algebraic system must be solved at each time step and other considerations such as accuracy and efficiency are less straightforward to achieve. Depending on the target application, the process limiting the maximum allowed time-step is generally different. In this paper, we introduce offline diagnostics to predict stability limits associated with internal gravity waves, advection, diffusion, and rotation. This suite of diagnostics is applied to a set of global, regional and coastal numerical simulations with several horizontal/vertical resolutions and different numerical models. We show that, for resolutions finer that 1/2◦, models with an eulerian vertical coordinate are generally constrained by vertical advection in a few hot spots and that numerics must be extremely robust to changes in Courant number. Based on those results, we review the stability and accuracy of existing numerical kernels in vogue in primitive equations oceanic models with a focus on advective processes and the dynamics of internal waves. We emphasize the additional value of studying the numerical kernel of oceanic models in the light of coupled space-time approaches instead of studying the time schemes independently from spatial discretizations. From this study, we suggest some guidelines for the development of temporal schemes in future generation multi-purpose oceanic models.
Liste complète des métadonnées

Littérature citée [69 références]  Voir  Masquer  Télécharger
Contributeur : Florian Lemarié <>
Soumis le : vendredi 26 juin 2015 - 17:26:18
Dernière modification le : mardi 29 mai 2018 - 12:51:07
Document(s) archivé(s) le : vendredi 9 octobre 2015 - 19:00:36


Fichiers produits par l'(les) auteur(s)




Florian Lemarié, Laurent Debreu, Gurvan Madec, Jérémie Demange, Jean-Marc Molines, et al.. Stability constraints for oceanic numerical models: implications for the formulation of time and space discretizations. Ocean Modelling, Elsevier, 2015, 92, pp.124-148. 〈10.1016/j.ocemod.2015.06.006〉. 〈hal-01065979v2〉



Consultations de la notice


Téléchargements de fichiers