Le contexte général de cette thèse est la simulation numérique de l'injection de carburant dans un moteur à combustion interne, afin d'améliorer son rendement et de limiter la production de polluants. De manière plus générale, ce travail s'applique à tout système industriel mettant en jeu un écoulement multiphasique constitué d'un carburant liquide injecté dans une chambre occupée initialement par du gaz, comme par exemple les moteurs automobiles ou aéronautiques, ou les turbomachines. Intrinsèquement, il est possible de simuler l'ensemble de l'écoulement avec les équations classiques de la dynamique des fluides sans avoir recours à des outils de modélisation supplémentaires liés au caractère diphasique. Mais, les tailles des structures générées pendant l'injection (gouttes de diamètre inférieur à 10 μm) conduisent à des temps de calculs prohibitifs pour une application industrielle. C'est pourquoi il est nécessaire d'introduire une modélisation diphasique. C'est dans ce contexte que deux régions sont formellement distinguées: le coeur liquide dense proche de l'injecteur, appelé écoulement à phases séparées, et le spray constitué d'une population de gouttes polydisperse (c'est-à-dire de tailles différentes) générées après le processus d'atomisation en aval de l'injecteur. Ce travail de thèse étudie les modèles Eulériens pour la description de spray évaporants et polydisperses, en vue d'applications industrielles. Ils représentent une alternative potentielle aux modèles Lagrangiens qui sont majoritairement utilisés en industrie mais présentant des inconvénients majeurs. Ainsi, le modèle multi-fluide est étudié dans un premier temps. Bien que prometteur, deux difficultés sont soulignées: le coût requis pour une description précise de la polydispersion, et son incapacité à décrire les croisements de gouttes (particle trajectory crossing, PTC, en anglais). La thèse propose des solutions à ces deux limitations. Ces solutions reposent chacune sur des méthodes de moments. Premièrement, le modèle appelé Eulerian Size Multi Size Moment (EMSM) permet de résoudre des sprays évaporants et polydisperses de manière bien plus efficace que le modèle multi-fluide. Des outils mathématiques sont utilisés pour fermer le système d'équations associé au modèle, et combinés à des schémas de types volumes finis appelés schémas cinétiques, afin de préserver la réalisabilité du vecteur de moments, pour le transport et l'évaporation. Une réponse à la seconde limitation est apportée avec le modèle appelé Eulerian Multi Velocity Moment (EMVM) basé sur le transport de moments en vitesse d'ordre élevé. Une distribution bimodale peut être localement reconstruite à partir des moments en utilisant une méthode de quadrature de moments (Quadrature Method of Moment, QMOM en anglais) en une ou plusieurs dimensions d'espace. De la même manière que précédemment, l'utilisation de schémas cinétiques permet de préserver la réalisabilité du vecteur de moment. De plus, une étude mathématique approfondie de la dynamique du système en une dimension d'espace en révèle toute la complexité et représente une étape indispensable en vue de l'élaboration de schémas de transport d'ordre élevé (supérieur ou égal à 2).Afin de les tester, ces deux modèles ainsi que les outils numériques associés sont implémentés dans MUSES3D, un code académique de simulation numérique directe (Direct Numerical Simulation DNS en anglais) dédié à l'évaluation des modèles de spray. Des résultats de grande qualité démontrent le potentiel des modèles. L'extension du modèle EMSM dans un contexte industriel est ensuite considérée, avec son implémentation dans IFP-C3D, un code résolvant des écoulements réactifs sur des maillages non structurés et mobiles (dû au mouvement du piston) dans un formalisme RANS (Reynolds Averaged Navier Stokes) en présence de sprays. Le formalisme ALE (Arbitrary Lagrangian Eulerian en anglais) est utilisé et le modèle EMSM réécrit dans ce formalisme afin de mener des calculs en maillage mobile. De plus, une étude numérique a permis d'étendre les propriétés de précision et de stabilité obtenues en maillage fixe. La robustesse du modèle EMSM est alors démontrée avec succès dans IFP-C3D sur un cas impliquant un mouvement de piston, ainsi que dans le cadre d'une comparaison avec le code MUSES3D. Enfin, des résultats très encourageants prouvent la faisabilité d'un calcul d'injection dans une chambre de combustion d'un spray polydisperse avec le modèle EMSM.

Le contexte général de cette thèse est la simulation numérique de l'injection de carburant dans un moteur à combustion interne, afin d'améliorer son rendement et de limiter la production de polluants. De manière plus générale, ce travail s'applique à tout système industriel mettant en jeu un écoulement multiphasique constitué d'un carburant liquide injecté dans une chambre occupée initialement par du gaz, comme par exemple les moteurs automobiles ou aéronautiques, ou les turbomachines. Intrinsèquement, il est possible de simuler l'ensemble de l'écoulement avec les équations classiques de la dynamique des fluides sans avoir recours à des outils de modélisation supplémentaires liés au caractère diphasique. Mais, les tailles des structures générées pendant l'injection (gouttes de diamètre inférieur à 10 μm) conduisent à des temps de calculs prohibitifs pour une application industrielle. C'est pourquoi il est nécessaire d'introduire une modélisation diphasique. C'est dans ce contexte que deux régions sont formellement distinguées: le coeur liquide dense proche de l'injecteur, appelé écoulement à phases séparées, et le spray constitué d'une population de gouttes polydisperse (c'est-à-dire de tailles différentes) générées après le processus d'atomisation en aval de l'injecteur. Ce travail de thèse étudie les modèles Eulériens pour la description de spray évaporants et polydisperses, en vue d'applications industrielles. Ils représentent une alternative potentielle aux modèles Lagrangiens qui sont majoritairement utilisés en industrie mais présentant des inconvénients majeurs. Ainsi, le modèle multi-fluide est étudié dans un premier temps. Bien que prometteur, deux difficultés sont soulignées: le coût requis pour une description précise de la polydispersion, et son incapacité à décrire les croisements de gouttes (particle trajectory crossing, PTC, en anglais). La thèse propose des solutions à ces deux limitations. Ces solutions reposent chacune sur des méthodes de moments. Premièrement, le modèle appelé Eulerian Size Multi Size Moment (EMSM) permet de résoudre des sprays évaporants et polydisperses de manière bien plus efficace que le modèle multi-fluide. Des outils mathématiques sont utilisés pour fermer le système d'équations associé au modèle, et combinés à des schémas de types volumes finis appelés schémas cinétiques, afin de préserver la réalisabilité du vecteur de moments, pour le transport et l'évaporation. Une réponse à la seconde limitation est apportée avec le modèle appelé Eulerian Multi Velocity Moment (EMVM) basé sur le transport de moments en vitesse d'ordre élevé. Une distribution bimodale peut être localement reconstruite à partir des moments en utilisant une méthode de quadrature de moments (Quadrature Method of Moment, QMOM en anglais) en une ou plusieurs dimensions d'espace. De la même manière que précédemment, l'utilisation de schémas cinétiques permet de préserver la réalisabilité du vecteur de moment. De plus, une étude mathématique approfondie de la dynamique du système en une dimension d'espace en révèle toute la complexité et représente une étape indispensable en vue de l'élaboration de schémas de transport d'ordre élevé (supérieur ou égal à 2).Afin de les tester, ces deux modèles ainsi que les outils numériques associés sont implémentés dans MUSES3D, un code académique de simulation numérique directe (Direct Numerical Simulation DNS en anglais) dédié à l'évaluation des modèles de spray. Des résultats de grande qualité démontrent le potentiel des modèles. L'extension du modèle EMSM dans un contexte industriel est ensuite considérée, avec son implémentation dans IFP-C3D, un code résolvant des écoulements réactifs sur des maillages non structurés et mobiles (dû au mouvement du piston) dans un formalisme RANS (Reynolds Averaged Navier Stokes) en présence de sprays. Le formalisme ALE (Arbitrary Lagrangian Eulerian en anglais) est utilisé et le modèle EMSM réécrit dans ce formalisme afin de mener des calculs en maillage mobile. De plus, une étude numérique a permis d'étendre les propriétés de précision et de stabilité obtenues en maillage fixe. La robustesse du modèle EMSM est alors démontrée avec succès dans IFP-C3D sur un cas impliquant un mouvement de piston, ainsi que dans le cadre d'une comparaison avec le code MUSES3D. Enfin, des résultats très encourageants prouvent la faisabilité d'un calcul d'injection dans une chambre de combustion d'un spray polydisperse avec le modèle EMSM.

Dans le cadre de la simulation de procédés de mise à forme par la méthode des éléments finis, le formalisme Arbitraire Lagrangien Eulérien (ALE) permet de découpler le mouvement du maillage et de la matière. Pour de très grandes déformations, la qualité du maillage peut être ainsi améliorée sans avoir recours à une procédure de remaillage complexe et coûteuse. Un second domaine d'application du formalisme ALE est la simulation de procédés stationnaires pour lesquels le maillage peut rester fixe dans la direction de l'écoulement de matière. Ce type de maillage quasi eulérien permet de diminuer le nombre d'éléments finis du modèle numérique par rapport à une simulation lagrangienne classique. En conséquence, le temps de calcul est également réduit. Bien que le formalisme ALE ne soit pas nouveau, il est rarement utilisé en pratique. D'une part, les techniques de repositionnement de noeuds tridimensionnelles ne sont pas évidentes à mettre en oeuvre et, d'autre part, il n'existe pas de schéma de transfert de données précis adapté aux éléments finis à plus d'un point de Gauss tels que ceux utilisés en dynamique implicite. Cette thèse tente de combler ces deux lacunes: une méthode de repositionnement de noeud sur les surfaces libres du solide est présentée. Elle est très robuste et permet de conserver la forme des surfaces malgré le mouvement arbitraire du maillage. Concernant le transfert des données, un schéma de convection 3D précis au second ordre et utilisable pour des éléments finis à plusieurs points de Gauss est déduit d'un opérateur de projection. L'efficacité et la généralité de l'algorithme ALE qui en résulte sont ensuite démontrées sur une série d'applications de plus en plus complexes: impact d'une barre de Taylor, traction d'une barre d'Hopkinson, double extrusion, laminage, planage et profilage à froid. Dans chaque cas, les résultats ALE sont comparés avec des résultats lagrangiens traditionnels ainsi que des mesures expérimentales lorsque celles-ci sont disponibles.

Dans le cadre de la simulation de procédés de mise à forme par la méthode des éléments finis, le formalisme Arbitraire Lagrangien Eulérien (ALE) permet de découpler le mouvement du maillage et de la matière. Pour de très grandes déformations, la qualité du maillage peut être ainsi améliorée sans avoir recours à une procédure de remaillage complexe et coûteuse. Un second domaine d'application du formalisme ALE est la simulation de procédés stationnaires pour lesquels le maillage peut rester fixe dans la direction de l'écoulement de matière. Ce type de maillage quasi eulérien permet de diminuer le nombre d'éléments finis du modèle numérique par rapport à une simulation lagrangienne classique. En conséquence, le temps de calcul est également réduit. Bien que le formalisme ALE ne soit pas nouveau, il est rarement utilisé en pratique. D'une part, les techniques de repositionnement de noeuds tridimensionnelles ne sont pas évidentes à mettre en oeuvre et, d'autre part, il n'existe pas de schéma de transfert de données précis adapté aux éléments finis à plus d'un point de Gauss tels que ceux utilisés en dynamique implicite. Cette thèse tente de combler ces deux lacunes: une méthode de repositionnement de noeud sur les surfaces libres du solide est présentée. Elle est très robuste et permet de conserver la forme des surfaces malgré le mouvement arbitraire du maillage. Concernant le transfert des données, un schéma de convection 3D précis au second ordre et utilisable pour des éléments finis à plusieurs points de Gauss est déduit d'un opérateur de projection. L'efficacité et la généralité de l'algorithme ALE qui en résulte sont ensuite démontrées sur une série d'applications de plus en plus complexes: impact d'une barre de Taylor, traction d'une barre d'Hopkinson, double extrusion, laminage, planage et profilage à froid. Dans chaque cas, les résultats ALE sont comparés avec des résultats lagrangiens traditionnels ainsi que des mesures expérimentales lorsque celles-ci sont disponibles.

Le contexte général de cette thèse est la simulation numérique de l’injection de carburant dans un moteur à combustion interne, afin d’améliorer son rendement et de limiter la production de polluants. Intrinsèquement, il est possible de simuler l’ensemble de l’écoulement avec les équations classiques de la dynamique des fluides sans avoir recours à des outils de modélisation supplémentaires liés au caractère diphasique. Mais, les tailles des structures générées pendant l’injection (gouttes de diamètre < à 10 μm) conduisent à des temps de calculs prohibitifs pour une application industrielle. C’est pourquoi il est nécessaire d’introduire une modélisation diphasique. C’est dans ce contexte que deux régions sont formellement distinguées: le cœur liquide dense proche de l’injecteur, appelé écoulement à phases séparées, et le spray constitué d’une population de gouttes polydisperse générées après le processus d’atomisation en aval de l’injecteur. Ce travail de thèse étudie les modèles Eulériens pour la description de spray évaporants et polydisperses, en vue d’applications industrielles. Ils représentent une alternative potentielle aux modèles Lagrangiens qui sont majoritairement utilisés en industrie mais présentant des inconvénients majeurs. Ainsi, le modèle multi-fluide est étudié dans un premier temps. Bien que prometteur, deux difficultés sont soulignées: le coût requis pour une description précise de la polydispersion, et son incapacité à décrire les croisements de gouttes (particle trajectory crossing, PTC). La thèse propose des solutions à ces deux limitations. Elles reposent sur des méthodes de moments. Premièrement, le modèle appelé Eulerian Size Multi Size Moment (EMSM) permet de résoudre des sprays évaporants et polydisperses de manière bien plus efficace que le modèle multi-fluide. Des outils mathématiques sont utilisés pour fermer le système d’équations associé au modèle, et combinés à des schémas de types volumes finis appelés schémas cinétiques, afin de préserver la réalisabilité du vecteur de moments, pour le transport et l’évaporation. Une réponse à la seconde limitation est apportée avec le modèle appelé Eulerian Multi Velocity Moment (EMVM) basé sur le transport de moments en vitesse d’ordre élevé. Une distribution bimodale peut être localement reconstruite à partir des moments en utilisant une méthode de quadrature de moments ( QMOM) en une ou plusieurs dimensions d’espace. De la même manière, l’utilisation de schémas cinétiques permet de préserver la réalisabilité du vecteur de moment. De plus, une étude mathématique approfondie de la dynamique du système en une dimension d’espace en révèle toute la complexité et représente une étape indispensable en vue de l’élaboration de schémas de transport d’ordre élevé (supérieur ou égal à 2).Afin de les tester, ces deux modèles ainsi que les outils numériques associés sont implémentés dans MUSES3D, un code académique de simulation numérique directe (DNS) dédié à l’évaluation des modèles de spray. Des résultats de grande qualité démontrent le potentiel des modèles. L’extension du modèle EMSM dans un contexte industriel est ensuite considérée, avec son implémentation dans IFP-C3D, un code résolvant des écoulements réactifs sur des maillages non structurés et mobiles dans un formalisme RANS (Reynolds Averaged Navier Stokes) en présence de sprays. Le formalisme ALE (Arbitrary Lagrangian Eulerian) est utilisé et le modèle EMSM réécrit dans ce formalisme afin de mener des calculs en maillage mobile. De plus, une étude numérique a permis d’étendre les propriétés de précision et de stabilité obtenues en maillage fixe. La robustesse du modèle EMSM est alors démontrée avec succès dans IFP-C3D sur un cas impliquant un mouvement de piston, ainsi que dans le cadre d’une comparaison avec le code MUSES3D. Enfin, des résultats très encourageants prouvent la faisabilité d’un calcul d’injection dans une chambre de combustion d’un spray polydisperse avec le modèle EMSM. / The general context of the PhD is the simulation of fuel injection in an internal combustion engine, in order to improve its thermal and ecological efficiency. This work more generally concerns any industrial device involving a multiphase flow made of liquid fuel injected in a chamber filled with gaz: automotive or aircraft engines, or turbo machines. In and of itself, it is possible to simulate this flow without any modeling. However the small structures created during injection (droplets of diameter until 10 μm or less) lead to a prohibitive computational cost for any industrial application. Therefore modeling is necessary. In this context, two areas are formally distinguished: the dense liquid core close to the injector called separate-phase flow, and the spray made of a polydisperse droplet population (i.e. droplets with different sizes) generated after the atomization processes downstream of the injector. This PhD work investigates Eulerian models for the description of polydisperse evaporating sprays, for industrial computations. They represent a potential alternative to Lagrangian models, widely used at present, yet suffering from major drawbacks. Thus, the Multi-Fluid model is assessed. Although it is very promising, two difficulties are highlighted: its cost for a precise description of polydispersity, and its inability to describe particle trajectory crossing (PTC). Solutions to these two limitations are considered. Both rely on high order moment methods. First, the Eulerian Multi Size Moment (EMSM) proposes a much more efficient resolution of polydisperse evaporating sprays than the Multi-Fluid model does. Mathematical tools are used to close the model and combined with original finite volume kinetic-based schemes in order to preserve the moment-set integrity, for evaporation and advection. An answer to the second limitation is provided with the Eulerian Multi Velocity Moment (EMVM) based on high order velocity moments. A bimodal velocity distribution can be locally reconstructed for the moments using the quadrature method of moments (QMOM), in one or multi-dimensions. Here also, finite volume kinetic-based schemes are studied in order to preserve the moment set integrity. Moreover, a mathematical study of the one-dimensional dynamic system highlights its peculiarity and constitutes a necessary basis for the design of high order numerical schemes. In order to assess them, both the models and their numerical tools are implemented in the MUSES3D code, an academic DNS solver that provides a framework devoted to spray method evaluation. Achievements of the EMSM and the EMVM models are presented. The extension of the EMSM model to an industrial context is then considered, with its implementation in the IFP-C3D code, a 3D unstructured reactive flow solver with spray. In order to perform computations within a moving domain (due to the piston movement) the Arbitrary Lagrangian Eulerian (ALE) formalism is used. A numerical study has been achieved, in order to extent to this formalism the properties of accuracy and stability of the EMSM model, which already induces strong stability condition in an Eulerian approach. The robustness of the EMSM model in the IFP-C3D code has been successfully demonstrated on a case involving a moving piston, and also on a comparison with the MUSES3D code. Moreover, very encouraging results demonstrate the feasibility of the EMSM model for spray injection.