LEGI - UMR 5519

Financement : LabEx Tec21

Cette thèse a pour objectif la simulation numérique du phénomène d’ébullition sur maillages non structurés. L’ébullition est le changement de phase des particules fluides de la phase liquide vers la phase vapeur sous l’action des flux thermiques à l’interface séparant les deux phases. Il s’agit donc d’un phénomène rencontré au sein d’écoulements diphasiques et piloté par le taux de transfert de masse à l’interface. Ce taux de transfert de masse est calculé à partir des flux thermiques de part et d’autre de l’interface. Cela implique donc la nécessité d’adopter une méthode de suivi d’interface très précise pour localiser l’interface à tout instant de la simulation. Les équations de Navier-Stokes sont alors couplées à l’équation de la chaleur par l’intermédiaire du taux de transfert de masse à l’interface. De telles simulations ont été menées par Tanguy et al. (J. Comput. Phys., 2014) sur des maillages cartésiens axisymétriques en deux dimensions. Dans cette thèse, nous étendons cette méthodologie à des maillages non structurés en trois dimensions (maillages composés de tétraèdres non réguliers utiles pour décrire des géométries complexes). Pour ce faire, nous avons développé un solveur spécifique dans le code YALES2 (code diphasique basé sur la méthode des volumes finis pour des maillages 3D non structurés). Le suivi de l’interface est assuré par la méthode Level Set. Le changement de phase engendre des discontinuités de vitesse et de pression à l’interface qui dépendent notamment du taux de transfert de masse. Ces discontinuités sont prises en compte par la méthode Ghost Fluid à l’aide de deux champs de vitesse et deux champs de température. Cette méthodologie étant déjà bien établie pour des maillages structurés cartésiens, l’apport de cette thèse est la possibilité de simuler le changement de phase par ébullition sur des maillages non structurés en trois dimensions. Les spécificités des maillages non structurés ont nécessité de nombreux développements pour la réinitialisation de la fonction Level Set après advection, ainsi que l’utilisation d’opérateurs d’ordres élevés pour le calcul du taux de transfert de masse à l’interface. L’ensemble des développements proposés est finalement validé sur maillages non structurés à l’aide du cas-test analytique d’une bulle 3D en expansion dans un liquide surchauffé au repos.

Publications

Currently in Preparation or Submitted

2019

Sahut, G., Ghigliotti, G., Marty, P., & Balarac, G. (2019). Evaluation of Level Set reinitialization algorithms for phase change simulation on unstructured grids.

  

Peer-reviewed Publications

2023

Berthelon, T., Sahut, G., Leparoux, J., Balarac, G., Lartigue, G., Bernard, M., et al. (2023). Toward the use of LES for industrial complex geometries. Part II: Reduce the time-to-solution by using a linearised implicit time advancement. Journal of Turbulence, 24(6-7), 311–329.

  

2021

Sahut, G., Ghigliotti, G., Balarac, G., Bernard, M., Moureau, V., & Marty, P. (2021). Numerical simulation of boiling on unstructured grids. Journal of Computational Physics, , 110161.

  

Conference Proceedings

2021

Atmani, Y., Pecquery, F., Cailler, M., Moureau, V., Lartigue, G., Mercier, R., et al. (2021). Consistent scalar transport with front capturing methods: application to two-phase heat transfer. In ICLASS 2021, 15th Triennial International Conference on Liquid Atomization and Spray Systems. Edinburgh, United Kingdom.

  

2019

Sahut, G., Ghigliotti, G., Bégou, P., Marty, P., & Balarac, G. (2019). Numerical simulation of boiling on 3D unstructured grids. In ICMFHT’19 – 4th International Conference on Multiphase Flow and Heat Transfer. Rome, Italy.

  

2018

Sahut, G., Ghigliotti, G., Marty, P., & Balarac, G. (2018). Numerical Simulation of Boiling. In DTPF – Dispersed Two-Phase Flows. Toulouse, France.

  

Ph.D. Theses

2019

Sahut, G. (2019). Numerical simulation of boiling on unstructured grids. Ph.D. thesis, Université Grenoble Alpes, .