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.

CHAPITRE 3.MODÈLE STOCHASTIQUE <strong>DE</strong> LANGEVINLa vitesse Lagrangienne <strong>de</strong>s petites échelles v > i (t|X 0,t 0 ) doit être déterminée. Onva appliquer les équations <strong>de</strong> Langevin et le modèle à une particule aux fluctuations <strong>de</strong>vitesse <strong>de</strong>s petites échelles.3.7.1 Formulation du modèleSelon les équations <strong>de</strong> Langevin en formulation discrète, on peut écrire le mouvementet la vitesse <strong>de</strong>s petites échelles <strong>de</strong>s particules sous la forme:⎧⎪⎨⎪⎩X > iv > in+1 = v> in+1 =n + v i>2n+1(1 − ∆t )TL> v i>∆tn + σ v >√∆tT > L(2 − ∆t )(3.56)TL> ξ n+1où v i> n et v i> n+1 sont respectivement la vitesse <strong>de</strong> particule Lagrangienne <strong>de</strong>s petiteséchelles à l’instant n∆t et à l’instant prochain (n+1)∆t . L’échelle <strong>de</strong> temps Lagrangienne<strong>de</strong> sous-maille est donnée par:T > L = 2σ2 v >C 0 ǫ > (3.57)où σ u> est l’écart type <strong>de</strong> vitesse <strong>de</strong>s petites échelles, et ǫ > la dissipation réalisée dans lazone <strong>de</strong>s petites échelles.3.7.2 Détermination <strong>de</strong>s gran<strong>de</strong>urs <strong>de</strong> sous-mailleEnergie moyenne <strong>de</strong>s petites échelles E sgEvaluation globale dans l’espace spectralPour déterminer le paramètre σ v > , on peut d’abord calculer l’énergie résiduelle <strong>de</strong>sous-maille E sg (en anglais, c’est subgrid energy). Naturellement, dans la simulation <strong>de</strong>sgran<strong>de</strong>s échelles, l’énergie totale peut être décomposée en <strong>de</strong>ux parties: l’une pour lesgran<strong>de</strong>s échelles, et l’autre pour les petites ( figure 3.4).E totale = E gran<strong>de</strong> + E sg (3.58)Selon la loi <strong>de</strong> Kolmogorov <strong>de</strong> 1941 ( Lesieur [15]), le spectre <strong>de</strong> l’énergie cinétiqued’une turbulence homogène isotrope à grand nombre <strong>de</strong> Reynolds, décroît en suivant une103

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

Saved successfully!

Ooh no, something went wrong!