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 ...
SIMULATIONS DES GRANDES ÉCHELLES. LE CODE AVBP Modélisation du flux de chaleur de sous-maille De manière analogue aux tensions de Reynolds de sous-maille, le flux de chaleur de sous-maille est représenté de la même façon que le flux de chaleur laminaire : avec q i t = −λ t ∂ ˜T ∂x i , (3.35) λ t = ρ ν tC p P r t . (3.36) P r t est le nombre de Prandtl de sous-maille. Moin et al. (1991) ont réalisé des simulations numériques directes de THI, de cisaillement homogène et de canal turbulent pour étendre le modèle de Smagorinsky dynamique de Germano et al. (1991) à des cas compressibles. Ils calculent le nombre de Prandtl de sousmaille dans ces trois configurations et montrent qu’il dépend fortement de l’écoulement considéré : les valeurs sont globalement entre 0.4 et 0.6 pour les écoulements de THI et de cisaillement homogène. Dans le canal turbulent, le Prandtl de sous-maille varie en fonction de la distance à la paroi. Au centre du canal, il est d’environ 0.5 ; au mur, il atteint des valeurs de l’ordre de 1 pour un nombre de Prandtl laminaire de 0.7. Ainsi, une procédure dynamique serait recommandée pour évaluer P r t . En pratique, P r t est souvent considéré comme étant constant avec des valeurs comprises entre 0.6 et 1.1 selon les applications. Dans AVBP, de même que les constantes C S et C w , il est constant, et λ t ne dépend donc que du choix du modèle pour ν t . Sa valeur est fixée à 0.9. Si cette valeur est discutable pour les écoulements libres, nous nous intéressons à l’écoulement proche d’une paroi. En conséquence, nous avons gardé cette valeur de P r t = 0.9. En outre, le raffinement utilisé dans nos calculs limite l’impact des modèles de sous-maille. 3.3.6 Filtrage Pour une description des différents types de filtrage et une discussion sur leur importance, le lecteur pourra s’appuyer sur les ouvrages traitant spécifiquement de SGE (par exemple Sagaut, 2002). Dans AVBP, le filtrage est implicite : c’est le maillage qui l’impose. L’échelle de coupure dans la représentation du spectre d’énergie turbulente est donc fixée par la taille du maillage. Dans les expressions des modèles de sous-maille, la longueur caractéristique du filtrage △ est définie à partir du volume de la maille V maille : △ = V 1/3 maille . 78
Chapitre 4 Simulations numériques de l’écoulement dans les parois multi-perforées Ce chapitre présente les simulations réalisées pour étudier l’écoulement autour de parois multiperforées. Les méthodes expérimentales sont confrontées à des difficultés importantes pour fournir les données nécessaires à la modélisation de l’écoulement dans des configurations réalistes : tailles caractéristiques, conditions opératoires extrêmes, difficulté des mesures proche paroi... On souhaite donc effectuer des simulations numériques —qui ne présentent pas ce genre de difficultés— pour obtenir ces données. Les simulations effectuées en RANS ont toutes cherché à caractériser le développement spatial d’un film de refroidissement le long de plusieurs rangées de perforations. Les simulations RANS sont cependant confrontées au problème de la modélisation de la turbulence. On sait que le modèle choisi influe fortement sur les résultats (voir chapitre 2). Ces simulations pourraient être effectuées en SGE ou en SND. Mais les coûts de calcul s’avèrent alors extrêmement élevés. Avant de procéder à ce genre de simulations, il convient donc de se demander si elles sont vraiment indispensables. C’est la question à laquelle nous allons répondre dans la première partie de ce chapitre (4.1). Nous verrons que dans cette thèse, le choix est fait de simuler une configuration particulière, dans laquelle la plaque perforée est considérée comme étant infinie. Dans ce cas, le domaine de calcul peut être réduit à une boîte fermée par des conditions de périodicité. Le cas périodique demande la mise en place d’une méthode de simulation particulière : la partie 4.2 détaille et compare deux méthodes permettant d’effectuer des simulations périodiques autour d’une paroi multi-perforée. L’une d’entre elles sera retenue. Dans le reste du chapitre (partie 4.3) nous étudions la sensibilité du calcul à différents paramètres (schémas numériques, maillages, taille du domaine de calcul, etc.). Ces tests sont effectués sur un cas isotherme, correspondant à la configuration de Miron (2005). L’extension de la méthode aux calculs de refroidissement à proprement parler (écoulement secondaire froid injecté dans écoulement principal chaud) sera présentée au chapitre 7.
- 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 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: 3.3 Le code de calcul AVBP de sa st
- 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
- Page 99 and 100: 2 S. Mendez and F. Nicoud COMBUSTIO
- Page 101 and 102: 4 S. Mendez and F. Nicoud (c) The i
- Page 103 and 104: 6 S. Mendez and F. Nicoud CALCULATI
- Page 105 and 106: 8 S. Mendez and F. Nicoud Name Numb
- Page 107 and 108: 10 S. Mendez and F. Nicoud x/d z/d
- Page 109 and 110: 12 S. Mendez and F. Nicoud main, it
- Page 111 and 112: 14 S. Mendez and F. Nicoud 1.0 a 1.
- Page 113 and 114: 16 S. Mendez and F. Nicoud ROWS 1 3
- Page 115 and 116: 18 S. Mendez and F. Nicoud Figure 1
- Page 117 and 118: 20 S. Mendez and F. Nicoud 30 ◦ 3
- Page 119 and 120: 22 S. Mendez and F. Nicoud & Mahesh
- Page 121 and 122: 24 S. Mendez and F. Nicoud Figure 1
- Page 123 and 124: 26 S. Mendez and F. Nicoud Figure 1
- Page 125 and 126: 28 S. Mendez and F. Nicoud film, th
- Page 127 and 128: 30 S. Mendez and F. Nicoud Region t
Chapitre 4<br />
Simulations numériques <strong>de</strong> l’écoulement<br />
dans les parois multi-perforées<br />
Ce chapitre présente les <strong>simulation</strong>s réalisées pour étudier l’écoulement <strong>autour</strong> <strong>de</strong> parois multiperforées.<br />
Les métho<strong>de</strong>s expérimentales sont confrontées à <strong>de</strong>s difficultés importantes pour fournir les<br />
données nécessaires à la modélisation <strong>de</strong> l’écoulement dans <strong>de</strong>s configurations réalistes : tailles caractéristiques,<br />
conditions opératoires extrêmes, difficulté <strong>de</strong>s mesures proche paroi... On souhaite donc<br />
effectuer <strong>de</strong>s <strong>simulation</strong>s numériques —qui ne présentent pas ce genre <strong>de</strong> difficultés— pour obtenir ces<br />
données.<br />
Les <strong>simulation</strong>s effectuées en RANS ont toutes cherché à caractériser le développement spatial d’un<br />
film <strong>de</strong> refroidissement le long <strong>de</strong> plusieurs rangées <strong>de</strong> perforations. Les <strong>simulation</strong>s RANS sont cependant<br />
confrontées au problème <strong>de</strong> la modélisation <strong>de</strong> la turbulence. On sait que le modèle choisi influe<br />
fortement sur les résultats (voir chapitre 2). Ces <strong>simulation</strong>s pourraient être effectuées en SGE ou en<br />
SND. Mais les coûts <strong>de</strong> calcul s’avèrent alors extrêmement élevés.<br />
Avant <strong>de</strong> procé<strong>de</strong>r à ce genre <strong>de</strong> <strong>simulation</strong>s, il convient donc <strong>de</strong> se <strong>de</strong>man<strong>de</strong>r si elles sont vraiment<br />
indispensables. C’est la question à laquelle nous allons répondre dans la première partie <strong>de</strong> ce chapitre<br />
(4.1). Nous verrons que dans c<strong>et</strong>te thèse, le choix est fait <strong>de</strong> simuler une configuration particulière, dans<br />
laquelle la plaque perforée est considérée comme étant infinie. Dans ce cas, le domaine <strong>de</strong> calcul peut<br />
être réduit à une boîte fermée par <strong>de</strong>s conditions <strong>de</strong> périodicité. Le cas périodique <strong>de</strong>man<strong>de</strong> la mise<br />
en place d’une métho<strong>de</strong> <strong>de</strong> <strong>simulation</strong> particulière : la partie 4.2 détaille <strong>et</strong> compare <strong>de</strong>ux métho<strong>de</strong>s<br />
perm<strong>et</strong>tant d’effectuer <strong>de</strong>s <strong>simulation</strong>s périodiques <strong>autour</strong> d’une paroi multi-perforée. L’une d’entre elles<br />
sera r<strong>et</strong>enue. Dans le reste du chapitre (partie 4.3) nous étudions la sensibilité du calcul à différents<br />
paramètres (schémas numériques, maillages, taille du domaine <strong>de</strong> calcul, <strong>et</strong>c.).<br />
Ces tests sont effectués sur un cas isotherme, correspondant à la configuration <strong>de</strong> Miron (2005).<br />
L’extension <strong>de</strong> la métho<strong>de</strong> aux calculs <strong>de</strong> refroidissement à proprement parler (écoulement secondaire<br />
froid injecté dans écoulement principal chaud) sera présentée au chapitre 7.