Le contexte général de cette thèse s'inscrit dans le cadre de la modélisation eulérienne instationnaire des écoulements turbulents anisothermes gaz - particules. La modélisation de ces écoulements est cruciale pour de nombreuses applications industrielles et pour la prédiction de certains phénomènes naturels. Par exemple, la combustion diphasique dans les moteurs automobiles et aéronautiques est précédée par l'injection et la dispersion de carburant liquide dans la chambre de combustion. Les phénomènes mis en jeu exigent alors une prédiction locale tenant compte du caractère instationnaire de l'écoulement turbulent et de la présence de géométries complexes. De plus, de nombreuses études expérimentales et numériques récentes ont mis en évidence le rôle prépondérant de l'inertie des particules sur les mécanismes de dispersion et de concentration préférentielle en écoulement turbulent. Ceci rend donc indispensable la prise en compte de ces mécanismes dans la modélisation diphasique. Au cours de ce travail de thèse, une approche eulérienne locale et instantanée a été développée pour prédire les écoulements gaz-particules anisothermes et turbulents. Elle est basée sur l'approche statistique du Formalisme Eulérien Mésoscopique (MEF) introduite par Février et al. (JFM, 2005). Cette approche a été ici étendue aux variables thermiques pour la prise en compte du caractère anisotherme de l'écoulement. Cette approche a été ensuite utilisée dans le cadre de la méthode des moments (Kaufmann et al., JCP, 2008), et un système d'équations locales et instantanées pour la phase dispersée a été proposé. La modélisation au premier ordre exige la fermeture des moments de second ordre apparaissant dans les équations de la quantité de mouvement et de l'énergie. La proposition de telles relations constitutives fait l'objet d'une partie de la thèse. Afin de fournir une méthode capable de prédire le comportement local, instantané et anisotherme de la phase dispersée dans des configurations `a une échelle réaliste, les équations pour la phase dispersée ont été filtrées et une modélisation aux grandes échelles (LES) est effectuée. Cette modélisation étends, par la prise en compte des variables thermiques, le travail de Moreau et al. (FtaC, 2010) sur l'approche LES Euler-Euler en conditions isothermes. L'approche complète est enfin appliquée aux résultats de simulation numérique d'un jet plan turbulent gazeux froid, chargé en particules, dans une turbulence homogène isotrope chaude monophasique. / The aim of this thesis is to provide an Eulerian modeling for the dispersed phase interacting with unsteady non-isothermal turbulent flows. The modeling of these flows is crucial for several industrial applications and for predictions of natural events. Examples are the combustion chambers of areo engines where the combustion is preceded by the injection and dispersion of liquid fuel. The prediction of such phenomena involves a local modeling of the mixture for taking into account the unsteady behavior of the turbulent flow and the presence of complex geometries. Moreover, many experimental and numerical studies have recently highlighted the significant role of the particle inertia on the mechanisms of dispersion and preferential concentration. Accounting for such mechanisms is therefore essential for modeling the particle-laden turbulent flows. In this thesis, a local and instantaneous Eulerian approach able to describe and to predict the local behavior of inertial particles interacting with non-isothermal turbulent flows has been developed. It is based on the statistical approach known as Mesoscopic Eulerian formalism (MEF) introduced by Février et al. (JFM, 2005). The statistical approach has been extended to the thermal quantities in order to account for the non-isothermal conditions into the modeling. This formalism is then used in the framework of the moment approach (Kaufmann et al., JCP, 2008) and a system of local and instantaneous equations for the non-isothermal dispersed phase has been suggested. The first order modeling requires to close second-order moments appearing in momentum and energy equations. The proposal of such constitutive relations makes the object of a part of this study. In order to provide an Eulerian approach usable in real configurations at industrial scale, the equations of the dispersed phase are filtered and the approach developed in the framework of the Large-Eddy Simulations. From the work of Moreau et al. (FTaC, 2010), the Eulerian-Eulerian LES approach is then extended to non-isothermal conditions. The whole modeling is then a priori tested against numerical simulations of a cold planar turbulent particle-laden jet crossing a homogeneous isotropic decaying hot turbulence.

De nombreuses applications en océanographie et en météorologie côtière et opérationnelle nécessitent la mise en place de modèles locaux haute-résolution pour lesquels les interactions océan-atmosphère doivent être correctement représentées. Dans ce cas, la principale difficulté est de raccorder de manière consistante la solution des deux modèles à l'interface air-mer. Le travail présenté ici vise à adapter les méthodes de décomposition de domaine de type algorithmes de Schwarz pour résoudre ce problème d'une manière mathématiquement satisfaisante. Une difficulté importante provient notamment du caractère fondamentalement turbulent des couches limites près de l'interface air-mer. Pour aborder cet aspect, nous présentons dans un premier temps une synthèse des divers schémas de paramétrisation utiles. Puis nous étudions une formulation idéalisée du problème, sous la forme d'un couplage simplifié de deux équations de diffusion modélisant le mélange turbulent dans les couches limites, les coefficients de diffusion turbulente étant fournis par des paramétrisations usuelles. Afin d'assurer une convergence rapide de la méthode, nous recherchons des conditions de transmission optimisées. Nous déterminons celles-ci de manière analytique dans le cas de coefficients de diffusion constants et discontinus à l'interface, puis nous étendons ces résultats dans le cas de coefficients variables à l'aide d'une approche nouvelle. Dans une dernière partie, nous montrons comment les méthodes usuelles de couplage océan-atmosphère peuvent être décrites dans le formalisme des méthodes de Schwarz. Puis nous proposons une première application réaliste (formation et propagation d'un cyclone tropical) et présentons des résultats numériques préliminaires, obtenus avec une méthode non optimisée.

Nous présentons une approche alternative aux diagonalisations exactes de calculs en modèles en couches pour l'étude microscopique de la structure nucléaire. Nous avons tout d'abord minimisé l'énergie totale du système par la résolution des équations Hartree-Fock dans la base définie par l'espace de valence du modèle en couches. Nous avons finalisé la mise au point de ce programme par comparaison des résultats avec les solutions de problèmes solubles analytiquement comme ceux d'un hamiltonien SU(3) ou par comparaison avec des diagonalisations exactes du modèle en couches. Ensuite, nous avons mis au point des procédures pour l'obtention de ces mêmes calculs avec contraintes pour une description complète de la surface d'énergie potentielle d'un noyau donné en fonction des degrés de liberté Q20 et Q22. La restauration du moment angulaire par méthode de projections a permis d'obtenir les spectres des noyaux rotationnels dans le cas axial et dans le cas triaxial. Enfin, la méthode des coordonnées génératrices a permis le mélange de ces déterminants de Slater non orthogonaux. Parmi les différentes applications, nous avons pu décrire les noyaux déformés le long de la ligne N=Z autour du80Zr ou un nouvel îlot d'inversion à N=50. / We present an alternative approach to shell diagonalizations for microscopic description of nuclear structure. First we minimized the total system energy solving the Hartree-Fock equations within the shell model valence space. The results are compared with exact shell model diagonalization and an exact soluble SU(3) hamiltonian. Then, we developed procedures to obtain the same of type calculations with constrained conditions on the quadrupole degrees of freedom to obtain the full potential energy surfaces. The angular momentum restauration was obtained through projection method to generate rotionnal spectra of nuclei in both axial and triaxial cases. Finally the generate coordinate method was applied to mix several of these non-orthogonal Slater determinants. Among several applications we managed to describe deformed nuclei along the N=Z line around 80Zr or a new island of deformation at N=50. Finally, a first application in the superheavy region predicts a spherical gap for the Z=114, N=184 isotope

De nos jours, la simulation des écoulements diphasiques a de plus en plus d’importance dans les chambres de combustion aéronautiques en tant qu’un des éléments requis pour analyser et maîtriser le processus complet de combustion, afin d’améliorer la performance du moteur et de mieux prédire les émissions polluantes. Dans les applications industrielles, la modélisation du combustible liquide trouvé en aval de l’injecteur sous forme de brouillard de gouttes polydisperse, appelé spray, est de préférence faite à l’aide de méthodes Eulériennes. Ce choix s’explique par les avantages qu’offrent ces méthodes par rapport aux méthodes Lagrangiennes, notamment la convergence statistique intrinsèque, le couplage aisé avec la phase gazeuse ainsi que l’efficacité pour le calcul haute performance. Dans la présente thèse, on utilise une approche Eulérienne basée sur une fermeture au niveau cinétique de type distribution Gaussienne Anisotrope (AG). L’AG résout des moments de vitesse jusqu’au deuxième ordre et permet de capter les croisements des trajectoires (PTC) à petite échelle de manière statistique. Le système d’équations obtenu est hyperbolique, le problème est bien-posé et satisfait les conditions de réalisabilité. L’AG est comparé au modèle monocinétique (MK) d’ordre 1 en vitesse. Il est approprié pour la description des particules faiblement inertielles. Il mène à un système faiblement hyperbolique qui peut générer des singularités. Plusieurs schémas numériques, utilisés pour résoudre les systèmes hyperboliques et faible- ment hyperboliques, sont évalués. Ces schémas sont classifiés selon leur capacité à traiter les singularités naturellement présentes dans les modèles Eulériens, sans perdre l’ordre global de la méthode ni rompre les conditions de réalisabilité. L’AG est testé sur un champ turbulent 3D chargé de particules dans des simulations numériques directes. Le code ASPHODELE est utilisé pour la phase gazeuse et l’AG est implémenté dans le code MUSES3D pour le spray. Les résultats sont comparés aux de simulations Lagrangiennes de référence et aux modèle MK. L’AG est validé pour des gouttes modérément inertielles à travers des résultats qualitatifs et quantitatifs. Il s’avère prometteur pour les applications complexes comprenant des PTC à petite échelle. Finalement, l’AG est étendu à la simulation aux grandes échelles nécessaire dans les cas réels turbulents dans le domaine industriel en se basant sur un filtrage au niveau cinétique. Cette stratégie aide à garantir les conditions de réalisabilités. Des résultats préliminaires sont évalués en 2D pour tester la sensibilité des résultats LES sur les paramètres des modèles de fermetures de sous mailles. / In aeronautical combustion chambers, the ability to simulate two-phase flows gains increasing importance nowadays since it is one of the elements needed for the full understanding and prediction of the combustion process. This matter is motivated by the objective of improving the engine performance and better predicting the pollutant emissions. On the industrial scale, the description of the fuel spray found downstream of the injector is preferably done through Eulerian methods. This is due to the intrinsic statistical convergence of these methods, their natural coupling to the gas phase and their efficiency in terms of High Performance Computing compared to Lagrangian methods. In this thesis, the use of Kinetic-Based Moment Method with an Anisotropic Gaussian (AG) closure is investigated. By solving all velocity moments up to second order, this model reproduces statistically the main features of small scale Particles Trajectories Crossing (PTC). The resulting hyperbolic system of equations is mathematically well-posed and satisfies the realizability properties. This model is compared to the first order model in the KBMM hierarchy, the monokinetic model MK which is suitable of low inertia particles. The latter leads to a weakly hyperbolic system that can generate δ-shocks. Several schemes are compared for the resolution of the hyperbolic and weakly hyperbolic system of equations. These methods are assessed based on their ability to handle the naturally en- countered singularities due to the moment closures, especially without globally degenerating to lower order or violating the realizability constraints. The AG is evaluated for the Direct Numerical Simulation of 3D turbulent particle-laden flows by using ASPHODELE solver for the gas phase, and MUSES3D solver for the Eulerian spray in which the new model is implemented. The results are compared to the reference Lagrangian simulation as well as the MK results. Through the qualitative and quantitative results, the AG is found to be a predictive method for the description of moderately inertial particles and is a good candidate for complex simulations in realistic configurations where small scale PTC occurs. Finally, within the framework of industrial turbulence simulations a fully kinetic Large Eddy Simulation formalism is derived based on the AG model. This strategy of directly applying the filter on the kinetic level is helpful to devise realizability conditions. Preliminary results for the AG-LES model are evaluated in 2D, in order to investigate the sensitivity of the LES result on the subgrid closures.