1

5GM - MECFD – TP2 StarCCM+

Ecoulement subsonique autour d'un profil d'aile

I. INTRODUCTION

Dans le cadre de ce TP, vous allez effectuer le calcul d'un écoulement subsonique autour d’un profil
d'aile. Il s'agit d'un profil d'aile supercritique de chez Boeing, ici en conditions de vol subsonique, qui
sera étudié en conditions de vol transsonique (pour lesquelles il est en fait conçu) dans un TP ultérieur.

Téléchargez le fichier "boeing_sub_2019.sim" sur Moodle. Le calcul est quasiment prêt à être lancé.
Vous allez auparavant expliquer les options qui ont été choisies, pour ensuite mener le calcul et
analyser la physique de l'écoulement. Pour cela, répondez aux questions qui vous sont posées sur la
feuille ci-dessous.

II. PARAMÈTRES DE SIMULATION

II.1 . Paramètres géométriques

Ouvrir la scène de géométrie, prenez connaissance de la géométrie utilisée puis répondez aux
questions suivantes :
- quelles sont les dimensions de l'aile (corde L= 1m , épaisseur maximale e =0.1m ) ?
- quelles sont les dimensions du domaine de calcul : longueur = 10.5 ; largeur = 8m ?
Vous semblent-elles adaptées ? Pourquoi ? Comment vérifier qu'elles le sont effectivement ?
Il faut être à un distance suffisant de l’entrée ~10 fois

- quelle est l'orientation du profil (corde) par rapport aux frontières du domaine de calcul?
horizontale
- Quelles sont les solutions pour pouvoir étudier plusieurs incidences ?

2
L'incidence peut être entrée directement dans les conditions initiales/aux limites, en projetant de façon
appropriée la vitesse sur les axes. Ici, on a préféré créer un repère relatif correspondant à la direction
du vent incident ; il suffit alors de modifier l'angle de ce repère par rapport au repère absolu (de la boite
de calcul et du profil) pour modifier l'incidence (se répercutant sur tout le calcul, notamment sur les axes
de projection pour la portance et la traînée). Pour vérifier l'orientation et localisation de ce repère par
rapport au repère global, faire "Tools>coordinate Systems> Laboratory >Local Coordinate Systems>
Cartesian 1". A quelle incidence en degréscorrespondent les cosinus directeurs actuels: = °

II.2 . Physique du problème : fluide considéré

On cherche à simuler les conditions de vol stationnaire dans de l'air sec à 5500m d'altitude avec un
nombre de Mach M=0,45.

- quel type de fluide considérer ? Pourquoi ?

Air compressible, M=0.45>0.3

Est-ce le cas dans le fichier de simulation? .

Pour le vérifier, suivre le cheminement ci-dessous, puis renseigner le tableau des caractéristiques
effectivement utilisées pour ce fluide dans le calcul (sans omettre les unités):

Chaleur spécifique à pression constante 1003.62

Masse molaire moléculaire 28.96
Constante du gaz r 287
Coefficient adiabatique  1.4
Viscosité dynamique  1,855.10e-5
(voir valeurs en annexe)

II.3 . Physique du problème : conditions aux limites

3
Sur le dessin ci-dessus, sans regarder ce qui a été effectivement choisi pour la simulation, indiquer le
type de conditions aux limites que vous adopteriez pour les différentes limites de la zone de calcul.
Repérez ensuite avec le cheminement ci-dessus le choix des différentes zones géométriques
correspondant aux différentes conditions aux limites. Dans StarCCM, on distingue le "type" de
condition, auquel correspondent différentes icônes, le choix de paramètres pour les décrire ( physics
conditions) et les valeurs de ces paramètres (physics values). Cela correspond-il à votre choix?

Dans le tableau ci-dessous, justifiez le type de conditions aux limites retenues et précisez les conditions
mathématiques sur les termes de pression p et de composantes de vitesse.
Limite Type de condition Justification et condition mathématique correspondante
(nom StarCCM+)
V donnée (nombre de Mach)
Direction de fluide donnée
Amont Free stream T
P

P=0
T
Aval Pressure outlet

Intrados Imperméable => composant de vitesse vertical =0

Extrados No slip
& Wall
Bord
d'attaque

"Dessous"

"Dessus"

Indiquer alors les conditions de vol retenues (conditions amont, c.a.d de l'écoulement non perturbé) :

Nombre de Mach M 0.45

Température T (K) 253
Pression p (Pa) 50500
Masse volumique  (kg/m3) 0.69
4
Vitesse U∞ (m/s) 143.48
Incidence (en °) 3 (0 dans le repère relatif)
II.4. Maillage

Sans regarder ce qui a été effectivement choisi pour la simulation, indiquer si certaines zones de
l'écoulement vont requérir des maillages plus fins que d'autres, et pourquoi.
Aux parois

Ouvrez la scène "maillage". Quel type de mailles est utilisé autour de la paroi?
Polyhedral, le maillage est composé de hexagone sauf qu’autour de la paroi où on a des mailles
rectangulaires de faible épaisseur dans le sens des gradients normaux à la paroi
Quelle est la taille de ces mailles? 0.0005 m
Vous semble-t-elle justifiée? ? .
Pourquoi (justifier par des calculs, en précisant bien les hypothèses de votre calcul)?
1 2 ¿2
τ = ρ U infini C f = ρU
2

U ¿ =U infini
√ Cf
2
Or C f dépend de l’écoulement dans la couche limite
1
Laminaire : C f = 0 ,5
R ex
1
Turbulent : C f = 0 ,2
R ex
ρU infini L 6
R e L =R e x (L)= =5.10 avec L=longeur de la corde
μ
 Turbulent sur 90% du profil
On prend la corrélation en turbulent et à mi-corde, pour avoir un C f moyen
6 0.0592 −3
R e L =2 ,5.10 → C f = =3 , 06.10
2
ReL
2

¿ −1
u =5 , 6 m. s
Si on veut utiliser les fonctions de paroi, on vise y +¿=50 ¿
¿
y∗u ∗ρ
Or y +¿= μ
¿

Donc y=2, 5.10−4

Or ce serait la valeur au centre de la maille -> correspondant
On dispose de corrélation pour évaluer y+ en assimilant le profil à une plaque plqne

5
On doit d’abord déterminer si la couche limite sera laminaire ou turbulente, pour ce la il faut calculer Re
basé sur la longueur (mi corde ici) parcourue : ReL/2 =(r0UL/2)

En conclusion, le maillage vous semble-t-il bien adapté au problème ?

Au cas de décollement le maillage ne sera pas adapté (ici petit incidence donc pas de décollement)

Faudrait-il le modifier en cas d'augmentation (typiquement jusqu'à 0,8) du nombre de Mach amont?
Pourquoi?
Oui, parce qu’il y aura un risque de l’onde de choc risque de faire décoller la couche limite

6
II.5 . Modèles et équations à résoudre

Sans regarder ce qui a été effectivement choisi pour la simulation, indiquer les équations à résoudre,
les conditions de résolution et le type de modèle :
- équations à résoudre : nombre, lesquelles, pourquoi ?
Les variables pour décrire l’écoulement :

- turbulent / laminaire ? Pourquoi ?

- stationnaire (steady) ou instationnaire (unsteady)? Pourquoi?

Vérifier, en suivant le cheminement ci-dessous, les équations et modèles effectivement utilisés.

Correspondent-ils à vos choix ci-dessus? Commentez les choix de modèles.

7
II.6 . Pour lancer le calcul : conditions initiales
Quelle est la distinction entre conditions aux limites et conditions initiales?

Sans regarder ce qui a été effectivement choisi pour la simulation, quelles conditions initiales seraient
envisageables quant à la vitesse ? à la pression ? Pour chaque cas proposé, indiquer si elles seraient
compatibles avec les conditions aux limites.

A l'aide du cheminement ci-contre, indiquer les conditions initiales choisies

Commentaire par rapport à vos propositions:

Pour affecter les conditions initiales, il est nécessaire d' "initialiser" le calcul : dans la barre de menu
"Solution  Initialize solution". L'avez-vous fait ?

8
 III. CALCUL

III.1. Préparer l'exploitation des résultats et le suivi du calcul : choix de valeurs de référence

Il est intéressant (voire nécessaire!) de mettre les résultats sous forme adimensionnelle, afin de pouvoir
critiquer les résultats du calcul par rapport à des valeurs de la littérature, de pouvoir les étendre à
d'autres configurations, etc. On s'intéresse ici aux caractéristiques aérodynamiques du profil : forces de
traînée et de portance, moment aérodynamique.
Le logiciel calcule-t-il des forces ou des coefficients de force?

Comment ces forces sont-elles calculées par le logiciel? Quelle composante de force est calculée ?

Rappeler la relation des forces et moment avec leurs coefficients adimensionnels.

Sur quel axe la résultante aérodynamique doit-elle être projetée pour obtenir la force de portance?
Même question pour la force de traînée. Vérifier que c'est bien le cas dans la simulation : cela est défini
dans les "Reports".

Donner les valeurs des paramètres utilisées ( dans les "Reports") pour adimensionner forces et
moments, et les critiquer (les justifier).

9
III.2. Suivi et convergence du calcul

Quels sont les différents critères pour juger de la convergence du calcul?

Des critères d'arrêt sont définis dans la simulation. Les préciser. Ces critères d'arrêt automatique sont-
ils suffisants (répondre en lien avec la question précédente)?

Des résultats vont aussi s'afficher dans la fenêtre de dialogue au cours du calcul : les résidus – de
façon automatique – et aussi les coefficients aérodynamique en raison des "Monitors" définis à partir
des "Reports". Afin de juger de manière plus visuelle de la convergence (et de l'évolution des résultats),
il est intéressant d'afficher ces résultats et résidus sous forme de graphes, les "plots" (graphes avec les
itérations en abscisse, traduisant de façon graphique les "monitors"). Ouvrez ceux qui sont déjà
préparés: résidus (s'ouvre tout seul sinon), les 3 coefficients aéro. pour qu'ils s'affichent en cours de
calcul.
Les avez-vous ouverts? : .

Un autre "Plot" est défini : le coefficient de pression, en fonction de la position le long de la corde. Il
permet de visualiser l'évolution de ce coefficient avec les itérations. Comment est défini le coefficient de
pression?

Vérifier la valeur des coefficients (chercher "Pressure coefficient" dans "ToolsField Functions",
cheminement ci-dessus) utilisés. Sont-ils conformes à ceux observés en III.1. ?

Enfin, avant de lancer le calcul, ouvrez, en double-cliquant sur les scènes, "Mach" et "pressure" que
vous pourrez ainsi visualiser en cours de calcul.
Sont-elles ouvertes ? . .

III.3. Calcul
Vous avez déjà initialisé le calcul (question II.6). Lancez-le (barre de menu "solutionRun" ou icône
avec le bonhomme qui court).

Observez les différents "plots" et scènes au cours du calcul.

10
IV. RESULTATS
IV.1. Convergence du calcul

Sur quel(s) critère(s) le calcul s'est-il arrêté ou avez-vous choisi de l'arrêter ?

Tous les critères définis en III.2 sont-ils respectés?

D'autres tests seraient ils nécessaires – éventuellement demandant des modifications profondes du
calcul, hors du cadre de ce TP – pour s’assurer de la validité du calcul ?

IV.2. Maillage et domaine ?

Il est nécessaire de faire certaines vérifications a posteriori.

Ouvrir le plot "Wall y+" qui vous renseigne sur la valeur de y+ au centre des mailles entourant le profil.
Les valeurs sont-elles conformes au(x) modèles(s) choisi ?

Vérifier les gradients près des conditions aux limites pour avoir une indication sur l'adéquation de la
taille du domaine.

11
IV.3. Aérodynamique

Coefficient de pression :
- dessiner ci-dessous l’évolution du coefficient de pression Cp en fonction de x autour du profil
- à quoi correspond x sur ce graphe ? . .
Cp

- y-a-t-il des valeurs remarquables de Cp ? Où et à quoi correspondent-elles ?

- commenter les variations de ce coefficient sur l’extrados et l’intrados.

Coefficients aérodynamiques et finesse :

- donner la valeur du coefficient de portance C L. Commenter cette valeur, par rapport aux ordres de
grandeur donnés en cours.

- donner la valeur du coefficient de traînée CD. Commenter de même.

- donner la valeur de la finesse f. Commenter de même.

12
 V. POURSUITES DU CALCUL

V.1. Evolution des coefficients aérodynamique avec l'incidence – polaire

Recommencez le calcul avec l'angle d'incidence  donné (à votre demande!) par l'enseignant. Les
résultats de tout le groupe permettront de tracer les graphes C L=f() et CD=f() ainsi que la polaire du
profil dans ces conditions de vol.

V.2. Influence de la compressibilité

S'il reste du temps, enregistrez la simulation sous un nouveau nom. Modifiez les "Models" et les
conditions aux limites de façon à lancer la simulation en supposant le gaz incompressible. Relancer le
calcul. Commenter les différences observées et conclure.

V.3. Influence de la vitesse/ du nombre de Mach

S'il reste du temps, relancer le calcul (avec le gaz compressible) avec des nombres de Mach de 0,8 et
2. Observez les résultats puis discutez avec l'enseignant.

13