Turbulence

Universit di Lecce --- Corso di
Fluidodinamica Numerica --- Ing.

Paola CINNELLA
Modlisation et calcul
dcoulements turbulents en
Arodynamique compressible
Arts et MtiersParisTech
Cours Arodynamique et Mthodes numriques
avances
Arts et MtiersParisTech Arodynamique et Mthodes Numriques Avances Paola CINNELLA
Plan de lexpos : 1re partie
Approches statistiques
Rappels sur les coulements turbulents
Equations moyennes et contraintes de
Reynolds
Problme de fermeture
Modles RANS
Modles bass sur lapproximation de
Boussinesq
Modles aux contraintes de Reynolds
Conditions aux limites
Loi-log
Termes bas Reynolds
Applications
Plan de lexpos : 2me partie
LES et modles hybrides
Limites des modles RANS
Large Eddy Simulation
Equations de NS filtres
Modles de paroi (WMLES)
Modles hybrides RANS/LES
Motivations
Modles zonaux
Modles deux couches (two-layer models)
Autres modles zonaux
Modles globaux
Detached eddy simulation (DES)
Mthodes autres que DNS (LNS, FCM)
Conditions dentre et aux interfaces
Modlisations de la turbulence et incertitudes
Rappels sur les coulements
turbulents
TURBULENCE Pas de dfinitions rigoureuses.
Description phnomnologique.
Irrgularit
La turbulence est un tat de mouvement irrgulier, (en
apparence) alatoire
Diffusivit
Intensification des changes de masse, q. mvt et nergie
Fluctuations 3D de la vorticit
Un coulement turbulent est toujours rotationnel
Tourbillons cohrents (eddies), mais localiss de faon alatoire
Lnergie est transfre des grandes aux petites chelles de
lcoulements par effet des interactions entre structures
tourbillonnaires (mcanisme de cascade)
= v
turbulents
Augmentation significative de la diffusivit
Augmentation des contraintes de cisaillement
Exemple : abaque de Moody
turbulents
The Onset of Two-Dimensional Grid Generated Turbulence in Flowing Soap Films
Maarten A. Rutgers, Xiao-lun Wu, and Walter I. Goldberg
Visualisations dcoulements turbulents
turbulents
Mme en rgime turbulent, lcoulement est rgi
par les quations de NS (ici, en formulation
incompressible)

+ conditions initiales et aux limites
Paradoxe : quations dterministesSolutions
alatoires??
2
p
t

=

+ = +

u 0
u
u u u
Equations de Lorentz (1963)
Certains systmes non linaires sont tellement sensibles aux
perturbations, meme infinitsimales, quelles se donnent lieu
des solutions (localement) compltement diffrentes,
quoiquelles conservent certaines proprits communes
(0) 0.5, (0) 0.1, (0) 0.3 x y z = = =
(0) 0.5, (0) 0.1 1 5, (0) 0.3 x y e z = = + =
Solutions pour des conditions initiales
perturbes
Ecoulements rgis par les qs de NS
Perturbations exprimentales
Impurets
Petites fluctuations de la pression, vitesse, temprature, etc.
Petites variations de la gomtrie
Perturbations numriques
Erreurs darrondi et de troncature, erreurs sur les conditions
initiales et aux limites
La turbulence nest pas un phnomne compltement
stochastique
Les approches statistiques peuvent savrer inadquates
Cascade de lnergie
Energie injecte par lcoulement
moyen dans les grandes chelles
Energie transfre des grandes aux
petites chelles (fluctuations)
Energie dissipe lchelle
molculaire
Instabilit des tourbillons
Simulation numrique de la turbulence
Simulation directe (DNS): les quations de NS sont
rsolues pour toutes les chelles
Applications limits des coulements bas Re
Large Eddy Simulation (LES): on rsout les grandes
chelles et on modlise les petites, qui ne sont pas
captures par le maillage
Ecoulements Re moyennement lev
Approche pratique (la plus suivie jusqu prsent)
Approche statistique
On rsout des quations pour le champ moyen
Simulation numrique
de la turbulence
Simulation directe impossible pour des
coulements grand Reynolds
Rsolution spatiale Re
9/4

Large eddy simulation (LES)
prometteuse, mais encore trop couteuse
pour des applications complexes
Cot Re
0.5
(coulements sans parois)

Cot Re
2.4
(coulements avec parois)

Approche statistique
Equations de Navier-Stokes
moyennes et compltes par un
modle de turbulence
14
Simulation numrique directe (DNS)
Le calcul doit rsoudre toutes les chelles spatiales et
temporelles de lcoulement. Les simulations les plus
avances ont atteint une rsolution spatiale de 4096
3

9
3
4
Re N
Large Eddy Simulation (LES)
La simulation ne doit rsoudre que les plus grandes chelles
Dans les couches limites toutefois, ces chelles ont une taille
trs petite (tourbillons en fer de cheval par exemple)

3 2.4
Re N
Reynolds Averaged Navier-Stokes (RANS)
Les quations ne dcrivent que le champ moyen, elles
ignorent les structures cohrentes

Robustes et peu couteuses, mais elles chouent pour des
coulements dans lesquels leffet des structures cohrentes
joue un rle crucial pour la physique
Dcrochage, tremblement
Mthodes hybrides LES/ RANS cherchent runir
les meilleures proprits des deux approches

Equations moyennes la Reynolds
Description statistique dun coulement turbulent
Equations pour le champ moyen
Plusieurs dfinitions de moyenne possibles
Densemble, spatiale, temporelle, de phase
Moyenne temporelle
Bien adapte des coulements statistiquement stationnaires

On dcompose les grandeurs instantanes en une partie moyenne et
une partie fluctuante

Proprits des oprateurs de moyenne

' + =
' 0; ' ';
i i

x x

= = + =

Equations de NS incompressibles
On dcompose vitesse et pression en une partie moyenne et une
partie fluctuante

On injecte dans les quations de NS
Aprs moyenne, on aboutit aux quations de Reynolds
contraintes visqueuses
0
( ) ' ' contraintes turbulentes
1 1
( ' ' )
(tenseur de Reynolds)
ij
i
i
ij
j
i
R
i i
i j i j
i j
j i j j
u
u
x
x
u u u u
u p u
u u
t x x x x
= =
=
= =

+ = +
Systme dquations formellement identique celui de NS mais ouvert

Equations moyennes la Favre
Equations de NS compressibles

Masse volumique fluctuante
Lapplication de la dcomposition de Reynolds aux quations NS
compressible conduit un systme dquations pour le champ moyen
qui nest pas formellement identique aux quations du dpart
Cest pourquoi on adopte la dcomposition de Favre, base sur des
moyennes pondres par rapport la masse volumique, e.g. :

avec

En introduisant des dcompositions semblables pour toutes les
variables, on aboutit aux quations de NS compressibles moyennes :

Systme dquations formellement identique celui de NS mais ouvert
Fermeture des qs RANS compressibles
Analysons les nouveaux termes
Tenseur des contraintes de Reynolds
Energie cintique turbulente
Flux de chaleur turbulent
Diffusion molculaire et transport turbulent dnergie

Comme lon verra par la suite, ces deux termes sont souvent ngligs
Tenseur des contraintes de Reynolds
Souvent modlis laide de lapproximation de Boussinesq

Le dernier terme permet davoir une trace gale -2k.
Utilisation dune loi de Fourier turbulente

introduction dun nombre de Prandtl turbulent constant (0.9 pour lair)

Diffusion molculaire et transport turbulent dnergie
Ngligs pour des coulements subsoniques et dans tous les modles
algbriques, ainsi que le terme

dans lexpression du tenseur de Reynolds car

Pour des coulements supersoniques, il faut assurer la conservation exacte
de lnergie moyenne
Faute de mieux, les deux termes sont modliss ensemble par une
expression linaire du type

Spcificits des coulements
compressibles turbulents
Variations de masse volumique et de temprature
Ecoulements Mach lev changes thermiques intenses, rchauffement par
frottement transfert dnergie aux fluctuations turbulentes
Rchauffement par frottement viscosit leve prs des parois sous-
couche visqueuse plus paisse
Tout de mme, pas de variation significatives de la structure de la turbulence
pour des coulements Mach < 5, pourvu quil ny ait pas de chocs!
Pour des coulements avec transferts thermiques, pas de variations
significatives si Tw/Te < 6
Tw=temprature de la paroi; Te=temprature lextrieur de la couche limite
Hypothse de Morkovin : la turbulence inaffecte par les effets de
compressibilit si <<1
Non vrifie en cas de combustion, couches de cisaillement haute vitesse (Mach>5)
Lhypothse simplificatrice dun nombre de Prandtl turbulent constant
partout et gal 0.9 (pour lair) est trs grossire
En ralit ce nombre varie en fonction de lcoulement considr
(admis quil soit constant pour un coulement donn!)
Exemples: 0.7 pour un jet, 0.5 pour une couche de mlange,

Entre dair
Gouverne
Profil
Exemples de dcollement de CL induit par la prsence de chocs dans lcoulement
Simulation
dun cycle de tremblement
sur le profil OAT15A
M
= 0.73
Re
c
= 4.2 10
6

= 3.5
Interactions choc/couche limite
M>1.3 dcollement de la couche limite (Dlery)
Epaississement de la couche limite
Augmentation de la longueur dinfluence amont
Dviation du choc vers lamont par rapport sa position en coulement idal

Validit des approximations introduites pour des coulements
compressibles en Arodynamique
En arodynamique externe, typiquement on a Mach < 3
Les hypothses sont assez bien respectes, sauf au voisinage des chocs
En cas dinteraction choc/couche limite, la plupart des modles chouent

Modles de turbulence pour les
quations RANS compressibles
Modles algbriques
Longueur de mlange (Prandtl, 1925)
Viscosit tourbillonnaire
Baldwin-Lomax, Cebeci-Smith, Michel
Modles une quation de transport
Energie cintique turbulente
Viscosit turbulente
Spalart-Allmaras, Baldwin-Barth
Modles deux quations de transport
K-, k-, k-l, k-, q-,
Modles aux tensions de Reynolds
Equations de transport pour les 6 composantes du tenseur
de Reynolds + une chelle de longueur
Launder-Reece-Rodi, Wilcox-Rubesin,
Approx.
Boussinesq
Modles algbriques
Aussi dits 0 quations
La viscosit tourbillonnaire est calcule en fonction dune chelle
de longueur opportunment choisie (longueur de mlange)
Cette longueur est cense reprsenter lanalogue du libre parcours
moyen des particules, qui influence la viscosit molculaire
Raisonnement dimensionnel

Modle de longueur de mlange (Prandtl, 1925) :
lt =lm = longueur de mlange (donne par une relation algbrique)

t t
t
t
l u C = =

y
u
l u
m t
=
y
u
l
m t
=
2

Modles algbriques
La longueur de mlange dpend de lcoulement
considr
Diffrence fondamentale avec le libre parcours moyen qui, lui, ne
dpend que des conditions thermodynamiques du fluide
Cette longueur doit tre spcifie au cas par cas, en fonction des
chelles caractristiques de lcoulement moyen
Typiquement, elle est reprsente comme une fonction de la
distance la paroi la plus proche
Casse-tte en cas de gomtrie complexes, coulements dcolls,
solveurs non structurs
De faon gnrale, la validit de lapproximation de
Boussinesq est discutable dans beaucoup de situations
Hypothse de la longueur de mlange de Prandtl
Analogie avec entre le transport turbulent et le transport
molculaire de quantit de mouvement
Modles algbriques
y
U(y)
Molecular transport
Turbulent transport
Prandtl Mixing Length Model
Lhypothse de longueur de mlange permet de fermer le
systme dquations

Constantes de proportionnalit c
1
et c
2
spcifier
Lquation pour la viscosit turbulente fait partie de la solution,
puisquelle dpend du gradient de vitesse moyen

Evolutions du modle de base

Van Driest (1956) wall damping
Clauser (1956) defect layer modification
Corrsin and Kistler (1954) intermittency
modification

Modle de Baldwin-Lomax

Clauser
Van Driest
Corrsin and Kistler
Validit des modles algbriques
La viscosit tourbillonnaires ne dpend que des
proprits locales de lcoulement
Lcoulement ne doit pas prsenter de variations abruptes
Pas de forts gradients de pression
Pas de dcollements
Faible courbure des lignes de courant
Ils donnent de bons rsultats pour des coulements
simple, pour lesquels ils ont t calibrs
Rapides et robustes
Pas trs gnraux
Mal adapts aux solveurs non structurs
Modles quations de transport
Les modles algbriques ne tiennent pas
compte des effets dhistoire
La turbulence ne sadapte pas instantanment aux
changements abruptes!!
Ide : crire des quations de transport pour
certaines proprits reprsentatives de la
turbulence
Typiquement, quation de transport pour lnergie
cintique des fluctuations turbulentes, k
Equation de transport pour k
1. On considre lq de la q. mvt en variables primitives

2. On multiplie par u
i

3. On moyenne
4. Le rsultat de ces oprations est
ji
i i
j
j i j
u u p
u
t x x x

+ = +

Diffusion de pression
Transport turbulent de k
Diffusion moleculaire de
k

=
Termes sans
quivalent en
incompressible
souvent
ngligs
Modles une ou deux quations de transport
Toujours bass sur lapproximation de Boussinesq

Vitesse de mlange dtermine partir de lnergie cintique
turbulente

Echelle de longueur calcule par une loi algbrique de type
longueur de mlange ou bien partir dune deuxime quation
de transport, par exemple

Modles une ou deux quations de transport
Modles de Spalart-Allmaras, Baldwin-Barth
Lquation de transport est formule directement en fonction
de la viscosit tourbillonnaire ( lextrieur de la c.l.)
Spalart-Allmaras

Seconde quation

Modle k-
Viscosit tourbillonnaire

Equation de transport pour k

Equation de transport pour le taux de dissipation

Constantes du modle
Modle k-omega SST
Cherche rassembler les meilleures proprits du
k-eps (indpendance des valeurs linfini) et du k-
omega de Wilcox (robustesse, comportement
asymptotique correct lapproche des parois)

Fait recours un limiteur sur la viscosit
tourbillonnaire pour viter une dissipation
excessive dans les zones dcolles
Modle k-omega SST

Modle k-omega SST

Fonctions de relaxation F
1
et F
2

Constantes de fermeture

Points faibles de
lapproximation de Boussinesq
Ecoulements avec variations brutales du taux de
dformation moyen

Points faibles de
lapproximation de Boussinesq
Ecoulements avec forte courbure
Ecoulements dans des conduites non
axisymtriques
Ecoulements en rotation et stratifis
Ecoulements 3D
Ecoulements dcolls
Ecoulements avec chocs
On peut driver une quation exacte pour les tensions de Reynolds
(ici, en incompressible)

22 nouvelles inconnues!!!

Modles aux tensions de Reynolds:
en 3D, 7 equations additionnelles (6 contraintes + chelle de longueur)
manque de robustesse (pas de viscosit tourbillonnaire)
Modles aux tensions de Reynolds algbriques (ASM)
On ne rsout que des quations pour lnergie cintique turbulente et pour
lchelle de longueur
fonction des chelles charactristiques et du gradient de vitesse moyen
relations algbriques implicites pour
couteux en temps de calcul, parfois instable
Modles aux tensions de Reynolds algbriques explicites (EASM)
Les relations prcdantes sont rcrites de faon explicite sous une
hypothse dquilibre local
On obtient des relations constitutives ou le premier terme est de type
Boussinesq, suivi par des termes qui sont des fonctions quadratiques du
gradient de vitesse moyen
R

Gatski and Speziale (1993) EASM:

with:

Equation pour k and pour , ou autre variable permettant de
construire une chelle de longueur.
Boussinesq
approx.
Nonlinear
terms
( )
ij ij ij ij
t
W W S S
C
k
C
2
3
2
2
2
2
2 2 2 2
1
2
* *
;
6 6 3
1 3
;

= =
+ + +
+
=
|
.
|
\
|
=
( )
|
.
|
\
|
+ +
|
.
|
\
|
=
ij kl kl kj ik t ki jk kj ik t
ij ij kk ij t
R
S S S S W S W S
k S S
ij

3
1
3
2
3
1
2
5 4
58
MODELE VITESSE AXIALE VITESSE
TRANSVERSE
Spalart - Allmaras

EASM k l

Conditions aux limites
Conditions au loin
Intensit turbulente, chelle de longueur

Conditions de paroi
Sous-couche visqueuse non rsolue (loi-
log) (30 < y
+
< 100)
Sous-couche visqueuse rsolue (y
+
< 2)
Conditions aux limites, Loi-Log
Localement, on suppose un coulement de type Couette 1D
stationnaire et sans gradient de pression, avec contrainte de
cisaillement la paroi constante

Corrections prvoir pour les couches limites compressibles
( partir de Mach 5 environ)

Les contraintes visqueuses dans les ps de la q. mvt pour les
cellules proches de la paroi sont remplaces par une expression
drive partir de la loi-log
Conditions aux limites, Loi-Log
On assume une condition dquilibre production=dissipation
prs de la paroi et les termes de production et dissipation
peuvent tre calculs partir de la loi-log par:

Conditions aux limites, termes bas Reynolds
Les quations de transport pour les variables turbulentes sont
dtermines sous lhypothse de grand nombre de Reynolds
Il faut assurer un comportement consistant des quations
lapproche des parois
Conditions aux limites, termes bas Reynolds
Pour obtenir le comportement asymptotique correct on introduit des
termes damortissement visqueux dans la sous-couche

Apparition de termes sources raides, manque de robustesse,
ncessit de maillages trs fins prs des parois.
Le modle k- ne ncessite pas de modifications prs
des parois, car il prsente un comportement
asymptotique assez correct.
Les conditions aux limites utilises sont

Mme en sa version bas-Reynolds, le modle est
robuste (pas de termes-sources raides)

Modle k-omega bas-Reynolds
Ecoulement transsonique faible incidence autour dun NACA64A010

Rsultats, modle de Baldwin-Lomax
Applications
Ecoulement transsonique forte incidence autour dun NACA64A010
Memes Mach et Re, alpha=4

Applications
Exprience,
J ohnson & Bachalo (1980)
EASM k-l
k-l
Exprience
Memes Mach et Re, alpha=4
Modle Baldwin-Lomax, effect de la constante Cwake

Cwake=0.25 Cwake=1

Applications
Mmes Mach et Re, alpha=4
Modle Baldwin-Lomax, effet de la constante Cwake

Cwake=0.25 Cwake=1

Applications
Coefficient de pression la paroi
Profils de vitesse/cisaillement pour x/c=0.67 (extrados)
Model CFL Total cost
B-L 2 0.433
S-A 2 0.533
k- 2 1
k-l 2 0.833
EASM k- 1 2.33
EASM k-l 1 1.70
Cot de calcul
Convergence level: 5 orders of magnitude
Equations de NS filtres : dcomposition des grandeurs en une partie
rsolue (filtre) et une partie de sous-maille (SGS)

La vitesse filtre dpend
du temps : ce ne sont que
les plus petites chelles
qui sont filtres et doivent
tre modlises

Dfinition des variables filtres

Plusieurs filtrages possibles. Un choix typique est le box filter, qui
consiste tout simplement filtrer toutes les chelles plus petites des
cellules du maillage
La valeur locale de la proprit filtre correspond alors sa moyenne
sur une cellule
La plupart des filtres commute avec les drives en espace
De plus, ils respectent les proprits :

Linarit :

Conservation dune constante :

Attention! En gnral
u v u v + = +
Soit a cte a a = =
u u
La plupart des filtres commute avec les drives en espace
De plus, ils respectent les proprits :

Linarit :

Conservation dune constante :

Attention! En gnral
u v u v + = +
Soit a cte a a = =
u u
Filtrage des quations de Navier-Stokes
Tout comme pour les quations moyennes, le filtrage des termes
non linaires donne lieu des termes supplmentaires

Tenseur de Lonard fluctuations des interactions entre chelles
rsolues =0 pour RANS
Tenseur des contraintes croises interactions directes entre les
chelles rsolues et celles non rsolues =0 pour RANS
Tenseur de Reynolds de sous-maille effet des chelles non
rsolues sur le champ rsolu

Leonard stresses Cross stresses Reynolds stresses
avec
' ' ' '
i j i j ij ij ij
ij i j i j ij i j i j ij i j
u u u u L C R
L u u u u C u u u u R u u
= + + +
= = =
Les contraintes de Lonard (L
ij
) sont de lordre de (x)
2
et sont
souvent ngliges
Les contraintes croises (C
ij
) sont typiquement modeles avec celles
de Reynolds

Le premier modle de sous-maille (qui reste encore parmi les plus
utiliss) est le modle de Smagorinsky, qui utilise une modlisation de
type viscosit tourbillonnaire pour le tenseur de sous-maille

SGS
ij ij ij
C R = +
Modle de Smagorinsky
LES dcoulements avec parois
Les modles SGS ncessitent damortissement dans la sous-couche visqueuse, du
genre fonction de van Driest, tout comme les modles RANS
La rsolution des structures cohrentes fortement anisotropes dans la rgion de
proche paroi maillages extrmement fins
limitation sur le nombre de Reynolds
Lois de paroi ou bien couplage avec modles RANS 0, une ou deux quations
de transport dans la rgion de proche paroi

Profil de vitesse reconstruit
partir de la loi-log
Plusieurs versions possibles
LES dcoulements avec parois
Les lois de parois ne fonctionnent que pour des coulements en
quilibre
En cas dcoulements dcolls, transition, etc., recours des
modles plus avancs
URANS dans la couche proche de la paroi, LES ailleurs (WMLES)

LES dcoulements compressibles
Mach de sous-maille gnralement petit (<0.4)
effets de compressibilit ngligeables vis--vis du modle de sous-maille
la majorit des effets de compressibilit ne concernent que les chelles rsolues
Difficults lies la ncessit de capturer les chocs
En compressible les schmas numriques doivent possder une certaine dissipation
intrinsque, afin de garantir la stabilit du calcul, prserver la monotonicit des
solutions, la positivit de certaines variables, les conditions dentropie, etc.
une dissipation numrique excessive peut galer, voir dpasser celle introduite
par le modle de sous-maille
recours des schmas dordre lev, avec une faible dissipation intrinsque
approche MILES (Monotonically Integrated LES)
Combination de RANS et de LES
Motivations
Cot LES : besoin de modles RANS prs des parois
(URANS enboite dans LES)
Cot RANS : besoin de prendre en compte des
coulements complexes, souvent instationnaires et dcolls
Approche dite URANS mais Cest quoi, au juste?
Lapproche Unsteady RANS (URANS)
Dans lapproche RANS, les instationnarit sont ignores
demble (moyenne temporelle)
Problme: dans certain cas, les calcul RANS narrivent pas
converger, cycle limite sur les rsidus observ
Si calcule correctement le transitoire (schma time accurate),
on trouve en effet un champ moyen dpendant du temps et
non stationnaire
La moyenne en temps de cette solution amliore les prdictions
RANS stationnaires
Les instationnarits captures par URANS sont cense
reprsenter les (trs) grosses structures cohrentes
Problmes
Lvolution en temps est souvent incorrecte
Par exemple, URANS tend produire toujours un
comportement priodique meme quand lcoulement nest pas
vraiment priodique.
Toutefois solution moyenne plus proche des expriences que
les rsultats stationnaires
Hypothses sous-jacentes
Large sparation dchelles
Frquence de forage (ventuelle) faible devant celles des
grosses structures

Deux catgories
Approches zonales
Equations diffrentes dans diffrentes zones du maillage
Permet de mieux optimiser le temps de calcul
Permet de mieux cibler les zones traiter avec lune ou lautre
approche
Traitement des interfaces dlicat
Mal adapt aux gomtries complexes et/ou solveurs non structurs
Approches globales
On rsout les memes quations partout
Le modle bascule de faon automatique de RANS LES
Gnralit, robustesse
Controle des zones dinterface plus difficile
Deux principales stratgies dhybridisation
Fonctions de relaxation
On calcule la viscosit tourbillonnaire par un
modle RANS et un modle LES
Fonction de relaxation permettant de basculer dun
modle lautre
DES (Detached Eddy simulation)
Meme modle dans les rgions RANS et LES
Le modle bascule automatiquement du mode
RANS au mode LES par modification dune chelle
de longueur
Labsence de structures turbulentes linterface
de passage RANS/LES introduit des erreurs
Contrainte de cisaillement pas assez rsolue
Viscosit tourbillonnaire trop faible
Logarithmic layer mismatch (LLM) erreur sur Cf
Keating and Piomelli, J oT, 2006
Beaucoup de recherches visent rsoudre
ce problme ainsi que dautres lis la
transition RANS/LES
Forage stochastique
Fonction de relaxation permettant de basculer
dun modle lautre
Autres approches
Fluctuation Correction Method (Labourasse &
Sagaut 2002, Salvetti et al. 2007)
URANS partout
+ quation de transport pour la partie rsolue
+ fonction de relaxation
Filtrage hybride (Germano 2004)
termes supplmentaire de forage, actifs
dans la rgion dinterface
Beaucoup de recherches visent rsoudre
ce problme et dautre lis au rgions de
transition RANS/LES
Forage stochastique
Fonction de relaxation permettant de basculer dun
modle lautre
DES (Detached Eddy simulation)
Meme modle dans les rgions RANS et LES
Le modle bascule automatiquement du mode
RANS au mode LES par modification dune chelle
de longueur
Labsence de structures turbulentes linterface
de passage RANS/LES introduit des erreurs
Contrainte de cisaillement pas assez rsolue
Viscosit tourbillonnaire trop faible
Logarithmic layer mismatch (LLM) erreur sur Cf
Costantinescu & Squire 2005
Ecoulement autour dune sphre
Re=10000
Dcrochage massif
RANS
DES
Ecoulements turbulents et incertitudes de
modlisation
Quel modle choisir?
Cela dpend du contexte, sensibilit et expertise de lutilisateur
Pour un modle donn, deux sources dincertitudes:
Constantes du modle (calibr sur des expriences, eux m^emes
sujettes des incertitudes, cas de calibrations choisis par le
dveloppeur, recalibration pour des cas diffrents)
Forme mathmatique du modle
Besoin dtabli des intervalles de confiance sur les
rsultats des simulations!
Utilisation de ces donne pour loptimisation robuste vis--vis
des incertitudes de modlisation

Turbulence

Transféré par

Droits d'auteur :

Formats disponibles

Turbulence

Transféré par

Informations du document

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Turbulence

Transféré par

Droits d'auteur :

Formats disponibles

Universit di Lecce --- Corso di

Fluidodinamica Numerica --- Ing.

Systme dquations formellement identique celui de NS mais ouvert

Modles aux tensions de Reynolds

Vous aimerez peut-être aussi