Author: Sara PavanPublisher:
ISBN:
Category :
Languages : fr
Pages : 0
Book Description
Cette thèse a pour objectif la construction de schémas d'ordre élevé et peu diffusifs pour le transport d'un scalaire dans les écoulements à surface libre, en deux ou trois dimensions. On souhaite en particulier obtenir des schémas robustes, qui gardent au niveau discret les propriétés mathématiques de l'équation de transport avec une faible diffusion numérique, et les utiliser sur des cas industriels. Dans ce travail deux méthodes numériques sont envisagées : une méthode aux volumes finis (VF) et une méthode aux résidus distribués (RD). Dans les deux cas, l'équation de transport est résolue avec une approche découplée, qui est la solution la plus avantageuse en termes de précision et de coûts de calcul. Pour ce qui concerne la méthode aux volumes finis, les équations de Saint-Venant couplées à l'équation du transport sont d'abord résolues avec un schéma dit vertex-centred où le flux numérique est approximé avec un solveur de Riemann appelé Harten-Lax-Van Leer-Contact [135]. A partir de cette approche, une formulation découplée est proposée. Cette dernière permet de résoudre l'équation du transport avec un pas de temps plus grand que celui de la formulation couplée. Cette idée a été d'abord proposée pour d'autres schémas dans [13]. Pour augmenter l'ordre de précision en espace, la technique MUSCL [89] est utilisée en combinaison avec l'approche découplée. Finalement, la problématique des zones sèches est abordée. Dans le cas de la méthode aux résidus distribués, les équations de Saint-Venant sont résolues avec une méthode éléments finis, et la méthode RD est utilisée seulement pour discrétiser l'équation du transport, en focalisant l'attention sur les problèmes non stationnaires. L'équation de continuité du fluide discrétisée est employée pour garantir la conservation de la masse et le principe du maximum. Pour obtenir des schémas d'ordre deux dans les problèmes non stationnaires, un schéma prédicteur-correcteur [112] est utilisé, en l'adaptant au cas de concentration moyennée sur la verticale. Une version d'ordre 1 mais peu diffusive, est aussi présentée dans ce travail. De plus, un schéma localement implicite, complètement nouveau, est aussi formulé pour pouvoir traiter le problème des bancs découvrant. Les deux techniques sont validées d'abord sur des cas simples, pour évaluer l'ordre de précision des schémas et ensuite sur des cas plus complexes pour vérifier aussi les autres propriétés numériques. Les résultats montrent que les nouveaux schémas sont à la fois précis et conservatifs, tout en gardant la monotonie comme le prévoient les démonstrations. Un cas d'application industriel est aussi présenté en conclusion. Le schéma prédicteur-correcteur RD est adapté aussi au cas 3D, sans aucun problème théorique nouveau, par rapport au cas 2D. Les propriétés de base des schémas sont validées sur des cas test préliminaires.
