Ce travail de thèse est consacré au développement d'une approche numérique permettant de conduire des simulations "low Mach number" d'écoulements réactifs. L'algorithme d'intégration retenu pour procéder à la résolution des équations de transport repose sur une méthode implicite de prédiction-correction (méthode de projection). Une contrainte physique est retenue pour garantir que le champ de vitesse est résolue correctement. Le code de calcul est soumis à plusieurs séries de vérifications préliminaires basées sur l'emploi de la méthode solution manufacturées pour des conditions incompressibles d'abord puis à masse volumique variable qui permettent de statuer quant à la bonne implémentation des schémas numériques retenus. Les performances de l'outil numérique en terme de stabilité et de robustesse sont elles-aussi analysées dans des situations simples: couche de mélange à densité variable en développement spatial et temporel. Le modèle numérique final repose sur l'emploi d'une méthode hybride LES / PDF. Pour ce qui concerne la représentation de la turbulence, deux fermetures sont implémentées pour représenter l'effet des fluctuations de vitesse non résolues. Il s'agit du modèle de Smagorinsky dans sa version dynamique ou non. La spécification de conditions aux limites turbulentes réalistes est elle-aussi analysée en détail et trois approches différentes sont considérées. Pour ce qui concerne la combustion, l'influence des fluctuations de composition aux échelles non résolues est pris en compte par le biais d'une résolution de la PDF scalaire de sous maille. Le modèle de PDF correspondant repose sur l'emploi d'une méthode de Monte Carlo. Des équations différentielles stochastiques, équivalentes aux équations de Fokker-Planck, sont résolues pour la variable de progrès de la réaction chimique. L'objectif final est aussi de pouvoir procéder, à moyen terme, à des simulations LES en géométries complexes et l'emploi du calcul distribué est essentiel. De ce point de vue, la méthode de décomposition de domaine retenue dans ce travail montre des niveau de performances relativement satisfaisants. Les capacités du modèle numérique résultant de ces développements sont illustrées sur deux configurations expérimentales. La première géométrie correspond à un écoulement très fortement turbulent de réactifs pré-mélangés dans un canal bidimensionnel. La seconde correspond à un jet rapide et non confiné de réactifs pré-mélangés.
Identifer | oai:union.ndltd.org:CCSD/oai:tel.archives-ouvertes.fr:tel-00665800 |
Date | 18 November 2011 |
Creators | Vedovoto, João, Marcelo |
Source Sets | CCSD theses-EN-ligne, France |
Language | English |
Detected Language | French |
Type | PhD thesis |
Page generated in 0.0017 seconds