Résolution des équations de Saint Venant avec transport de polluant par un schéma éléments finis et un schéma volumes finis

Abstract or Introduction

Dans ce travail nous avons modélisé numériquement les équations de Saint Venant qui régissent les écoulements en eau peu profonde, couplées à l'équation de diffusion d'un polluant. Nous avons utilisé deux techniques de calcul numérique parmi les plus efficaces, la méthode des éléments finis et la méthode des volumes finis.

En éléments finis nous utilisons le schéma de Taylor-Galerkin qui est un schéma centré de deuxième ordre dans l'espace et dans le temps. La solution est calculée par ce schéma en deux étapes, la première utilise une approximation constante par élément et la seconde une approximation linéaire. Nous obtenons un schéma conservatif, précis et simple à programmer ; néanmoins il génère des oscillations aux voisinages des gradients importants. Cet handicap est résolu par l'introduction d'un modèle de diffusion numérique du à MacCormak. Cette dissipation a pour caractéristique de n'agir que dans les régions à forts gradients et d'être inactive dans les régions lisses.

Le deuxième schéma est issu de la méthode des volumes finis et il est dû à Roe. On y utilise un maillage de type éléments finis non-structuré (cell-centered) ou les grandeurs physiques sont évaluées aux centres des cercles circonscrits aux triangles. La partie hyperbolique du système d'équations est approximée par un solveur de Riemann, puis le flux numérique est corrigé de façon à garantir l'entropie de la solution. Afin d'augmenter la précision du schéma dans l'espace à l'ordre deux, la technique MUSCL est utilisée. Quant à la discrétisation du temps, elle est effectuée par le schéma de Runge Kutta à l'ordre deux.

La validition des deux méthodes est démontrée sur le test du ressaut hydraulique dont la solution théorique est connue. Ensuite, des comparaisons sont effectuées sur une série d'exemples d'écoulements à surface libre. Les deux schémas se révèlent être puissants, précis et efficaces. Ils présentent des performances comparables quant à leur capacité de capter les ondes de choc. Néanmoins, le schéma de Taylor-Galerkin nécessite le réglage par l'utilisateur du coefficient de la diffusion numérique tandis que le shéma de Roe est plus gourmand en temps CPU.