Ce travail de thèse s'inscrit dans une problématique visant à étudier numériquement le bruit d'origine aérodynamique généré par les écoulements turbulents autour des véhicules en utilisant la méthode Boltzmann sur Réseau (LBM). Les objectifs de cette thèse sont l'étude des capacités aéroacoustiques de la LBM ainsi que l'élaboration d'un code de calcul tridimensionnel et parallèle.<br />Dans un premier temps, les élements historiques et théoriques de la LBM sont présentés ainsi que le développement permettant de passer de l'équation de Boltzmann aux équations de Navier-Stokes. La construction des modèles à vitesses discrètes est également décrite. Deux modèles basés sur des opérateurs de collision différents sont présentés : le modèle LBM-BGK et le modèle LBM-MRT. Pour l'étude des capacités aéroacoustiques de la LBM, une analyse de von Neumann est réalisée pour les modèles LBM-BGK et LBM-MRT ainsi que pour l'équation de Boltzmann à vitesse discrète (DVBE). Une comparaison avec les schémas Navier-Stokes d'ordre élevé est alors menée. Pour remédier aux instabilités numériques de la méthode Boltzmann sur Réseau intervenant lors de la propagation dans des directions particulières à M>0.1, des filtres sélectifs sont utilisés et leur effet sur la dissipation est étudié.<br />Dans un second temps, le code de calcul L-BEAM est présenté. La structure générale et les différentes techniques de calculs sont décrites. Un algorithme de transition de résolution est développé. La modélisation de la turbulence est abordée et le modèle de Meyers-Sagaut est implémenté dans le code. Enfin, des cas tests numériques sont utilisés pour valider le code et la simulation d'un écoulement turbulent complexe est réalisée.

Les calculs d'écoulements stationnaires peuvent être considérés comme efficace si l?état stationnaire est atteint pour un temps CPU réduit mais aussi si la place mémoire utilisée reste faible; cette dernière exigence devient primordiale pour les applications industrielles où le nombre de points de calcul est très important. Ceci vaut également pour les écoulements instationnaires, désormais classiquement résolus via une approche pas-de-temps dual pour laquelle les états physiques successifs sont vus comme des états stationnaires vis-à-vis d'un temps fictif. Le besoin crucial de méthodes implicites à faible encombrement mémoire a conduit au développement de traitements sans matrice. Pour les applications qui intéressent le CEA, à savoir la simulation d'écoulements réactifs multi-espèces à l'intérieur d'une enceinte de réacteur nucléaire à eau pressurisée, les méthodes doivent être assez versatiles pour traiter la gamme d'écoulements allant du quasi-incompressible au fortement compressible. Le préconditionnement bas-Mach des équations de Navier-Stokes permet d'appliquer en régime incompressible les schémas initialement conçus pour la simulation des écoulements compressibles. Le présent travail montre comment obtenir un traitement implicite sans matrice pour tout régime d'écoulement lorsque la phase implicite contient une matrice de préconditionnement; l'efficacité intrinsèque du schéma implicite sans matrice couplé à une technique de relaxation de type Jacobi par point (PJ) ou Symmetric Gauss-Seidel (SGS) est étudiée grâce à une analyse de Von Neumann; puis des comparaisons avec des méthodes implicites blocs standards sont effectuées. La méthode implicite sans matrice est finalement implémentée au sein du code non-structuré CAST3M et elle est appliquée à la modélisation d'un Té de mélange à faible nombre de Mach. Le schéma implicite sans matrice constitue une alternative compétitive pour la simulation des écoulements compressibles et faiblement compressibles en maillages non-structurés.