these simulation numerique et modelisation de l'ecoulement autour ...
these simulation numerique et modelisation de l'ecoulement autour ... these simulation numerique et modelisation de l'ecoulement autour ...
L’ÉCOULEMENT AUTOUR D’UNE PAROI MULTI-PERFORÉE U 1 ρ 1 V j ρ j ρ2 U 2 FIG. 2.6 - Définition des vitesses et masses volumiques caractéristiques de l’écoulement. rapport de flux de quantité de mouvement noté τ 2 : τ 1 = ρ jV j ρ 1 U 1 (2.2) τ 2 = ρ jV j 2 ρ 1 U 1 2 (2.3) où les vitesses U 1 , U 2 et V j sont des vitesses débitantes, ρ la masse volumique moyenne. Les quantités indicées j se rapportent au jet (voir la figure 2.6). Ces rapports sont parfois remplacés par le rapport de vitesses V j et le rapport de densité ρ j . Notons que suivant les auteurs, les vitesses de référence U 1 ρ 1 sont parfois des vitesses maximales. On peut ajouter à ces paramètres de soufflage les caractéristiques de chaque écoulement (nombre de Reynolds, nombre de Mach, épaisseur de quantité de mouvement, taux de turbulence ou encore gradient de pression longitudinal) ou encore des paramètres d’aspiration (rapports entre quantités caractéristiques de jet et quantités caractéristiques de l’écoulement secondaire). 2.1.4 Un modèle pour la multi-perforation Avant d’entamer la revue des travaux existants sur la multi-perforation, il semble indispensable de préciser quels types d’informations sont utilisables pour la modélisation que l’on souhaite effectuer. On veut remplacer, dans les calculs de chambre de combustion (et de chambre de combustion avec contournement), la paroi multi-perforée par un modèle. Un travail de modélisation doit donc être réalisé pour chacun des deux côtés de la paroi perforée pour permettre d’effectuer des simulations incluant le contournement. L’objectif de la modélisation est double. Il s’agit premièrement de remplacer une paroi composée de parties solides et des zones perforées par une paroi homogène, dont l’effet serait la combinaison des deux. La plaque multi-perforée doit donc être homogénéisée. Cette homogénéisation permet de s’affranchir des contraintes de tailles des perforations pour le maillage proche paroi. L’autre aspect du travail est la modélisation pure de l’écoulement. Pour mieux comprendre ce besoin de modélisation, on peut s’appuyer sur l’exemple de la paroi imperméable classique. Dans le cas d’une 36
2.2 Perte de charge au passage d’une paroi multi-perforée simple paroi solide, la paroi est homogène, ce qui ôte une difficulté par rapport à la paroi perforée. La condition physique de non-glissement au mur se traduit par une condition de vitesse nulle. Sur un maillage grossier, incapable de représenter finement l’écoulement proche paroi, la condition de nonglissement impose un gradient normal de vitesse à la paroi faux, comme schématisé sur la figure 2.7. Par exemple, l’utilisation d’un maillage trop grossier dans la direction normale à la paroi dans un écoulement en conduite mènerait à sous-estimer complètement la perte de charge. Pour résoudre ce problème, on utilise un modèle qui permet de mimer le véritable effet de la paroi sur l’écoulement : la condition de non-glissement peut alors être abandonnée et remplacée par une condition de frottement imposé à la paroi. L’estimation de ce frottement est possible grâce à la connaissance de la structure de l’écoulement près de la paroi (von Karman, 1934; Schlichting & Gersten, 2000). C’est la loi de paroi qui permet de déterminer le frottement pariétal à partir de la connaissance de l’écoulement au premier point au-dessus de la paroi. PREMIER POINT DU MAILLAGE GRADIENT FAUX GRADIENT AVEC MODELE PAROI FIG. 2.7 - Schéma d’une simulation mal résolue à la paroi. Estimation du gradient avec ou sans modèle. Le raisonnement est le même pour la multi-perforation. Une idée simple pour homogénéiser la multiperforation consiste à utiliser une condition en vitesse qui serait la moyenne pondérée des parties solides et perforées. σ est la porosité de la paroi, définie par l’équation 2.1. La condition de vitesse homogénéisée impose alors ⃗ V = ⃗ V j .σ, où ⃗ V j est le vecteur vitesse caractéristique de l’écoulement dans les perforations. Comme nous le verrons en chapitre 6, cette condition n’est pas représentative de l’écoulement proche paroi. Notamment, elle ne reproduit pas du tout les flux pariétaux, qui vont caractériser l’effet de la paroi perforée sur son environnement. Il convient donc de savoir le plus précisément possible ce que sont ces flux à la paroi perforée, et d’être capable ensuite de les modéliser. Les informations les plus pertinentes pour la modélisation sont donc : la connaissance des flux à la paroi perforée et les moyens de les estimer. Cela vaut pour les deux côtés de la paroi. 2.2 Perte de charge au passage d’une paroi multi-perforée En passant au travers d’une paroi perforée (d’une fente ou de perforations cylindriques par exemple), l’air de refroidissement subit une importante perte de charge. Cette perte de charge va déterminer quelle 37
- Page 1: UNIVERSITE MONTPELLIER II SCIENCES
- Page 5 and 6: Table des matières Remerciements 7
- Page 7 and 8: Remerciements Cette thèse a été
- Page 9 and 10: Liste des symboles Lettres romaines
- Page 11 and 12: Introduction générale La volonté
- Page 13 and 14: Chapitre 1 Contexte industriel et s
- Page 15 and 16: 1.1 Les turbines à gaz pour les tu
- Page 17 and 18: 1.1 Les turbines à gaz a) b) GAZ C
- Page 19 and 20: 1.2 La simulation numérique des é
- Page 21 and 22: 1.2 La simulation numérique des é
- Page 23 and 24: 1.2 La simulation numérique des é
- Page 25 and 26: 1.2 La simulation numérique des é
- Page 27 and 28: 1.3 Objectifs et plan de la thèse
- Page 29 and 30: Chapitre 2 L’écoulement autour d
- Page 31 and 32: 2.1 Présentation de la configurati
- Page 33 and 34: 2.1 Présentation de la configurati
- Page 35: 2.1 Présentation de la configurati
- Page 39 and 40: 2.3 La multi-perforation : aspects
- Page 41 and 42: 2.3 La multi-perforation : aspects
- Page 43 and 44: 2.3 La multi-perforation : aspects
- Page 45 and 46: 2.3 La multi-perforation : aspects
- Page 47 and 48: 2.4 Caractérisation aérodynamique
- Page 49 and 50: 2.4 Caractérisation aérodynamique
- Page 51 and 52: 2.4 Caractérisation aérodynamique
- Page 53 and 54: 2.4 Caractérisation aérodynamique
- Page 55 and 56: 2.4 Caractérisation aérodynamique
- Page 57 and 58: 2.4 Caractérisation aérodynamique
- Page 59 and 60: 2.4 Caractérisation aérodynamique
- Page 61 and 62: 2.4 Caractérisation aérodynamique
- Page 63 and 64: 2.4 Caractérisation aérodynamique
- Page 65 and 66: 2.4 Caractérisation aérodynamique
- Page 67 and 68: 2.5 Modélisation de la multi-perfo
- Page 69 and 70: Chapitre 3 Simulations des grandes
- Page 71 and 72: 3.2 Simulations des Grandes Echelle
- Page 73 and 74: 3.3 Le code de calcul AVBP 3.3.1 As
- Page 75 and 76: 3.3 Le code de calcul AVBP dans le
- Page 77 and 78: 3.3 Le code de calcul AVBP de sa st
- Page 79 and 80: Chapitre 4 Simulations numériques
- Page 81 and 82: 4.1 Quelles simulations pour la mod
- Page 83 and 84: 4.2 Choix de la méthode de simulat
- Page 85 and 86: 4.2 Choix de la méthode de simulat
2.2 Perte <strong>de</strong> charge au passage d’une paroi multi-perforée<br />
simple paroi soli<strong>de</strong>, la paroi est homogène, ce qui ôte une difficulté par rapport à la paroi perforée.<br />
La condition physique <strong>de</strong> non-glissement au mur se traduit par une condition <strong>de</strong> vitesse nulle. Sur un<br />
maillage grossier, incapable <strong>de</strong> représenter finement l’écoulement proche paroi, la condition <strong>de</strong> nonglissement<br />
impose un gradient normal <strong>de</strong> vitesse à la paroi faux, comme schématisé sur la figure 2.7.<br />
Par exemple, l’utilisation d’un maillage trop grossier dans la direction normale à la paroi dans un<br />
écoulement en conduite mènerait à sous-estimer complètement la perte <strong>de</strong> charge. Pour résoudre ce<br />
problème, on utilise un modèle qui perm<strong>et</strong> <strong>de</strong> mimer le véritable eff<strong>et</strong> <strong>de</strong> la paroi sur l’écoulement : la<br />
condition <strong>de</strong> non-glissement peut alors être abandonnée <strong>et</strong> remplacée par une condition <strong>de</strong> frottement<br />
imposé à la paroi. L’estimation <strong>de</strong> ce frottement est possible grâce à la connaissance <strong>de</strong> la structure <strong>de</strong><br />
l’écoulement près <strong>de</strong> la paroi (von Karman, 1934; Schlichting & Gersten, 2000). C’est la loi <strong>de</strong> paroi qui<br />
perm<strong>et</strong> <strong>de</strong> déterminer le frottement pariétal à partir <strong>de</strong> la connaissance <strong>de</strong> l’écoulement au premier point<br />
au-<strong>de</strong>ssus <strong>de</strong> la paroi.<br />
PREMIER POINT<br />
DU MAILLAGE<br />
GRADIENT<br />
FAUX<br />
GRADIENT<br />
AVEC MODELE<br />
PAROI<br />
FIG. 2.7 - Schéma d’une <strong>simulation</strong> mal résolue à la paroi. Estimation du gradient avec ou sans modèle.<br />
Le raisonnement est le même pour la multi-perforation. Une idée simple pour homogénéiser la multiperforation<br />
consiste à utiliser une condition en vitesse qui serait la moyenne pondérée <strong>de</strong>s parties soli<strong>de</strong>s<br />
<strong>et</strong> perforées. σ est la porosité <strong>de</strong> la paroi, définie par l’équation 2.1. La condition <strong>de</strong> vitesse homogénéisée<br />
impose alors ⃗ V = ⃗ V j .σ, où ⃗ V j est le vecteur vitesse caractéristique <strong>de</strong> l’écoulement dans les perforations.<br />
Comme nous le verrons en chapitre 6, c<strong>et</strong>te condition n’est pas représentative <strong>de</strong> l’écoulement proche<br />
paroi. Notamment, elle ne reproduit pas du tout les flux pariétaux, qui vont caractériser l’eff<strong>et</strong> <strong>de</strong> la paroi<br />
perforée sur son environnement. Il convient donc <strong>de</strong> savoir le plus précisément possible ce que sont ces<br />
flux à la paroi perforée, <strong>et</strong> d’être capable ensuite <strong>de</strong> les modéliser. Les informations les plus pertinentes<br />
pour la modélisation sont donc : la connaissance <strong>de</strong>s flux à la paroi perforée <strong>et</strong> les moyens <strong>de</strong> les estimer.<br />
Cela vaut pour les <strong>de</strong>ux côtés <strong>de</strong> la paroi.<br />
2.2 Perte <strong>de</strong> charge au passage d’une paroi multi-perforée<br />
En passant au travers d’une paroi perforée (d’une fente ou <strong>de</strong> perforations cylindriques par exemple),<br />
l’air <strong>de</strong> refroidissement subit une importante perte <strong>de</strong> charge. C<strong>et</strong>te perte <strong>de</strong> charge va déterminer quelle<br />
37