Schémas numériques d'advection et de propagation d'ondes de gravité dans les modèles de circulation océanique
Auteur / Autrice : | Jérémie Demange |
Direction : | Éric Blayo, Laurent Debreu, Patrick Marchesiello |
Type : | Thèse de doctorat |
Discipline(s) : | Mathématiques appliquées |
Date : | Soutenance le 21/10/2014 |
Etablissement(s) : | Grenoble |
Ecole(s) doctorale(s) : | École doctorale Mathématiques, sciences et technologies de l'information, informatique (Grenoble ; 1995-....) |
Partenaire(s) de recherche : | Laboratoire : Laboratoire Jean Kuntzmann (Grenoble) |
Jury : | Président / Présidente : Emmanuel Maitre |
Examinateurs / Examinatrices : Virginie Daru, Emmanuel Audusse | |
Rapporteurs / Rapporteuses : Daniel Le Roux, Gurvan Madec |
Résumé
Les modèles numériques d'océans régionaux tridimensionnels sont basés sur la résolution des équations primitives et utilisent pour la plupart des méthodes de résolution eulérienne de type différences finies sur des grilles décalées. Ces modèles doivent représenter fidèlement les transports et transferts d'énergie. L'amélioration de ces modèles numériques exige donc (i) l'identification des processus prépondérants, notamment en terme de dissipation, dans ces transferts et (ii) la construction de méthodes numériques respectant un certain nombre d'équilibres. La première partie du travail se concentre sur la propagation des ondes externes et internes de gravité. Nous nous intéresserons en premier lieu à la stabilité de la séparation en mode rapide (barotrope) et lents (baroclines) et montrons qu'elle peut être ameliorée en levant certaines hypothèses traditionnellement effectuées. Dans un second temps, nous étudions l'impact de la discrétisation (ordre des schémas, grilles décalées ou non) sur la propagation des ondes internes de gravité provenant du couplage vitesse pression. Une décomposition en modes verticaux nous permet également de proposer un schéma espace temps très efficace. La seconde partie étudie en détail les schémas d'advection de quantité de mouvement et de traceurs, tout particulièrement dans l'objectif d'une réduction de la diffusion diapycnale (diffusion dans les directions orthogonales aux couches de densité constante). Ce travail nous amène tout d'abord à porter notre attention sur les schémas d'advection verticaux souvent négligés au regard de la dimension horizontale. Les bonnes propriétés d'un schéma compact (et de ses variantes espace temps et monotones) sont mises en avant. Enfin nous analysons le comportement multidimensionnel de ces schémas d'advection.