Chapitre 3 : Modélisation de la

turbulence

Y. EL KHCHINE Mod-Syst-Energ 1 / 26
Phénomène de la turbulence
La turbulence est un régime d'écoulement qui commence à naître au-delà
d'un seuil d'instabilité de l'écoulement laminaire, essentiellement en raison de
la faible importance des eets visqueux comparativement aux autres eets
d'inertie, et après que cette instabilité devient irréversible, plus développée et
chaotique.

Y. EL KHCHINE Mod-Syst-Energ 2 / 26
Phénomène de la turbulence
Enregistrements temporels de la vitesse en un point donné pour la même
expérience répétée 4 fois :

Y. EL KHCHINE Mod-Syst-Energ 3 / 26
Phénomène de la turbulence

Y. EL KHCHINE Mod-Syst-Energ 4 / 26
Décomposition de Reynolds
On peut donc dénir la vitesse par sa valeur moyenne ū et sa partie
uctuante u dont la moyenne temporelle sera nulle.
′

Y. EL KHCHINE Mod-Syst-Energ 5 / 26
Décomposition de Reynolds
On peut donc dénir la vitesse par sa valeur moyenne ū et sa partie
uctuante u dont la moyenne temporelle sera nulle.
′

Lorsque l'on retranche la valeur moyenne ū des enregistrements on

obtient les uctuations de u . Le problème de la turbulence consiste à
′

comprendre et modéliser l'eet de ces uctuations sur l'écoulement

moyen.

Y. EL KHCHINE Mod-Syst-Energ 5 / 26
Décomposition de Reynolds
On peut donc dénir la vitesse par sa valeur moyenne ū et sa partie
uctuante u dont la moyenne temporelle sera nulle.
′

Lorsque l'on retranche la valeur moyenne ū des enregistrements on

obtient les uctuations de u . Le problème de la turbulence consiste à
′

comprendre et modéliser l'eet de ces uctuations sur l'écoulement

moyen.
La méthode proposée par Reynolds consiste à décomposer chaque
variable (de vitesse, de pression, de température, etc...) en une partie
moyenne ϕ et une partie uctuante ϕ de moyenne temporelle nulle.
′

Y. EL KHCHINE Mod-Syst-Energ 5 / 26
Décomposition de Reynolds
′
ui = ūi + ui
où ūi est la valeur moyenne de la vitesse et ui est la uctuation de la vitesse
′

par rapport à la valeur moyenne ūi (avec u¯i′ = 0).

Y. EL KHCHINE Mod-Syst-Energ 6 / 26
Décomposition de Reynolds
′
ui = ūi + ui
où ūi est la valeur moyenne de la vitesse et ui est la uctuation de la vitesse
′

par rapport à la valeur moyenne ūi (avec u¯i′ = 0).

Cet opérateur de moyenne vérie les propriétés suivantes :
f + g = f¯ + ḡ
αf = αf¯
∂f ∂ f¯
∂x = ∂x
g f¯ = f¯ḡ
fg = f¯ḡ + f ′ g ′
αf + βg = αf¯ + β ḡ
f¯′ = g¯′ = 0
avec α et β sont des constantes.
Y. EL KHCHINE Mod-Syst-Energ 6 / 26
Décomposition de Reynolds

Régime laminaire

Pour un uide newtonien, écoulement laminaire, incompressible et pas

d'échange thermique, on a :

Y. EL KHCHINE Mod-Syst-Energ 7 / 26
Décomposition de Reynolds

Régime laminaire

Pour un uide newtonien, écoulement laminaire, incompressible et pas

d'échange thermique, on a :
=0
 ∂vi

 ∂xi
 
∂vi ∂vi ∂p
 ρ
 ∂t + vj ∂xj
= − ∂x i
+ ρgi + µ∆vi

Y. EL KHCHINE Mod-Syst-Energ 7 / 26
Décomposition de Reynolds

Régime laminaire

Pour un uide newtonien, écoulement laminaire, incompressible et pas

d'échange thermique, on a :
=0
 ∂vi

 ∂xi
 
∂vi ∂vi ∂p
 ρ
 ∂t + vj ∂xj
= − ∂x i
+ ρgi + µ∆vi

Ici, 4 équations scalaires (une équation de continuité, 3 équations de

quantité de mouvement) avec 4 inconnus (3 pour le champ de vitesse
et la pression). Donc, le système complètement résolu.

Y. EL KHCHINE Mod-Syst-Energ 7 / 26
Décomposition de Reynolds

Régime laminaire

Pour un uide newtonien, écoulement laminaire, incompressible et pas

d'échange thermique, on a :
=0
 ∂vi

 ∂xi
 
∂vi ∂vi ∂p
 ρ
 ∂t + vj ∂xj
= − ∂x i
+ ρgi + µ∆vi

Ici, 4 équations scalaires (une équation de continuité, 3 équations de

quantité de mouvement) avec 4 inconnus (3 pour le champ de vitesse
et la pression). Donc, le système complètement résolu.
Régime turbulent

En appliquant la décomposition de Reynolds :

Y. EL KHCHINE Mod-Syst-Energ 7 / 26
Décomposition de Reynolds
Equation de continuité

En introduisant la décomposition de Reynolds et en prenant la moyenne

d'ensemble, ces équations deviennent :

Y. EL KHCHINE Mod-Syst-Energ 8 / 26
Décomposition de Reynolds
Equation de continuité

En introduisant la décomposition de Reynolds et en prenant la moyenne

d'ensemble, ces équations deviennent :
 ′
  ′
  ′

∂ui ∂ u+u ∂ v +v ∂ w +w
=0⇒ + + =0
∂xi ∂x ∂y ∂z

Y. EL KHCHINE Mod-Syst-Energ 8 / 26
Décomposition de Reynolds
Equation de continuité

En introduisant la décomposition de Reynolds et en prenant la moyenne

d'ensemble, ces équations deviennent :
 ′
  ′
  ′

∂ui ∂ u+u ∂ v +v ∂ w +w
=0⇒ + + =0
∂xi ∂x ∂y ∂z
Ce qui donne
∂u ∂v ∂w ∂u ′ ∂v ′ ∂w ′
+ + + + + =0
∂x ∂y ∂z ∂x ∂y ∂z

Y. EL KHCHINE Mod-Syst-Energ 8 / 26
Décomposition de Reynolds
Equation de continuité

En introduisant la décomposition de Reynolds et en prenant la moyenne

d'ensemble, ces équations deviennent :
 ′
  ′
  ′

∂ui ∂ u+u ∂ v +v ∂ w +w
=0⇒ + + =0
∂xi ∂x ∂y ∂z
Ce qui donne
∂u ∂v ∂w ∂u ′ ∂v ′ ∂w ′
+ + + + + =0
∂x ∂y ∂z ∂x ∂y ∂z
et on a
∂u ′ ∂v ′ ∂w ′
= = =0
∂x ∂y ∂z

Y. EL KHCHINE Mod-Syst-Energ 8 / 26
Décomposition de Reynolds
Equation de continuité

En introduisant la décomposition de Reynolds et en prenant la moyenne

d'ensemble, ces équations deviennent :
 ′
  ′
  ′

∂ui ∂ u+u ∂ v +v ∂ w +w
=0⇒ + + =0
∂xi ∂x ∂y ∂z
Ce qui donne
∂u ∂v ∂w ∂u ′ ∂v ′ ∂w ′
+ + + + + =0
∂x ∂y ∂z ∂x ∂y ∂z
et on a
∂u ′ ∂v ′ ∂w ′
= = =0
∂x ∂y ∂z
Finalement
∂u ∂v ∂w
+ + =0
∂x ∂y ∂z
Y. EL KHCHINE Mod-Syst-Energ 8 / 26
Décomposition de Reynolds

En introduisant la moyenne de cette équation, on obtient :

∂u ∂v ∂w ∂u ∂v ∂w
+ + =0⇒ + + =0
∂x ∂y ∂z ∂x ∂y ∂z

Y. EL KHCHINE Mod-Syst-Energ 9 / 26
Décomposition de Reynolds

En introduisant la moyenne de cette équation, on obtient :

∂u ∂v ∂w ∂u ∂v ∂w
+ + =0⇒ + + =0
∂x ∂y ∂z ∂x ∂y ∂z

C'est l'équation de continuité moyennée

Y. EL KHCHINE Mod-Syst-Energ 9 / 26
Décomposition de Reynolds

En introduisant la moyenne de cette équation, on obtient :

∂u ∂v ∂w ∂u ∂v ∂w
+ + =0⇒ + + =0
∂x ∂y ∂z ∂x ∂y ∂z

C'est l'équation de continuité moyennée

Equations de quantité de mouvement
Les équations de quantité de mouvement moyennées :
 
∂ ūi ∂ ūi ∂ p̄ ∂ ′ ′
ρ + u¯j =− + ρgi + µ∆ūi −ρ uu
∂t ∂xj ∂xi ∂xi i j

Y. EL KHCHINE Mod-Syst-Energ 9 / 26
Décomposition de Reynolds

En introduisant la moyenne de cette équation, on obtient :

∂u ∂v ∂w ∂u ∂v ∂w
+ + =0⇒ + + =0
∂x ∂y ∂z ∂x ∂y ∂z

C'est l'équation de continuité moyennée

Equations de quantité de mouvement
Les équations de quantité de mouvement moyennées :
 
∂ ūi ∂ ūi ∂ p̄ ∂ ′ ′
ρ + u¯j =− + ρgi + µ∆ūi −ρ uu
∂t ∂xj ∂xi ∂xi i j

On pose Rij = −ρui′ uj′ c'est le tenseur des contraintes de Reynolds.

Y. EL KHCHINE Mod-Syst-Energ 9 / 26
Décomposition de Reynolds
Démonstration
∂ (u i + ui′ )
∂ (u i + ui′ ) 1 ∂ (p + p ′ ) ∂ 2 (u i + ui′ )
+ u j + uj′ =− + ḡi + ν
∂t ∂xj ρ ∂xi ∂xj ∂xj
∂ (ūi + ui′ )
∂ (ui + ui′ ) 1 ∂ (p + p ′ ) ∂ 2 (u i + ui′ )
+ uj + uj′ =− + ḡi + ν
∂t ∂xj ρ ∂xi ∂xj ∂xj

Y. EL KHCHINE Mod-Syst-Energ 10 / 26
Décomposition de Reynolds
Démonstration
∂ (u i + ui′ )
∂ (u i + ui′ ) 1 ∂ (p + p ′ ) ∂ 2 (u i + ui′ )
+ u j + uj′ =− + ḡi + ν
∂t ∂xj ρ ∂xi ∂xj ∂xj
∂ (ūi + ui′ )
∂ (ui + ui′ ) 1 ∂ (p + p ′ ) ∂ 2 (u i + ui′ )
+ uj + uj′ =− + ḡi + ν
∂t ∂xj ρ ∂xi ∂xj ∂xj
Le terme non linéaire (ui uj ) est simplié en :


ui uj = (ui + ui′ ) uj + uj′ = u i u j + u i uj′ + ui′ u j + ui′ uj′ = ui uj + ui′ uj′

Y. EL KHCHINE Mod-Syst-Energ 10 / 26
Décomposition de Reynolds
Démonstration
∂ (u i + ui′ )
∂ (u i + ui′ ) 1 ∂ (p + p ′ ) ∂ 2 (u i + ui′ )
+ u j + uj′ =− + ḡi + ν
∂t ∂xj ρ ∂xi ∂xj ∂xj
∂ (ūi + ui′ )
∂ (ui + ui′ ) 1 ∂ (p + p ′ ) ∂ 2 (u i + ui′ )
+ uj + uj′ =− + ḡi + ν
∂t ∂xj ρ ∂xi ∂xj ∂xj
Le terme non linéaire (ui uj ) est simplié en :


ui uj = (ui + ui′ ) uj + uj′ = u i u j + u i uj′ + ui′ u j + ui′ uj′ = ui uj + ui′ uj′
d'où :

∂ 2 ui ∂ui′ uj′
 
∂u i ∂u i ∂p
ρ + uj =− + ρgi + µ −ρ
∂t ∂xj ∂xi ∂xj ∂xj ∂xj

Y. EL KHCHINE Mod-Syst-Energ 10 / 26
Décomposition de Reynolds
Démonstration
∂ (u i + ui′ )
∂ (u i + ui′ ) 1 ∂ (p + p ′ ) ∂ 2 (u i + ui′ )
+ u j + uj′ =− + ḡi + ν
∂t ∂xj ρ ∂xi ∂xj ∂xj
∂ (ūi + ui′ )
∂ (ui + ui′ ) 1 ∂ (p + p ′ ) ∂ 2 (u i + ui′ )
+ uj + uj′ =− + ḡi + ν
∂t ∂xj ρ ∂xi ∂xj ∂xj
Le terme non linéaire (ui uj ) est simplié en :


ui uj = (ui + ui′ ) uj + uj′ = u i u j + u i uj′ + ui′ u j + ui′ uj′ = ui uj + ui′ uj′
d'où :

∂ 2 ui ∂ui′ uj′
 
∂u i ∂u i ∂p
ρ + uj =− + ρgi + µ −ρ
∂t ∂xj ∂xi ∂xj ∂xj ∂xj
Alors la forme nale des équations de quantité de mouvement est :
 
∂u i
ρ ∂t
+ u j ∂u
∂xj
i ∂p
= − ∂xi
+ ρgi + ∂
∂xj
(τij + Rij )

Ces équations appelées équations de Reynolds Navier-Stokes moyennées (équations

RANS).
Y. EL KHCHINE Mod-Syst-Energ 10 / 26
Equations de Reynolds Navier-Stokes moyennées
(RANS)
La contrainte totale appliquée sur une particule uide est
 
∂ ūi ∂ u¯j ′ ′
τijtot = µ + −ρui uj
∂xj ∂xi | {z }
| {z } contraintes turbulentes
contraintes visqueuses

Y. EL KHCHINE Mod-Syst-Energ 11 / 26
Equations de Reynolds Navier-Stokes moyennées
(RANS)
La contrainte totale appliquée sur une particule uide est
 
∂ ūi ∂ u¯j ′ ′
τijtot = µ + −ρui uj
∂xj ∂xi | {z }
| {z } contraintes turbulentes
contraintes visqueuses
Si on écrit la contrainte de Reynolds sous forme matricielle
 
u′ u′ u′ v ′ u′ w ′
 ′ ′
Rij = −ρ  u v v′v′ v′w′ 

u′ w ′ v′w′ w′w′

Y. EL KHCHINE Mod-Syst-Energ 11 / 26
 Equations de Reynolds Navier-Stokes moyennées
(RANS)
La contrainte totale appliquée sur une particule uide est
 
∂ ūi ∂ u¯j ′ ′
τijtot = µ + −ρui uj
∂xj ∂xi | {z }
| {z } contraintes turbulentes
contraintes visqueuses
Si on écrit la contrainte de Reynolds sous forme matricielle
 
u′ u′ u′ v ′ u′ w ′
 ′ ′
Rij = −ρ  u v v′v′ v′w′ 

u′ w ′ v′w′ w′w′
Ce tenseur symétrique possède 6 composantes indépendantes. Si on fait le
bilan des inconnues de notre système toujours composé de quatre équations,
on a cette fois-ci : ū, v̄ , w̄ , p̄ mais aussi u ′ u ′ , u ′ v ′ , u ′ w ′ , v ′ v ′ , w ′ w ′ et v ′ w ′ .

Y. EL KHCHINE Mod-Syst-Energ 11 / 26
 Equations de Reynolds Navier-Stokes moyennées
(RANS)
La contrainte totale appliquée sur une particule uide est
 
∂ ūi ∂ u¯j ′ ′
τijtot = µ + −ρui uj
∂xj ∂xi | {z }
| {z } contraintes turbulentes
contraintes visqueuses
Si on écrit la contrainte de Reynolds sous forme matricielle
 
u′ u′ u′ v ′ u′ w ′
 ′ ′
Rij = −ρ  u v v′v′ v′w′ 

u′ w ′ v′w′ w′w′
Ce tenseur symétrique possède 6 composantes indépendantes. Si on fait le
bilan des inconnues de notre système toujours composé de quatre équations,
on a cette fois-ci : ū, v̄ , w̄ , p̄ mais aussi u ′ u ′ , u ′ v ′ , u ′ w ′ , v ′ v ′ , w ′ w ′ et v ′ w ′ .
Soit au total 10 inconnues pour seulement quatre équations. Il est nécessaire
d'introduire d'autres équations pour fermer le système⇒ équations de modèle
de turbulence.
Y. EL KHCHINE Mod-Syst-Energ 11 / 26
Classication des modèles de turbulence
Modèles de premier ordre ou modèles à viscosité turbulente
• Modèle à zéro équation : longueur de mélange
• Modèle à une équation : Spalart-Allmaras
• Modèle à deux équations : k − ε Standard ; k − ε RNG ; k − ε Realizable ;
k − ω ; k − ω SST

Y. EL KHCHINE Mod-Syst-Energ 12 / 26
Classication des modèles de turbulence
Modèles de premier ordre ou modèles à viscosité turbulente
• Modèle à zéro équation : longueur de mélange
• Modèle à une équation : Spalart-Allmaras
• Modèle à deux équations : k − ε Standard ; k − ε RNG ; k − ε Realizable ;
k − ω ; k − ω SST

Modèles de seconde ordre

• Algebric Stress Model
• Reynolds Stress model

Y. EL KHCHINE Mod-Syst-Energ 12 / 26
Classication des modèles de turbulence
Modèles de premier ordre ou modèles à viscosité turbulente
• Modèle à zéro équation : longueur de mélange
• Modèle à une équation : Spalart-Allmaras
• Modèle à deux équations : k − ε Standard ; k − ε RNG ; k − ε Realizable ;
k − ω ; k − ω SST

Modèles de seconde ordre

• Algebric Stress Model
• Reynolds Stress model

Autres méthodes de résolution

• Simulation des grandes échelles (Large Eddy Simulation (LES))
• Simulation des tourbillons détachés (Detached Eddy Simulation (DES))
• Simulation numérique directe (Direct Numerical Simulation (DNS))

Y. EL KHCHINE Mod-Syst-Energ 12 / 26
Classication des modèles de turbulence

Y. EL KHCHINE Mod-Syst-Energ 13 / 26
Modélisation de la turbulence

Les modèles du premier ordre ou modèles à viscosité turbulente basés sur

l'hypothèse, dite de Boussinesq, qui consiste à modéliser directement les
contraintes de Reynolds à l'aide du concept de viscosité turbulente.

Y. EL KHCHINE Mod-Syst-Energ 14 / 26
Modélisation de la turbulence

Les modèles du premier ordre ou modèles à viscosité turbulente basés sur

l'hypothèse, dite de Boussinesq, qui consiste à modéliser directement les
contraintes de Reynolds à l'aide du concept de viscosité turbulente.
Par analogie à la contrainte visqueuse qui s'écrit comme suit :
 
∂ ūi ∂ ūj
τl = µ +
∂xj ∂xi

Y. EL KHCHINE Mod-Syst-Energ 14 / 26
Modélisation de la turbulence

Les modèles du premier ordre ou modèles à viscosité turbulente basés sur

l'hypothèse, dite de Boussinesq, qui consiste à modéliser directement les
contraintes de Reynolds à l'aide du concept de viscosité turbulente.
Par analogie à la contrainte visqueuse qui s'écrit comme suit :
 
∂ ūi ∂ ūj
τl = µ +
∂xj ∂xi

De la même manière, on dénit la contrainte turbulente

 
′ ′ ∂ ūi ∂ ūj
τt = −ρu v = µt +
∂xj ∂xi

µ propriété du uide et µt propriété de l'écoulement.

Y. EL KHCHINE Mod-Syst-Energ 14 / 26
Modélisation de la turbulence
Considérons un écoulement parallèle à l'axe x ci-dessous.

Les composantes temporelles moyennes : ū = u (y

¯ ) ; v̄ = w̄ = 0.

Y. EL KHCHINE Mod-Syst-Energ 15 / 26
Modélisation de la turbulence
Considérons un écoulement parallèle à l'axe x ci-dessous.

Les composantes temporelles moyennes : ū = u (y ¯ ) ; v̄ = w̄ = 0.

Ce tourbillon va prélever de la quantité de mouvement dans la zone supérieure
A où la vitesse vaut ū + ℓ ∂y
∂ ū
pour la transporter vers la zone B de plus basse
vitesse où la vitesse ne vaut que ū − ℓ ∂y∂ ū
(eet de mélange).
Les expériences de laboratoire montrent que u et v sont du même ordre de
′ ′

grandeur :
′ ′ ∂ ū
u ≈v ≈ℓ
∂y
Y. EL KHCHINE Mod-Syst-Energ 15 / 26
Modèles à 0 équation de transport
Modèle de longueur de mélange

A l'origine, les modèles à zéro équation déterminent algébriquement la viscosité

de turbulence à partir d'une échelle de longueur caractéristique de l'agitation.

Y. EL KHCHINE Mod-Syst-Energ 16 / 26
Modèles à 0 équation de transport
Modèle de longueur de mélange

A l'origine, les modèles à zéro équation déterminent algébriquement la viscosité

de turbulence à partir d'une échelle de longueur caractéristique de l'agitation.
Principalement conçu pour les écoulements cisaillés simples, le schéma de base
s'écrit :  
∂ ū ∂ ū
avec µt = ρℓm u = ρℓ2m
′
−ρu ′ v ′ = µt
∂y ∂y

Y. EL KHCHINE Mod-Syst-Energ 16 / 26
Modèles à 0 équation de transport
Modèle de longueur de mélange

A l'origine, les modèles à zéro équation déterminent algébriquement la viscosité

de turbulence à partir d'une échelle de longueur caractéristique de l'agitation.
Principalement conçu pour les écoulements cisaillés simples, le schéma de base
s'écrit :  
∂ ū ∂ ū
avec µt = ρℓm u = ρℓ2m
′
−ρu ′ v ′ = µt
∂y ∂y

La longueur de mélange ℓm est prescrite algébriquement.

Y. EL KHCHINE Mod-Syst-Energ 16 / 26
Modèles à 0 équation de transport
Modèle de longueur de mélange

A l'origine, les modèles à zéro équation déterminent algébriquement la viscosité

de turbulence à partir d'une échelle de longueur caractéristique de l'agitation.
Principalement conçu pour les écoulements cisaillés simples, le schéma de base
s'écrit :  
∂ ū ∂ ū
avec µt = ρℓm u = ρℓ2m
′
−ρu ′ v ′ = µt
∂y ∂y

La longueur de mélange ℓm est prescrite algébriquement.

La longueur de mélange est une propriété de l'écoulement.

Y. EL KHCHINE Mod-Syst-Energ 16 / 26
Modèles à 0 équation de transport
Modèle de longueur de mélange

A l'origine, les modèles à zéro équation déterminent algébriquement la viscosité

de turbulence à partir d'une échelle de longueur caractéristique de l'agitation.
Principalement conçu pour les écoulements cisaillés simples, le schéma de base
s'écrit :  
∂ ū ∂ ū
avec µt = ρℓm u = ρℓ2m
′
−ρu ′ v ′ = µt
∂y ∂y

La longueur de mélange ℓm est prescrite algébriquement.

La longueur de mélange est une propriété de l'écoulement.
Sa valeur reète la nature de l'agitation propre à chaque conguration.
• Jet plan : ℓm (x) = 0, 09δ0,5 (x), δ0,5 épaisseur de vitesse moitié.

Y. EL KHCHINE Mod-Syst-Energ 16 / 26
Modèles à 0 équation de transport
Modèle de longueur de mélange

A l'origine, les modèles à zéro équation déterminent algébriquement la viscosité

de turbulence à partir d'une échelle de longueur caractéristique de l'agitation.
Principalement conçu pour les écoulements cisaillés simples, le schéma de base
s'écrit :  
∂ ū ∂ ū
avec µt = ρℓm u = ρℓ2m
′
−ρu ′ v ′ = µt
∂y ∂y

La longueur de mélange ℓm est prescrite algébriquement.

La longueur de mélange est une propriété de l'écoulement.
Sa valeur reète la nature de l'agitation propre à chaque conguration.
• Jet plan : ℓm (x) = 0, 09δ0,5 (x), δ0,5 épaisseur de vitesse moitié.
• Jet rond : ℓm (x) = 0, 075δ0,5 (x).

Y. EL KHCHINE Mod-Syst-Energ 16 / 26
Modèles à 0 équation de transport
Modèle de longueur de mélange

A l'origine, les modèles à zéro équation déterminent algébriquement la viscosité

de turbulence à partir d'une échelle de longueur caractéristique de l'agitation.
Principalement conçu pour les écoulements cisaillés simples, le schéma de base
s'écrit :  
∂ ū ∂ ū
avec µt = ρℓm u = ρℓ2m
′
−ρu ′ v ′ = µt
∂y ∂y

La longueur de mélange ℓm est prescrite algébriquement.

La longueur de mélange est une propriété de l'écoulement.
Sa valeur reète la nature de l'agitation propre à chaque conguration.
• Jet plan : ℓm (x) = 0, 09δ0,5 (x), δ0,5 épaisseur de vitesse moitié.
• Jet rond : ℓm (x) = 0, 075δ0,5 (x).

2
4
• Ecoulement dans une conduite : ℓRm = 0, 14 − 0, 08 ℓ − Ry − 0, 06 ℓ − Ry ,
avec R : rayon de la conduite.
Y. EL KHCHINE Mod-Syst-Energ 16 / 26
Modèles à 0 équation de transport

Avantages :

• Facile à mettre en ÷uvre

• Temps de calcul rapides
• Bonnes prévisions pour les écoulements simples pour lesquels il existe
des corrélations expérimentales pour la longueur de mélange.

Y. EL KHCHINE Mod-Syst-Energ 17 / 26
Modèles à 0 équation de transport

Avantages :

• Facile à mettre en ÷uvre

• Temps de calcul rapides
• Bonnes prévisions pour les écoulements simples pour lesquels il existe
des corrélations expérimentales pour la longueur de mélange.
Inconvénients :

• Totalement incapable de décrire des écoulements où l'échelle de

longueur turbulente varie : tout élément avec séparation ou circulation
• Calcule uniquement les propriétés d'écoulement moyen et la contrainte
de cisaillement turbulente.

Y. EL KHCHINE Mod-Syst-Energ 17 / 26
Modèles à 0 équation de transport

Avantages :

• Facile à mettre en ÷uvre

• Temps de calcul rapides
• Bonnes prévisions pour les écoulements simples pour lesquels il existe
des corrélations expérimentales pour la longueur de mélange.
Inconvénients :

• Totalement incapable de décrire des écoulements où l'échelle de

longueur turbulente varie : tout élément avec séparation ou circulation
• Calcule uniquement les propriétés d'écoulement moyen et la contrainte
de cisaillement turbulente.
Utilisation :

• Parfois utilisé pour de simples écoulements aérodynamiques externes

• À peu près complètement ignorés dans les programmes commerciaux
de CFD aujourd'hui.

Y. EL KHCHINE Mod-Syst-Energ 17 / 26
Modèles à 1 équation de transport
Modèle Spalart-Allmaras

But : Relier la viscosité de turbulence à une propriété de l'agitation variable

dans le champ de l'écoulement.

Y. EL KHCHINE Mod-Syst-Energ 18 / 26
Modèles à 1 équation de transport
Modèle Spalart-Allmaras

But : Relier la viscosité de turbulence à une propriété de l'agitation variable

dans le champ de l'écoulement.
Formulation : Schéma de Prandtl-Kolmogorov
p
µt = Cµ ℓm k̄

Y. EL KHCHINE Mod-Syst-Energ 18 / 26
Modèles à 1 équation de transport
Modèle Spalart-Allmaras

But : Relier la viscosité de turbulence à une propriété de l'agitation variable

dans le champ de l'écoulement.
Formulation : Schéma de Prandtl-Kolmogorov
p
µt = Cµ ℓm k̄

avec
• k̄ = 21 ui ui , énergie cinétique (cinématique) locale de l'agitation turbulente,

Y. EL KHCHINE Mod-Syst-Energ 18 / 26
Modèles à 1 équation de transport
Modèle Spalart-Allmaras

But : Relier la viscosité de turbulence à une propriété de l'agitation variable

dans le champ de l'écoulement.
Formulation : Schéma de Prandtl-Kolmogorov
p
µt = Cµ ℓm k̄

avec
• k̄ = 21 ui ui , énergie cinétique (cinématique) locale de l'agitation turbulente,
• ℓ echelle de longueur caractéristique de l'agitation turbulente.

Y. EL KHCHINE Mod-Syst-Energ 18 / 26
Modèles à 1 équation de transport
Modèle Spalart-Allmaras

But : Relier la viscosité de turbulence à une propriété de l'agitation variable

dans le champ de l'écoulement.
Formulation : Schéma de Prandtl-Kolmogorov
p
µt = Cµ ℓm k̄

avec
• k̄ = 21 ui ui , énergie cinétique (cinématique) locale de l'agitation turbulente,
• ℓ echelle de longueur caractéristique de l'agitation turbulente.
Fermeture
• La valeur de k̄ s'obtient par résolution d'une équation de transport modélisée,
• La référence de longueur ℓm reste prescrite algébriquement.

Y. EL KHCHINE Mod-Syst-Energ 18 / 26
Modèles à 1 équation de transport
Modèle Spalart-Allmaras

But : Relier la viscosité de turbulence à une propriété de l'agitation variable

dans le champ de l'écoulement.
Formulation : Schéma de Prandtl-Kolmogorov
p
µt = Cµ ℓm k̄

avec
• k̄ = 21 ui ui , énergie cinétique (cinématique) locale de l'agitation turbulente,
• ℓ echelle de longueur caractéristique de l'agitation turbulente.
Fermeture
• La valeur de k̄ s'obtient par résolution d'une équation de transport modélisée,
• La référence de longueur ℓm reste prescrite algébriquement.

2
1 ∂ ∂ 2 k̄

∂ k̄ ∂ k̄ ∂ ūi ∂
ν ∂ui
+ u¯j = −ui uj − kuj − − (pui ) + ν
∂t x¯j ∂xj ∂xj 2 ∂xj ρ ∂xi ∂xi ∂xi
Y. EL KHCHINE Mod-Syst-Energ 18 / 26
Modèles à une équation de transport

Avantages :

•Ecoulements attachés au paroi

•Ecoulements avec séparation et recirculation douces

Y. EL KHCHINE Mod-Syst-Energ 19 / 26
Modèles à une équation de transport

Avantages :

•Ecoulements attachés au paroi

•Ecoulements avec séparation et recirculation douces
Inconvénients :

•Ecoulements massivement séparés

•Ecoulements de cisaillement libre.

Y. EL KHCHINE Mod-Syst-Energ 19 / 26
Modèles à 2 équation de transport
Modèle k −ε Standard

But : Eviter la prescription algébrique d'une échelle de longueur.

Y. EL KHCHINE Mod-Syst-Energ 20 / 26
Modèles à 2 équation de transport
Modèle k −ε Standard

But : Eviter la prescription algébrique d'une échelle de longueur.

La viscosité de turbulence est donnée localement par :

k̄ 2
µt = Cµ
ε̄

Y. EL KHCHINE Mod-Syst-Energ 20 / 26
Modèles à 2 équation de transport
Modèle k −ε Standard

But : Eviter la prescription algébrique d'une échelle de longueur.

La viscosité de turbulence est donnée localement par :

k̄ 2
µt = Cµ
ε̄
• k̄ energie cinétique d'agitation,
• ε̄ taux de dissipation,
Cµ coecient constant.

Y. EL KHCHINE Mod-Syst-Energ 20 / 26
Modèles à 2 équation de transport
Modèle k −ε Standard

But : Eviter la prescription algébrique d'une échelle de longueur.

La viscosité de turbulence est donnée localement par :

k̄ 2
µt = Cµ
ε̄
• k̄ energie cinétique d'agitation,
• ε̄ taux de dissipation,
Cµ coecient constant.
Les valeurs de k̄ et ε̄ s'obtiennent en résolvant des équations de transport
modélisées, qui s'écrivent à grand nombre de Reynolds de turbulence :

Y. EL KHCHINE Mod-Syst-Energ 20 / 26
Modèles à 2 équation de transport
Modèle k −ε Standard

But : Eviter la prescription algébrique d'une échelle de longueur.

La viscosité de turbulence est donnée localement par :

k̄ 2
µt = Cµ
ε̄
• k̄ energie cinétique d'agitation,
• ε̄ taux de dissipation,
Cµ coecient constant.
Les valeurs de k̄ et ε̄ s'obtiennent en résolvant des équations de transport
modélisées, qui s'écrivent à grand nombre de Reynolds de turbulence :
 
∂ k̄ ∂ k̄ ∂ ūi ∂ νt k̄
+ u¯j = −ui uj + − ε̄
∂t x¯j ∂xj ∂xj σk x¯j

ε̄2
 
∂ ε̄ ∂ ε̄ k̄ ∂ ūi ∂ νt ε̄
+ u¯j = −Cε1 ui uj + − Cε2
∂t x¯j ε̄ ∂xj ∂xj σε x¯j k̄
Y. EL KHCHINE Mod-Syst-Energ 20 / 26
Modèles à 2 équation de transport
Coecients Cµ = 0, 09 ; Cε1 = 1, 44 ; Cε2 = 1, 92 ; σk = 1 ; σε = 1, 3

Y. EL KHCHINE Mod-Syst-Energ 21 / 26
Modèles à 2 équation de transport
Coecients Cµ = 0, 09 ; Cε1 = 1, 44 ; Cε2 = 1, 92 ; σk = 1 ; σε = 1, 3
Avantages :

• Relativement simple à mettre en ÷uvre

• Conduit à des calculs stables qui convergent facilement
• Des prédictions raisonnables pour de nombreux écoulements
□ Couche limite 2D attachée,
□ Ecoulement en conduit à section circulaire,
□ Certains écoulements libres et à recirculation

Y. EL KHCHINE Mod-Syst-Energ 21 / 26
Modèles à 2 équation de transport
Coecients Cµ = 0, 09 ; Cε1 = 1, 44 ; Cε2 = 1, 92 ; σk = 1 ; σε = 1, 3
Avantages :

• Relativement simple à mettre en ÷uvre

• Conduit à des calculs stables qui convergent facilement
• Des prédictions raisonnables pour de nombreux écoulements
□ Couche limite 2D attachée,
□ Ecoulement en conduit à section circulaire,
□ Certains écoulements libres et à recirculation
Inconvénients :

Mauvaises prévisions pour :

• Ecoulement tourbillonnaire et tournant
• Ecoulement avec séparation forte
• Proche paroi (on ajoute loi de paroi)
• Certains écoulements non connés
• Ecoulements pleinement développés dans des conduits non circulaires.
Y. EL KHCHINE Mod-Syst-Energ 21 / 26
Modèles à 2 équation de transport

Modèle k −ε RNG
Ce modèle est une amélioration du modèle k − ε standard par l'utilisation de
techniques basées sur la théorie des groupes de renormalisation ( utilisé en
mécanique quantique ).

Avantages :

Prévisions améliorées pour :

• Courbure aérodynamique et vitesse de déformation élevées
• Ecoulements transitoires
• Transfert de chaleur et de masse dans les parois.
Inconvénients :

• Mais toujours ne prévoit pas la propagation d'un jet rond correctement.

Y. EL KHCHINE Mod-Syst-Energ 22 / 26
Modèles à 2 équation de transport

Modèle k −ε Realizable
Améliorations
• Partage la même équation d'énergie cinétique turbulente que le modèle
k − ε standard
• Amélioration de l'équation pour ε
• Cµ variable au lieu de constante
• Amélioration des performances pour les écoulements impliquants :
• Jets plans et ronds
• Couche limite sous forte gradients de pression ou séparation
• Rotation, recirculation
• Lignes de courant courbées

Y. EL KHCHINE Mod-Syst-Energ 23 / 26
Modèles à 2 équation de transport
Modèle k −ω
La viscosité de turbulence est donnée par :
k̄
µt =
ω
• k̄ énergie cinétique d'agitation,
• ω fréquence caractéristique des mouvements d'agitation k̄ε̄ .
Les équations de transport modélisées de k̄ et ω s'écrivent :
 
∂ k̄ ∂ k̄ ω ∂ ūi ∂ ∂ω
+ u¯j = −α ui uj + (ν + σνt ) − βω 2
∂t ∂xj k̄ ∂xj ∂xj ∂xj
 
∂ω ∂ω ∂ ūi ∂ ∗ ∂ k̄
+ u¯j = −ui uj + (ν + σ νt ) − β ∗ k̄ω
∂t ∂xj ∂xj ∂xj ∂xj
Coecients : α = 5/9 ; β = 3/40 ; β ∗ = 9/100 ; σ = 1/2 et σ ∗ = 1/2
Y. EL KHCHINE Mod-Syst-Energ 24 / 26
Modèles à 2 équation de transport

Avantages :

• Des performances bonnes pour les écoulements transitoires et ceux

avec recirculations
• Très stable numériquement et converge rapidement par rapport au
modèle k − ε
• Moins compliquée et plus économique en temps, notamment sa version
low-Re, pour le traitement du voisinage de paroi, par rapport au modèle
k − ε avec lois de paroi étant donné qu'elle ne fait pas appel au calcul de
distance à la paroi y + et d'autres grandeurs.

Y. EL KHCHINE Mod-Syst-Energ 25 / 26
Modèles à 2 équation de transport

Avantages :

• Des performances bonnes pour les écoulements transitoires et ceux

avec recirculations
• Très stable numériquement et converge rapidement par rapport au
modèle k − ε
• Moins compliquée et plus économique en temps, notamment sa version
low-Re, pour le traitement du voisinage de paroi, par rapport au modèle
k − ε avec lois de paroi étant donné qu'elle ne fait pas appel au calcul de
distance à la paroi y + et d'autres grandeurs.
Inconvénients :

• Très sensible à la condition aux limites sur ω sur les frontières libre
dans le cas des écoulements cisaillés libres.

Y. EL KHCHINE Mod-Syst-Energ 25 / 26
Modèles à 2 équation de transport

Modèle k −ω SST
Ce modèle comporte les avantages de : k − ω en proche paroi et de k − ε loin
du paroi.
Il comporte également un terme supplémentaire de diusion dans l'équation
de transport de la dissipation pour que la prédiction dans les zones transitoires
soit bonne.

Y. EL KHCHINE Mod-Syst-Energ 26 / 26