Résumé:

La caractérisation d'écoulements d'air intérieur est un enjeu d'importance dans un contexte social et réglementaire visant à optimiser la consommation énergétique et le confort thermique dans le bâtiment. Parallèlement, la qualité de l'air intérieur est responsable de milliers de décès dans le monde et suscite l'intérêt croissant de la communauté scientifique et les institutions; la crise épidémique liée au virus SARS-COVID 19 souligne l'utilité des études de milieux confinés.Dans ce but, la simulation est un outil efficace pouvant reproduire de manière fiable l'écoulement considéré et de manière moins coûteuse comparée à l'expérimental. Cette thèse se focalise sur le développement d'un outil numérique pour les simulations aérauliques à l'échelle locale (Computational Fluid Dynamics, CFD). Ce dernier a pour but d'être appliquée dans des études de milieux intérieurs allant du résidentiel (aide à la maîtrise d'oeuvre, étude de qualité d'air intérieur) à l'industriel (sûreté des centrales nucléaires, ventilation d'enceintes sportives).Après avoir posé le contexte et les enjeux liés à l'aéraulique dans le premier chapitre, l'identification des phénomènes physiques caractérisant l'écoulement d'air intérieur et le choix des équations à utiliser sont présentés dans un second chapitre.Au troisième chapitre, pour répondre aux enjeux de modélisation, un schéma volumes finis d'ordre 2 en temps pour les écoulements à densité variable est proposé, pour des solutions régulières et discontinues. Développé pour l'air sec tout d'abord, le schéma est implémenté dans le logiciel CFD sous licence libre code_saturne. Ce dernier rentre dans la famille des theta-schéma de type prédiction correction. Le second ordre en temps est atteint grâce à une localisation de variables décalées. De plus, l'énergie totale est conservée grâce à la résolution de l'équation de l'énergie interne complétée par un terme source dérivé de l'équation discrète de l'énergie cinétique. Enfin, la variation de pression est prise en compte en linéarisant l'équation d'état, ce qui conduit à une équation d'Helmholtz à résoudre pour la pression. Les termes liés à cette dernière sont implicités, menant à des calculs plus rapides tout en évitant toute contrainte de type CFL liée aux ondes acoustiques. Après une analyse numérique menant à des nouvelles contraintes de stabilité, le schéma en temps est vérifié et validé par de nombreux cas allant de zéro à trois dimensions et du régime incompressible au compressible, représentatifs des enjeux de modélisation aérauliques. Des simulations turbulentes d'ordre un et deux (RANS, LES) sont également réalisées.Le quatrième chapitre présente des développements complémentaires dans permettant d'étendre le schéma à l'air humide avec changement de phase. L'équilibre thermodynamique est considéré et la fraction massique d'eau, sous forme gazeuse et liquide, est transportée. Dans le but d'utiliser les équations choisies auparavant, le changement de phase est traité en utilisant la méthode de Newton à partir de l'énergie interne résolue. De la même manière, une analyse numérique et des vérifications sont réalisées. Des nouvelles conditions CFL sont également étudiées. Enfin, dans le cinquième chapitre, l'outil numérique est appliqué pour pour caractériser l'écoulement d'air intérieur au sein du Stade Pierre de Coubertin dans le cadre des jeux olympiques 2024 de Paris. Le maillage numérique est créé à partir d'un nuage de points issu de mesures de scanners 3D. Des premières simulations mènent à l'identification de zones d'intérêt dynamiques et thermiques et à la création d'un protocole expérimental pour une campagne de mesures. S'en suit une validation du schéma dans le but de reproduire l'évolution de la concentration de particules lors de la finale de la ligue de Handball française, lors de laquelle un pic de PM10 a été mesuré.

Lien Teams