Le relief Alpin modifie de manière importante la circulation de l'atmosphère. Grâce à l'étude de différentes configurations, en période hivernale et estivale, et différentes régions des Alpes, ce travail expose comment se réalisent le mélange et le transport des polluants gazeux et particulaires depuis leurs sources d'émissions jusque dans la troposphère libre. Grâce à des résultats d'observation réalisés dans le Y grenoblois ces dernières années et des simulations numériques provenant d'une chaîne de modèles eulériens, ce travail explicite les mécanismes mis en jeu dans un cas de pollution photochimique estivale et un cas de pollution hivernale par les particules fines et par le benzène. Après la validation des modèles en utilisant des données expérimentales des campagnes POVA, GRENOPHOT 99 et ESCOMPTE, la chaîne de modèles est utilisée pour simuler deux périodes très polluées dans la région de Grenoble : la canicule d'août 2003 et un long épisode d'inversion thermique en février 2005. Après coup, les incertitudes sur les résultats du modèle sont calculées. En appliquant la chaîne de modèles sur un événement très particulier faisant intervenir de l'ozone provenant de la stratosphère on montre comment les différentes échelles de temps et d'espace de l'atmosphère sont intimement couplées dans le relief complexe des Alpes.

Cette thèse s'inscrit dans un projet de modélisation numérique de la dispersion atmosphérique de polluants à travers le code de mécanique des fluides numérique Code_Saturne. L'objectif est de pouvoir simuler la dispersion atmosphérique de polluants en environnement complexe, c'est-à-dire autour de centrales, sites industriels ou en milieu urbain. Dans ce contexte, nous nous concentrons sur la modélisation de la dispersion des polluants à micro-échelle, c'est-à-dire pour des distances de l'ordre de quelques mètres à quelques kilomètres et correspondant à des échelles de temps de l'ordre de quelques dizaines de secondes à quelques dizaines de minutes : on parle de modélisation en champ proche. L’approche suivie dans ces travaux de recherche suit une formulation hybride eulérienne/lagrangienne, où les champs dynamiques moyens relatifs au fluide porteur (pression, vitesse, température, turbulence) sont calculés via une approche eulérienne et sont ensuite fournis au solveur lagrangien. Ce type de formulation est couramment utilisé dans la littérature atmosphérique pour son efficacité numérique. Le modèle lagrangien stochastique considéré dans nos travaux est le Simplified Langevin Model (SLM), développé par Pope (1985,2000). Ce modèle appartient aux méthodes communément appelées méthodes PDF (Probability Density Function), et, à notre connaissance, n'a pas été exploité auparavant dans le contexte de la dispersion atmosphérique. Premièrement, nous montrons que le SLM respecte le critère dit de mélange homogène (Thomson, 1987). Ce critère, essentiel pour juger de la bonne qualité d'un modèle lagrangien stochastique, correspond au fait que si des particules sont initialement uniformément réparties dans un fluide incompressible, alors elles doivent le rester. Nous vérifions le bon respect du critère de mélange homogène pour trois cas de turbulence inhomogène représentatifs d'une large gamme d'applications pratiques : une couche de mélange, un canal plan infini, ainsi qu'un cas de type atmosphérique mettant en jeu un obstacle au sein d'une couche limite neutre. Nous montrons que le bon respect du critère de mélange homogène réside simplement en la bonne introduction du terme de gradient de pression en tant que terme de dérive moyen dans le modèle de Langevin (Pope, 1987; Minier et al., 2014; Bahlali et al., 2018c). Nous discutons parallèlement de l'importance de la consistance entre champs eulériens et lagrangiens dans le cadre de telles formulations hybrides eulériennes/lagrangiennes. Ensuite, nous validons le modèle dans le cas d'un rejet de polluant ponctuel et continu, en conditions de vent uniforme et turbulence homogène. Dans ces conditions, nous disposons en effet d'une solution analytique nous permettant une vérification précise. Nous observons que dans ce cas, le modèle lagrangien discrimine bien les deux différents régimes de diffusion de champ proche et champ lointain, ce qui n'est pas le cas d'un modèle eulérien à viscosité turbulente (Bahlali et al., 2018b).Enfin, nous travaillons sur la validation du modèle sur plusieurs campagnes expérimentales en atmosphère réelle, en tenant compte de la stratification thermique de l'atmosphère et de la présence de bâtiments. Le premier programme expérimental considéré dans nos travaux concerne le site du SIRTA (Site Instrumental de Recherche par Télédétection Atmosphérique), dans la banlieue sud de Paris, et met en jeu une stratification stable de la couche limite atmosphérique. La seconde campagne étudiée est l'expérience MUST (Mock Urban Setting Test). Réalisée aux Etats-Unis, dans le désert de l'Utah, cette expérience a pour but de représenter une ville idéalisée, au travers d'un ensemble de lignées de conteneurs. Deux rejets ont été simulés et analysés, respectivement en conditions d'atmosphère neutre et stable (Bahlali et al., 2018a) / This Ph.D. thesis is part of a project that aims at modeling pollutant atmospheric dispersion with the Computational Fluid Dynamics code Code_Saturne. The objective is to simulate atmospheric dispersion of pollutants in a complex environment, that is to say around power plants, industrial sites or in urban areas. In this context, the focus is on modeling the dispersion at micro-scale, that is for distances of the order of a few meters to a few kilometers and corresponding to time scales of the order of a few tens of seconds to a few tens of minutes: this is also called the near field area. The approach followed in this thesis follows a hybrid Eulerian/Lagrangian formulation, where the mean dynamical fields relative to the carrier fluid (pressure, velocity, temperature, turbulence) are calculated through an Eulerian approach and are then provided to the Lagrangian solver. This type of formulation is commonly used in the atmospheric literature for its numerical efficiency. The Lagrangian stochastic model considered in our work is the Simplified Langevin Model (SLM), developed by Pope (1985,2000). This model belongs to the methods commonly referred to as PDF (Probability Density Function) methods, and, to our knowledge, has not been used before in the context of atmospheric dispersion. First, we show that the SLM meets the so-called well-mixed criterion (Thomson, 1987). This criterion, essential for any Lagrangian stochastic model to be regarded as acceptable, corresponds to the fact that if particles are initially uniformly distributed in an incompressible fluid, then they must remain so. We check the good respect of the well-mixed criterion for three cases of inhomogeneous turbulence representative of a wide range of practical applications: a mixing layer, an infinite plane channel, and an atmospheric-like case involving an obstacle within a neutral boundary layer. We show that the good respect of the well-mixed criterion lies simply in the good introduction of the pressure gradient term as the mean drift term in the Langevin model (Pope, 1987; Minier et al., 2014; Bahlali et al., 2018c). Also, we discuss the importance of consistency between Eulerian and Lagrangian fields in the framework of such Eulerian/Lagrangian hybrid formulations. Then, we validate the model in the case of continuous point source pollutant dispersion, under uniform wind and homogeneous turbulence. In these conditions, there is an analytical solution allowing a precise verification. We observe that in this case, the Lagrangian model discriminates well the two different near- and far-field diffusion regimes, which is not the case for an Eulerian model based on the eddy-viscosity hypothesis (Bahlali et al., 2018b).Finally, we work on the validation of the model on several experimental campaigns in real atmosphere, taking into account atmospheric thermal stratification and the presence of buildings. The first experimental program considered in our work has been conducted on the `SIRTA' site (Site Instrumental de Recherche par Télédétection Atmosphérique), in the southern suburb of Paris, and involves a stably stratified surface layer. The second campaign studied is the MUST (Mock Urban Setting Test) experiment. Conducted in the United States, in Utah's desert, this experiment aims at representing an idealized city, through several ranges of containers. Two cases are simulated and analyzed, respectively corresponding to neutral and stable atmospheric stratifications (Bahlali et al., 2018a)

Cette thèse porte sur l'étude de modèles et de méthodes numériques pour les écoulements diphasiques compressibles à faible nombre de Mach. Toutes les méthodes numériques développées dans cette étude sont basées sur une formulation de type volumes finis en maillages non structurés. La première partie de cette thèse propose une analyse du comportement des schémas décentrés de type Godunov dans la limite des faibles nombres de Mach. Nous expliquons de manière rigoureuse les raisons pour lesquelles ces schémas aboutissent à des approximations imprécises lorsque les écoulements sont très proches de l'incompressible. Nous développons alors des méthodes de préconditionnement adaptées qui permettent de retrouver de bonnes approximations. Ce premier travail complète un certain nombre de travaux récents sur l'analyse des schémas décentrés dans la limite des faibles nombres de Mach. Le deuxième point abordé dans cette thèse est un travail de modélisation où nous développons à partir d'un modèle bi-fluides à sept équations de type Baer-Nunziato, un modèle réduit à cinq équations pour les écoulements diphasiques. Ce travail présente une méthode originale de réduction de systèmes hyperboliques avec termes sources raides. Nous développons pour ce modèle un schéma numérique implicite et suivant la stratégie utilisée dans la première partie de cette thèse, une technique de préconditionnement adaptée aux écoulements à faible vitesse. Les expériences numériques réalisées montrent que ce modèle est bien adapté au calcul d'ondes de détonations ainsi qu'à la simulation d'interfaces entre fluides compressibles. Enfin la dernière partie de cette thèse porte sur l'étude d'un modèle à sept équations pour le calcul d'écoulements diphasiques à faible nombre de Mach. On développe des méthodes numériques implicites basées sur des solveurs de Riemann approchés, permettant de réduire les coûts de calcul pour ce type de régime.

La pollution de l'environnement par les pesticides reste une préoccupation sociale et environnementale importante. Lors de l'application, une partie des pesticides peut contaminer l'environnement (dérive). La diminution des pollutions repose sur la maîtrise de la taille et la vitesse des gouttes en sortie de buse. L'objectif principal de cette thèse est de développer une approche Eulérienne afin d'estimer les caractéristiques initiales des gouttes produites telles que la taille et la vitesse à la sortie de buse. Le modèle Eulérien d'atomisation considère l'écoulement diphasique d'un liquide et d'un gaz comme un écoulement turbulent d'un seul fluide avec une masse volumique variable, variant entre celle du gaz et celle du liquide pulvérisé. Une équation de transport pour la fraction massique liquide moyenne permet de décrire la dispersion du liquide dans la phase gazeuse. La turbulence est modélisée par une approche aux tensions de Reynolds (RSM) en résolvant les équations de transport de chacune des six composantes du tenseur de Reynolds. Par ailleurs, une équation de transport de l'interface volumique liquide/gaz est considérée. La production de l'interface volumique qui exprime la création des gouttes est fonction, d'une part, à grande échelle, du gradient de vitesse moyenne et, d'autre part, à petite échelle, de la turbulence. La destruction de l'interface volumique qui exprime la coalescence des gouttes, quant à elle, prend en compte la tension de surface qui s'oppose à la désintégration de la surface liquide. L'évolution spatiale des rayons et des vitesses des gouttes produites en sortie de buse est décrite en couplant l'équation de la fraction massique liquide moyenne et celle de l'interface volumique à celles de la conservation de la masse, de la quantité de mouvement et de la turbulence. L'étude est faite avec le code de calculs FLUENT V.12 en utilisant les Fonctions Définies par l'Utilisateur (UDF "User Defined Function") pour adapter le code à nos conditions d'injection très particulières puisque le rapport des masses volumiques liquide et gaz est de l'ordre de 800. Les résultats obtenus par le modèle indiquent la formation d'une nappe conique creuse constituée de grosses gouttelettes et la présence d'une zone de recirculation près de l'axe du spray constituée de gouttelettes plus petites, conformément aux expérimentations. La comparaison du Diamètre Moyen de Sauter calculé par le modèle et mesuré expérimentalement par l'Anémométrie Phase Doppler montre un bon accord.