27.06.2013 Views

Simulation numérique de l'essorage et du refroidissement d'un film ...

Simulation numérique de l'essorage et du refroidissement d'un film ...

Simulation numérique de l'essorage et du refroidissement d'un film ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

2.5 Modélisations <strong>de</strong> la turbulence 51<br />

Le problème rési<strong>de</strong> dans la détermination <strong>de</strong> la longueur caractéristique correcte pour<br />

utiliser le modèle dans le cas d’un filtre anisotrope. Deux approches peuvent être utilisées :<br />

– La première consiste à définir une seule échelle <strong>de</strong> longueur pour représenter le<br />

filtre. Les modèles restent alors les mêmes que ceux définis pour le cas isotrope,<br />

en utilisant par exemple les viscosités <strong>de</strong> sous-maille scalaires pour représenter le<br />

procédé <strong>de</strong> casca<strong>de</strong> inverse. Ceci implique une légère modification <strong>de</strong>s modèles <strong>de</strong><br />

sous-maille puisque seul le calcul <strong>de</strong> l’échelle <strong>de</strong> coupure caractéristique est modifié.<br />

C<strong>et</strong>te approche est valable en théorie seulement pour les cas d’anisotropie faible,<br />

pour laquelle les longueurs <strong>de</strong> coupure sont <strong>du</strong> même ordre <strong>de</strong> magnitu<strong>de</strong>.<br />

– La secon<strong>de</strong> approche est basée sur l’intro<strong>du</strong>ction <strong>de</strong> différentes échelles <strong>de</strong> longueur<br />

dans le modèle. Ceci implique parfois <strong>de</strong>s modifications majeures dans le modèle<br />

isotrope, telles que la définition <strong>de</strong>s viscosités tensorielles <strong>de</strong> sous-maille pour re-<br />

présenter le procédé <strong>de</strong> casca<strong>de</strong> inverse. En théorie, c<strong>et</strong>te approche prend mieux en<br />

compte l’anisotropie <strong>du</strong> filtre, mais rend plus difficile l’étape <strong>de</strong> modélisation.<br />

Certains auteurs [Fureby 99] pensent possible <strong>de</strong> se dispenser <strong>de</strong> modèle <strong>de</strong> sous-maille<br />

pour faire tourner un calcul LES en utilisant simplement un schéma <strong>du</strong> premier ordre pour<br />

compenser la casca<strong>de</strong> d’énergie <strong>et</strong> conserver la régularité <strong>de</strong> la solution. Cependant un tel<br />

calcul ne peut traiter les régions <strong>de</strong> proche paroi <strong>de</strong>s couches limites.<br />

Diffusivité <strong>de</strong> sous-maille<br />

La simulation <strong>de</strong>s gran<strong>de</strong>s échelles appliquée aux écoulements turbulents anisothermes<br />

implique l’intro<strong>du</strong>ction d’une diffusivité <strong>de</strong> sous-maille calculée à partir <strong>de</strong> la viscosité <strong>de</strong><br />

sous-maille <strong>du</strong> modèle hydrodynamique.<br />

Les composantes <strong>du</strong> vecteur <strong>de</strong> flux <strong>de</strong> sous-maille ΠjT vu à la relation (2.101) peuvent<br />

s’écrire en fonction <strong>du</strong> gradient <strong>de</strong> température résolue <strong>et</strong> d’une diffusivité <strong>de</strong> sous-maille<br />

asm :<br />

∂T<br />

ΠjT = −asm<br />

(2.144)<br />

∂xj<br />

C<strong>et</strong>te diffusivité asm est souvent obtenue à partir <strong>de</strong> l’analogie <strong>de</strong> Reynolds en définis-<br />

sant un nombre <strong>de</strong> Prandtl <strong>de</strong> sous-maille tel que :<br />

asm = νsm · P rsm<br />

(2.145)

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

Saved successfully!

Ooh no, something went wrong!