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 ...

21.04.2014 Views

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

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!