Le modèle multi-fluide permet de décrire par une approche Eulérienne les sprays polydispersés et apparaît donc comme une méthode indiquée pour les applications de combustion diphasique. Sa pertinence pour la simulation à l'échelle d'applications industrielles est évaluée dans ce travail, par sa mise en oeuvre dans des configurations bi-dimensionnelle et tri-dimensionnelle plus représentatives de ce type de simulations. Cette évaluation couple une étude de faisabilité en terme de coût de calcul avec une analyse de la précision obtenue, par des comparaisons avec les résultats de méthodes de références pour la description des sprays. Afin de définir une telle référence, une hiérarchisation des modèles de spray est proposée dans ce travail, soulignant les niveaux de modélisation associée aux diverses méthodes. Une première configuration d'écoulements tourbillonnaires est utilisée pour caractériser la méthode multi-fluide. L'étude de la structure mathématique du système de lois de conservation permet d'analyser la formation de singularités et de fournir les outils permettant d'évaluer leur impact sur la modélisation. Cette étude permet également de dériver un schéma numérique robuste et efficace pour des configurations bi- et tri-dimensionnelle. La description des dynamiques de gouttes conditionnées par la taille est évaluée dans ces configurations tourbillonnaires au moyen de comparaisons quantitatives, sur des champs instantanés, où le multi-fluide est confronté à une méthode Lagrangienne, ainsi qu'à des résultats expérimentaux. Afin d'évaluer le comportement de la méthode multi-fluide dans des configurations plus représentatives des problématiques industrielles, le solveur MUSES3D est développé, permettant, entre autres, une évaluation fine des méthodes de résolution des sprays. Une implémentation originale de la méthode multi-fluide, conciliant généricité et efficacité pour le calcul parallèle, est réalisée. Le couplage de ce solveur avec le code ASPHODELE, développé au CORIA, permet d'effectuer une évaluation opérationnelle des approches Euler/Lagrange et Euler/Euler pour la description des écoulements diphasiques à inclusions dispersées. Finalement, le comportement de la méthode multi-fluide dans des jets bi-dimensionnels et dans une turbulence homogène isotrope tri-dimensionnelle permet de montrer sa précision pour la description de la dynamique de sprays évaporant dans des configurations plus complexes. La résolution de la polydispersion du spray permet de décrire précisément la fraction massique de combustible en phase vapeur, un élément clé pour les applications de combustion. De plus, l'efficacité du calcul parallèle par décomposition de domaine avec la méthode multi-fluide permet d'envisager son utilisation à l'échelle d'applications industrielles.

Dans divers domaines de la physique, certains phénomènes sont modélisés par des systèmes hyperboliques non-conservatifs. En particulier, dans le domaine de la physique des plasmas, dont l'un des champs d'application majeur est la Fusion par Confinement Inertiel, le système d'Euler bi-température, modélisant les phénomènes de transport de particules chargées, en est un exemple. La difficulté de l'étude de ces systèmes réside dans la présence de termes non-conservatifs, qui empêchent la définition classique des solutions faibles. Pour parvenir à une définition de ce type de solutions, on a recours à l'emploi de systèmes cinétiques sous-jacents. Dans ce manuscrit, on s'intéresse à l'étude numérique de ces systèmes cinétiques pour la résolution du système d'Euler bi-température.Ce manuscrit se divise en deux parties. La première partie contient l'étude numérique du système d'Euler bi-température. Dans un premier chapitre, on résout numériquement les équations en dimension 1 d'espace par le biais d'un système sous-jacent issu de la physique des plasmas: le système de Vlasov-BGK-Ampère. On présente une méthode numérique préservant l'asymptotique pour ce système sous-jacent et on montre, par des simulations numériques, que le schéma limite obtenu donne des résultats consistants avec Euler bi-température. Dans un second chapitre, on résout le même modèle en dimension 2 d'espace par un système sous-jacent de type BGK discret. On démontre une inégalité d'entropie pour les solutions issues du modèle sous-jacent, ainsi qu'une inégalité discrète de dissipation d'entropie pour le schéma.Dans la deuxième partie de ce manuscrit, on s'intéresse au développement de méthodes numériques pour quelques modèles cinétiques. On considère ici le cas des écoulements raréfiés de mélanges de gaz, dans l'optique d'une application aux cas des plasmas. Premièrement, on présente un schéma cinétique adaptatif et dynamique en vitesse pour les gaz inertes. Par l'emploi de lois de conservation discrètes, la solution est approchée sur un ensemble de vitesses discrètes local et dynamique. Dans un second temps, on propose une extension de cette méthode visant à améliorer les performances de celle-ci. Puis, ces deux versions de la méthode sont comparées à la méthode classique sur grille fixe uniforme sur une série de cas tests.Enfin, dans le dernier chapitre, on propose une méthode numérique pour la résolution d'une extension de ces équations, prenant en compte la présence de réactions chimiques au sein du mélange. Le contexte considéré est celui des réactions chimiques bi-moléculaires réversibles lentes. La méthode proposée, de type implicite-explicite, est linéaire, stable et conservative. / In various domains of physics, several phenomena can be modeled via the use of nonconservative hyperbolic systems. In particular, in plasma physics, in the process of developping and understanding the phenomena leading to Inertial Confinement Fusion, the bi-temperature Euler sytem can be used to model particle transport phenomena in a plasma. The difficulty of the mathematical study of such systems dwells in the presence of so-called non-conservative products, which prevent the classical definition of weak solutions via distribution theory. To attempt to define these quantities, it is useful to supplement the hyperbolic system with an underlying kinetic model. In this work, the objective is the numerical study of such kinetic systems in order to solve the bi-temperature Euler system.This manuscript is split in two parts. The first one contains the study of the bi-temperature Euler system. In the first chapter, this system in dimension 1 is solved by the use of an underlying kinetic model sprung from plasma physics: the Vlasov-BGK-Ampère system. An asymptotic-preserving numerical method is introduced, and it is shown that the scheme obtained in the limit is consistant with a scheme for teh bi-temperature Euler system. In the following chapter, the same hyperbolic model in dimension 2 is studied, this time via a discrete-BGK type underlying model. An entropy inequality is proved for solutions coming from the kinetic model, as well as a discrete entropy dissipation inequality.In the second part of the manuscript, we are interested in the development of numerical schemes for gas mixture rarefied flows. Firstly, an adaptive kinetic scheme is introduced for inert gas mixtures. By the use of discrete conservation laws, the solution is approximated on a set of discrete velocities that depends on space, time and species. Secondly, an extension of the method is proposed in order to improve the efficiency of the first method. Finally, the two methods are compared to the classical fixed grid method on a series of test cases.In the last chapter, a numerical method is proposed for rarefied flows of reacting mixtures. The setting considered is the case of slow bimolecular reversible chemical reactions. The method introduced is an explicit-implicit treatment of the relaxation operator, which is shown to be stable, linear and conservative.

Le modèle multi-fluide permet de décrire par une approche Eulérienne les sprays polydispersés et apparaît donc comme une méthode indiquée pour les applications de combustion diphasique. Sa pertinence pour la simulation à l’échelle d’applications industrielles est évaluée dans ce travail, par sa mise en oeuvre dans des configurations bi-dimensionnelle et tri-dimensionnelle plus représentatives de ce type de simulations. Cette évaluation couple une étude de faisabilité en terme de coût de calcul avec une analyse de la précision obtenue, par des comparaisons avec les résultats de méthodes de références pour la description des sprays. Afin de définir une telle référence, une hiérarchisation des modèles de spray est proposée dans ce travail, soulignant les niveaux de modélisation associée aux diverses méthodes. Une première configuration d’écoulements tourbillonnaires est utilisée pour caractériser la méthode multi-fluide. L’étude de la structure mathématique du système de lois de conservation permet d’analyser la formation de singularités et de fournir les outils permettant d’évaluer leur impact sur la modélisation. Cette étude permet également de dériver un schéma numérique robuste et efficace pour des configurations bi- et tri-dimensionnelle. La description des dynamiques de gouttes conditionnées par la taille est évaluée dans ces configurations tourbillonnaires au moyen de comparaisons quantitatives, sur des champs instantanés, où le multi-fluide est confronté à une méthode Lagrangienne, ainsi qu’à des résultats expérimentaux. Afin d’évaluer le comportement de la méthode multi-fluide dans des configurations plus représentatives des problématiques industrielles, le solveur MUSES3D est développé, permettant, entre autres, une évaluation fine des méthodes de résolution des sprays. Une implémentation originale de la méthode multi-fluide, conciliant généricité et efficacité pour le calcul parallèle, est réalisée. Le couplage de ce solveur avec le code ASPHODELE, développé au CORIA, permet d’effectuer une évaluation opérationnelle des approches Euler/Lagrange et Euler/Euler pour la description des écoulements diphasiques à inclusions dispersées. Finalement, le comportement de la méthode multi-fluide dans des jets bi-dimensionnels et dans une turbulence homogène isotrope tri-dimensionnelle permet de montrer sa précision pour la description de la dynamique de sprays évaporant dans des configurations plus complexes. La résolution de la polydispersion du spray permet de décrire précisément la fraction massique de combustible en phase vapeur, un élément clé pour les applications de combustion. De plus, l’efficacité du calcul parallèle par décomposition de domaine avec la méthode multi-fluide permet d’envisager son utilisation à l’échelle d’applications industrielles. / The multi-fluid model, providing a Eulerian description of polydisperse sprays, appears as an interesting method for two-phase combustion applications. Its relevance as a numerical tool for industrial device simulations is evaluated in this work. This evaluation assesses the feasibility of multi-fluid simulations in terms of computational cost and analyzes their precision through comparisons with reference methods for spray resolution. In order to define such a reference, the link between the available methods for spray resolution is provided, highlighting their corresponding level of modeling. A first framework of 2-D vortical flows is used to assess the mathematical structure of the multi-fluid model governing system of equations. The link between the mathematical peculiarities and the physical modeling is provided, and a robust numerical scheme efficient for 2-D/3-D configurations is designed. This framework is also used to evaluate the multi-fluid description of evaporating spray sizeconditioned dynamics through quantitative, time-resolved, comparisons with a Lagrangian reference and with experimental data. In order to assess the multi-fluid efficiency in configurations more representative of industrial devices, a numerical solver is designed, providing a framework devoted to spray method evaluation. An original implementation of the multifluid method, combining genericity and efficiency in a parallel framework, is achieved. The coupling with a Eulerian/Lagrangian solver for dispersed two-phase flows, developed at CORIA, is conducted. It allows a precise evaluation of Euler/Lagrange versus Euler/Euler approaches, in terms of precision and computational cost. Finally, the behavior of the multi-fluid model is assessed in 2D-jets and 3-D Homogeneous Isotropic Turbulence. It illustrates the ability of the method to capture evaporating spray dynamics in more complex configurations. The method is shown to describe accurately the fuel vapor mass fraction, a key issue for combustion applications. Furthermore, the method is shown to be efficient in domain decomposition parallel computing framework, a key issue for simulations at the scale of industrial devices.