13.07.2015 Views

FLUIDES EN ÉCOULEMENT Méthodes et modèles Jacques PADET

FLUIDES EN ÉCOULEMENT Méthodes et modèles Jacques PADET

FLUIDES EN ÉCOULEMENT Méthodes et modèles Jacques PADET

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

3.4.3.1. – MODÈLES À BAS NOMBRE DE REYNOLDSUn défaut du modèle k - ε standard est qu’il surévalue le frottement <strong>et</strong> le flux dechaleur dans les écoulements de couche limite avec gradient de pression adverse.Les <strong>modèles</strong> dits (improprement) « à bas Reynolds » essaient de corriger c<strong>et</strong>inconvénient. Ils consistent à introduire des termes correctifs pour mieux tenir compte deseff<strong>et</strong>s de parois. Ceux-ci prennent généralement la forme de fonctions d’amortissement,notées f µ <strong>et</strong> f 2 .La première, f µ , est appliquée à la viscosité turbulente, qui devient :2kν t = Cµf µ(3.71)εPar ricoch<strong>et</strong>, dans l’équation en ε, le facteur du premier terme C ε1 C µ k, qui s’écrivaitaussi en raison de (3.60) :εCε 1 Cµk = Cε1 ν tkdevient de ce fait : C f kC ε 1 µ µLa seconde, f 2 , est appliquée au dernier terme de l’équation en ε.Formellement, dans (3.68) l’équation en k est donc inchangée tandis que l’équation enε est corrigée comme suit :Vi∂ε∂xi= Cε 1Cµfµk⎛⎜∂V⎝∂xij+∂V∂xij⎞⎟⎠∂V∂xij+∂∂xi⎪⎧⎛⎨⎜ν+⎪⎩ ⎝ν t ⎞ ∂ε⎪⎫− Cσ⎟ ⎬ε ⎠ ∂xi⎪⎭ε 2f22εk(3.72)Quant à ces fonctions d’amortissement f µ <strong>et</strong> f 2 , elles ont été ajustées (pour ne pas direbricolées) à l’aide d’un nouveau paramètre, le « nombre de Reynolds turbulent » R t défini par :2kR t = (3.73)ν εL’appellation inattendue de « Reynolds turbulent » trouve une justification partiellegrâce aux échelles de turbulence (§ 3.6.4) mais elle est malencontreuse, car R t n’a rien à voiravec le nombre de Reynolds Re de l’écoulement.Parmi les diverses fonctions d’amortissement qui ont été proposées, citons seulementcomme exemples :−2,5( 1 + Rt / 50f = e) 2µ (Jones <strong>et</strong> Launder) (3.74a)f22Rt−36= 1 − 0,22 e(Hanjalic) (3.74b)Mais pourquoi appelle–t–on ces <strong>modèles</strong> « bas Reynolds » ? C’est que, d’après (3.60)<strong>et</strong> (3.73) :

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!