Introduction à la CFD

Bouhelal Abdelhamid, Arezki Smaili

To cite this version:

Bouhelal Abdelhamid, Arezki Smaili. Introduction à la CFD. École d’ingénieur. Introduction to CFD
(Computational Fluid Dynamics), Algérie. 2022, pp.106. �hal-04427690�

HAL Id: hal-04427690

https://hal.science/hal-04427690
Submitted on 31 Jan 2024

HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est

archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents
entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non,
lished or not. The documents may come from émanant des établissements d’enseignement et de
teaching and research institutions in France or recherche français ou étrangers, des laboratoires
abroad, or from public or private research centers. publics ou privés.
Copyright
 Ecole Nationale Polytechnique (ENP)
Département de Génie Mécanique

Cours
Introduction à la CFD
(Computational Fluid Dynamics)

Dr. BOUHELAL Abdelhamid │ Pr. SMAILI Arezki

Ecole Nationale Polytechnique, Alger
 Qu’est ce que la CFD ?

Introduction à la CFD - ENP @2021-2022 2

 Qu’est ce que la CFD I.1 Introduction

▪ La CFD (Computational Fluid Dynamics) ou la dynamique des fluides numérique est

la science permettant de prédire les caractéristiques des écoulements des fluides, le
transfert thermique, le transfert de masse, les réactions chimiques …. en résolvant les
équations différentielles qui régissent ces processus à partir d’un codage
numérique.
▪ L'écoulement des fluides est couramment étudié de l'une des trois manières suivantes:
• Dynamique des fluides expérimentale.
• Dynamique des fluides théorique.
• Dynamique des fluides numérique (CFD).
▪ La CFD est principalement basée sur la résolution des
équations de Navier-Stokes numériquement à l’aide d’un
ordinateur.

Introduction à la CFD - ENP @2021-2022 3

 Pourquoi utiliser la CFD ?

Introduction à la CFD - ENP @2021-2022 4

 Pourquoi utiliser la CFD ? I.1 Introduction
▪ Avantages de l’utilisation la CFD par rapport l’expérimental:
• CFD fournit une base de données haute fidélité pour diagnostiquer le champ
d'écoulement.
• Simulation de phénomènes physiques difficiles à expérimenter (simulations grandeur
nature: e.g., navires et avions; Effets environnementaux: e.g., éoliennes, météo,… ;
phénomènes dangereux: e.g., combustion, explosion, etc.; ).
Mesure expérimentale Simulation CFD
Coûteuse Moins chère
Lente Rapide
Séquentielle Parallèle
Mono usage Multi-usages
▪ La CFD, en pratique, est utilisée pour:
• Études de conception et optimisation de (nouveaux) produits;
• Études des performances de systèmes/installations;
• Fournir des données complémentaires aux données théoriques et expérimentales.
Introduction à la CFD - ENP @2021-2022 5
 Applications de la CFD

Introduction à la CFD - ENP @2021-2022 6

 Applications de la CFD I.1 Introduction

▪ Quelques exemples d’applications:

• Aérodynamique;
• Hydrodynamique;
• Energies renouvelables;
• Turbomachines;
• Électroniques;
• Génie des procédés chimiques;
• Génie maritime;
• Combustion;
• Génie de l'environnement;
• Sports; Champ de vitesse autour d’une maquette de voiture.
• Génie biomédical, etc.

Introduction à la CFD - ENP @2021-2022 7

 Bref historique

Introduction à la CFD - ENP @2021-2022 8

 Bref historique I.1 Introduction

▪ Premiers travaux CFD en 1922.

▪ Richardson a résolu manuellement les équations de Laplace en
utilisant la méthode des différences finies: Application prévision
météorologique.
▪ Il a introduit le principe de calcul
parallèle: Richardson à proposer une
solution qu'il a appelée "l'usine de
prévisions".
o L'« usine » aurait rempli un vaste
stade de 64 000 personnes. Lewis Fry Richardson
o Chacun, armé d'une calculatrice (1881 - 1953)
mécanique, effectuerait une partie
du calcul.

Introduction à la CFD - ENP @2021-2022 9

 Bref historique I.1 Introduction

▪ A.Thom (1933): Première solution numérique pour

un écoulement autour un cylindre.
A.Thom (The Flow Past Circular Cylinders at Low Speeds),
Proc. Royal Society, A141, pp. 651-666, London, 1933.

▪ Au cours des années 1970, un groupe travaillant sous la direction

de Pr. Brian Spalding, (London), développe:
• Codes d’écoulements paraboliques (GENMIX).
• Codes basés sur la Vorticité-fonction de courant.
• Le principe de la méthode des volumes finis.
• L'algorithme SIMPLE (couplage vitesse-pression).
• La forme des équations 𝑘 − 𝜖 (utilisées aujourd’hui).
▪ En 1980, Suhas V. Patankar (encadré par Spalding) a publié
Numerical Heat Transfer and Fluid Flow, probablement le livre le Brian Spalding
plus influent sur la CFD à ce jour. (La base de la méthode des (1923 - 2016)
volumes finis).
Introduction à la CFD - ENP @2021-2022 10
 Bref historique I.1 Introduction

▪ Au cours des années 1980-2000, la CFD était utilisée principalement pour la recherche à
partir de codes “maison”: quand on voulait effectuer un calcul CFD, il fallait écrire un
programme.
▪ Il s'agit de la période au cours de laquelle les codes ou les logiciels CFD commerciaux se
sont développés dont la plupart sont disponibles aujourd'hui:
o Fluent (UK and US).
o CFX (UK and Canada).
o Polyflow (Belgium).
o Phoenix (UK).
o Star CCM+ (UK).
o Flow 3d (US).
o ESI/CFDRC (US).
o SCRYU (Japan).
Voir le site: www.cfdreview.com

Introduction à la CFD - ENP @2021-2022 11

 Bref historique I.1 Introduction
Logiciels CFD aujourd’hui
Commercial codes:
▪ Ansys-Fluent ( https://www.ansys.com/products/fluids/ansys-fluent );
▪ Ansys-CFX ( https://www.ansys.com/products/fluids/ansys-cfx );
▪ COMSOL Multiphysics ( https://www.comsol.com/ );
▪ EasyCFD (http://www.easycfd.net/);
▪ SimFlow ( https://sim-flow.com/ );
▪ CONVERGE ( https://convergecfd.com/ ).

Open source:
▪ OpenFlower ( http://openflower.sourceforge.net/ );
▪ OpenFoam (https://www.openfoam.com/ );
▪ Code_Saturne (https://www.code-saturne.org/cms/web/ );
▪ OpenLB (https://www.openlb.net/ );
▪ SLFCFD (http://slfcfd.sourceforge.net/ ).
Introduction à la CFD - ENP @2021-2022 12
 Bref historique I.1 Introduction

ANSYS INC:
o Création en 1970 à Canonsburg (États unis)
2006: acquisition de Fluent Inc.
o Fluent est un logiciel CFD avec interface graphique.
o Le plus célèbre au monde.

Introduction à la CFD - ENP @2021-2022 13

 Equations gouvernantes

Introduction à la CFD - ENP @2021-2022 14

 Equations gouvernantes I.1 Introduction
▪ Les équations gouvernantes le mouvement des fluides sont:
o L’équation de continuité qui traduit le principe de conservation de masse;
o Les équations de quantité de mouvement: la variation de quantité de mouvement est
égale à la somme des forces agissant sur le fluide (deuxième loi de Newton);
o L’équation d’énergie: l’énergie est conservée d’après le premier principe de la
thermodynamique.
▪ Pour les problèmes impliquant du transfert de chaleur, l’équation d’énergie nous permettra
d'obtenir la température du fluide.
▪ Si l’écoulement est, en plus, compressible, nous aurons également besoin d'une équation reliant
la masse volumique à la température (et à la pression). Cela se présente sous la forme d'une
équation d'état thermodynamique (e.g. équation des gaz parfaits).
▪ Dans certains cas, nous pouvons supposer que le fluide est isotherme (sa température est
constante), et donc l'équation de conservation de l'énergie ne sera pas nécessaire.
▪ Pour dériver les équations gouvernantes le mouvement de fluide (et le transfert de chaleur), il est
nécessaire de considérer le fluide comme un milieu continu (approche macroscopique: les
mouvements moléculaires sont ignorés).
▪ Une particule fluide est le plus petit élément possible de fluide dont les propriétés macroscopiques
ne sont pas influencées par des molécules individuelles (longueur macroscopique de 1 µm et plus).

Introduction à la CFD - ENP @2021-2022 15

 Equations gouvernantes I.1 Introduction
Comment peut-on décrire le mouvement des fluides ?
Description Lagrangienne Description Eulérienne
Suivre chaque particule individuellement: Suivre le mouvement des particules dans une
repère fixé à la particule fluide. région fixe: repère fixé à la région d’intérêt.

Champ de vitesse dans la région

Trajectoire d’une particule: d’étude (volume de contrôle):

Variation = dérivée 𝑫𝑽 𝛛𝑽
𝒅𝑽𝒑 (𝒕) totale (particulaire): 𝒂= = + 𝑽. 𝜵 . 𝑽
Variation = dérivée temporelle: 𝒂= 𝑫𝒕 𝝏𝒕
𝒅𝒕
dérivée locale dérivée convective
Introduction à la CFD - ENP @2021-2022 16
 Equations gouvernantes I.1 Introduction

Comment pouvons-nous dériver les équations régissant

l'écoulement du fluide ?

Technique du Bilan sur Volume de Contrôle: Appliquer le

bilan (de masse, de quantité de mouvement, d'énergie) sur
un volume de control infinitésimal (description Eulérienne).
Technique de Théorème de Transport de Reynolds (TTR):
Appliquer des lois de conservation à une parcelle de fluide
lagrangien. Ensuite, appliquer le TTR pour transformer ces
lois en une formulation eulérienne.

Introduction à la CFD - ENP @2021-2022 17

 Equations gouvernantes I.1 Introduction
Théorie du Transport de Reynolds (TTR)
▪ Soit 𝑭 une propriété extensive associée à la parcelle fluide:

𝒇 peut être une fonction

scalaire ou vectorielle

▪ Le théorème de transport de Reynolds relie la dérivée lagrangienne d'une propriété avec les
changements de cette propriété dans le volume de contrôle eulérien fixe comme suit:

Variation de 𝑭 dans le VC Débit net de 𝑭 au travers les SC

Introduction à la CFD - ENP @2021-2022 18

 Equations gouvernantes Equation de continuité
Equation de continuité: La masse est conservée dans un système fluide en mouvement.
▪ Le bilan de masse implique que le taux d'augmentation de la masse dans un élément fluide
(variation de masse pendant le temps) est égal la somme des masses qui entre et sort.
▪ On considère un élément de volume fluide dans un repère catésien: 𝒅ϑ = 𝒅𝒙 𝒅𝒚 𝒅𝒛
𝝏𝒎 𝝏𝝆
▪ Pendant le temps 𝒅𝒕, la variation de sa masse s’écrit: 𝐝𝒎 = 𝒅𝒕 = dϑ 𝒅𝒕
𝝏𝒕 𝝏𝒕
▪ Cette variation de masse est égale à la somme des masses qui entre et sort.
▪ Suivant l’axe 𝑥, le fluide entre avec la vitesse 𝒖𝒙 et sort avec la vitesse 𝒖𝒙+𝒅𝒙.
o La masse entrant pendant 𝒅𝒕 est: 𝝆𝒖 𝐝𝒚 𝒅𝒛 𝒅𝒕 𝒙
o La masse sortant pendant 𝒅𝒕 est: 𝝆𝒖 𝐝𝒚 𝒅𝒛 𝒅𝒕 𝒙+𝒅𝒙
 Le bilan de masse sur 𝑥 donne: 𝝆𝒖 𝒙 − 𝝆𝒖 𝒙+𝐝𝒙 𝐝𝒚 𝒅𝒛 𝒅𝒕
▪ Faisons un développement du 1er ordre: 𝝆𝒖 𝒙+𝒅𝒙 ≈ 𝝆𝒖 𝒙 + 𝝏 𝝆𝒖 /𝝏𝒙 𝒅𝒙
▪ Il reste alors sur l’axe 𝑥: − 𝝏 𝝆𝒖 /𝝏𝒙 𝒅ϑ 𝒅𝒕
▪ Par analogie, sur y: − 𝝏 𝝆𝒗 /𝝏𝒚 𝒅ϑ 𝒅𝒕 , et sur z: − 𝝏 𝝆𝒘 /𝝏𝒛 𝒅ϑ 𝒅𝒕
𝝏𝝆 𝝏 𝝆𝒖 𝝏 𝝆𝒗 𝝏 𝝆𝒘
▪ Au final: 𝐝𝒎 = 𝒅ϑ 𝐝𝒕 = − + + 𝒅ϑ 𝐝𝒕
𝝏𝒕 𝝏𝒙 𝝏𝒚 𝝏𝒛
𝝏𝝆 𝝏 𝝆𝒖 𝝏 𝝆𝒗 𝝏 𝝆𝒘
 + + + =𝟎
𝝏𝒕 𝝏𝒙 𝝏𝒚 𝝏𝒛 𝝏𝝆
Ou: +𝜵 𝝆𝑽 =𝟎 Equation de continuité
𝝏𝒕
Introduction à la CFD - ENP @2021-2022 19
 Equations gouvernantes Equation de continuité
Equation de Continuité: Dérivation à l'aide du théorème de transport de Reynolds
▪ D’après le TTR, la forme intégrale de l'équation de conservation de la masse s’écrit:
𝝏𝝆
න 𝒅ϑ + න 𝝆(𝑽 ⋅ 𝒏)𝒅𝑨 = 𝟎
𝝏𝒕
𝑽𝑪 𝑺𝑪

▪ En appliquant le théorème de divergence:

න 𝝆 𝑽 ⋅ 𝒏 𝒅𝑨 = න 𝜵 𝝆𝑽 𝒅ϑ
𝑺𝑪 𝑽𝑪

𝝏𝝆 𝛛𝝆
 න 𝒅ϑ + න 𝜵 𝛒𝑽 𝐝ϑ = 𝟎 ⇔ න + 𝜵 𝛒𝑽 𝒅ϑ = 𝟎
𝝏𝒕 𝝏𝒕
𝑽𝑪 𝐕𝐂 𝑽𝑪

▪ L'intégrale est satisfaite pour un volume arbitraire; elle doit être nulle:

𝝏𝝆
+𝜵 𝝆𝑽 =𝟎 On retrouve l’équation de continuité
𝝏𝒕
Introduction à la CFD - ENP @2021-2022 20
 Equations gouvernantes Equation de quantité de mouvement
Equation de quantité de mouvement: On applique le principe fondamental de la dynamique à
un élément de volume fluide:
𝑫𝑽
෍ 𝐝𝑭 = 𝐝𝑭𝒔 + 𝐝𝑭𝐛 = 𝝆𝒅ϑ
𝑫𝒕
▪ Les forces appliquées sont les forces de surface d𝐹Ԧ𝑠 et les forces de volume d𝐹Ԧb .
❑ Forces de volume 𝐝𝐅Ԧ𝐛 :
Forces gravitationnelles, électriques, magnétiques, flottaisons. La force de volume la plus courante
est la force due à l'accélération gravitationnelle 𝑔Ԧ : 𝐝𝑭𝒃 = 𝝆𝒅ϑ 𝒈
❑ Forces de surface 𝐝𝑭𝒔
Forces de pression (normales à la surface) et forces de de viscosité (frottements).
La force totale de surface (due à la pression et à la viscosité) est
exprimée comme: 𝒅𝑭𝒔 = 𝝈 ന 𝒏𝒅𝑨 = 𝜵ന 𝝈𝒅ϑ
𝝈𝒙𝒙 𝝈𝒙𝒚 𝝈𝒙𝒛
ന = 𝝈𝒚𝒙 𝝈𝒚𝒚 𝝈𝒚𝒛 est le tenseur des contraintes.
Où: 𝝈
𝝈𝒛𝒙 𝝈𝒛𝒚 𝝈𝒛𝒛
𝑫𝑽 𝝏𝑽
 𝝆 =𝝆 + 𝑽. 𝜵 𝑽 = 𝜵ന 𝝈 + 𝝆𝒈
𝑫𝒕 𝝏𝒕
Equation de Cauchy
Introduction à la CFD - ENP @2021-2022 21
 Equations gouvernantes Equation de quantité de mouvement
Dérivation à l'aide du théorème de transport de Reynolds
▪ D’après le TTR, la forme intégrale de l'équation de quantité de mouvement s’écrit:
𝛛 𝛒𝑽
න 𝒅ϑ + න 𝛒𝑽 𝑽 ⋅ 𝒏 𝐝𝐀 = 𝑭𝒔 + 𝑭𝒃
𝝏𝒕
𝑽𝑪 𝑺𝑪

❑ 𝑭𝒔 : forces de surface: 𝑭𝒔 = න 𝝈
ന 𝒏𝒅𝑨 ❑ 𝑭𝒃 : forces de volume: 𝑭𝒃 = න 𝝆𝒈𝒅ϑ
𝑺𝑪 𝐕𝑪
▪ On applique le théorème de divergence:

න 𝛒𝑽 𝑽 ⋅ 𝒏 𝐝𝐀 = න 𝜵 𝛒𝑽 𝑽 𝒅ϑ 𝑭𝒔 = න 𝝈 𝝈𝒅ϑ
ന 𝒏𝒅𝑨 = න 𝜵ന
𝐒𝐂 𝑽𝑪 𝑺𝑪 𝐕𝑪

𝛛 𝛒𝑽 Le volume de contrôle
 න + 𝜵 𝛒𝑽 𝑽 𝒅ϑ = න 𝜵ന
𝝈 + 𝝆𝒈 𝒅ϑ
𝝏𝒕 est arbitraire.
𝑽𝑪 𝐕𝑪

𝝏 𝝆𝑽 Equation de Cauchy
Donc: + 𝜵 𝝆𝑽 𝑽 = 𝜵ന
𝝈 + 𝝆𝒈 (Forme alternative)
𝝏𝒕
Introduction à la CFD - ENP @2021-2022 22
 Equations gouvernantes Equation de quantité de mouvement
Dérivation à l'aide du théorème de transport de Reynolds
▪ En développant, on obtient:
𝜕 𝜌𝑉
+ 𝛻 𝜌𝑉 𝑉 = 𝛻𝜎ധ + 𝜌𝑔Ԧ
𝜕𝑡
𝜕 𝜌𝑉 𝜕𝑉 𝜕𝜌
=𝜌 +𝑉 Et 𝛻 𝜌𝑉 𝑉 = 𝜌𝛻 𝑉𝑉 + 𝑉𝛻 𝜌𝑉
𝜕𝑡 𝜕𝑡 𝜕𝑡

Avec: 𝜌𝛻 𝑉𝑉 = 𝜌 𝑉𝛻 𝑉  𝛻 𝜌𝑉 𝑉 = 𝜌 𝑉𝛻 𝑉 + 𝑉𝛻 𝜌𝑉

𝜕𝑉 𝜕𝜌
 𝜌 +𝑉 + 𝜌 𝑉𝛻 𝑉 + 𝑉𝛻 𝜌𝑉 = 𝛻𝜎ധ + 𝜌𝑔Ԧ
𝜕𝑡 𝜕𝑡
𝜕𝑉 𝜕𝜌
 𝜌 +𝑉 + 𝛻 𝜌𝑉 + 𝜌 𝑉𝛻 𝑉 = 𝛻𝜎ധ + 𝜌𝑔Ԧ
𝜕𝑡 𝜕𝑡

Donc:
𝝏𝑽 On retrouve l’équation
𝝆 + 𝑽𝜵 𝑽 = 𝜵ന
𝝈 + 𝝆𝒈
𝝏𝒕 de Cauchy
Introduction à la CFD - ENP @2021-2022 23
 Equations gouvernantes Equation de quantité de mouvement
▪ Contrairement à l’équation de continuité, l’équation de Cauchy est une équation vectorielle.
▪ L’équation de Cauchy représente ainsi trois équations scalaires, une pour chaque axe de
coordonnées dans les problèmes 3D.
▪ Pour un repère Cartésien, les équations de Cauchy s'écrivent:

𝐷𝑉 𝜕𝑉
𝜌 =𝜌 + 𝑉𝛻 𝑉 = 𝛻𝜎ധ + 𝜌𝑔Ԧ
𝐷𝑡 𝜕𝑡
𝐷𝑢 𝜕𝑢 𝜕𝑢 𝜕𝑢 𝜕𝑢 𝜕𝜎xx 𝜕𝜎yx 𝜕𝜎𝑧x
𝜌 =𝜌 +𝑢 +𝑣 +𝑤 = + + + 𝜌𝑔𝑥
𝐷𝑡 𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑥 𝜕𝑦 𝜕𝑧

𝐷𝑣 𝜕𝑣 𝜕𝑣 𝜕𝑣 𝜕𝑣 𝜕𝜎x𝑦 𝜕𝜎y𝑦 𝜕𝜎𝑧𝑦

𝜌 =𝜌 +𝑢 +𝑣 +𝑤 = + + + 𝜌𝑔y
𝐷𝑡 𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑥 𝜕𝑦 𝜕𝑧

𝐷𝑤 𝜕𝑤 𝜕𝑤 𝜕𝑤 𝜕𝑤 𝜕𝜎x𝑧 𝜕𝜎y𝑧 𝜕𝜎𝑧𝑧

𝜌 =𝜌 +𝑢 +𝑣 +𝑤 = + + + 𝜌𝑔z
𝐷𝑡 𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑥 𝜕𝑦 𝜕𝑧

On ne peut pas résoudre aucun problème de mécanique des fluides en utilisant l'équation de Cauchy !

Introduction à la CFD - ENP @2021-2022 24

 Equations gouvernantes Equation
Équations
de quantité
de quantité
de mouvement
de mouvement
▪ Le tenseur des contraintes 𝝈ന représente les contraintes dues à la pression et les contraintes dues à
la viscosité.
▪ La pression agit toujours vers l'intérieur, et normale à la surface, et les contraintes visqueuses
agissent vers l’extérieur normales et tangentielles à la surface.
▪ On peut, donc, séparer la pression des contraintes dues à la viscosité:
𝜎xx 𝜎x𝑦 𝜎x𝑧 −𝑝 0 0 𝜏xx 𝜏x𝑦 𝜏x𝑧
𝜎ധ = 𝜎yx 𝜎yy 𝜎yz = 0 −𝑝 0 + 𝜏yx 𝜏yy 𝜏yz
𝜎zx 𝜎zy 𝜎zz 0 0 −𝑝 𝜏zx 𝜏zy 𝜏zz

Tenseur de contraintes visqueuses

▪ 𝝈𝒅ϑ = 𝜵 −𝒑ധ𝑰 + 𝝉ധ 𝒅𝐯 = −𝜵𝒑 + 𝜵ധ𝝉 𝒅ϑ
Les forces de surface peuvent être exprimées par: 𝒅𝑭𝒔 = 𝜵ന
▪ L’équation générale du mouvement peut être reformuler comme:
𝑫𝑽 𝝏𝑽
𝝆 =𝝆 + 𝑽. 𝜵 𝑽 = −𝜵𝒑 + 𝜵ധ𝝉 + 𝝆𝒈
𝑫𝒕 𝝏𝒕
Cas particuliers: (i) Fluide au repos, 𝑽 =0: 𝜵𝒑 = 𝝆𝒈 (Equation Fondamentale de la Statique)
𝑫𝑽
(ii) Fluide non-visqueux (parfait), 𝝉ധ =0: 𝝆 = −𝜵𝒑 + 𝝆𝒈 (Equation d’Euler)
𝑫𝒕
Introduction à la CFD - ENP @2021-2022 25
 Equations gouvernantes Equation de quantité de mouvement
Fluides Newtoniens et les équations de Navier-Stokes
▪ Des observations expérimentales ont montré que pour les fluides Newtoniens, il existe une
relation linéaire entre les contraintes visqueuses et le taux de déformations pures.
▪ Navier et Stokes ont proposé un modèle de tenseur des contraintes visqueuses pour les fluides
Newtoniens.
▪ Les hypothèses fondamentales proposées par Navier et Stokes lors de la dérivation des
composantes du tenseur des contraintes (𝜏)Ӗ sont:
o La rotation pure d'un élément fluide n'engendre aucune contrainte visqueuse.
o L'écoulement du fluide est isotrope, ce qui signifie que ses propriétés sont les mêmes
lorsqu'elles sont mesurées dans différentes directions.
o Le tenseur des contraintes visqueuses (𝜏)Ӗ et le tenseur des taux de déformations pures
𝜖 Ӗ sont liés par l’expression suivante (relation de Navier et Stokes):

𝝉ധ = 𝟐𝝁ധ𝝐 + 𝝀 𝜵𝑽ധ𝑰

Où 𝝀 = −𝟐/𝟑 𝝁 (hypothèse de Stokes).

Introduction à la CFD - ENP @2021-2022 26

 Equations gouvernantes Equation de quantité de mouvement
Fluides Newtoniens et les équations de Navier-Stokes
Rappel: La déformation totale du fluide en mouvement est caractérisé par le tenseur des taux
de déformations (𝑫
ന ). Ce tenseur se décompose en deux tenseurs:

𝝏𝒖 𝝏𝒖 𝝏𝒖 𝝏𝒖 𝟏 𝝏𝒖 𝝏𝒗 𝟏 𝝏𝒖 𝝏𝒘
+ + Tenseur
𝝏𝒙 𝝏𝒚 𝝏𝒛 𝝏𝒙 𝟐 𝝏𝒚 𝝏𝒙 𝟐 𝝏𝒛 𝝏𝒙
symétrique
𝝏𝒗 𝝏𝒗 𝝏𝒗 𝟏 𝝏𝒗 𝝏𝒖 𝝏𝒗 𝟏 𝝏𝒗 𝝏𝒘
ന=
𝑫 = + + 𝝐ധ
𝝏𝒙 𝝏𝒚 𝝏𝒛 𝟐 𝝏𝒙 𝝏𝒚 𝝏𝒚 𝟐 𝝏𝒛 𝝏𝒚
𝝏𝒘 𝝏𝒘 𝝏𝒘 𝟏 𝝏𝒘 𝝏𝒖 𝟏 𝝏𝒘 𝝏𝒗 𝝏𝒘
𝝏𝒙 𝝏𝒚 𝝏𝒛 + +
𝟐 𝝏𝒙 𝝏𝒛 𝟐 𝝏𝒚 𝝏𝒛 𝝏𝒛
𝟏 𝝏𝒖 𝝏𝒗 𝟏 𝝏𝒖 𝝏𝒘
𝟎 −
𝟐 𝝏𝒚 𝝏𝒙
−
𝟐 𝝏𝒛 𝝏𝒙 Tenseur
𝟏 𝝏𝒖 𝝏𝒗 𝟏 𝝏𝒗 𝝏𝒘 antisymétrique
+ − + 𝟎 − ന
𝝎
𝟐 𝝏𝒚 𝝏𝒙 𝟐 𝝏𝒛 𝝏𝒚
𝟏 𝝏𝒖 𝝏𝒘 𝟏 𝝏𝒗 𝝏𝒘
− + − + 𝟎
𝟐 𝝏𝒛 𝝏𝒙 𝟐 𝝏𝒛 𝝏𝒚
𝝐ധ ∶ Tenseur des taux de déformations pures (déformation angulaire et volumétrique).
𝝎ന : Tenseur des taux de rotations pures.
Introduction à la CFD - ENP @2021-2022 27
 Equations gouvernantes Equation de quantité de mouvement
Fluides Newtoniens et les équations de Navier-Stokes
▪ En introduisant le modèle de Navier-Stokes dans l’équation de qté de mvt, on obtient:
𝑫𝑽 𝝏𝑽 Modèle de Navier-Stokes:
qté de mvt: 𝝆 =𝝆 + 𝑽𝜵 𝑽 = −𝜵𝒑 + 𝜵ധ𝝉 + 𝝆𝒈
𝑫𝒕 𝝏𝒕 𝟐
𝝉ധ = 𝟐𝝁ധ𝝐 − 𝝁 𝜵𝑽ധ𝑰
𝟑
𝑫𝑽 𝝏𝑽 𝟐
 𝝆 =𝝆 + 𝑽𝜵 𝑽 = −𝜵𝒑 + 𝜵 𝟐𝝁ധ𝝐 − 𝝁 𝜵𝑽ധ𝑰 + 𝝆𝒈
𝑫𝒕 𝝏𝒕 𝟑

 𝝆 𝑫𝑽 = 𝝆 𝝏𝑽 + 𝑽𝜵 𝑽 = −𝜵𝒑 + 𝟐𝝁𝜵ധ 𝟐
𝝐 − 𝝁𝜵 𝜵𝑽 + 𝝆𝒈
𝑫𝒕 𝝏𝒕 𝟑
▪ Explicitons le terme 2𝜇𝛻𝜖 Ӗ (on considère la projection suivant l’axe x):
𝜕 𝜕𝑢 1 𝜕 𝜕𝑢 𝜕𝑣 1 𝜕 𝜕𝑢 𝜕𝑤 𝜕2𝑢 𝜕2𝑢 𝜕2𝑣 𝜕2𝑢 𝜕2𝑤
𝟐𝝁𝜵ധ𝝐 = 2𝜇 + + + + =𝜇 2 2+ + + +
𝐱 𝜕𝑥 𝜕𝑥 2 𝜕𝑦 𝜕𝑦 𝜕𝑥 2 𝜕𝑧 𝜕𝑧 𝜕𝑥 𝜕𝑥 𝜕𝑦 2 𝜕𝑦𝜕𝑥 𝜕𝑧 2 𝜕𝑧𝜕𝑥
𝜕2𝑢 𝜕2𝑢 𝜕2𝑣 𝜕2𝑢 𝜕2𝑤 𝜕 2 𝑢 𝜕 2 𝑢 𝜕 2 𝑢 𝜕 𝜕𝑢 𝜕𝑣 𝜕𝑤
=𝜇 2 2+ + + + =𝜇 + + + + +
𝜕𝑥 𝜕𝑦 2 𝜕𝑥𝜕𝑦 𝜕𝑧 2 𝜕𝑥𝜕𝑧 𝜕𝑥 2 𝜕𝑦 2 𝜕𝑧 2 𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑧
𝜕2𝑢 𝜕2𝑢 𝜕2𝑢 𝝏
=𝜇 2
+ 2
+ 2
= 𝝁𝜟𝒖 + 𝝁 𝜵𝑽  𝟐𝝁𝜵ധ𝝐 = 𝝁𝜟𝑽 + 𝝁𝜵 𝜵𝑽
𝜕𝑥 𝜕𝑦 𝜕𝑧 𝝏𝒙
Introduction à la CFD - ENP @2021-2022 28
 Equations gouvernantes Equation de quantité de mouvement
Fluides Newtoniens et les équations de Navier-Stokes
L'équation de quantité de mouvement devient donc:
𝑫𝑽 𝝏𝑽 𝟐
𝝆 =𝝆 + 𝑽𝜵 𝑽 = −𝜵𝒑 + 𝟐𝝁𝜵ധ𝝐 − 𝝁𝜵 𝜵𝑽 + 𝝆𝒈 Avec: 𝟐𝝁𝜵ധ𝝐 = 𝝁𝜟𝑽 + 𝝁𝜵 𝜵𝑽
𝑫𝒕 𝝏𝒕 𝟑
𝑫𝑽 𝝏𝑽 𝟐
 𝝆 =𝝆 + 𝑽𝜵 𝑽 = −𝜵𝒑 + 𝝁𝜟𝑽 + 𝝁𝜵 𝜵𝑽 − 𝝁𝜵 𝜵𝑽 + 𝝆𝒈
𝑫𝒕 𝝏𝒕 𝟑

𝑫𝑽 𝝏𝑽 𝟏 Equation de Navier-Stokes
𝝆 =𝝆 + 𝑽𝜵 𝑽 = −𝜵𝒑 + 𝝁𝜟𝑽 + 𝝁𝜵 𝜵𝑽 + 𝝆𝒈 (cas compressible)
𝑫𝒕 𝝏𝒕 𝟑

Cas particulier: écoulement incompressible:

Pour un écoulement incompressible (𝛻𝑉 = 0), l’équation de Navier-Stokes devient:

𝑫𝑽 𝝏𝑽 Equation de Navier-Stokes
𝝆 =𝝆 + 𝑽𝜵 𝑽 = −𝜵𝒑 + 𝝁𝜟𝑽 + 𝝆𝒈 (cas incompressible)
𝑫𝒕 𝝏𝒕

Introduction à la CFD - ENP @2021-2022 29

 Equations gouvernantes Equation d’énergie
Le changement total d'énergie dans le système est
Equation d’énergie: égal à la différence entre la chaleur transférée au
1er principe de la thermodynamique: système et le travail effectué par le système sur
son environnement. (si 𝑊ሶ est effectué sur le
système  (dE/dt = Qሶ + 𝑊) ሶ
L'énergie par unité de masse contenue dans un système se compose:

Avec:

▪ Pour fermé l'équation de l'énergie, nous avons besoin:

(i) d’une équation d’état: (e.g. gaz parfaits: )
(ii) des relations thermodynamiques reliant les variables énergétiques aux propriétés de base
telles que la température et la pression.

Gaz parfaits:

Introduction à la CFD - ENP @2021-2022 30

 Equations gouvernantes Equation d’énergie
▪ Le principal mécanisme de transfert de chaleur 𝑸ሶ vers le système fluide est la conduction
thermique. Puisque la chaleur est transférée vers/depuis le système à travers la surface limite,
nous pouvons la représenter comme une intégrale de surface du flux de chaleur (W/m2):

Loi de Fourier

▪ En thermodynamique, il existe deux principaux mécanismes de

travail effectués par un système sur son environnement: le
travail de pression et le travail de cisaillement visqueux.

Travail de pression:

Travail des contraintes visqueuses

(ou dissipassions visqueuses) :

TTR:

Introduction à la CFD - ENP @2021-2022 31

 Equations gouvernantes Equation d’énergie

▪ On peut placer le terme de travail de pression à l'intérieur du terme de convection et

introduire la loi de Fourier et simplifier le membre de droite pour obtenir:

▪ On applique le théorème de divergence aux intégrales de surface, on obtient:

Enfin, on obtient la forme différentielle de

l'équation de conservation de l'énergie:
Introduction à la CFD - ENP @2021-2022 32
 Equations gouvernantes Equation d’énergie
▪ En utilisant le modèle de Navier-Stokes des contraintes visqueuses, l'équation de l'énergie
devient:

▪ L’équation d’énergie, après le réarrangement, s’écrit:

𝑫𝑬 𝝏𝑬
𝝆 =𝝆 + 𝑽𝜵 𝑬 = 𝜵 𝐤𝜵𝐓 − 𝒑𝜵𝑽 + 𝚽
𝑫𝒕 𝝏𝒕

𝟐 𝟐 𝟐 𝟐 𝟐 𝟐
𝝏𝒖 𝝏𝒗 𝝏𝒘 𝝏𝒖 𝝏𝒗 𝝏𝒖 𝝏𝒘 𝝏𝒗 𝝏𝒘 𝟐
Avec: 𝚽=𝝁 𝟐 + + + + + + + + − 𝝁 𝜵𝑽
𝟐
𝝏𝒙 𝝏𝒚 𝝏𝒛 𝝏𝒚 𝝏𝒙 𝝏𝒛 𝝏𝒙 𝝏𝒛 𝝏𝒚 𝟑

Introduction à la CFD - ENP @2021-2022 33

 Résumé
Equations gouvernantes –
équations CFD

Introduction à la CFD - ENP @2021-2022 34

 Equations gouvernantes FluidesI.2 Équation
quantitéde
Newtoniens
Équations de decontinuité
compressibles
mouvement
𝝏𝝆
▪ Equation de continuité (Conservation de masse): +𝜵 𝝆𝑽 =𝟎
𝝏𝒕

𝑫𝑽 𝝏𝑽 𝟏
▪ Equations de quantité du mouvement: 𝝆 =𝝆 + 𝑽𝜵 𝑽 = −𝜵𝒑 + 𝝁𝜟𝑽 + 𝝁𝜵 𝜵𝑽 + 𝝆𝒇
𝑫𝒕 𝝏𝒕 𝟑

𝑫𝑬 𝝏𝑬
▪ Equation d’énergie: 𝝆 =𝝆 + 𝑽𝜵 𝑬 = 𝜵 𝐤𝜵𝐓 − 𝒑𝜵𝑽 + 𝚽
𝑫𝒕 𝝏𝒕

𝟐 𝟐 𝟐 𝟐 𝟐 𝟐
𝝏𝒖 𝝏𝒗 𝝏𝒘 𝝏𝒖 𝝏𝒗 𝝏𝒖 𝝏𝒘 𝝏𝒗 𝝏𝒘 𝟐 𝟐
Avec: 𝚽 = 𝝁 𝟐 + + + + + + + + − 𝝁 𝜵𝑽
𝝏𝒙 𝝏𝒚 𝝏𝒛 𝝏𝒚 𝝏𝒙 𝝏𝒛 𝝏𝒙 𝝏𝒛 𝝏𝒚 𝟑

Energie interne: 𝑬 = 𝑪𝒗 𝑻 Enthalpie: 𝑯 = 𝑪𝒑 𝑻

𝒑
▪ Equation d’état: 𝑭 𝝆, 𝑻, 𝒑 = 𝟎 Gaz parfait: = 𝒓𝑻
𝝆

Introduction à la CFD - ENP @2021-2022 35

 Equations gouvernantes
▪ Résumé: équations de Navier-Stokes incompressible (dans un repère cartésien):

𝛛𝒖 𝛛𝒗 𝛛𝒘
+ + =𝟎
𝛛𝒙 𝛛𝒚 𝛛𝒛
𝝏𝒖 𝝏𝒖 𝝏𝒖 𝝏𝒖 𝝏𝒑 𝝏𝟐 𝒖 𝝏𝟐 𝒖 𝝏𝟐 𝒖
𝝆 +𝒖 +𝒗 +𝒘 =− +𝝁 𝟐
+ 𝟐 + 𝟐 + 𝝆𝒈𝒙
𝝏𝒕 𝝏𝒙 𝝏𝒚 𝝏𝒛 𝝏𝒙 𝝏𝒙 𝝏𝒚 𝝏𝒛
𝝏𝒗 𝝏𝒗 𝝏𝒗 𝝏𝒗 𝝏𝒑 𝝏𝟐 𝒗 𝝏𝟐 𝒗 𝝏𝟐 𝒗
𝝆 +𝒖 +𝒗 +𝒘 =− +𝝁 + + + 𝝆𝒈𝒚
𝝏𝒕 𝝏𝒙 𝝏𝒚 𝝏𝒛 𝝏𝒚 𝝏𝒙𝟐 𝝏𝒚𝟐 𝝏𝒛𝟐
𝝏𝒘 𝝏𝒘 𝝏𝒘 𝝏𝒘 𝝏𝒑 𝝏𝟐 𝒘 𝝏𝟐 𝒘 𝝏𝟐 𝒘
𝝆 +𝒖 +𝒗 +𝒘 =− +𝝁 + + + 𝝆𝒈𝒛
𝝏𝒕 𝝏𝒙 𝝏𝒚 𝝏𝒛 𝝏𝒛 𝝏𝒙𝟐 𝝏𝒚𝟐 𝝏𝒛𝟐
𝝏𝒖 𝝏𝑻 𝝏𝑻 𝝏𝑻 𝝏 𝝏𝑻 𝝏 𝝏𝑻 𝝏 𝝏𝑻
𝝆𝑪𝒑 +𝒖 +𝒗 +𝒘 = 𝐤 + 𝐤 + 𝐤 +𝚽
𝝏𝒕 𝝏𝒙 𝝏𝒚 𝝏𝒛 𝝏𝒙 𝝏𝒙 𝝏𝒚 𝝏𝒚 𝝏𝒛 𝝏𝒛

Introduction à la CFD - ENP @2021-2022 36

 Equations gouvernantes Equation générale de transport
Equation générale de transport (Equations CFD)

▪ En utilisant une variable générale 𝝓, la forme conservatrice de toutes les équations CFD peut
utilement être écrite sous la forme suivante:

𝝏(𝝆𝝓)
+ 𝜵. (𝝆𝑽𝝓) = 𝜵 ⋅ (𝜞 𝜵𝝓) + 𝑺𝝓
𝝏𝒕
Taux d’augmentation Taux net de 𝝓 Taux Taux d’augmentation
de 𝝓 dans l’élément en sortie de d’augmentation de 𝝓 du à la source
l’élément fluide de 𝝓 du à la
fluide
(convection)
diffusion

𝝏(𝝆𝝓)
+ 𝒅𝒊𝒗(𝝆𝑽𝝓) = 𝒅𝒊𝒗(𝜞 𝒈𝒓𝒂𝒅 𝝓) + 𝑺𝝓
𝝏𝒕

Introduction à la CFD - ENP @2021-2022 37

 Equations gouvernantes Equation générale de transport
Equation générale de transport (Equations CFD)
𝝏(𝝆𝝓)
+ 𝜵. (𝝆𝑽𝝓) = 𝜵 ⋅ (𝜞 𝜵𝝓) + 𝑺𝝓
𝝏𝒕

Continuity 𝝓=𝟏 𝜞=𝟎 𝑺𝝓 = 𝟎

Flow
Eqts
Momentum 𝝓=𝑽 𝜞=𝝁 𝑺𝝓 = −𝛁𝒑 + 𝝆𝒈

Energy 𝝓 = 𝑪𝒑 𝑻 𝜞 = 𝒌/𝑪𝒑 𝑺𝝓 = 𝚽 + 𝒒ሶ

+ Turbulence Concentration Combustion Other eqts

Introduction à la CFD - ENP @2021-2022 38

 CFD – Comment cela fonctionne

Introduction à la CFD - ENP @2021-2022 39

 CFD – Comment cela fonctionne I.1 Introduction
Étapes de la CFD
Pre-processing Simulation Post-processing

Géométrie (CAO) Maillage

ANSYS ICEM CFD

Codes of CFD : TECPLOT
GRIDGEN
FLUENT MATLAB
GMSH
CFX FIELDVIEW
SOLIDWORKS STAR CCM+
CATIA Paraview
CODES MAISONS
PRO-ENGINEER
OpenFOAM
UNIGRAPHICS
FreeCAD

Introduction à la CFD - ENP @2021-2022 40

 CFD – Comment cela fonctionne I.1 Introduction

• L'analyse commence par un modèle

mathématique d'un problème physique. Goulot de
remplissage
• Lois de conservation de la masse, de la
quantité de mouvement et de l’énergie Bouteille
doivent être satisfaites dans toute la région
d'intérêt.
• Propriétés des fluides sont modélisées
empiriquement.
• Hypothèses simplificatrices sont posées pour
que le problème posé soit facile à résoudre
(e.g., état stationnaire, incompressible, non
visqueux, 2D…)
• Définir des conditions initiales et aux limites
Exemple: Domaine pour le
remplissage de la bouteille
appropriées pour le problème.

Introduction à la CFD - ENP @2021-2022 41

 CFD – Comment cela fonctionne I.1 Introduction

Discrétisation: La CFD applique des méthodes numériques pour

développer des approximations des équations dans la région de
fluide d'intérêt.
Les équations différentielles sont transformées en un système
d’équations algébriques.
• On définit un quadrillage formé de cellules: Le maillage.
• Le système d’équations algébriques est résolu
numériquement (sur un ordinateur) pour les toutes
variables (V, P, T…) à chaque nœud ou cellule.
• Système d'équations sont résolues simultanément afin de
fournir la solution.
Maillage pour le
remplissage de la bouteille

Introduction à la CFD - ENP @2021-2022 42

 CFD – Comment cela fonctionne I.1 Introduction
Maillage
Types de maillage
Pour les géométries simples,
les maillages quad/hexa
peuvent fournir des solutions
triangle quadrilateral
de haute qualité avec moins de
cellules qu'un maillage tri/tetra.

tetrahedron hexahedron

Pour les géométries complexes,

les maillages quad/hexa sont
difficiles, on utilise un maillage
pyramid prism/wedge
tri/tetra ou un maillage par blocks.

On doit raffiner la zone proche paroi !

Introduction à la CFD - ENP @2021-2022 43

 CFD – Comment cela fonctionne I.1 Introduction
Maillage

Maillage structuré: Maillage non structuré:

Formé par des lignes continues. Formé par des lignes discontinues.
Maillage idéal.

Introduction à la CFD - ENP @2021-2022 44

 CFD – Comment cela fonctionne I.1 Introduction
Maillage
Exemples des maillages structurés: Exemples des maillages non structurés:

Introduction à la CFD - ENP @2021-2022 45

 CFD – Comment cela fonctionne I.1 Introduction
Méthodes numériques utilisées en logiciels CFD
• La méthode des volumes finis (~80%).
• La méthode des éléments finis (~15%).
• La méthode des différences finies (~5%).

▪ Méthode des différences finies

• Historiquement, la plus ancienne des
trois. (a été utilisée par Richardson en
1910).
Avantage: facile à mettre en œuvre.
Inconvénients: limité aux maillages
simples, et ne conserve pas la quantité Principe: Basée sur l’approximation des
de mouvement, l’énergie et la masse sur dérivées: Les dérivées première et seconde
des maillages grossiers. sont approchées par des développements en
série de Taylor. Le système algébrique
résultant est résolu numériquement.
Introduction à la CFD - ENP @2021-2022 46
 CFD – Comment cela fonctionne I.1 Introduction
Méthodes numériques utilisées en logiciels CFD
• La méthode des volumes finis (~80%).
• La méthode des éléments finis (~15%).
• La méthode des différences finies (~5%).
▪ Méthode des éléments finis
• La première utilisation: Courant (1943) pour
résoudre un problème de torsion.
• Clough (1960): a donné à la méthode son
nom.
Avantage: précision maximale sur les maillages
grossiers. Excellente pour les problèmes
dominés par la diffusion: e.g., écoulements
visqueux à surface libre. Principe: Basée sur l’approximation
Inconvénients: lent pour les problèmes de des solutions: Chercher des solutions
grandes dimensions et peu adapté aux approchées sur des domaines décrets.
écoulements turbulents.
Introduction à la CFD - ENP @2021-2022 47
 CFD – Comment cela fonctionne I.1 Introduction
Méthodes numériques utilisées en logiciels CFD
• La méthode des volumes finis (~80%).
• La méthode des éléments finis (~15%). node
• La méthode des différences finies (~5%).
edge
▪ Méthode des volumes finis face
cell

• La première utilisation bien documentée a été cell

node

faite par Evans et Harlow (1957) center

• Patankar et Spalding (1972-1980): applications face
pour les problèmes des écoulements.
Avantage: masse, quantité de mouvement, énergie cell
conservée même sur des maillages grossiers. Principe: Intégration des équations
Utilisation pratique pour des maillages non sur des VCs. Intégration la forme
structurés; conservative des équations sur chaque
Inconvénients: mal-estimation de la diffusion volume de contrôle pour obtenir les
lorsque des schémas simples sont utilisés. équations discrétisées qui conservent
toutes les grandeurs physiques.
Introduction à la CFD - ENP @2021-2022 48
 Introduction à la modélisation
de la turbulence

Introduction à la CFD - ENP @2021-2022 49

 Introduction I.1 Introduction
La plupart des écoulements industriels et naturels sont turbulents !

Introduction à la CFD - ENP @2021-2022 50

 Introduction I.1 Introduction
Écoulements laminaires vs écoulements turbulents

Ecoulement laminaire Ecoulement turbulent

• Mouvement ordonné et régulier. • Mouvement irrégulier et aléatoire.
• Stable: absence de fluctuations. • Instable: présence de fluctuations.
• Une seule échelle de longueur. • Rotationnels (tourbillons non nuls).
• Faibles nombres de Reynolds. • Large gamme d’échelles de longueur.
• Dissipation d’énergie.
La plupart des écoulements en pratique • Diffusivité: grande capacité de mélange.
sont turbulents ! • Très sensible aux conditions initiales.
• Nombres de Reynolds élevés.
Introduction à la CFD - ENP @2021-2022 51
 Introduction I.1 Introduction
Transition du régime laminaire au régime turbulent
Nombre de Reynolds
𝜌𝑈𝐿∗
𝜕𝑢𝑖 𝜕𝑢𝑖 𝜕𝑝 𝜕 𝜕𝑢𝑖 𝑅𝑒𝐿∗ =
𝜌 + 𝑢𝑗 =− + 𝜇 𝜇
𝜕𝑡 𝜕𝑥𝑗 𝜕𝑥𝑖 𝜕𝑥𝑗 𝜕𝑥𝑗 𝐿∗ = 𝑥, 𝑑, 𝑑ℎ , 𝑒𝑡𝑐.
Forces d′inertie Forces de viscosité
𝜕𝑢𝑖
𝜌𝑢𝑗
Forces d′ inertie 𝜕𝑥𝑗 𝜌𝑈 2 /𝐿∗ 𝜌𝑈𝐿∗
𝑅𝑒𝐿∗ = = ~ 2 =
Forces de viscosité 𝜕 𝜕𝑢𝑖 𝜇𝑈/𝐿 ∗ 𝜇
𝜇
𝜕𝑥𝑗 𝜕𝑥𝑗
Écoulements externes
𝑅𝑒𝑥 ≥ 5 × 105 le long d'une surface D'autres facteurs tels que la rugosité de
𝑅𝑒d ≥ 2 × 104 autour d'un obstacle surface, le soufflage, l'aspiration et d'autres
sources de perturbation peuvent provoquer
Écoulements internes une transition vers la turbulence à des
𝑅𝑒dh ≥ 2300 nombres de Reynolds inférieurs.
Introduction à la CFD - ENP @2021-2022 52
 Description physique de la
turbulence

Introduction à la CFD - ENP @2021-2022 53

 Description physique de la turbulence I.1 Introduction
Cascade d’énergie – Théorème de Richardson (1922)
• L'écoulement turbulent est constitué de tourbillons de différentes tailles.
• Les grands tourbillons extraient l'énergie de l'écoulement moyen (production d'énergie).
• L'énergie est transférée des tourbillons plus grands aux tourbillons plus petits.
• Les tourbillons plus grands contiennent la majeure partie de l'énergie.
• Dans les plus petits tourbillons, l'énergie turbulente est convertie en énergie interne par
dissipation visqueuse (dissipation d'énergie).

Production
l

Dissipation

Introduction à la CFD - ENP @2021-2022 54

 Description physique de la turbulence I.1 Introduction
Echelles intégrales (grandes échelles de turbulence)
Considérons un tourbillon de la plus grande taille caractérisé par: 𝑢l
• l : longueur caractéristique. l
• 𝑢l : vitesse caractéristique.
• 𝜏l : temps caractéristique (durée de vie).i.e. (𝜏l ≡ l/𝑢l ).
Soit un tourbillon de la plus petite taille ayant:
• 𝜂 : longueur caractéristique.
• 𝑢𝜂 : vitesse caractéristique. 𝑢η 𝜂
• 𝜏η : temps caractéristique (durée de vie).i.e. (𝜏𝜂 ≡ 𝜂 /𝑢𝜂 ).
Les plus grands tourbillons ont une énergie cinétique d'ordre 𝑘 ~ 𝑢l2 . Il s’en suit:

l Longueur intégrale (de grands tourbillons).

1 ′
Echelles intégrales 𝑘 = 𝑢𝑖 𝑢𝑖 ′
ne dépendent que 𝑢l = 𝑘1/2 Échelle intégrale de vitesse. 2
de 𝑘 et de l. Énergie cinétique
𝜏l = l.𝑘 −1/2 Échelle intégrale de temps. turbulente
Introduction à la CFD - ENP @2021-2022 55
 Description physique de la turbulence I.1 Introduction
Petites échelles de turbulence – théorie de Kolmogorov (1941)
• Le taux de dissipation (dissipation par unité de masse et de temps) est:  ~ 𝑢𝜂2 /𝜏η ~ 𝑢𝜂3 /𝜂
• Nombre de Reynolds basé sur la longueur de petite échelle 𝑅𝑒𝜂 = 𝑢𝜂 𝜂/𝜈 = 1.0 (i.e., effets
visqueux sont comparables aux effets inertiels).
• Donc 𝑅𝑒𝜂 = 𝑢𝜂 𝜂/𝜈 = 1.0 ⇒ 𝑢𝜂 = 𝜈/𝜂 et  ~ 𝑢𝜂3 /𝜂 ⇒  ~ 𝜈 3 /𝜂 4 ⇒ 𝜂 ~ 𝜈 3/4 /1/4
• Echelle de vitesse: 𝑢𝜂 = 𝜈/𝜂 avec 𝜂 ~ 𝜈 3/4 /1/4 ⇒ 𝑢𝜂 ~ 1/4 𝜈 1/4
• Echelle du temps 𝜏𝜂 ~ 𝜂 /𝑢𝜂 avec 𝑢𝜂 ~ 1/4 𝜈 1/4 et 𝜂 ~ 𝜈 3/4 /1/4 ⇒ 𝜏𝜂 ~ 𝜂 /𝑢𝜂 ~ 𝜈 1/4 /1/2
Petites échelles (échelles de Kolmogorov): l'échelle à laquelle l'énergie est dissipée.

3 1/4
Echelles de Kolmogorov ne 𝜈 𝜈 1/2
dépendent que de  et de 𝜈 .
𝜂~ 𝑢𝜂 ~ 𝜈 1/4
𝜏𝜂 ~
 
Echelle de longueur Echelle de vitesse Echelle de temps
Toutes les informations sur la géométrie sont perdues au niveau de petits tourbillons lorsque
l'énergie passe dans la cascade (car 𝜂 ≠𝑓(l)). Donc, le Mvt de petites tourbillons est universel.
Introduction à la CFD - ENP @2021-2022 56
 Description physique de la turbulence I.1 Introduction
Rapport entre grandes et petites échelles
• Selon la théorie de Kolmogorov, les taux d'énergie mise en cascade sont égaux:
⇒ Le taux de production (𝑢l /l ) égale le taux de dissipation ( ~ 𝑢𝜂3 /𝜂 ). Donc:  ~ 𝑢l /l
3 3

• Estimation du taux de dissipation: 𝑘 ~ 𝑢l2 ⇒  ~ (𝑘 1/2 )3 /l ⇒  ~𝑘 3/2 /l ⇒ l ~𝑘 3/2 /

• Rapports entre les échelles de Kolmogorov et les échelles intégrales:
3 1/4 −3/4
• D’après 𝜂 ~ 𝜈 3 / 1/4 et  ~ 𝑢l /l ⇒ 𝜂 ~
3
𝜈 l/𝑢
3
l ⇒ 𝜂 4 /l 4 ~𝜈 3 /𝑢l3 l 3 ⇒ 𝜂/l ~ 𝑅𝑒l
1/4 4 −1/4
• D’après 𝑢𝜂 ~ 𝜈 1/4 et  ~ 𝑢l /l ⇒ 𝑢𝜂 ~ 𝜈𝑢l /l
3 3
⇒ 𝑢𝜂 /𝑢l ~ 𝜈/ l 𝑢l ⇒ 𝑢𝜂 /𝑢l ~ 𝑅𝑒l
3 1/2 2 1/2
• D’après 𝜏𝜂 ~ 𝜈/ 1/2 et  ~ 𝑢l /l ⇒ 𝜏𝜂 ~ 𝜈l/𝑢l
3
⇒ 𝜏𝜂 ~ 𝜈 𝜏l /𝑢l ⇒ 𝜏𝜂2 /𝜏l ~𝜈/𝑢l2
2 1/2 −1/2
⇒ 𝜏𝜂2 /𝜏l 𝜏l ~ 𝜈/𝜏l 𝑢l ⇒ 𝜏𝜂 /𝜏l ~ 𝜈𝑢l /l𝑢l
2
⇒ 𝜏𝜂 /𝜏l ~ 𝜈/l𝑢l 1/2 ⇒ 𝜏𝜂 /𝜏l ~𝑅𝑒l

𝜂 −3/4 𝑢𝜂 −1/4 𝜏𝜂 −1/2

~ 𝑅𝑒l ~ 𝑅𝑒l ~ 𝑅𝑒l
l 𝑢l 𝜏l
Rapport d'échelle de longueur Rapport d'échelle de vitesse Rapport d'échelle de temps
Introduction à la CFD - ENP @2021-2022 57
 Simulation Numérique Directe
(Direct Numerical Simulation -DNS)

DNS RANS

Introduction à la CFD - ENP @2021-2022 58

 Simulation numérique directe I.1 Introduction
Simulation numérique directe
• La simulation numérique directe (DNS) est basée sur la résolution des équations
instantanées de Navier-Stokes, en résolvant toutes les échelles, jusqu'aux échelles de
dissipation de Kolmogorov, sans utiliser de modèles de turbulence.
−3/4
• Pas spatial nécessaire: Δ𝑥 = 𝜂 ~ 𝑅𝑒l
3/4
3
• Nombre des points nécessaire pour un volume 𝐿 : Δ𝑥 = 𝐿 /𝜂 ~ 𝑅𝑒l 3
⇒ 𝑁3𝐷 = 𝑵𝟑 ~𝑹𝒆𝟗/𝟒
l
−3/4
• Pas de temps nécessaire: Δt < Δ𝑥/𝑢l ⇒ Δt ~ 𝑅𝑒l
𝟑/𝟒
• Nombre de pas de temps nécessaire: M = 𝑇𝑒𝑚𝑝𝑠_𝑟é𝑒𝑙_𝑠𝑖𝑚𝑢𝑙𝑎𝑡𝑖𝑜𝑛 /Δt ⇒ 𝐌 ~ 𝑹𝒆l
⇒ 𝑵𝟑 𝑴 ~𝑹𝒆𝟑
𝟗/𝟒 𝟑/𝟒
• Nombre d’opérations en virgules flottantes nécessaire: 𝑵 𝑴 ~𝑹𝒆
𝟑
𝑹𝒆l l
l
• Temps nécessaire (en jours) pour effectuer une simulation DNS en supposant 1000
opérations par mod par étape pour un taux de calcul 1 gigafloop:

103 𝑁 3 𝑀 Exemple: pour 𝑅𝑒l = 105 , le temps

𝑡𝐷𝑁𝑆 = [𝐷𝑎𝑦𝑠] nécessaire de calcul est 32 ans !
109 × 60 × 60 × 24
On a besoin de modéliser la turbulence
Introduction à la CFD - ENP @2021-2022 59
 Modélisation de la turbulence

Introduction à la CFD - ENP @2021-2022 60

 Approches de modélisation de la turbulence I.1 Introduction

RANS
DNS LES
(Reynolds Averaged Navier-
(Direct Numerical Simulation) (Large Eddy Simulation)
Stokes Simulation)

▪ Résoudre les équations de

▪ Résoudre les équations
N-S moyennées dans le
▪ Résoudre les équations instationnaires filtrées.
temps (champs moyen).
instationnaires de N-S. ▪ Les échelles de turbulence
▪ Moins cher que les autres
▪ Résoudre toutes les importantes sont
approches.
échelles de turbulence. directement résolues.
▪ L'approche la plus utilisée
▪ Aucune modélisation n'est ▪ Les petits tourbillons sont
en pratique.
requise. modélisés.
▪ Manque d’information
▪ Le coût est trop prohibitif ! ▪ Moins cher que la DNS
(champ fluctuant).
▪ Pas pratique pour les ▪ Les ressources nécessaires
▪ Modèles de turbulence sont
écoulements industriels ! sont encore trop
requises pour fermer les
importantes.
équations RANS.

Introduction à la CFD - ENP @2021-2022 61

 Approches de modélisation de la turbulence I.1 Introduction
Spectre d'énergie

• Longueur d’ondes: 𝜆 = 2𝜋/l

∞
• 𝑘 = ‫׬‬0 𝐸 𝜆 𝑑𝜆; 𝐸: densité d'énergie cinétique turbulente.

Introduction à la CFD - ENP @2021-2022 62

 Description statistique de la
turbulence
& Approche RANS

Introduction à la CFD - ENP @2021-2022 63

 Description statistique de la turbulence I.1 Introduction
Description statistique
• En raison de la nature aléatoire des fluctuations turbulentes, il
est utile de les décrire à l'aide de description statistique.
• La moyenne temporelle pour une variable aléatoire 𝑢 à
𝑢ത
travers d’une période d'échantillonnage T est définie par:
𝑡+𝑇
𝑢′
1
𝑢ത = lim න 𝑢(𝑡)𝑑𝑡
𝑇→∞ T 𝑢′
𝑡
• Chaque variable (𝑢, 𝑣, 𝑤, 𝑝, 𝑇, etc.) est décomposée en une 𝑢ത (𝑡)
partie moyenne et une partie fluctuante:
𝑢= 𝑢ത + 𝑢′
Moyenne Fluctuante
• La période d'échantillonnage T est supposée être suffisamment longue pour que les
fluctuations aléatoires soient en moyenne nulles: 𝑡+𝑇
ഥ 1
𝑢′ = lim න 𝑢′ 𝑡 𝑑𝑡 = 0
𝑇→∞ T
𝑡
Introduction à la CFD - ENP @2021-2022 64
 Description statistique de la turbulence I.1 Introduction
Règles de moyenne temporelle
• Pour deux variables 𝜙1 et 𝜙2 , on peut montrer que les règles suivantes s'appliquent:

Notez que 𝜙1 ′ 𝜙2 ′ ≠ 0

• Les moyennes temporelles des dérivées spatiales sont les dérivées spatiales des variables
moyennes temporelles:

Introduction à la CFD - ENP @2021-2022 65

 Description statistique de la turbulence I.1 Introduction
Equations RANS
• Equations de N-S pour un écoulement incompressible isotherme sans forces de volume:
𝜕𝑢𝑗
=0 Continuité
𝜕𝑥𝑗
𝜕𝑢𝑖 𝜕𝑢𝑖 𝜕𝑝 𝜕 𝜕𝑢𝑖
𝜌 + 𝑢𝑗 =− + 𝜇 qté de Mvt
𝜕𝑡 𝜕𝑥𝑗 𝜕𝑥𝑖 𝜕𝑥𝑗 𝜕𝑥𝑗

• On introduisant la moyenne temporelle, l’équation de continuité s’écrit:

𝜕𝑢𝑗 𝜕 (ത
𝑢𝑗 + 𝑢𝑗 ′) 𝜕 𝑢ഥ𝑗 𝜕 𝑢𝑗 ′ 𝜕 𝑢ഥ𝑗
=0 ⇒ =0 ⇒ + =0 ⇒ =0
𝜕𝑥𝑗 𝜕𝑥𝑗 𝜕𝑥𝑗 𝜕𝑥𝑗 𝜕𝑥𝑗

𝜕𝑢𝑗 Eqt de continuité moyennée

𝜕𝑢𝑗 ′
= 0 ⇒ 𝜕 (ത
𝑢𝑗 + 𝑢𝑗 ′)
=0 ⇒ =0
𝜕𝑥𝑗 𝜕𝑥𝑗 𝜕𝑥𝑗

Le champ moyen et le champ fluctuant vérifient séparément l’équation de continuité !

Introduction à la CFD - ENP @2021-2022 66
 Description statistique de la turbulence I.1 Introduction
Equations RANS
• Equations de quantité de mouvement:
𝜕𝑢𝑖 𝜕𝑢𝑖 𝜕𝑝 𝜕 𝜕𝑢𝑖 𝜕𝑢𝑖 𝜕𝑢𝑖 𝑢𝑗 𝜕𝑝 𝜕 𝜕𝑢𝑖
𝜌 + 𝑢𝑗 =− + 𝜇 ≡ 𝜌 + =− + 𝜇
𝜕𝑡 𝜕𝑥𝑗 𝜕𝑥𝑖 𝜕𝑥𝑗 𝜕𝑥𝑗 𝜕𝑡 𝜕𝑥𝑗 𝜕𝑥𝑖 𝜕𝑥𝑗 𝜕𝑥𝑗

Car: 𝜕(𝑢𝑖 𝑢𝑗 ) = 𝑢𝑖 𝜕𝑢𝑗 + 𝑢𝑗 𝜕𝑢𝑖 = 𝑢𝑗 𝜕𝑢𝑖

𝜕𝑥𝑗 𝜕𝑥𝑗 𝜕𝑥𝑗 𝜕𝑥𝑗
continuité
• Equations de quantité de mouvement moyennées:

𝜕𝑢𝑖 𝜕𝑢𝑖 𝑢𝑗 𝜕𝑝 𝜕 𝜕𝑢𝑖 𝜕𝑢ത 𝑖 𝑢𝑗 𝑢𝑗 + 𝑢𝑗 ′𝑢𝑗 ′ 𝜕 𝑝ҧ 𝜕 𝜕𝑢ത 𝑖

𝜌
𝜕𝑡
+
𝜕𝑥𝑗
=− +
𝜕𝑥𝑖 𝜕𝑥𝑗
𝜇
𝜕𝑥𝑗
⇒ 𝜌 𝜕𝑡 + 𝜕𝑥𝑗
=− +
𝜕𝑥𝑖 𝜕𝑥𝑗
𝜇
𝜕𝑥𝑗

𝜕𝑢ത 𝑖 𝜕𝑢ത 𝑖 𝜕𝑝ҧ 𝜕 𝜕𝑢ത 𝑖

⇒ 𝜌 𝜕𝑡 + 𝑢ത𝑗 𝜕𝑥 = − 𝜕𝑥 + 𝜕𝑥 𝜇 𝜕𝑥 − 𝜌𝑢𝑖′ 𝑢𝑗′ Eqt de qté de Mvt moyennée
𝑗 𝑖 𝑗 𝑗

Introduction à la CFD - ENP @2021-2022 67

 Description statistique de la turbulence I.1 Introduction
Equations RANS
• Equations de Navier-Stokes moyennées RANS:
𝜕𝑢ത 𝑖 𝜕𝑢ത 𝑖 𝜕𝑢ത 𝑖 𝜕𝑝ҧ 𝜕 𝜕𝑢ത 𝑖
=0 𝜌 + 𝑢ത𝑗 =− + 𝜇 − 𝜌𝑢𝑖′ 𝑢𝑗′
𝜕𝑥𝑖 𝜕𝑡 𝜕𝑥𝑗 𝜕𝑥𝑖 𝜕𝑥𝑗 𝜕𝑥𝑗
𝑢′ 𝑢′ 𝑢′ 𝑣′ 𝑢′ 𝑤′
• Contraintes turbulentes (tenseur de Reynolds): 𝑅𝑖𝑗 = −𝜌𝑢𝑖′ 𝑢𝑗′ = 𝑣 ′ 𝑢′ 𝑣 ′ 𝑣′ 𝑣 ′ 𝑤′
𝑤 ′ 𝑢′ 𝑤 ′ 𝑣′ 𝑤 ′ 𝑤′
𝑅𝑖𝑗 est un tenseur symétrique: −𝜌𝑢𝑗′ 𝑢𝑖′ = −𝜌𝑢𝑖′ 𝑢𝑗′

• Nombre des inconnues = 10 (ത 𝑢𝑖 , 𝑝ҧ et les 6 composantes de 𝑢𝑖′ 𝑢𝑗′ );

• Nombre d’équations = 4 (continuité et 3 éqts de qté de Mvt).
• Le nombre d’équations < nombre d’inconnues;
• Problème de fermeture.
Modèles de turbulence

Introduction à la CFD - ENP @2021-2022 68

 Modèles de turbulence I.1 Introduction
Modèles de turbulence – approche RANS
𝜕𝑢ത 𝑖
=0
𝜕𝑥𝑖 Reynolds-stress tensor
𝜕𝑢ത 𝑖 𝜕𝑢ത 𝑖 𝜕𝑝ҧ 𝜕 𝜕𝑢ത 𝑖
𝜌 + 𝑢ത𝑗 =− + 𝜇 − 𝜌𝑢𝑖′ 𝑢𝑗′
𝜕𝑡 𝜕𝑥𝑗 𝜕𝑥𝑖 𝜕𝑥𝑗 𝜕𝑥𝑗
Boussinesq hypotheses
Reynolds
stress
model

Eddy viscosity models

Algebraic Spalart– 𝑘–𝜖 𝑘–𝜖 𝑘–𝜔

𝑘–𝜖 𝑘–𝜔 𝑘–𝜔
models Allmaras Standard
RNG Realizable Standard BSL SST
𝛾-𝑅𝑒𝜃 𝑅𝑆𝑀

0 eqt 1 eqt 4 eqts 6 eqts

2 eqts
Introduction à la CFD - ENP @2021-2022 69
 Modèles de turbulence I.1 Introduction
Hypothèse de Boussinesq
• On peut décomposer le tenseur de Reynolds comme suit:
𝑢′2 + 𝑣 ′2 + 𝑤 ′2
𝑅𝑖𝑗 = −𝜌𝑢𝑖′ 𝑢𝑗′ = 𝑅𝑖𝑠𝑜𝑡𝑟𝑜𝑝𝑒 + 𝑅𝑎𝑛𝑖𝑠𝑜𝑡𝑟𝑜𝑝𝑒 = −𝜌 + 𝑅𝑎𝑛𝑖𝑠𝑜𝑡𝑟𝑜𝑝𝑒
3
𝑢𝑖 ′𝑢𝑖 ′ 𝑢′2 + 𝑣 ′2 + 𝑤 ′2
• Puisque: 𝑘 = = ⇒ 𝑢𝑖 ′𝑢𝑖 ′ = 2𝑘
2 2 Champ moyen + hypothèse
2 de fluide Newtonien
• Donc: 𝑅𝑖𝑗 = −𝜌𝑢𝑖 𝑢𝑗 = −𝜌 𝑘𝛿𝑖𝑗 + 𝑅𝑎𝑛𝑖𝑠𝑜𝑡𝑟𝑜𝑝𝑒
′ ′
3
𝜕𝑢ത 𝑖 𝜕𝑢ത𝑗 2
Modèle de Boussinesq: 𝑅𝑖𝑗 = −𝜌𝑢𝑖′ 𝑢𝑗′ = 𝜇𝑡 + − 𝜌 𝑘𝛿𝑖𝑗
𝜕𝑥𝑗 𝜕𝑥𝑖 3

viscosité turbulente ou viscosité tourbillonnaire (eddy viscosity)

• La viscosité turbulente est fonction de l'espace et du temps. Ce n'est pas une propriété de fluide !
• Le nombre des inconnues est réduit de 6 à 2: (𝑘 et 𝜇𝑡 ).
• Dans certains modèles 𝑘 est négligée (hypothèse: contraintes de cisaillement dominent).
Introduction à la CFD - ENP @2021-2022 70
 Modèles de turbulences I.1 Introduction
Modèle algébrique (zéro équation) – longueur de mélange (Prandtl 1925)
• Analyse dimensionnelle de la viscosité turbulence:
𝜇𝑡 𝑚2 1 2 𝑑𝑈
= 𝜈𝑡 = ×𝑚 ≡ × 𝑙𝑜𝑛𝑔𝑢𝑒𝑢𝑟 2
𝜌 𝑠 𝑠 𝑑𝑦

• Prandtl (1925) a proposé le 1er modèle de turbulence:

𝑑 𝑢ത
𝜇𝑡 = 𝜌l2𝑚
𝑑𝑦

Avec: l𝒎 = 𝜿𝒚 la longueur de mélange. 𝜅 = 0.41: constante de Von Karman.

• Le modèle de Prandtl ne cause pas aucun problème de convergence !

• Il est valide pour des écoulements 2D attachés (sans séparation) et 𝜕𝑝/𝜕𝑥𝑖 faibles.
• Différents modèles ont été proposés pour améliorer la précision de modèle de
Prandtl: e.g. Van Driest (1956), Cebeci Smith (1967) et Baldwin Lomax (1978).
Introduction à la CFD - ENP @2021-2022 71
 Modèles de turbulences I.1 Introduction
Modèle algébrique (zéro équation) – Loi de paroi (Prandtl 1925)
• Soit 𝑢𝜏 la vitesse de frottement local à la paroi: 𝜏𝑤 ~𝜌𝑢𝑖′ 𝑢𝑗′ ~ 𝜌𝑢𝜏2 ⇒ 𝑢𝜏 = 𝜏𝑤 /𝜌
• Soit 𝑦𝜏 une distance proche à la paroi de taille 𝜂: 𝑦𝜏 ~ 𝜂 : 𝑅𝑒𝑦𝜏 = 𝑦𝜏 𝑢𝜏 /𝜈=1.0 ⇒ 𝑦𝜏 = 𝜈/𝑢𝜏
• On définit une vitesse 𝑢+ et une distance 𝑦 + adimensionnelle, telles que:
𝑢
ത 𝑦 +
= 𝑦/𝑦 𝑦𝑢𝜏
𝑢+ = et 𝜏 ⇒ +
𝑦 =
𝑢𝜏 𝜈
• Dans la zone très proche à la paroi, la viscosité moléculaire domine (𝜏𝑡𝑢𝑟𝑏 = 𝜌𝑢𝑖′ 𝑢𝑗′ = 0):
𝑢ത (𝑢+ 𝑢𝜏 )
𝜏𝑤 = 𝜏𝑙𝑎𝑚 = 𝜇 ⇒ 𝜏𝑤 = 𝜇 + ⇒ 𝑢+ = 𝑦 +
𝑦 (𝑦 𝜈/u𝜏 )
• Dans la zone sufisamant loin à la paroi: (𝜏𝑡𝑢𝑟𝑏 ≠ 0), on utilisant le modèle de Prandtl:
2 2 2
′ ′ 2
𝑑 𝑢ത 2
𝑑 𝑢ത 𝜏𝑤 2 2
𝑑 𝑢ത
𝜏𝑤 = 𝜏𝑡𝑢𝑟𝑏 = −𝜌𝑢𝑖 𝑢𝑗 = 𝜌l𝑚 = 𝜌 𝜅𝑦 ⇒ = 𝑢𝜏 = 𝜅𝑦
𝑑𝑦 𝑑𝑦 𝜌 𝑑𝑦
𝑑 𝑢ത 𝜈𝑦 + 𝑑 𝑢+ 𝑢𝜏 +
𝑑𝑢 +
1 𝐶 = 5.0 (Exp.
⇒ 𝑢𝜏 = 𝜅𝑦 ⇒ 𝑢𝜏 = 𝜅 ⇒ 1 = 𝜅𝑦 ⇒ + +
𝑢 = ln 𝑦 + 𝐶 parois lisses)
𝑑𝑦 𝑢𝜏 𝑑(𝑦 + 𝜈/𝑢𝜏 ) 𝑑𝑦 + 𝜅
Introduction à la CFD - ENP @2021-2022 72
 Modèles de turbulences I.1 Introduction
Loi de paroi

Inner layer
Couche interne

Outer layer
Zone linéaire (Couche externe) y est la distance normale à la paroi.
𝑢 𝜌𝑢𝜏 𝑦
=
𝑢𝜏 𝜇
La limite supérieure de la région logarithmique
dépend du nombre de Reynolds.
Couche Région entièrement
tampon turbulente (région
Sous-couche
logarithmique)
visqueuse

𝒚+ = 𝟓𝟎𝟎

Introduction à la CFD - ENP @2021-2022 73

 Modèles de turbulences I.1 Introduction
Loi de paroi

Introduction à la CFD - ENP @2021-2022 74

 Modèles de turbulences I.1 Introduction
Traitement proche paroi
Moins de nœuds sont nécessaires normales à la paroi lorsque la loi logarithmique est
utilisée (par rapport à la résolution de la sous-couche visqueuse)

High Reynolds approach Low Reynolds approach

30 ≤ y + ≤400

𝑦+ = 1
Loi logarithmique de paroi est Approche de résolution de sous-
utilisée pour résoudre la couche couche visqueuse utilisée pour
limite résoudre la couche limite

Introduction à la CFD - ENP @2021-2022 75

 Modèles de turbulences I.1 Introduction
Traitement proche paroi
Estimation la hauteur du premier élément (ou nœud) près de la paroi 𝑦𝑝 .
Estimer le coefficient de frottement:

Formule de Schlichting: 𝐶𝑓 = 2 log10 𝑅𝑒 − 0.65 −2.3 𝑓𝑜𝑟 𝑅𝑒 < 109 (Ecoulements externes)

Formule de Colebrook (ou diagramme de Moody) (Ecoulements internes)

Calculer 𝑦𝑝 en utilisant la relation inverse de 𝑦 + : Exercice:

Estimer la hauteur de la première maille
1 2
𝜏𝑤 𝜌𝑈∞ 𝑅 proche paroi nécessaire pour simuler un
𝜏𝑤 = 𝐶𝑓 𝜌𝑈∞ 𝑢𝜏 = 𝑅𝑒 =
2 𝜌 𝜇 écoulement autour d’un profile d’une éolienne
pour les modèles de turbulence suivants:
(i) Spalart-Allmaras (haut Reynolds), et
𝑦+𝜇 (ii) k − ω SST (faible Reynolds).
hauteur du 1er nœud proche de paroi: 𝑦𝑝 =
𝑢𝜏 𝜌 Données: corde = 1𝑚 , 𝑈∞ = 10 𝑚/𝑠 , 𝜇𝑎𝑖𝑟
= 1.5 × 10−5 𝑃𝑎. 𝑠, 𝜌𝑎𝑖𝑟 = 1.19 𝑘𝑔/𝑚3 .

Introduction à la modélisation de la turbulence 76

 Modèles de turbulences I.1 Introduction
Modèle à une équation – Spalart-Allmaras (1992)
𝜒3 𝜈ǁ Conçu spécialement pour les
Viscosité turbulente: 𝜇𝑡 = 𝜌𝜈𝑓
ǁ 𝜈1 𝑓𝜈1 = 3 3 ; 𝜒=
𝜒 + 𝐶𝑣1 𝜈 applications aérospatiales
2
𝜕 𝜕 1 𝜕 𝜕𝜈ǁ 𝜕𝜈ǁ
Equation de transport: 𝜌𝜈ǁ + 𝜌𝜈𝑢
ǁ 𝑖 = 𝐺𝜈 + 𝜇 + 𝜌𝜈ǁ + 𝐺𝑏2 𝜌 − 𝑌𝜈
𝜕𝑡 𝜕𝑥𝑖 𝜎𝜈෥ 𝜕𝑥𝑖 𝜕𝑥𝑗 𝜕𝑥𝑗

𝜈ǁ 𝜒 1 𝜕𝑢𝑖 𝜕𝑢𝑗
Production: 𝐺𝜈 = 𝐶𝑏1 𝜌𝜈ǁ 𝑆 + 2 2 𝑓𝜈2 ; 𝑓𝜈2 = 1 − ;𝑆 = 2𝛺𝑖𝑗 𝛺𝑖𝑗 ; 𝛺𝑖𝑗 = −
𝜅 𝑑 1 + 𝜒𝑓𝑣1 2 𝜕𝑥𝑗 𝜕𝑥𝑖

6 1Τ6 2
1 + 𝐶𝑤3 𝜈ǁ 𝜈ǁ
Destruction: 𝑌𝜈 = 𝐶𝑤1 𝜌 𝑟 + 𝐶𝑤2 𝑟6 −𝑟 6 ; 𝑟=
𝑆𝜅 2 𝑑 2 + 𝜈𝑓
ǁ 𝜈2
𝑟 + 𝐶𝑤2 𝑟 6 − 𝑟 6 + 𝐶𝑤3 𝑑

2
Constantes du modèle: 𝐶𝑏1 = 0.1335, 𝐶𝑏2 = 0.622, 𝜎𝜈෥ = , 𝐶𝜈1 = 7.1, 𝐶𝑤1 = 3.1945
3
𝐶𝑤2 = 0.3, 𝐶𝑤3 = 2.0, 𝜅 = 0.4187
Introduction à la CFD - ENP @2021-2022 77
 Modèles de turbulences I.1 Introduction
Modèle à deux équations 𝑘 − 
1
• Une équation pour l’énergie cinétique turbulente 𝑘 = 𝑢𝑖′ 𝑢𝑖′ peut être obtenue en
2
prenant la trace de l’opération suivante sur les équations de Navier-Stokes:
′ ′
𝜕𝑢𝑖 𝜕𝑢𝑖 𝜕𝑝 𝜕 𝜕𝑢𝑖
𝑢𝑖 𝑁 𝑢𝑗 + 𝑢𝑗 𝑁 𝑢𝑖 = 0 𝑁 𝑢𝑖 = 𝜌 + 𝑢𝑘 + − 𝜇 = 0 Navier-Stokes
𝜕𝑡 𝜕𝑥𝑗 𝜕𝑥𝑖 𝜕𝑥𝑘 𝜕𝑥𝑘
• Équation pour 𝑘:

Dissipation visqueuse de 𝑘: Les termes du côté droit sont complexes

et nécessitent une modélisation.
Introduction à la CFD - ENP @2021-2022 78
 Modèles de turbulences I.1 Introduction
Modèle à deux équations: 𝑘 −  (𝑆𝑡𝑎𝑛𝑑𝑎𝑟𝑑) – Launder- Spalding (1974)
~𝜌𝑘1/2 l 𝑘2
• Viscosité turbulente: l ~𝑘 3/2
/ 𝑒𝑡 𝜇 𝑇 ⇒ 𝜇𝑡 = 𝜌𝐶𝜇


𝜕 𝜌𝑘 𝜕 𝜌𝑘 𝑢ഥ𝑖 𝜕 𝜇𝑡 𝜕𝑘
• Eqt pour 𝑘: + = 𝜇+ + 𝐺𝑘 − 𝜌𝜖
𝜕𝑡 𝜕𝑥𝑖 𝜕𝑥𝑗 𝜎𝑘 𝜕𝑥𝑗

𝜕 𝜌 𝜕 𝜌𝑢ഥ𝑖 𝜕 𝜇𝑡 𝜕  2
• Eqt pour : + = 𝜇+ + 𝐶1𝜖 𝐺𝑘 − 𝐶2𝜖 𝜌
𝜕𝑡 𝜕𝑥𝑖 𝜕𝑥𝑗 𝜎𝜖 𝜕𝑥𝑗 𝑘 𝑘

𝜕𝑢ഥ𝑗 𝜕𝑢ത 𝑖 𝜕𝑢ത𝑗 Le modèle de 𝑘 −  Standard a été

La production: 𝐺𝑘 = −𝜌𝑢𝑖′ 𝑢𝑗′ = 𝜇𝑡 +
𝜕𝑥𝑖 𝜕𝑥𝑗 𝜕𝑥𝑖 conçu pour les écoulements confinés

• Constantes du modèle: 𝐶1𝜖 = 1.44, 𝐶2𝜖 = 1.92, 𝐶𝜇 = 0.09, 𝜎𝑘 = 1.0, 𝜎𝜖 = 1.3

Introduction à la CFD - ENP @2021-2022 79

 Modèles de turbulences I.1 Introduction
Autres modèles à deux équations 𝑘 − 
• Modifications ont été introduites pour améliorer la précision du modèle 𝑘 −  Standard:

𝐶𝜇 𝜂3 1 − 𝜂Τ4.38 1Τ2
• Modèle 𝑘 −  RNG: 𝐶2𝜖 = 1.68 + 𝜂 = 𝑘 2𝑆𝑖𝑗 𝑆𝑖𝑗 ൗ
1 + 𝛽𝜂3
(Yakhot et al 1992)
Développé pour tenir compte l’effet de
rotation.

1 𝛺෨ 𝑖𝑗 = 𝛺𝑖𝑗 − 2𝜖𝑖𝑗𝑘 𝑤𝑘
• Modèle 𝑘 −  Realizable: 𝐶𝜇 = 𝑘𝑈 ∗ 𝑈∗ = 𝑆𝑖𝑗 𝑆𝑖𝑗 + 𝛺෨ 𝑖𝑗 𝛺෨ 𝑖𝑗
𝐴0 + 𝐴𝑠 𝛺𝑖𝑗 = 𝛺ത𝑖𝑗 − 2𝜖𝑖𝑗𝑘 𝑤𝑘
(Shih et al 1995) 
Constante de viscosité Le modèle Realizable est développé
turbulente modifée. pour tenir compte la séparation de BL.
1 𝑆𝑖𝑗 𝑆𝑗𝑘 𝑆𝑘𝑖
Constantes: 𝐴0 = 4.04, 𝐴𝑠 = 6 cos 𝜙 , 𝜙 = cos −1 6𝑊 , 𝑊 = , ሚ=
𝑆 𝑆𝑖𝑗 𝑆𝑖𝑗
3 𝑆ሚ 3
Introduction à la CFD - ENP @2021-2022 80
 Modèles de turbulences I.1 Introduction
Modèles à deux équations 𝑘 − 𝜔
• Autre classe des modèles à 2 éqts est basée sur deux équations de transport: 𝑘 et ω.
ω: le taux de dissipation spécifique, le rapport de  à 𝑘: i.e. 𝜔 = /𝑘
𝑘
• Viscosité turbulente: 𝜇𝑡 = 𝜌𝛼 ∗ Développé pour les écoulements externes.
ω
∗
∗
𝛼 0 + 𝑅 𝑒𝑡 Τ𝑅𝑘 𝜌𝑘
• 𝒌 − 𝝎 Standard (Wilcox 1993): 𝛼 = 1 + 𝑅 𝑒 Τ𝑅 𝑅𝑒𝑡 =
𝜇𝜔
𝑡 𝑘

𝜕 𝜌𝑘 𝜕 𝜌𝑘𝑢𝑖 𝜕 𝜇𝑡 𝜕𝑘
• Eqt pour 𝑘: + = 𝜇+ + 𝐺𝑘 − 𝑌𝑘
𝜕𝑡 𝜕𝑥𝑖 𝜕𝑥𝑗 𝜎𝑘 𝜕𝑥𝑗

𝜕 𝜌𝜔 𝜕 𝜌𝜔𝑢𝑖 𝜕 𝜇𝑡 𝜕𝜔
• Eqt pour 𝜔: + = 𝜇+ + 𝐶𝜔 − 𝑌𝜔
𝜕𝑡 𝜕𝑥𝑖 𝜕𝑥𝑗 𝜎𝜔 𝜕𝑥𝑗

𝛼𝜔 𝛼∞ 𝛼0 + 𝑅 𝑒𝑡 Τ𝑅𝑤 1
𝐶𝜔 = 𝐺 𝛼= ∗ 𝛼∞ = 0.52, 𝛼0 = , 𝜎𝑘 = 𝜎𝜔 = 2.0
𝑘 𝑘 𝛼 1 + 𝑅 𝑒𝑡 Τ𝑅𝑤 9
Introduction à la CFD - ENP @2021-2022 81
 Modèles de turbulences I.1 Introduction
Modèles à deux équations: 𝑘 − 𝜔 𝑆𝑆𝑇 (Menter 1994)
• Modifications ont été introduites pour améliorer la précision du modèle 𝑘 −  de Wilcox:

∗
1 𝑘 500𝜇
𝛼 = 𝐹 = 𝑡𝑎𝑛ℎ 𝜙22 𝜙2 = 𝑚𝑎𝑥 2 ,
1 𝑆𝐹 0.09𝜔𝑦 𝜌𝑦 2 𝜔
𝑚𝑎𝑥 ∗ ,
𝛼 0.31𝜔

Amélioration du transport de Modèle k – ω transformé à partir du

Couche
contraintes turbulentes de externe modèle standard k – 
cisaillement (Shear Stress
Transport) Couche interne
(sous-couche
visqueuse, Modèle de Wilcox k – ω modifié
couche log)
Wall 𝑦+ = 1

Introduction à la CFD - ENP @2021-2022 82

 Modèles de turbulences I.1 Introduction
Reynolds stresses models (RSM)- modèles du second ordre
• Une équation exacte pour le transport du tenseur de Reynolds 𝑅𝑖𝑗 = −𝜌𝑢𝑖′ 𝑢𝑗′ peut être
obtenue en effectuant l’opération suivante : 𝑢′ 𝑁 𝑢 + 𝑢′ 𝑁 𝑢 = 0
𝑖 𝑗 𝑗 𝑖
Équations RSM:
𝜕 𝜕 𝜕
𝜌𝑢𝑖′ 𝑢𝑗′ + ′ ′
𝜌𝑢𝑘 𝑢𝑖 𝑢𝑗 = − 𝜌𝑢𝑖′ 𝑢𝑗′ 𝑢𝑘′ + 𝑝′ 𝛿𝑘𝑗 𝑢𝑖′ + 𝛿𝑖𝑘 𝑢𝑗′
𝜕𝑡 𝜕𝑥𝑘 𝜕𝑥𝑘
Local Time Derivative 𝐶𝑖𝑗 ≡ 𝐶𝑜𝑛𝑣𝑒𝑐𝑡𝑖𝑜𝑛 𝐷𝑇,𝑖𝑗 ≡ Turbulence Diffision
𝜕 𝜕 𝜕𝑢𝑗 𝜕𝑢𝑖
+ 𝜇 𝑢𝑖′ 𝑢𝑗′ −𝜌 𝑢𝑖′ 𝑢𝑘′ ′ ′
+ 𝑢𝑗 𝑢𝑘 −𝜌𝛽 𝑔𝑖 𝑢𝑗′ 𝜃 + 𝑔𝑗 𝑢𝑗′ 𝜃
𝜕𝑥𝑘 𝜕𝑥𝑘 𝜕𝑥𝑘 𝜕𝑥𝑘
𝐺𝑖𝑗 ≡Buoyancy Production
𝐷𝐿,𝑖𝑗 ≡Molucular Diffision 𝑃𝑖𝑗 ≡ Stress Production

𝜕𝑢𝑖′ 𝜕𝑢𝑗′ 𝜕𝑢𝑖′ 𝜕𝑢𝑗′

+ 𝑝′ + − 2𝜇 + − 2𝜌Ω𝑘 𝑢𝑗′ 𝑢𝑚
′
𝜖𝑖𝑘𝑚 + 𝑢𝑖′ 𝑢𝑚
′
𝜖𝑗𝑘𝑚
𝜕𝑥𝑗 𝜕𝑥𝑖 𝜕𝑥𝑘 𝜕𝑥𝑘
𝐹𝑖𝑗 ≡ Production by System Rotation
𝜙𝑖𝑗 ≡Pressure Strain 𝜖𝑖𝑗 ≡ Dissipation

Les termes 𝐷𝑇,𝑖𝑗 ,𝐺𝑖𝑗 ,𝜙𝑖𝑗 et 𝜖𝑖𝑗 doivent être modélisé. RSM évite l'hypothèse de Boussinesq.
Introduction à la CFD - ENP @2021-2022 83
 Modèles de turbulences I.1 Introduction
Exemple #1 – Écoulement turbulent à travers un obstacle rectangulaire

• Prédiction de la séparation;
• Comparaison entre différents modèles de turbulence

xR
U0
ReD = 50,000

Recirculation zone Reattachment point

Introduction à la CFD - ENP @2021-2022 84

 Modèles de turbulences I.1 Introduction
Exemple #1 – Écoulement turbulent à travers un obstacle rectangulaire

Wall

Inlet Outlet

Wall
Symmetry

Introduction à la CFD - ENP @2021-2022 85

 Modèles de turbulences I.1 Introduction
Exemple #1 – Écoulement turbulent à travers un obstacle rectangulaire

Contours of Turbulent Kinetic Energy (m2/s2)

0.70

0.63

0.56
Standard k–ε RNG k–ε
0.49

0.42

0.35

0.28

0.21
Realizable k–ε Reynolds Stress
0.14

0.07

0.00

Introduction à la CFD - ENP @2021-2022 86

 Modèles de turbulences I.1 Introduction
Exemple #1 – Écoulement turbulent à travers un obstacle rectangulaire

Predicted separation bubble:

Skin
Standard k–ε (SKE) Friction
Coefficient
Cf × 1000

Realizable k–ε (RKE)

Distance Along
Plate, x / D
SKE severely underpredicts the size of Experimentally observed
reattachment point is at
the separation bubble, while RKE
x / D = 4.7
predicts the size exactly.

Introduction à la CFD - ENP @2021-2022 87

 Modèles de turbulences I.1 Introduction
Exemple #2 – Écoulement turbulent dans un mélangeur

0.1 m
0.12 m

Uin = 20 m/s

0.2 m

0.97 m

Introduction à la CFD - ENP @2021-2022 88

 Modèles de turbulences I.1 Introduction

Modèle Avantages / Inconvénients

Spalart – Economique (1 eqt). Convenable pour les écoulements externes/internes légèrement complexes (quasi-2D) et
Allmaras aux écoulements de couche limite sous gradient de pression (par exemple, profils aérodynamiques, ailes,
fuselages d'avions, missiles, coques de navires).
Mauvaise prédiction pour les écoulements 3D complexes, les écoulements à cisaillement libre, les
écoulements avec une forte séparation.
Standard k–ε Robuste. Largement utilisé pour les écoulements internes. Convenable pour calculer les itérations initiales et
aux études paramétriques.
Fonctionne mal pour les écoulements complexes impliquant un grand gradient de pression, une séparation,
une forte courbure de ligne de courant.
Standard k–ω Performances supérieures pour la couche limite délimitée par une paroi, le cisaillement libre et les
écoulements à faible nombre de Reynolds. Convient aux écoulements complexes de la couche limite sous
gradient de pression défavorable et séparation (aérodynamique externe et turbomachines).
Sousestimation la position de transition et surestimation de la separation.
SST k–ω Offre des avantages similaires à ceux de k–ω Standard. La dépendance à la distance du mur rend cela moins
adapté aux écoulements de cisaillement libres.
N’est pas convenable pour les écoulements de transition, compressibles, géométries complexes.
Reynolds Stress Physiquement le modèle RANS le plus solide. Évite l'hypothèse de viscosité turbulente isotrope (Boussinesq).
Convenable pour les écoulements 3D complexes avec une forte courbure aérodynamique, un fort
tourbillon/rotation (par exemple, un conduit incurvé, des passages d'écoulement rotatifs, des chambres de
combustion à tourbillon avec un très grand tourbillon d'entrée, des cyclones... ).
Plus de temps CPU et de mémoire requis. Plus difficile à converger en raison du couplage étroit des équations.

Introduction à la CFD - ENP @2021-2022 89

 Simulation à Grandes Échelles
(Large Eddy Simulation - LES)

Introduction à la CFD - ENP @2021-2022 90

 Large Eddy Simulation I.1 Introduction
Large Eddy Simulation

• L'approche LES est basée sur la résolution des équations de Navier-Stokes instantanées
filtrées par un filtre spécial et modéliser les tourbillons de sous-mailles qui ne sont pas
capturés par le maillage.

Introduction à la CFD - ENP @2021-2022 91

 Large Eddy Simulation I.1 Introduction
Large Eddy Simulation
Plusieurs modèles de filtrage sont disponible

Instantaneous Resolved Subgrid

component Scale Scale
Filter, Δ

Filtered N-S
equation

(Subgrid scale Turbulent stress)

• Smagorinsky sub-grid model

Plusieurs modèles de sous-mailles existent: • WALE sub-grid model
• Dynamic SGS model
• Autres …
Introduction à la CFD - ENP @2021-2022 92
 Large Eddy Simulation I.1 Introduction
Equations de Navier-Stokes filtrées

• Une variable filtrée (indiquée par une barre supérieure) est définie par:

𝜙ത = න 𝜙 𝑥 ′ 𝐺 𝑥, 𝑥 ′ 𝑑𝑥′
𝐷

où 𝑫 est le domaine fluide, et 𝑮 la fonction de filtre qui détermine l'échelle des tourbillons résolus.
• La discrétisation en volumes finis elle-même fournit implicitement l'opération de filtrage:

où 𝑽 est le volume d'une cellule de calcul.

Introduction à la CFD - ENP @2021-2022 93

 Large Eddy Simulation I.1 Introduction
Equations de Navier-Stokes filtrées

• N-S equation:

𝜕𝑢𝑖 𝜕𝑢𝑖 𝑢𝑗 1 𝜕𝑝 𝜕 𝜕𝑢𝑖

+ =− + 𝜈
𝜕𝑡 𝜕𝑥𝑗 𝜌 𝜕𝑥𝑖 𝜕𝑥𝑗 𝜕𝑥𝑗
• Filtered N-S equation:

𝜕𝑢lj 𝑖 𝜕𝑢lj 𝑖 𝑢lj 𝑗 1 𝜕𝑝lj 𝜕 𝜕𝑢lj 𝑖 𝜕𝜏𝑖𝑗

+ =− + 𝜈 −
𝜕𝑡 𝜕𝑥𝑗 𝜌 𝜕𝑥𝑖 𝜕𝑥𝑗 𝜕𝑥𝑗 𝜕𝑥𝑗

• Sub-grid scale (SGS) stress: 𝜏𝑖𝑗 ≡ 𝑢𝑖 𝑢𝑗 − 𝑢lj 𝑖 𝑢lj 𝑗

Needs modeling

Introduction à la CFD - ENP @2021-2022 94

 Large Eddy Simulation I.1 Introduction
Modéles de sous-mailles

1
• Subgrid stress: turbulent viscosity: 𝜏𝑖𝑗 − 𝜏𝑘𝑘 𝛿𝑖𝑗 = −2𝜌𝑣𝑡 𝑆ሜ𝑖𝑗
3

• Smagorinsky model (Smagorinsky, 1963): 𝑣𝑡 = 𝐶𝑠 Δሜ 𝑆ሜ

2

Need ad-hoc near wall damping

• Dynamic model (Germano et al., 1991): 𝑣𝑡 = 𝐶𝐷 Δሜ 2 𝑆ሜ

Local adaptation of the Smagorinsky constant
1/2
• Dynamic subgrid kinetic energy transport model (Kim & Menon 2001): 𝑣𝑡 = 𝐶𝑘 𝑘𝑠𝑔𝑠 Δሜ

3/2
𝜕𝑘𝑠𝑔𝑠 𝜕𝑢lj 𝑗 𝑘𝑠𝑔𝑠 𝜕𝑢lj 𝑖 𝑘𝑠𝑔𝑠 𝜕 𝜈𝑠𝑔𝑠 𝜕𝑘𝑠𝑔𝑠
+ = −𝜏𝑖𝑗 − 𝐶𝜀 +
𝜕𝑡 𝜕𝑥𝑗 𝜕𝑥𝑗 Δሜ 𝜕𝑥𝑗 𝜎𝑘 𝜕𝑥𝑗

Introduction à la CFD - ENP @2021-2022 95

 Modèles hybrides RANS/LES
Detached Eddy Simulation (DES)

Introduction à la CFD - ENP @2021-2022 96

 Modèles hybrides: LES-RANS I.1 Introduction
Detached Eddy Simulation (DES)
• Pour les écoulements délimités par une paroi à haut nombre de Reynolds, l’approche LES
devient prohibitive pour résoudre la région proche paroi.
• L'utilisation de RANS dans les régions proches parois réduirait considérablement l'exigence de
résolution de maillage.

• Exemple: modèle hybride RANS/LES basé sur le modèle de turbulence Spalart-Allmaras:

2
𝐷 𝜈ǁ 𝜈ǁ 1 𝜕 𝜕𝜈ǁ
= 𝐶𝑏1 𝑆ሚ𝜈ǁ − 𝐶𝑤1 𝑓𝑤 + 𝜇 + 𝜌𝜈ǁ +. . .
𝐷𝑡 𝑑ሜ 𝜎𝜈෥ 𝜕𝑥𝑗 𝜕𝑥𝑗

𝑑ሜ = min 𝑑𝑤 , 𝐶DES Δ

L’approche DES est une alternative pratique à LES pour les écoulements à nombre de
Reynolds élevé dans les applications aérodynamiques externes.

Introduction à la CFD - ENP @2021-2022 97

 Références utilisées I.1 Introduction

Introduction à la CFD - ENP @2021-2022 98

 Fin
Merci pour votre attention

Dr. BOUHELAL Abdelhamid│ Pr. SMAILI Arezki Introduction à la CFD - ENP @2021-2022
 Annexe:
Exemples d’application sur
ANSYS Fluent

Introduction à la CFD - ENP @2021-2022 100

 Exemples d’applications: CFD avec Ansys Fluent
Ex#1: 2D Steady-State Heat Conduction

Consider a 2D solid material of aluminum with constant

properties (k = 202.4 W/m.K, Cp= 871 J/kg.K, and 𝜌 =
2719 kg/m3) of dimensions LxW = 0.06 m x 0.09 m,
subjected to isothermal boundary conditions as
presented in the opposite figure.
• Simulate the heat transfer in the solid material using
the ANSYS Fluent software and compare the
temperature contours with the results of the finite
difference method and the analytical solutions.

The tutorial is available at this link: https://sites.google.com/g.enp.edu.dz/dr-abdelhamid-bouhelal

Introduction à la CFD - ENP @2021-2022 101
 Exemples d’applications: CFD avec Ansys Fluent
Ex#1: 2D Steady-State Heat Conduction

Fluent (Mesh 18x12) FDM (Mesh 18x12) Analytical

Introduction à la CFD - ENP @2021-2022 102

 Exemples d’applications: CFD avec Ansys Fluent
Ex#1: 2D Steady-State Heat Conduction

Temperature at the middle line (x=L/2)

Method Max absolute
error (%)
Fluent 0.3555
FDM 0.8405

Introduction à la CFD - ENP @2021-2022 103

 Exemples d’applications: CFD avec Ansys Fluent
Ex#2: Compressible subsonic flow through a converging nozzle

Consider an inviscid flow through a convergent nozzle of geometric configuration as shown

below:

Using Ansys Fluent, simulate this flow and quantitatively confirm

the theoretical variation laws shown in the opposite figure:

The tutorial is available at this link: https://sites.google.com/g.enp.edu.dz/dr-abdelhamid-bouhelal

Introduction à la CFD - ENP @2021-2022 104
 Exemples d’applications: CFD avec Ansys Fluent
Ex#2: Compressible subsonic flow through a converging nozzle

Pression Nombre de Mach

Masse volumique
Introduction à la CFD - ENP @2021-2022 105
 Exemples d’applications I.1 Introduction
Ex#3: 2D Turbulent Pipe Flow with Heat Transfer

Velocity

Temperature

❑ Non dimensional
outlet velocity
Turbulent Kinetic
energy, k
Turbulent
dissipation rate,
epsilon

Eddy viscosity

❑ Pressure at the
centerline
Relative total
pressure

This tutorial and more applications are available at this link: https://sites.google.com/g.enp.edu.dz/dr-abdelhamid-bouhelal
Introduction à la CFD - ENP @2021-2022 106