La description lagrangienne consiste à suivre dans le temps les particules fluides le long de leurs trajectoires. Elle est souvent utilisée pour résoudre les équations d'Euler en dynamique multi-matériaux, car elle permet de préserver exactement les interfaces des matériaux au niveau discret. Cependant, elle n'est pas toujours suffisante pour décrire la dynamique des gaz compressibles. C'est le cas si le rapport entre les forces d'inertie et visqueuse est localement faible (nombre de Reynolds petit). Le terme de diffusion moléculaire des équations de Navier-Stokes n'est alors plus négligeable et doit être pris en compte. Si le nombre de Reynolds est grand, les modélisations de la turbulence (Reynolds Averaged Navier-Stokes-RANS ou Large-Eddy-Simulation-LES) sont généralement basées sur l'approximation de Boussinesq. Celle-ci consiste à approcher l'effet de la turbulence non résolue sur l'écoulement moyen par un opérateur de diffusion, proportionnel au tenseur des contraintes. Le stage a pour objectif de proposer, analyser numériquement et tester différentes stratégies de discrétisation de l'opérateur de diffusion, et son couplage au solveur de Godunov acoustique pour les équations d'Euler [1,2]. Les stratégies retenues seront testées en 1D et 2D sur maillages non réguliers dans un code écrit en C++. Une bonne connaissance du lan-gage est souhaitée. [1] B. Després et C. Mazeran, Arch. Rat. Mech. Anal., 2005, 178, pages 327-372. [2] P.-H. Maire, R. Abgrall, J. Breil et J. Ovadia, Siam J. Sci. Comp., 2007, 29.