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

Create successful ePaper yourself

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

3.5.MODÈLE STOCHASTIQUE À UNE PARTICULE EN TURBULENCE HOMOGÈNE ISOTROPEAfin <strong>de</strong> déterminer les coefficients intervenant dans l’équation <strong>de</strong> Langevin (3.2), onva se positionner dans le cas le plus simple: une turbulence homogène isotrope et statistiquementstationnaire sans vitesse moyenne.⎧⎨u i (x,y,z,t) = 〈 u i (x,y,z,t) 〉 + u ′ i(x,y,z,t)⎩〈 u i (x,y,z,t) 〉 = 0(3.4)où l’on a fait la décomposition du champ <strong>de</strong> vitesse <strong>de</strong> particule u i en sa moyenne 〈 u i 〉et sa fluctuation u ′ i .Par analogie avec les travaux <strong>de</strong> Langevin, on développe une métho<strong>de</strong> pour rechercherles variations <strong>de</strong> vitesse <strong>de</strong>s particules avec le temps. Cette approche caractérisant latrajectoire aléatoire d’une particule dans un écoulement turbulent est issue d’une analogieavec le modèle <strong>de</strong> déplacement microscopique d’une molécule suivant la théorie dumouvement brownien. On peut écrire une équation <strong>de</strong> Langevin classique sous la forme:⎧⎪⎨⎪⎩dv i (t|X 0 ,t 0 )dt〈f i 〉 = 0= av i (t|X 0 ,t 0 ) + f i(3.5)On considère qu’il y a <strong>de</strong>ux forces qui contrôlent le mouvement <strong>de</strong> la particule, l’unest la force d’amortissement qui est proportionnelle à la vitesse <strong>de</strong> la particule, l’autre estla force aléatoire qui est nulle en moyenne. La difficulté est <strong>de</strong> comprendre cette forcealéatoire.3.5 Modèle stochastique à une particule en turbulencehomogène isotrope3.5.1 Approche discrèteOn suppose dans ce paragraphe que la turbulence est homogène, isotrope et statistiquementstationnaire. On se limite ici à l’analyse du cas idéal <strong>de</strong> la dispersion turbulentetridimensionnelle. Dans la suite, les variables Eulériennes seront indicées par la lettre Eet les variables Lagrangiennes par la lettre L .La particule issue <strong>de</strong> la source sera convectée par le champ fluctuant dans les troisdirections. Elle sera déterminée par sa position X i (t|X 0 ,t 0 ) et sa vitesse v i (t|X 0 ,t 0 ) . Ons’intéresse dès lors à la détermination du modèle qui va nous donner le champ fluctuant.Dans un écoulement où le nombre <strong>de</strong> Reynolds est suffisamment grand, on pourra faire88

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

Saved successfully!

Ooh no, something went wrong!