13.07.2015 Views

Th`ese de Doctorat L'ECOLE CENTRALE DE LYON Ecole Doctorale ...

Th`ese de Doctorat L'ECOLE CENTRALE DE LYON Ecole Doctorale ...

Th`ese de Doctorat L'ECOLE CENTRALE DE LYON Ecole Doctorale ...

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.5.MODÈLE STOCHASTIQUE À UNE PARTICULE EN TURBULENCE HOMOGÈNE ISOTROPEOn ne peut pas obtenir cette valeur par intégration <strong>de</strong> la fonction <strong>de</strong> corrélationLagrangienne pendant la simulation numérique, surtout pour le temps court. Il fauttoutefois modéliser ce paramètre important pour résoudre l’équation <strong>de</strong> Langevin. Dansle chapitre 5 , on va montrer que ce paramètre joue encore un rôle important dans lemélange turbulent et la dispersion d’un scalaire.En effet, pour une turbulence homogène isotrope et statistiquement stationnaire, ona un modèle classique pour ce temps Lagrangien:T L = 2σ2 vC 0 ǫ(3.24)où C 0 est une constante, σ v l’écart type <strong>de</strong> vitesse, et ǫ la moyenne d’ensemble du taux<strong>de</strong> dissipation <strong>de</strong> l’énergie cinétique turbulente.La constante C 0 varie en fait <strong>de</strong> 2.0 à 10.0 selon <strong>de</strong>s travaux différents. Par exemple:C 0 = 2.0 ∼ 3.0 dans la thèse <strong>de</strong> Michelot [50]; Sawford [52], en utilisant les résultats<strong>de</strong>s simulations numériques <strong>de</strong> Yeung [33], trouve que C 0 doit être égale à 7; Monin& Yaglom [13] proposent une valeur <strong>de</strong> 2.1; Yeung utilise C 0 est égale à 6.5 dans sonétu<strong>de</strong> [31]. Une comparaison menée par Du [53] entre les courbes <strong>de</strong> dispersion verticaled’une source ponctuelle dans une turbulence <strong>de</strong> grille obtenues par le modèle à <strong>de</strong>uxéchelles <strong>de</strong> temps <strong>de</strong> Sawford [52] et <strong>de</strong>s mesures expérimentales conduit à préférer unevaleur <strong>de</strong> 3.0 ± 0.5 . Il trouve qu’il y a une influence du nombre <strong>de</strong> Reynolds Re λ surcette constante C 0 .De plus, Sawford [54] a proposé une formule pour modifier ce paramètre:où ˜C 0 = 6.0 .C 0 =˜C 01 + 1.75 ˜C 02Re−1.64λ(3.25)<strong>de</strong>ux vérifications <strong>de</strong> C 0Nous reviendrons sur le problème posé par cette ’constante’ car, si c’est l’uniqueparamètre du modèle stochastique à une particule et une échelle <strong>de</strong> temps, sa valeurexacte en est d’autant plus importante pour la justesse <strong>de</strong>s résultats statistiques obtenuspar simulation numérique. Nous avons vérifié la valeur <strong>de</strong> cette constante par <strong>de</strong>ux approches:1. On peut déterminer cette constante C 0 par les résultats statistiques <strong>de</strong> la simulationDNS.C 0 = 2σ2 uT L ǫ = 4 k3T L ǫ(3.26)94

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

Saved successfully!

Ooh no, something went wrong!