Rapport P6 2021 35
Rapport P6 2021 35
Rapport P6 2021 35
STPI/P6/2021 – 35
Étudiants :
ANG Cécile
COSTE Lucien
KERRIEN Antoine
LE GAC Tanguy
NEEL Caroline
Objectifs du projet :
Mots-clés du projet :
❖ Écoulements turbulents
❖ Reynolds
❖ Navier-Stokes
❖ Simulation
STPI/P6/2021 - 35
3
STPI/P6/2021 - 35
4
1. INTRODUCTION
La mécanique des fluides est une branche de la physique qui consiste à étudier l’écoulement
des fluides, qui peuvent être en mouvement ou au repos, et qui sont généralement soumis à
des forces et des contraintes. Un fluide est un milieu matériel parfaitement déformable. Les
liquides, les gaz, voire les plasmas à très haute température sont qualifiés de fluides.
De plus, cette notion se divise en deux branches : d’une part, la statique des fluides pour les
fluides au repos, et d’autre part, la dynamique des fluides qui concerne les fluides en
mouvement. En effet, cette dernière consiste à comprendre, anticiper et calculer les
écoulements et mouvements des fluides. C’est cette notion qui nous intéressera dans le cadre
de notre projet.
Notre objectif est de modéliser et de simuler des fluides en écoulement turbulent à travers
différents volumes. La turbulence caractérise l’imprévisibilité d’un écoulement : plus celui-ci
sera turbulent, plus il sera difficile d’anticiper son évolution et plus il sera désordonné.
De surcroît, le domaine d’étude est vaste, et de nombreux théorèmes et lois peuvent être
appliqués pour étudier l’écoulement des fluides. Cependant, la théorie peut se heurter à
l’expérience car en effet, la modélisation des écoulements peut s’avérer invalide ou imprécise
comparée à l’étude expérimentale et pratique d’un fluide en mouvement.
L’écoulement des fluides est régi par de nombreuses équations, notamment les équations de
Navier-Stokes. Toutefois, il faut savoir que ces équations sont impossibles à résoudre
analytiquement. De plus, la viscosité, parfois négligée dans les cas idéaux, influence la
stabilité des fluides, et ils deviennent alors plus ou moins imprévisibles.
On comprend alors que théorie et expérience sont étroitement liées et complémentaires.
STPI/P6/2021 - 35
5
2. ORGANISATION DU TRAVAIL
La première partie de notre projet est constituée d’un travail de documentation sur la
mécanique des fluides, afin de se familiariser avec les notions principales du sujet
puisqu’auparavant nous n’avions que peu de connaissances en mécanique des fluides. Ainsi,
au début du projet, notre professeur référent, M. Safdari Shadloo, nous a fourni des sources
à consulter (voir bibliographie) ainsi que des cours vidéo à regarder, et de notre côté, nous
devions compléter ces ressources et approfondir nos connaissances par des recherches
personnelles. En effet, pour comprendre exactement en quoi consiste notre projet, il est
nécessaire de comprendre et d’assimiler les principes et théories de la dynamique des fluides.
Pour la partie documentation, nous avons donc tous fait des recherches chacun de notre côté
puis chaque semaine on essayait de mettre en commun nos différentes recherches, afin de
s’assurer que tout le monde comprenait bien les notions principales, puis on envoyait un
compte-rendu de ce que nous avions retenu à notre enseignant référent. Nous nous sommes
également regroupés tous ensemble chaque semaine avec M. Safdari Shadloo, afin
d’expliciter les points importants de notre étude, puis d’orienter notre réflexion et nos
recherches sur de nouvelles notions.
Nous sommes ensuite passés à la deuxième partie du sujet qui consiste à modéliser et simuler
des écoulements turbulents sur le logiciel Ansys Fluent. Il a également fallu se documenter
sur le logiciel pour comprendre son fonctionnement et apprendre à l’utiliser, puisqu’aucun
d’entre nous ne le connaissait. Pour cela, notre professeur nous a fourni des vidéos nous
montrant comment l’utiliser. Lorsque nous avions des modélisations à effectuer, nous
essayions de nous réunir à 3 ou 4 minimum en fonction des disponibilités de chacun et nous
rédigions un compte-rendu de la simulation pour en discuter avec M. Safdari Shadloo.
Concernant la répartition du travail, nous avons partagé équitablement les tâches entre les
différents membres du groupe, tout en s’assurant que chaque personne abordait au moins
une fois chaque aspect principal du projet, puisque nous estimons important de comprendre
et ainsi de pouvoir expliquer l’intégralité du projet pour chaque membre du groupe.
De plus, nous avons parfois dû faire face à des difficultés mais nous avons su les surmonter.
En effet, comme dit précédemment, n’ayant aucune connaissance du logiciel Ansys Fluent,
nous avons dû travailler tous ensemble afin bien manipuler le programme pour réussir nos
simulations.
STPI/P6/2021 - 35
6
3.1.1. Généralités
Tout d’abord, un fluide incompressible est un fluide dont on ne peut pas changer le volume,
on ne peut pas le comprimer dans un espace plus restreint, à l’inverse d’un fluide
compressible. Les liquides sont des fluides incompressibles.
De plus, un fluide newtonien est un fluide dont la viscosité ne dépend pas des contraintes qui
lui sont appliquées, c’est-à-dire que sa viscosité ne varie pas quand on le met en mouvement.
Cependant, la viscosité peut dépendre d’autres facteurs comme la température. En ce qui
concerne un fluide non-newtonien, elle dépend en partie de la contrainte de cisaillement qui
lui est appliquée. De surcroît, il faut savoir qu’il existe deux types de contraintes possibles sur
une face d’un matériau : tangentielle et normale. La contrainte de cisaillement est appliquée
de manière tangentielle sur la surface du fluide.
Dans un écoulement laminaire, le fluide s’écoule dans un tuyau comme s’il glissait
parallèlement aux parois qui le guident, et les lignes d’écoulement sont alors parallèles, tandis
que dans un écoulement turbulent, le fluide s’écoule avec des variations de vitesse brusques
et aléatoires en chaque point. De même, les lignes d’écoulement sont désordonnées et au
contraire de l'écoulement laminaire, elles ne sont pas linéaires.
Pour commencer, un solide peut résister à une contrainte de cisaillement par une déviation
statique, au contraire d’un fluide, qui va lui se déplacer et se déformer continuellement sous
l’effet de la contrainte de cisaillement qui lui est appliquée. Il existe deux types de fluide : les
liquides et les gaz. Un liquide est défini par ses molécules proches et à forte force de cohésion,
il tend à conserver son volume. Au contraire, le gaz avec ses molécules espacées et à force
de cohésion négligeable, est libre de se dilater ou de se contracter. Il n’a donc pas de volume
défini.
Cependant, il existe de nombreux cas limites où la matière peut avoir l’apparence d’un solide
et peut résister d’abord aux contraintes de cisaillement mais finir par se déformer lentement
en présentant un comportement fluide (bitume, plomb). D’autres substances résistent sous
l’effet d’un faible cisaillement, mais cèdent et s’écoulent comme des fluides particulièrement
lorsqu'elles subissent une contrainte plus importante (mélanges colloïdes et boueux).
Néanmoins, la distinction entre liquide et gaz est parfois floue. En effet, à des températures
et pressions supérieures au point critique d’un fluide, lorsqu’une seule phase existe
(ressemblant principalement à un gaz), on ne peut pas utiliser les modèles physiques tels que
la loi des gaz parfaits. Ainsi, on fera l’hypothèse dans notre étude que les phases liquides et
gazeuses sont clairement définies.
Dans un premier temps, la viscosité mesure la résistance à l’écoulement d'un fluide. En effet,
le coefficient de viscosité mesure l’aptitude d’un fluide à résister aux contraintes. La viscosité
des fluides newtoniens varie en fonction de la température et de la pression. Il faudra donc
dans notre étude bien prendre en compte ces notions. Effectivement, la température a un effet
important sur la viscosité, tandis que l’effet de la pression est plus modéré. Pour donner une
STPI/P6/2021 - 35
7
idée un peu plus imagée de la viscosité, nos mouvements sont ralentis lorsque l’on est
immergé dans l’eau, ceci s’explique par le fait que la viscosité de l’eau est nettement
supérieure à celle de l’air, ainsi l’eau va mieux résister à nos mouvements que l’air.
Dans un second temps, la notion de masse volumique est importante. Celle-ci représente la
𝑑𝑚
masse du fluide par unité de volume : 𝜌 = 𝑑𝑉 . Elle s’exprime en kg/m3. La densité est
également une notion élémentaire de la mécanique des fluides. Elle se calcule en faisant le
rapport entre la masse volumique d’un matériau et celle de l’eau.
Pour finir, les propriétés thermodynamiques sont aussi des propriétés importantes pour les
fluides avec la température, la pression, les capacités thermiques…
Tout d’abord, le champ de vitesse fait partie intégrante des propriétés d’un écoulement. Il est
principalement important pour les problèmes de débit d’un fluide. Le vecteur vitesse est alors
noté 𝑉 = 𝑖 ∗ 𝑢(𝑥, 𝑦, 𝑧, 𝑡) + 𝑗 ∗ 𝑣(𝑥, 𝑦, 𝑧, 𝑡) + 𝑘 ∗ 𝑤(𝑥, 𝑦, 𝑧, 𝑡). En ce qui concerne le vecteur
accélération du fluide, il est obtenu en dérivant le vecteur vitesse par rapport au temps ce qui
donne :
𝑑𝑉 𝑑𝑢 𝑑𝑣 𝑑𝑤 𝜕𝑉 𝜕𝑉 𝜕𝑉 𝜕𝑉 𝜕𝑉
𝑎= 𝑑𝑡
= 𝑖 𝑑𝑡 + 𝑗 𝑑𝑡 + 𝑘 𝑑𝑡
= 𝜕𝑡
+ (𝑢 𝜕𝑥 + 𝑣 𝜕𝑦 + 𝑤 𝜕𝑧 ) = 𝜕𝑡
+ (𝑉. 𝛻)𝑉,
𝜕 𝜕 𝜕
avec 𝛻 l’opérateur gradient 𝛻 = 𝑖 𝜕𝑥 + 𝑗 𝜕𝑦 + 𝑘 𝜕𝑧, et 𝑖, 𝑗 et 𝑘 les vecteurs unitaires du repère
utilisé.
STPI/P6/2021 - 35
8
STPI/P6/2021 - 35
9
expériences. On peut ajouter qu’elle est plus fiable pour les fluides en écoulement laminaire,
monophasé et chimiquement stable.
Il existe plusieurs conditions limites. La condition de non-glissement, qui existe pour un fluide
réel, et pour laquelle les particules fluides adhèrent à la paroi. La vitesse est donc nulle proche
de la paroi, on a alors : 𝑢 = 𝑣 = 𝑤 = 0. Pour un fluide non-visqueux, aussi appelé fluide parfait,
la condition s’atténue. Une autre condition à prendre en compte concerne la température.
Effectivement, si la température de surface est constante, on a T=T(surface). Cependant, si
on ne connaît pas cette température de surface, on a alors une condition limite liée au flux de
𝜕𝑇 −𝑞
chaleur entre la paroi et le fluide, sous forme de gradient à la paroi : (𝜕𝑛) = 𝑤
avec 𝑛 la
direction normale à la surface et 𝑞 la conductivité thermique. Ce sont des conditions spatiales
du domaine de calcul.
Malgré tout, pour résoudre les équations, il est nécessaire d’imposer des conditions à
l’écoulement. Les conditions aux limites cinématiques sont des conditions que doit vérifier le
champ de vitesse et les conditions aux limites dynamiques doivent être vérifiées par les
champs de contrainte et de vitesse aux frontières du domaine.
La résolution d’un problème de dynamique des fluides à l’aide de la méthode numérique est
constituée de 4 étapes successives :
• Identification du problème en 2 étapes : Identification du type de problème, du type
d’écoulement… Il faut ensuite choisir le modèle mathématique qui correspond le mieux
et définir le domaine d'étude ainsi que sa géométrie.
• Prétraitement :
o Génération du maillage : Définition de tous les paramètres liés
à la maille comme leur nombre, leur taille ou leur forme.
o Discrétisation de l’espace : approximer les dérivées spatiales
en fonction des volumes, des particules...
o Discrétisation du temps : approximer les dérivées temporelles.
• Simulation CFD :
o Logiciel de CFD : choix des algorithmes, des outils d’algèbre
linéaire ou encore du langage de programmation. Figure 2 : Exemple
o Lancement de la simulation : la qualité du résultat dépend de de maillages
tous les paramètres cités précédemment mais aussi du critère différents pour une
d'arrêt de la simulation. même géométrie
STPI/P6/2021 - 35
10
• Post-traitement :
o Analyse : Basée sur des calculs ou des visuels extraits de la simulation (à
l’aide de couleurs, en 1, 2 ou 3 dimensions), elle a pour but la visualisation et
l’analyse des résultats, et la validation du modèle utilisé.
o Vérification : il faut faire attention à bien faire la différence entre l’incertitude
et l’erreur. L’incertitude est la déviation potentielle. Au contraire, l’erreur est la
déviation reconnue (attention : une petite erreur locale peut entraîner une
erreur globale importante). Il existe différents types d’erreurs : erreurs connues
(modélisation physique, erreurs de discrétisation, convergence itérative,
erreurs d’arrondi) et erreurs inconnues (erreurs de programmation,
d’utilisation…). On vérifie aussi le code CFD et le modèle utilisé pour s'assurer
d’éliminer les éventuelles sources d’erreur. On peut également comparer les
résultats obtenus par simulation aux résultats expérimentaux si une expérience
a également été faite. Pour cela il faut examiner la convergence du résultat
obtenu, examiner la cohérence, examiner la grille.
Cette méthode est intéressante du point de vue rentabilité. En effet, elle est moins chère et
plus rapide qu'une étude expérimentale. Cependant, elle n’est pas fiable à 100% (imprécision,
modèle peut être inadéquat, limite de la puissance de calcul, mauvaise modélisation du
problème...), il faut donc bien vérifier les résultats obtenus et les hypothèses initiales. Pour
finir, elle est utilisée pour l'analyse et la conception, la connaissance et l'exploration de la
physique des fluides.
Il existe trois types de systèmes de coordonnées pouvant être utilisés pour modéliser des
fluides. Le premier est le système Eulérien, il consiste à suivre les caractéristiques de
l’écoulement dans un volume fixe par exemple, un écoulement à travers une boîte. Le second
est le système Lagrangien, il est basé sur les particules de fluide qui bougent dans le champ
d’écoulement. Pour finir, il existe le système mixte Eulérien-Lagrangien. Quel que soit le
modèle utilisé, les équations obtenues doivent être les mêmes.
Pour déterminer le type d’écoulement d’un fluide on utilise le nombre adimensionnel dit de
Reynolds d’équation :
𝜌𝑉𝐷
𝑅𝑒 =
𝜇
avec 𝜌 la densité, 𝑉 la vitesse, 𝐷 le diamètre du conduit dans lequel s’écoule le fluide et 𝜇 la
viscosité dynamique du fluide. Il est nécessaire pour la détermination de la nature de
l’écoulement étudié.
Si Re>4000 alors l’écoulement est turbulent (imprévisible, mélanges
rapides).
Si 2300<Re<4000 alors c’est un état de transition (c’est-à-dire pas
complètement désordonné, éclatement turbulent).
Si Re < 2300 alors l’écoulement est laminaire (prévisible, mélange lent).
Figure 3 : schéma
correspondant aux trois
types d'écoulement
La turbulence d’un fluide se caractérise par son irrégularité, sa diffusivité,
sa tendance rotative, son caractère dissipatif avec un besoin d’apport
d’énergie constant et de grandes échelles de temps et d’espace.
Justement, les échelles de Kolmogorov permettent d’obtenir une échelle de temps, de
longueur et de vitesse. Elles permettent donc de définir les écoulements turbulents.
STPI/P6/2021 - 35
11
𝜈3
• Echelle de longueur de Kolmogorov : 𝜂 = ( )1/4
𝜀
𝜈
• Echelle de temps de Kolmogorov : 𝜏𝜂 = (𝜀 )1/2
• Echelle de vitesse de Kolmogorov : 𝑢𝜂 = (𝜈𝜀) 1/4
avec 𝜀 qui représente le taux moyen de dissipation de l’énergie cinétique de turbulence par
unité de masse, et 𝑣 est la viscosité cinématique du fluide.
Premièrement, la méthode DNS n’est basée sur aucun modèle de turbulence. Elle doit donc
résoudre numériquement les échelles de Kolmogorov ainsi que les équations de Navier-
Stokes, ce qui mène à des calculs assez compliqués faisant intervenir de nombreuses
variables et des échelles très réduites.
C’est donc une méthode assez lente car elle demande beaucoup de calculs et donc,
beaucoup de ressources pour effectuer la simulation. Elle reste donc une méthode
relativement chère. Mais la méthode DNS est la plus précise des trois elle permet d’avoir
beaucoup de détails sur le système. La faisabilité de cette méthode repose sur le fait que les
plus petites échelles des écoulements turbulents ne sont pas infinies et peuvent donc être
modélisées. La simulation DNS est donc assez précise pour pouvoir simuler ces petites
échelles de temps et d’espace. Elle est basée sur une relation entre le nombre de Reynolds
3
𝐿 𝜈3
et les échelles de Kolmogorov : ≈ 𝑅𝑒 4 avec 𝜂 = ( )1/4 et 𝐿 une dimension caractéristique
𝜂 𝜀
de la géométrie du problème. L’exigence de stockage dans une simulation comme celle-ci
augmente rapidement avec le nombre de Reynolds.
Pour la plupart des applications industrielles, cette méthode requiert trop de ressources et est
donc inutilisable. Mais il est vrai que pour certains nombres de Reynolds assez peu élevés,
cette méthode est utilisable et les écoulements simulés présentent alors des caractéristiques
très proches de la réalité.
Cette méthode a été introduite par Smagorinsky en 1963. Le principe de ce modèle est le
suivant : en appliquant un “filtre passe-bas” aux équations de Navier-Stokes, on élimine les
plus petites échelles de longueurs (c’est à dire les petits tourbillons), nous laissant ainsi un
système d’équations simplifié et plus rapide à résoudre. Ce filtre peut être appliqué sur les
champs spatial et temporel, et a pour expression :
∞ ∞
̅̅̅̅̅̅̅̅̅
𝜙(𝑥, 𝑡) = ∫ ∫ 𝜙(𝑟, 𝜏)𝐺(𝑥 − 𝑟, 𝑡 − 𝜏)𝑑𝜏𝑑𝑟
−∞ −∞
STPI/P6/2021 - 35
12
𝐿 𝐿
DNS Δ= Δ=
32 16
Figure 4 : comparaison du champ de vitesse obtenu avec la méthode DNS puis avec deux filtres LES
Les conditions limites pour un problème traité grâce à la méthode LES jouent un rôle important
sur les résultats. En effet, elles doivent être assez simples pour être modélisées, apporter
suffisamment d’informations sur l’écoulement du fluide (on peut penser à la vitesse ou la
viscosité) et être solutions des équations de Navier-Stokes. Afin d’obtenir les meilleurs
résultats, les conditions limites devront être choisies avec le plus grand soin.
Pour étudier les mouvements des fluides, l’approche aux équations de Navier-Stokes
moyennées (RANS) est la plus répandue parmi les logiciels de simulation dans le secteur
industriel, car cette méthode présente un temps de calcul plus court que les méthodes LES et
DNS, et n’a pas besoin de recourir à des maillages extrêmement fins. Cependant, cette
méthode est très empirique et ne s’appuie pas sur de solides démonstrations théoriques.
De plus, la méthode RANS s’appuie principalement sur ce qu’on appelle la décomposition de
Reynolds, qui consiste, pour une grandeur 𝑓, à la décomposer de la façon suivante :
𝑓 = 𝑓 + 𝑓′, avec 𝑓 la valeur moyenne, dite moyenne de Reynolds et 𝑓′ qui représente les
fluctuations, c’est-à-dire l’écart par rapport à la moyenne.
La méthode RANS va utiliser cette décomposition pour toutes les variables de l’équation de
Navier-Stokes (𝑢, 𝑝, 𝑇). Ainsi, quand un écoulement devient turbulent, les équations de Navier-
STPI/P6/2021 - 35
13
Stokes, même si elles sont toujours valables, ne sont plus utilisables car il devient trop
compliqué de suivre tous les détails de cet écoulement. La décomposition de Reynolds permet
alors de simplifier les équations de Navier-Stokes pour donner les équations de Navier-Stokes
moyennées, qui prennent en compte la moyenne de Reynolds (Reynolds-averaged Navier-
Stokes equations - RANS).
Pour un fluide newtonien incompressible en écoulement stationnaire, les équations de Navier-
Stokes moyennées en coordonnées cartésiennes sont :
𝜕𝑢𝑖 𝜕 𝜕𝑢𝑖 𝜕𝑢𝑗
𝜌𝑢𝑗 = 𝜌𝑓𝑖 + [−𝑝𝛿𝑖𝑗 +𝜇( + ) − 𝜌𝑢𝑖′ 𝑢𝑗 ′].
𝜕𝑥𝑗 𝜕𝑥𝑗 𝜕𝑥𝑗 𝜕𝑥𝑖
On peut remplacer la viscosité 𝜇 par le tenseur des contraintes 𝜎𝑖𝑗 , ce qui donne l’équation de
Reynolds :
𝜕𝑢𝑖 𝜕𝑢𝑖 𝜕
𝜌( + ∑ 𝑢𝑗 ( )) = 𝜌𝑔𝑖 + ∑( )(𝜎𝑖,𝑗 − 𝜌𝑢𝑖′ 𝑢𝑗 ′)
𝜕𝑡 𝜕𝑥𝑗 𝜕𝑥𝑗
𝑗 𝑗
• Les modèles du second ordre dans lesquels les tensions de Reynolds sont
directement calculées.
Les modèles du premier ordre sont appelés modèles à viscosité turbulente. Ils traitent de la
dissipation de l’énergie cinétique et se concentrent sur l’énergie dissipée à un niveau
physiquement correct. Dans les modèles à viscosité turbulente, on néglige les petites cellules
de la maille tout en considérant qu’elles sont toujours présentes lors de la modélisation, ce
qui permet de ne pas avoir à les simuler. A noter qu’il existe plus de 200 modèles de ce type,
qui sont classés en fonction du nombre d’équations de transport résolues en plus des
équations de RANS : modèles à une équation (modèle de Prandtl, modèle de Baldwin-Barth,
modèle de Spalart-Allmaras…), modèles à deux équations (modèle k-𝜀, modèle k-𝜔, modèle
SST…).
Les modèles du second ordre sont appelés modèles de transport de la contrainte de
Reynolds. Ils résolvent les équations de transport des tensions de Reynolds 𝑢𝑖′ 𝑢𝑗 ′et des flux
de chaleur turbulents 𝑢𝑖′ 𝑇′, qui sont longues et compliquées à résoudre. Ces modèles du
second ordre sont en effet plus compliqués à mettre en œuvre que les modèles du premier
ordre. Effectivement, ils sont recommandés pour les écoulements turbulents complexes en
3D et réclament beaucoup de ressources.
STPI/P6/2021 - 35
14
• Modèle k-𝜀
Ce modèle est utilisé pour les flux turbulents. Celui-ci fonctionne sur un système de deux
équations dites de transport pour ainsi donner une description la plus complète de la
turbulence.
La première variable est l’énergie cinétique turbulente notée k :
𝜕(𝜌𝑘) 𝜕(𝜌𝑘𝑢𝑖 ) 𝜕 𝜇𝑡 𝜕𝑘
+ = [ ] + 2𝜇𝑡 𝐸𝑖𝑗 𝐸𝑖𝑗 − 𝜌𝜀
𝜕𝑡 𝜕𝑥𝑖 𝜕𝑥𝑗 𝜎𝑘 𝜕𝑥𝑗
La seconde variable est le taux de dissipation de l’énergie cinétique notée 𝜀 :
𝜕(𝜌𝜀) 𝜕(𝜌𝜀𝑢𝑖 ) 𝜕 𝜇𝑡 𝜕𝜀 𝜀 𝜀²
+ = [ ] + 𝐶1𝜀 2𝜇𝑡 𝐸𝑖𝑗 𝐸𝑖𝑗 − 𝐶2𝜀 𝜌
𝜕𝑡 𝜕𝑥𝑖 𝜕𝑥𝑗 𝜎𝜀 𝜕𝑥𝑗 𝑘 𝑘
Ce modèle est utile pour les couches de cisaillement planaire et les flux de recirculation. Il est
également très utile pour les fluides turbulents car il nécessite seulement de connaître les
conditions limites. C’est pour cela qu’il est utilisé dans des applications telles que les
écoulements industriels ou environnementaux. Cependant, il est inapproprié de l’utiliser dans
de nombreux cas comme : les flux non définis, les couches limites circulaires, les flux rotatifs...
Il est toutefois utile pour avoir une première représentation globale d’un écoulement.
• Modèle k-𝜔
Ce modèle cherche à prédire la turbulence grâce à deux équations pour deux variables,
𝜀
l’énergie cinétique de turbulence k et la dissipation spécifique 𝜔 telle que 𝜔 ≈ :
𝑘
𝜕(𝜌𝑘) 𝜕(𝜌𝑢𝑗 𝑘) 𝜕 𝜌𝑘 𝜕𝑘 𝜕𝑢𝑖
𝜕𝑡
+ 𝜕𝑥 = 𝜌𝑃 − 𝛽 ∗ 𝜌𝜔𝑘 + 𝜕𝑥 [(𝜇 + 𝜎𝑘 ) ] avec 𝑃
𝜔 𝜕𝑥𝑗
= 𝜏𝑖𝑗 𝜕𝑥𝑗
𝑗 𝑗
STPI/P6/2021 - 35
15
Tout d’abord, la discrétisation spatiale des solveurs Fluent est basée sur la méthode des
volumes finis qui exploite des approximations d’intégrales pour effectuer les calculs.
𝜕
∫ 𝜌𝛷𝑑𝑉 + ∮ 𝜌𝛷𝑉. 𝑑𝐴 = ∮ 𝛤𝛻𝛷. 𝑑𝐴 + ∫ 𝑆𝛷 𝑑𝑉
𝜕𝑡 𝑉 𝐴 𝐴 𝑉
Le maillage est alors constitué de volumes finis (ou surfaces finies en 2D et segments en 1D)
qui sont des petits volumes disjoints. L’équation générale du transport est appliquée à chaque
cellule et discrétisée. Les équations discrétisées nécessitent des informations sur les faces et
au centre de chaque cellule. Les équations discrétisées peuvent être exprimées de la manière
suivante :
𝑎𝑝 𝜙𝑝 + ∑ 𝑎𝑛𝑏 𝜙𝑛𝑏 = 𝑏𝑝
𝑛𝑏
Elles sont résolues itérativement. Les coefficients 𝑎𝑝 et 𝑎𝑛𝑏 sont des fonctions des variables
solutions. Ces coefficients sont écrits pour utiliser les variables solutions des itérations
précédentes.
A noter que les solutions peuvent être explicites ou implicites. Par exemple, pour la
linéarisation explicite, la valeur inconnue dans chaque cellule est calculée à partir de relations
incluant uniquement des valeurs connues. A contrario, pour la linéarisation implicite, 𝜙𝑝 et
𝜙𝑛𝑏 sont supposées inconnues et résolues en utilisant des techniques de résolution
d’équations linéaires. Les équations peuvent également être discrétisées grâce à des
méthodes d’interpolation. De plus, les procédés d’interpolation peuvent être utilisés pour
déterminer les termes convectifs ou les gradients des variables solutions.
En définitive, dans le but de coupler la pression et le champ de vitesse, Ansys Fluent propose
deux approches numériques : l’une qui utilise un solveur basé sur la pression tandis que l’autre
utilise un solveur basé sur la densité. Dans les deux cas, le champ de vitesse est obtenu à
partir de l’équation du moment.
Avant d’utiliser directement Ansys Fluent, nous sommes passés par l’outil Ansys Workbench.
En effet, cette interface est très utile pour effectuer une à plusieurs simulations, puisqu'elle
nous indique les principales étapes à suivre et les différents logiciels Ansys à utiliser pour
effectuer une modélisation complète en partant d’une page blanche.
Ainsi, pour un modéliser un écoulement turbulent avec Ansys Fluent, Les 5 étapes à suivre
sont :
• Création de la géométrie : cette étape consiste à créer la géométrie 3D de notre
problème en utilisant un programme appelé Ansys DesignModeler mais il est
également possible d’importer des géométries depuis des logiciels de CAO comme
SolidWorks ou Catia.
• Génération du maillage : on crée un maillage à partir de notre géométrie. Il est
important de trouver le bon compromis : un maillage précis donnera de meilleurs
résultats mais augmentera les temps de calcul ; à l’inverse, un maillage grossier
donnera des résultats approximatifs, mais plus rapidement.
STPI/P6/2021 - 35
16
Pour cette configuration, nous avons décidé de nous concentrer sur l’évolution de la
température en sortie de la maille. Les modèles de simulation que nous avons utilisés pour
cette géométrie sont les modèles RANS k-𝜀, k-𝜔 standard, k-𝜔 SST ainsi que le modèle LES.
En plus de faire varier le modèle utilisé, nous avons effectué plusieurs simulations pour
plusieurs rayons pour l’entrée dans le tube : 12,5 et 25 millimètres pour ainsi observer
l’influence du maillage sur les propriétés thermodynamiques.
• Modèle k-𝜀
o Rayon 12,5 mm
Figure 5 : Courbe de température et simulation du champ de température dans le tube avec le modèle k-ε, pour
un rayon de 12,5 mm pour la petite entrée
On obtient une courbe de température qui converge rapidement, après environ 125 itérations,
vers une température de 301 K et on remarque qu’en sortie de la maille la température est
plus chaude sur la partie droite du tube que sur la partie gauche, avec un écart de l’ordre de
6 K environ.
STPI/P6/2021 - 35
17
o Rayon 25 mm
Figure 6 : Courbe de température et simulation du champ de température dans le tube avec le modèle k-ε, pour
un rayon de 25 mm pour la petite entrée
Figure 7 : Courbe de température et simulation du champ de température dans le tube avec le modèle k-𝜔, pour
un rayon de 12,5 mm pour la petite entrée
On obtient quasiment les mêmes résultats qu’avec la méthode précédente, avec une
température qui converge vers 301 K aussi. Cependant, dans ce cas-ci, c’est avec un nombre
d’itérations plus important (150 itérations), soit 25 de plus qu’avec la méthode précédente que
la température converge. La méthode k-𝜔 est plus précise près des parois ce qui peut
expliquer les itérations supplémentaires avant la convergence, même s’il n’y a pas beaucoup
plus d’itérations par rapport à la méthode k-𝜀. De plus, le profil de température est très
ressemblant au précédent et on observe les mêmes écarts.
o Rayon 25 mm
Figure 8 : Courbe de température et simulation du champ de température dans le tube avec le modèle k-ω, pour
un rayon de 25 mm pour la petite entrée
STPI/P6/2021 - 35
18
La température converge vers 309,5 K, c’est très proche de ce qu’on a obtenu avec le modèle
k-𝜀. Les écarts de température en sortie entre la partie gauche et la partie droite sont encore
plus marqués pour le grand rayon (jusqu’à 16 K d’écart).
• Modèle standard de turbulence k-ω SST
o Rayon 12,5 mm
Figure 9 : Courbe de température et simulation du champ de température dans le tube avec le modèle k-ω SST,
pour un rayon de 12,5 mm pour la petite entrée
La température converge vers 300,5 K après 125 itérations, soit autant d’itérations que pour
le modèle k-𝜀. Effectivement, cela est dû au fait que le modèle k-𝜔 SST s’appuie sur le modèle
k-𝜀, qui est utile pour modéliser l’écoulement loin des limites de la géométrie. De plus, on
retrouve un écart pouvant aller jusqu’à 6 K en sortie du tube mais le profil est représenté plus
précisément qu’avec le modèle k-𝜀 au niveau de la paroi du tube. Ainsi, on retrouve bien les
caractéristiques du modèle k-𝜔 SST qui utilise les deux modèles vus au-dessus.
o Rayon 25 mm
Figure 10 : Courbe de température et simulation du champ de température dans le tube avec le modèle k-ω SST,
pour un rayon de 25 mm pour la petite entrée
Comme pour les modèles précédents avec le grand rayon, la température converge vers 309
K environ. Le profil de température est similaire à celui du k-𝜔 standard avec les mêmes écarts
allant jusqu’à 16 K.
STPI/P6/2021 - 35
19
• Modèle LES :
Nous avons fait une modélisation avec ce modèle qui n’utilise pas les équations moyennées
afin de voir si les différences sont importantes.
o Rayon 12,5 mm
Figure 11 : Courbe de température et simulation du champ de température dans le tube avec le modèle LES, pour
un rayon de 25 mm pour la petite entrée
Les résultats obtenus ne convergent pas et ne sont donc pas exploitables. Ce modèle ne
permet pas dans les conditions du problème d’étudier la température en sortie du tube avec
cet écoulement. Les exemples que nous avons traité avec notre professeur se concentrent
sur des modèles basés sur la méthode aux équations moyennées (RANS). De plus, il faut
admettre que nous maîtrisons moins cette méthode.
o Rayon 25 mm
Figure 12 : Courbe de température et simulation du champ de température dans le tube avec le modèle LES, pour
un rayon de 25 mm pour la petite entrée
Il n’y a toujours pas de convergence avec cette méthode, même en augmentant le rayon. Pour
conclure, on peut considérer que pour notre étude, la méthode LES n’est pas la plus adaptée
et ne permet pas de modéliser correctement cet écoulement.
Nous avons ensuite étudié un écoulement externe autour d’un profil aérodynamique en forme
d’aile d’avion en utilisant des modèles adaptés à cette étude. Les modèles de simulation que
nous avons utilisés pour cette géométrie sont les modèles RANS k-𝜀, k-𝜔 standard, k-𝜔 SST
ainsi que le modèle Spalart-Allmaras. Nous avons donc étudié l’impact du modèle sur la
simulation de l’écoulement. De plus, pour notre étude, nous nous concentrons sur les
coefficients de traînée et de portance du fluide en question, ici l’air.
Pour chaque modélisation, nous avons choisi de simuler 300 itérations afin d’observer
l’évolution des coefficients de traînée et de portance sur une étude assez conséquente. Afin
d’étudier la divergence des résultats et les différences entre les modélisations, nous avons
décidé de relever la valeur des coefficients de traînée et de portance au bout de la 300ème
itération. Ainsi, nous observons qu’au terme des différentes simulations, le coefficient de
portance converge toujours vers la même valeur : 0,72, tandis que le coefficient de trainée
STPI/P6/2021 - 35
20
converge quant à lui vers 0,05. Ainsi malgré le nombre d’itérations et les différents modèles
choisis, ces deux coefficients convergent toujours vers les mêmes valeurs. Nous pouvons
donc en déduire que le résultat obtenu est le même quel que soit le modèle utilisé pour un
écoulement externe.
De plus, nous avons aussi relevé les valeurs obtenues pour coefficient de quantité de
mouvement au terme des 300 itérations, selon les différentes modélisations. De nouveau,
nous retrouvons le même résultat, même si le modèle de turbulence n’est pas toujours le
même, la valeur du coefficient de quantité de mouvement converge vers la même valeur.
Comparaison des valeurs obtenues : Tableau récapitulatif des résultats obtenus avec les
différents modèles pour l’écoulement externe :
Il est intéressant de constater que pour cet écoulement, les coefficients de portance et de
traînée sont les paramètres les plus importants de la modélisation. Ils dépendent de
l’écoulement et ne peuvent être modifiés dans AnsysFluent. Si le maillage ou la géométrie de
l’écoulement étudié change, ces coefficients seront également modifiés en fonction des
nouveaux paramètres du fluide sur l’aile. Les caractéristiques de la simulation comme le
modèle utilisé n’ont donc que peu d’influence pour cet écoulement externe.
STPI/P6/2021 - 35
21
4. CONCLUSION ET PERSPECTIVES
D’une part, grâce à ce projet, nous avons étudié et approfondi nos connaissances en
mécanique des fluides, et plus particulièrement sur les fluides en écoulement turbulent.
D’autre part, cette étude nous a permis de nous familiariser avec un logiciel de modélisation,
Ansys Fluent, que plusieurs d’entre nous continuerons d’utiliser au cours de leur formation à
l’INSA. Nous nous sommes ainsi confrontés à des problèmes complexes et avons découvert
les caractéristiques des différents modèles de dynamique des fluides. Enfin, ce projet nous a
permis de développer notre esprit critique et nos capacités d’analyse grâce à la partie
recherche documentaire du projet puisqu’il faut toujours faire comparer et vérifier les
différentes sources d’information utilisées pour se documenter, particulièrement sur internet.
De plus, ce projet nous a permis de travailler en groupe, alors qu’on ne se connaissait pas
tous avant de commencer le projet, de s’organiser et de se répartir les tâches. C’est
intéressant de pouvoir travailler avec des personnes qu’on ne connaît pas puisque cela nous
permet de faire de nouvelles rencontres, mais surtout cela nous prépare au monde de
l’entreprise plus tard où on ne pourra pas choisir les personnes avec qui on travaille, il faut
donc faire preuve de flexibilité et s’adapter aux différentes méthodes de travail et personnalités
des personnes qui composent notre groupe de projet. En effet, malgré les conditions
sanitaires et l’organisation particulière des cours en semi-présentiel, nous avons su composer
avec les disponibilités de chacun pour se retrouver régulièrement et ainsi progresser dans
l’avancement du projet.
De surcroît, ce projet nous a permis de mieux comprendre certains éléments abordés en cours
mais aussi certains phénomènes physiques, en aéronautique notamment, comme
l’importance que le dimensionnement d’une pièce peut avoir sur la résistance aux turbulences
de celle-ci.
Enfin, la découverte d’un logiciel complexe comme Ansys Fluent nous a permis de réaliser
que la modélisation en dynamique des fluides est un domaine très vaste dont nous avons
seulement étudié une petite partie. La modélisation d’écoulements turbulents est un domaine
concret que l’on retrouve dans de nombreux domaines industriels, il est donc intéressant de
commencer à manipuler ce genre de logiciels dès maintenant puisque nous serons sûrement
amenés à les réutiliser plus tard dans notre cursus à l’INSA ou même dans le monde du travail.
Pour conclure, ce projet nous a permis d’avoir un aperçu des notions étudiées dans les filières
Energétique et Propulsion ou encore Mécanique, la majorité des membres du groupe étant
intéressée par ces départements de spécialité. Enfin, ce projet en groupe est une expérience
enrichissante qui pourra être utile et constructive pour notre avenir en tant qu’ingénieur, tant
d’un point de vue pratique que théorique.
STPI/P6/2021 - 35
22
6. BIBLIOGRAPHIE
Images :
[1] Page de couverture : https://boowiki.info/images/b540_20/310px-Airplane_vortex.jpg
Vidéos :
[2] 5 cours magistraux réalisés par M. Safdari Shadloo.
Livres :
[3] Bruce R. Munson, Donald F. Young, Theodore H. Okiishi et Wade W. Huebsch,
« Fundamentals of Fluid Mechanics 6th edition », John Wiley & Sons, 2009
[4] Yunus A. Çengel et John M. Cimbala, « Fluid Mechanics Fundamentals and
Applications 3rd edition», McGraw-Hill, 2014
[5] Franck M. White, « Fluid Mechanics 7th edition », McGraw-Hill, 2011
Liens internet :
[6] Page Wikipédia sur le modèle k-𝜀 ∶ https://en.wikipedia.org/wiki/K-
epsilon_turbulence_model
[7] Page Wikipedia sur le modèle Spalart-Allmaras :
https://en.wikipedia.org/wiki/Spalart%E2%80%93Allmaras_turbulence_model
[8] Page Wikipédia sur la modélisation de la turbulence :
https://fr.wikipedia.org/wiki/Mod%C3%A9lisation_de_la_turbulence
[9] Page qui présente les principales propriétés des fluides :
https://lesfluides8sciences.weebly.com/les-proprieacuteteacutes-des-fluides.html
[10] Page Wikipédia sur les équations de Navier-Stokes :
https://fr.wikipedia.org/wiki/%C3%89quations_de_Navier-Stokes
[11] Page Wikipédia sur la moyenne de Reynolds :
https://fr.wikipedia.org/wiki/Moyenne_de_Reynolds
[12] Page Wikipédia sur les échelles de Kolmogorov :
https://en.wikipedia.org/wiki/Kolmogorov_microscales
[13] Page Wikipédia sur la viscosité : https://fr.wikipedia.org/wiki/Viscosit%C3%A9
[14] Comparaison des modèles de turbulence :
http://hmf.enseeiht.fr/travaux/bei/beiep/content/g24/2-comparaison-differents-
modeles-turbulence
[15] Page Wikipédia sur le modèle k-𝜔 :
https://en.wikipedia.org/wiki/K%E2%80%93omega_turbulence_model
[16] Page Wikipédia sur le modèle k-𝜔 SST :
https://en.wikipedia.org/wiki/Menter%27s_Shear_Stress_Transport
[17] Article sur les méthodes DNS, LES et RANS pour les transferts de chaleur
turbulents dans les couches limites :
https://www.sciencedirect.com/science/article/abs/pii/S0142727X13000714
[18] Page du site ideal simulations sur les modèles de turbulence CFD :
https://www.idealsimulations.com/resources/turbulence-models-in-cfd/
[19] Page Wikipédia sur la méthode LES :
https://en.wikipedia.org/wiki/Large_eddy_simulation
[20] Page du site XCV Wiki sur la méthode DNS :
https://fr.xcv.wiki/wiki/Direct_numerical_simulation
[21] Page du site Simscale sur les méthodes k-𝜔 et k-𝜔 SST :
https://www.simscale.com/docs/simulation-setup/global-settings/k-omega-
sst/#:~:text=%20K-Omega%20and%20K-
Omega%20SST%20%201%20Overview.,describes%20the%20size%20of%20large
%20energy-containing...%20More%20
STPI/P6/2021 - 35
23
Documents :
[22] Tutoriel d’utilisation d’Ansys Fluent pour la géométrie coudée :
http://users.abo.fi/rzevenho/ansys%20fluent%2018%20tutorial%20guide.pdf
[23] Thèse « Simulation des Grandes Echelles d’écoulements turbulents avec
transferts de chaleur », A. Chatelain, 2004 : https://tel.archives-ouvertes.fr/tel-
00006845v2/document
[24] Thèse « Direct Numerical Simulation and Wall-Resolved Large Eddy
Simulation in Nuclear Thermal Hydraulics », I. Tiselj, C. Flageul, J. Oder, 2019 :
https://hal.archives-ouvertes.fr/hal-02154872/document
[25] Thèse « Simulation des grandes échelles (SGE) d’écoulements de fluides
incompressibles », B. Brangeon, 2012 : https://hal.archives-ouvertes.fr/hal-
00682561/document
[26] Thèse sur le « Développement de modèles de turbulence adaptés à la
simulation des écoulements de convection naturelle à haut nombre de Rayleigh », D.
Vanpouille, 2013 : https://tel.archives-ouvertes.fr/tel-01021502/document
STPI/P6/2021 - 35
24
7. ANNEXE
7.1. Annexe 1
7.2. Annexe 2
Figure 13 : Schéma des différentes procédures de résolution numérique disponibles sur Fluent
STPI/P6/2021 - 35
25
7.3. Annexe 3
STPI/P6/2021 - 35
26
7.5. Annexe 4
k-𝜀
k-𝜔 standard
k-𝜔 SST
Spalart-Allmaras
STPI/P6/2021 - 35