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 2.4.1 Structure aérodynamique de l’écoulement côté injection Dynamique générale Du côté chaud de la plaque, de l’air de refroidissement est injecté au travers de petites perforations régulièrement disposées en quinconce. Les données disponibles sur la structure fine de l’écoulement proviennent essentiellement d’expériences isothermes. Pour décrire la structure de cet écoulement, nous nous baserons tout d’abord sur les résultats de thèse de Miron (2005). Miron a cherché à se situer dans une configuration la plus réaliste possible en termes de géométrie et d’écoulement pour effectuer des mesures de vitesse longitudinale et verticale par Anémométrie Laser Doppler (LDA), du côté injection de la paroi. Pour effectuer des mesures précises de vitesse, le dispositif expérimental a été construit à l’échelle 10 par rapport à la géométrie typique dans les moteurs : la paroi fait 10 mm d’épaisseur au lieu de 1, les perforations 5 mm de diamètre au lieu de 0.5. La plaque est constituée de 12 rangées de 11 ou 12 trous en quinconce, inclinés à α = 30 ◦ avec l’écoulement, sans composante transverse (β = 0 ◦ ). L’espacement entre deux rangées consécutives est de X/d = 5.84. La distance transverse entre deux trous d’une même rangée est de Z/d = 6.74. La plaque sépare deux canaux de 24 d de haut. Deux écoulement principaux sont établis dans ces canaux. Une grille bouche en partie le côté aspiration pour faire monter la pression et injecter du fluide côté injection à travers la paroi. Les mesures de vitesse effectuées sur ce dispositif expérimental vont permettre de décrire l’écoulement côté injection. La description sera faite pour un point de fonctionnement particulier dont les caractéristiques sont les suivantes : U i est la vitesse moyenne au centre de chaque canal i = 1, 2 et le nombre de Reynolds Re i est basé sur U i et la demi-hauteur du canal. Côté injection (canal 1), on a U 1 = 4.5 m.s −1 et Re 1 = 17700. Côté aspiration (canal 2), on a U 2 = 2.26 m.s −1 et Re 2 = 8900. La différence de pression entre les deux canaux est ∆P = P 2 −P 1 = 42 Pa. Aucune mesure de débit n’est disponible. Aussi, les caractéristiques du jet ont été basées sur le maximum de module vitesse mesuré juste au-dessus de la perforation. Cette vitesse est de 8 m.s −1 . Basés sur cette vitesse, le taux de soufflage et le nombre de Reynolds du jet sont respectivement de 1.78 et 2600. La description détaillée des résultats expérimentaux de Miron (2005) répond à deux objectifs : a) le caractère isotherme et grande échelle mis à part, ce dispositif expérimental est celui qui correspond le mieux aux caractéristiques de la multi-perforation dans les turbines à gaz réelles. Ces résultats permettent donc de se représenter ce qu’est l’aérodynamique de l’écoulement près d’une paroi multi-perforée, du côté injection de la plaque. D’autre part, ces mesures ont servi de référence à une série de simulations effectuées dans cette thèse. Leur description permet donc de se familiariser avec des données qui seront réutilisées dans ce document à plusieurs reprises. Deux séries de profils seront présentées pour cette expérience. La première montre l’évolution des profils de vitesse immédiatement en aval des jets en fonction de la rangée considérée. Les profils sont mesurés côté injection, 2.92 d en aval du jet central des rangées 5, 7, 9 et 11, depuis la paroi (y = 0) jusqu’au centre du canal 1 (y = 12 d). Les positions de mesure sont reportées sur la figure 2.9. Les symboles utilisés sont ceux des profils tracés en figure 2.10. On est ici dans la zone dite « établie ». Les profils de vitesse longitudinale moyenne U (figure 2.10a) montrent une forme caractéristique des écoulements de multi-perforation. A environ y = d, un pic très net marque la présence du jet issu de la perforation située en amont. La position du jet ne bouge pas suivant les rangées. De même, de la 48
2.4 Caractérisation aérodynamique de l’écoulement ROWS 1 3 5 7 9 11 FIG. 2.9 - Localisation des mesures de la figure 2.10. rangée 5 à la rangée 11, la valeur de U ne varie pratiquement pas dans le coeur du jet. Pour le point de fonctionnement considéré, la valeur de ce pic se situe autour de U = 6.3 m.s −1 . Au-dessus d’une hauteur d’environ y = 8 d, tous les profils se rejoignent. Cela permet d’estimer la hauteur de film, distance de la plaque à partir de laquelle l’écoulement n’est pas influencé par l’injection. Cette épaisseur de film augmente légèrement avec le nombre de rangées. Au-dessus du premier pic, à environ y = 3 d, on observe la présence d’un autre pic de vitesse. Ce pic secondaire montre une évolution différente du pic de vitesse marquant le jet en amont. La vitesse de ce pic augmente nettement avec le nombre de rangées. Miron (2005) le qualifie de coeur du film. Il est principalement la signature du jet situé encore en amont (jet de la rangée 7 pour le profil mesuré rangée 9, par exemple). Le jet est plus ou moins mélangé an fonction des vitesses environnantes : la valeur de vitesse de ce pic augmente donc avec le nombre de rangées. C’est ce que nous appellerons l’effet d’accumulation. Les jets viennent nourrir petit à petit le film de refroidissement en y apportant de la masse et la quantité de mouvement associée. Il est à noter qu’aucune étude ne mentionne l’existence d’une limite pour la vitesse de ce pic. La dernière caractéristique notable de ces profils de U se situe près de la paroi (y = 0.5 d) : la présence d’un petit pic de vitesse peut être remarquée. Ce troisième pic est la manifestation du contournement du jet par l’écoulement principal. Comme l’expliquent Yavuzkurt et al. (1980a), un jet présentant un taux de soufflage τ 1 suffisamment fort (typiquement plus de 0.5) décolle à sa sortie de la perforation. La zone située sous ce jet, immédiatement en aval du trou, est une zone de dépression qui attire l’écoulement principal incident. Celui-ci contourne le jet en subissant un fort mouvement en direction de la paroi. L’écoulement contournant le jet ayant des vitesses de plus en plus élevées avec le nombre de rangées (voir le pic à y = 3 d), la valeur de vitesse du pic proche paroi a tendance à augmenter. Cet effet est appelé entrainment process, par Yavuzkurt et al. (1980a). Il constitue une caractéristique fondamentale des écoulements de refroidissement par multi-perforation, et permet d’expliquer plusieurs phénomènes 49
- 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 and 36: 2.1 Présentation de la configurati
- Page 37 and 38: 2.2 Perte de charge au passage d’
- 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: 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
- Page 87 and 88: 4.3 Sensibilité des simulations au
- Page 89 and 90: 4.3 Sensibilité des simulations au
- Page 91 and 92: 4.3 Sensibilité des simulations au
- Page 93 and 94: 4.3 Sensibilité des simulations au
- Page 95 and 96: 4.3 Sensibilité des simulations au
- Page 97 and 98: Chapitre 5 Simulations numériques
2.4 Caractérisation aérodynamique <strong>de</strong> l’écoulement<br />
ROWS<br />
1 3 5 7 9 11<br />
FIG. 2.9 - Localisation <strong>de</strong>s mesures <strong>de</strong> la figure 2.10.<br />
rangée 5 à la rangée 11, la valeur <strong>de</strong> U ne varie pratiquement pas dans le coeur du j<strong>et</strong>. Pour le point<br />
<strong>de</strong> fonctionnement considéré, la valeur <strong>de</strong> ce pic se situe <strong>autour</strong> <strong>de</strong> U = 6.3 m.s −1 . Au-<strong>de</strong>ssus d’une<br />
hauteur d’environ y = 8 d, tous les profils se rejoignent. Cela perm<strong>et</strong> d’estimer la hauteur <strong>de</strong> film,<br />
distance <strong>de</strong> la plaque à partir <strong>de</strong> laquelle l’écoulement n’est pas influencé par l’injection. C<strong>et</strong>te épaisseur<br />
<strong>de</strong> film augmente légèrement avec le nombre <strong>de</strong> rangées. Au-<strong>de</strong>ssus du premier pic, à environ y = 3 d,<br />
on observe la présence d’un autre pic <strong>de</strong> vitesse. Ce pic secondaire montre une évolution différente du<br />
pic <strong>de</strong> vitesse marquant le j<strong>et</strong> en amont. La vitesse <strong>de</strong> ce pic augmente n<strong>et</strong>tement avec le nombre <strong>de</strong><br />
rangées. Miron (2005) le qualifie <strong>de</strong> coeur du film. Il est principalement la signature du j<strong>et</strong> situé encore<br />
en amont (j<strong>et</strong> <strong>de</strong> la rangée 7 pour le profil mesuré rangée 9, par exemple). Le j<strong>et</strong> est plus ou moins<br />
mélangé an fonction <strong>de</strong>s vitesses environnantes : la valeur <strong>de</strong> vitesse <strong>de</strong> ce pic augmente donc avec le<br />
nombre <strong>de</strong> rangées. C’est ce que nous appellerons l’eff<strong>et</strong> d’accumulation. Les j<strong>et</strong>s viennent nourrir p<strong>et</strong>it<br />
à p<strong>et</strong>it le film <strong>de</strong> refroidissement en y apportant <strong>de</strong> la masse <strong>et</strong> la quantité <strong>de</strong> mouvement associée. Il est<br />
à noter qu’aucune étu<strong>de</strong> ne mentionne l’existence d’une limite pour la vitesse <strong>de</strong> ce pic. La <strong>de</strong>rnière<br />
caractéristique notable <strong>de</strong> ces profils <strong>de</strong> U se situe près <strong>de</strong> la paroi (y = 0.5 d) : la présence d’un<br />
p<strong>et</strong>it pic <strong>de</strong> vitesse peut être remarquée. Ce troisième pic est la manifestation du contournement du<br />
j<strong>et</strong> par l’écoulement principal. Comme l’expliquent Yavuzkurt <strong>et</strong> al. (1980a), un j<strong>et</strong> présentant un taux<br />
<strong>de</strong> soufflage τ 1 suffisamment fort (typiquement plus <strong>de</strong> 0.5) décolle à sa sortie <strong>de</strong> la perforation. La zone<br />
située sous ce j<strong>et</strong>, immédiatement en aval du trou, est une zone <strong>de</strong> dépression qui attire l’écoulement<br />
principal inci<strong>de</strong>nt. Celui-ci contourne le j<strong>et</strong> en subissant un fort mouvement en direction <strong>de</strong> la paroi.<br />
L’écoulement contournant le j<strong>et</strong> ayant <strong>de</strong>s vitesses <strong>de</strong> plus en plus élevées avec le nombre <strong>de</strong> rangées<br />
(voir le pic à y = 3 d), la valeur <strong>de</strong> vitesse du pic proche paroi a tendance à augmenter. C<strong>et</strong> eff<strong>et</strong> est<br />
appelé entrainment process, par Yavuzkurt <strong>et</strong> al. (1980a). Il constitue une caractéristique fondamentale<br />
<strong>de</strong>s écoulements <strong>de</strong> refroidissement par multi-perforation, <strong>et</strong> perm<strong>et</strong> d’expliquer plusieurs phénomènes<br />
49