Contexte
La modélisation eulérienne d’écoulements compressibles multiphasiques, comprenant des phases potentiellement solide et/ou fluide, est d’une grande importance pour la simulation numérique d’applications de l’ingénieur pour des milieux fluides (cavitation au voisinage de pâle, lithotripsie) comme pour des milieux solides (compaction de milieux hétérogènes [1], poreux). Elle nécessite entre autres de pouvoir suivre les interfaces au cours du temps et traiter des mélanges au sein d’un maillage de calcul fixe. Parmi les différentes approches existantes, la méthode des interfaces diffuses [2] permet de considérer les interfaces et les mélanges avec les mêmes modèle et schéma numérique dans tout le domaine de calcul, les interfaces étant considérées comme des zones de mélange.
Cependant, les interfaces sont autorisées à diffuser numériquement, et s’étalent donc sur plusieurs mailles au cours de la simulation. La simulation numérique de chocs sur ces interfaces étalées peut faire apparaître un phénomène de piégeage d’onde résultant d’une vitesse du son du mélange plus faible de celle de chaque phase, conduisant à des problèmes de convergence du calcul. Une approche suivie a consisté, après l’étape de convection des ondes, à effectuer une relaxation à vitesse finie des pressions de chaque phase dans une maille de calcul [2], ce qui permet de limiter le phénomène de piégeage d’onde et améliore la convergence des calculs.
On s’intéresse dans ce stage à une modélisation mathématique de la relaxation qui est décrite dans un cadre thermodynamiquement cohérent, avec des potentiels d’état et de dissipation. Cette modélisation peut se réécrire comme un principe variationnel, c’est-à-dire un problème d’optimisation d’une fonctionnelle, conduisant à un intégrateur variationnel discret précis à l’ordre un. Ce solveur a été implémenté dans un code de démonstration écrit en Python, et a montré son bon fonctionnement et sa robustesse. En revanche, la simulation numérique de scénarii d’écoulements multiphasiques multi-dimensionnels plus complexes requiert l’utilisation d’un langage et d’un code de calcul plus efficaces.Objectifs
L'objectif de ce stage est d’implémenter ce solveur variationnel de relaxation des pressions dans le code de calcul CFD open-source ECOGEN [3], développé en langage C++. ECOGEN (https://code-mphi.github.io/ECOGEN/docs/sphinx_docs/index.html). Ce code permet de résoudre des modèles d’écoulements multiphasiques avec la méthode des interfaces diffuses, utilise un schéma volumes finis avec divers types de grilles de calcul (cartésiennes, AMR et non structurées), et peut fonctionner sur plusieurs milliers de coeurs en utilisant MPI. Des cas de test simples permettront de valider l’implémentation, avant de tester des configurations multi-dimensionnelles plus riches.
Références
[1] N. Favrie, K. Schmidmayer, and J. Massoni. A multiphase irreversible-compaction model for granular-porous materials. Continuum Mechanics and Thermodynamics, pages 1–25, 2022.
[2] K. Schmidmayer, J. Cazé, F. Petitpas, E. Daniel, and N. Favrie. Modelling interactions between
waves and diffused interfaces. International Journal for Numerical Methods in Fluids, 95(2):215–241, 2023.
[3] K. Schmidmayer, F. Petitpas, S. Le Martelot, and É. Daniel. ECOGEN: An open-source tool for
multiphase, compressible, multiphysics flows. Computer Physics Communications, 251:107093, 2020.
