Nouvelles méthodes numériques pour les écoulements en eaux peu profondes
Creators: name: Beljadid, Abdelaziz
Publisher: name: Université d'Ottawa / University of Ottawa
Year: 2015
Subjects: Méthodes numériques
écoulements peu profonds
systèmes hyperboliques
méthode volumes finis
schémas semi-implicites
analyse de stabilité

description: Dans ce projet de recherche, on s’intéresse au développement et à l’évaluation de nouvelles méthodes numériques pour les écoulements peu profonds. De nouvelles techniques de discrétisation spatiales et temporelles des équations sont proposées. Une partie de la thèse est dédiée au développement d’une méthode des volumes finis explicite d’ordre élevé et d’une famille de schémas semi-implicites qui sont efficaces pour la modélisation des processus lents et rapides dans les écoulements océaniques et atmosphériques. La deuxième partie du projet de recherche concerne la construction d’un schéma numérique efficace sans solveur de Riemann pour les écoulements peu profonds avec une topographie variable sur un maillage non structuré. Dans cette partie de la thèse, une nouvelle approche est proposée pour l'analyse de stabilité des schémas numériques non structurés pour les équations en eaux peu profondes. Dans la troisième partie de la thèse, deux schémas de volumes finis sont développés pour les lois de conservation sur des surfaces courbes qui ont un large potentiel d’être appliqués aux écoulements peu profonds sur la sphère. Dans ces cas, les schémas numériques sont développés en adoptant la démarche suivie par Stanley Osher. Cette démarche consiste à utiliser des systèmes hyperboliques simples qui génèrent des phénomènes d'ondes complexes et des solutions qui ont différentes structures. Ces solutions sont très efficaces pour tester les méthodes numériques. Dans notre cas, nous avons utilisé les équations de Burgers qui ont joué un rôle très important dans le développement des schémas numériques à capture de chocs en mécanique des fluides. Dans le premier article, une nouvelle méthode des volumes finis décentrée explicite est proposée pour le système de Saint-Venant avec un terme source qui comprend le paramètre de Coriolis en utilisant un maillage non structuré. La plupart des schémas numériques décentrés, efficaces pour les ondes rapides (ondes de gravité), conduisent à un niveau d'amortissement élevé pour les ondes lentes (ondes de Rossby). La méthode proposée donne de bons résultats à la fois pour les ondes de gravité et les ondes de Rossby. Les techniques proposées sont suffisantes pour supprimer le bruit numérique des ondes courtes sans amortissement des ondes longues, telles que les ondes de Rossby qui sont essentielles dans le transport de l’énergie dans les océans et l'atmosphère. Dans le cas où le système comprend une large gamme de fréquences des ondes, ce qui est le cas des écoulements atmosphériques, il est important d’utiliser des méthodes semi-implicites afin d’opter pour un pas de temps optimal. La méthode semi-implicite semi-lagrangienne à deux niveaux (SETTLS) proposée par Hortal (2002) a une région de stabilité absolue indépendante du nombre de Courant-Friedrichs-Lewy (CFL). La plupart des modèles de prévision numérique atmosphérique utilisent cette méthode comme schéma temporel. Cependant, la méthode SETTLS peut générer des oscillations pour le traitement du terme non linéaire surtout pour le cas des solutions qui ont un caractère oscillatoire. Pour remédier à ce problème, dans le deuxième article, nous avons proposé une nouvelle classe de schémas semi-implicites semi-lagrangiens potentiellement applicables aux modèles atmosphériques. Cette classe de schémas numériques présente plusieurs avantages de stabilité, de précision et de convergence. De bons résultats sont obtenus en comparaison à d'autres schémas semi-implicites semi-lagrangiens et méthodes semi-implicites de type prédicteur-correcteur. Dans le troisième article, un nouveau schéma équilibre partiellement centré est développé pour la résolution numérique des équations de Saint-Venant avec une topographie variable sur un maillage non structuré. Cette méthode est stable et simple puisqu'elle ne fait pas appel à la résolution du problème de Riemann. La méthode proposée est précise pour le cas des solutions discontinues et peut être appliquée aux écoulements peu profonds avec une topographie variable et une géométrie complexe où l'utilisation des maillages non structurés est avantageuse. Motivé par de nombreuses applications en dynamique des fluides, dans le projet de thèse on s’intéresse également au développement de méthodes numériques dans le cas des surfaces courbes. L'objectif est de concevoir des méthodes numériques robustes et efficaces pour le cas des solutions discontinues et qui préservent la structure fondamentale des équations, notamment les propriétés liées à la géométrie. Pour développer ces méthodes, l'approche suivie par Stanley Osher est adoptée et les équations de Burgers sont utilisées vu leur importance pour le développement des schémas numériques à capture de chocs. Dans le quatrième article, une méthode des volumes finis satisfaisant la compatibilité géométrique est développée pour les lois de conservation sur la sphère. Cette méthode est basée sur la résolution du problème de Riemann généralisé et l'approche du «splitting» directionnel en latitude et en longitude sur la sphère. Les dimensions géométriques sont considérées de manière analytique et la forme discrète du schéma numérique proposé respecte la propriété de compatibilité géométrique. La méthode proposée est stable et précise pour le cas des solutions discontinues de grands chocs et amplitudes en comparaison avec des schémas numériques très connus. Une nouvelle classification des flux est proposée en introduisant les notions de flux feuilletés et de flux génériques. Le comportement asymptotique des solutions est étudié en fonction de la nature du flux et les propriétés des solutions discontinues sont analysées. Les résultats démontrent la capacité et le potentiel de la méthode proposée pour la résolution des lois de conservation sur la sphère dans le cas des solutions discontinues. Ce schéma numérique pourrait être étendu au cas des équations de Saint-Venant sur la sphère. Dans le cinquième article, on propose un schéma numérique efficace respectant la propriété de compatibilité géométrique pour les lois de conservation sur la sphère. La méthode proposée présente plusieurs avantages, notamment de bons résultats dans le cas des solutions discontinues avec des chocs d’amplitudes moyennes, une faible dissipation numérique et une simplicité puisqu'elle ne fait pas appel à la résolution du problème de Riemann. Cette méthode pourrait être étendue au cas des équations de Saint-Venant sur la sphère. Dans le sixième article, une nouvelle approche est proposée pour analyser la stabilité des schémas numériques appliqués aux écoulements peu profonds. Cette méthode utilise la notion du pseudo spectre des matrices. La méthode proposée est efficace en comparaison avec les méthodes couramment utilisées telles que la stabilité asymptotique et la stabilité de Lax-Richtmyer. Cette approche est utile pour le choix du type de maillage, des emplacements appropriés des variables primitives (hauteur et vitesses), et de la méthode de discrétisation la plus stable.