Développement d'un code de calcul multiphasique multiconstituants
Auteur / Autrice : | Raphaël Di Chiara Roupert |
Direction : | Gerhard Schäfer, Michel Quintard |
Type : | Thèse de doctorat |
Discipline(s) : | Mécanique des fluides |
Date : | Soutenance en 2009 |
Etablissement(s) : | Strasbourg |
Mots clés
Mots clés contrôlés
Résumé
Une méthode pour la simulation des écoulements triphasiques compressibles est proposé en tenant compte de la gravité et de la capillarité. Les équations sont ré-écrites en formulation fractionnaire avec une équation en pression globale et deux équations de saturation. La pression globale satisfaisant une condition dite de «différentiel totale» (DT) a été introduite par [Chavent et JafTré, 1986] afin de simplifier la formulation mathématique de l'écoulement. Une méthode d'optimisation avec contraintes est utilisée pour déterminer sur un diagramme ternaire les variables second aires à partir des saturations effectives. Une autre approche d'interpolation DT récemment mise au point par [Chavent, 2008] est également mis en oeuvre. Par rapport à la première formulation, la nouvelle approche est équivalente et plus efficace que la formulation classique. Nous discutons de la construction des variables à l'aide d'éléments finis composites spectraux, en particulier les conditions aux limites qui doivent tenir compte des données diphasiques connues les limites du diagramme ternaire [di Chiara Roupert et al. 2010]. Deux méthodes numériques performantes sont ensuite utilisées pour résoudre l'équation de pression globale et les deux équations de saturation en eau et en huile. Les éléments finis discontinus sont utilisés pour représenter le terme convectif des deux équations de saturation alors que la méthode des éléments finis mixtes est choisie de résoudre l'équation de pression globale et la partie diffusive des équations de saturation. Des résultats numériques sont donnés pour certains cas de tests analytiques et comparés aux autres approches.