Robust numerical schemes for Eulerian spray DNS and LES in two-phase turbulent flows

International audienceLarge Eddy Simulation (LES) and Direct numerical Simulation (DNS) of polydisperse evaporating sprays with Eulerian models are very promising tools for high performance computing of combustion applications. They are able to describe the turbulent dispersion and evaporation and properly predict the combustion regimes. However, the spray system of conservation equations has a convective part which is either similar to gas dynamics Euler equations with a real gas type state law or to the pressureless gas dynamics (PGD), depending on the local flow regime and droplet Stokes number; so, they usually involve singularities due to model closure assumptions and require dedicated numerical schemes. Besides, it is desirable to cope with exactly zero droplet density in some zones of the flow, especially near the injection zone, where droplets are injected in only some spatial locations. Even if the issue has been successfully tackled in de Chaisemartin (2009); Fréret et al. (2010) in the framework of PGD with the use of accurate kinetic schemes, it cannot be directly extended to general gas dynamics. The purpose of the present contribution is to introduce a new generation of numerical methods based on relaxation schemes which are able to treat both PGD and general gas dynamics, as well as to cope in a robust manner with vacuum zones and natural singularities of the resulting system of conservation equations. The proposed hybrid relaxation scheme and algorithms are validated through comparisons with analytical solutions and other numerical strategies on 1D and 2D configurations. They exhibit a very robust behavior and are a very promising candidate for more complex applications since they provide solutions to key numerical issues of the actual Eulerian spray DNS and LES models

A multi-Gaussian quadrature method of moments for simulating high Stokes number turbulent two-phase flows

With the great increase in computational resources, Large Eddy Simulation (LES) of industrial configurations is now an efficient and tractable tool. Numerous applications involve a liquid or solid disperse phase carried by a gaseous flow field (eg, fuel injection in automotive or aeronautical engines, fluidized beds, and alumina particles in rocket boosters). To simulate this kind of flow, one may resort to a Number Density Function (NDF), which satisfies a kinetic equation

Two-size moment Eulerian multi-fluid method describing the statistical trajectory crossing: modeling and numerical scheme

International audienceHigh fidelity modeling and simulation of turbulent dispersed two-phase flows is still a major challenge for many applications. Eulerian approaches are well suited for high performance computations of such flows. Recently, hybrid Eulerian methods that combine the multi-fluid approach-where the size is discretized-and the moment method were developed. On the one hand, in order to capture efficiently the size-polydispersion, two size moments were used on each interval of the size discretization (Two Size Moment method). On the other hand, the Anisotropic Gaussian (AG) velocity closure has been introduced as a relevant model to describe velocity dispersion occurring when the particles from the disperse phase have a significant inertia compared to the time scales of the flow, leading to particle trajectory crossings. The purpose of this contribution is to develop a model able to describe size and velocity dispersion, coupling the two-size moment Eulerian multi-fluid method and the anisotropic velocity closure. Adapted numerical schemes based on a relaxation method are provided. This new model (AG-TSM) is then evaluated on various test cases relevant to rocket propulsion and two-phase combustion

Simulation Numérique Directe des sprays dilués anisothermes avec le Formalisme Eulérien Mésoscopique

Le contexte général de cette thèse est la Simulation Numérique Directe des écoulements diphasiques dilués anisothermes. Un accent particulier est mis sur la détermination précise de la dispersion des particules et du transfert de chaleur entre la phase porteuse et dispersée. Cette dernière est décrite à l’aide d’une approche Eulérienne aux moments : le Formalisme Eulérien Mésoscopique (FEM) [41, 123], récemment étendu aux écoulements anisothermes [78]. Le principal objectif de ce travail est de déterminer si ce formalisme est capable de prendre en compte de manière précise l’inertie dynamique et thermique des particules dans un écoulement turbulent, et particulièrement dans une configuration avec un gradient moyen. Le code de calcul utilisé est AVBP. La simulation numérique d’un spray dilué avec une approche Eulerienne soulève des questions supplémentaires sur les méthodes numériques et les modèles employés. Ainsi, les méthodes numériques spécifiques aux écoulements diphasiques implémentées dans AVBP [69, 103, 109] ont été testées et revisitées. L’objectif est de proposer une stratégie numérique précise et robuste qui résiste aux forts gradients de fraction volumique de particule provoqués par la concentration préférentielle [132], tout en limitant la diffusion numérique. Ces stratégies numériques sont comparées sur une série de cas tests de complexité croissante et des diagnostics pertinents sont proposés. Par exemple, les dissipations dues à la physique et au numérique sont extraites des simulations et quantifiées. Le cas test du tourbillon en deux dimensions chargé en particules est suggéré comme une configuration simple pour mettre en évidence l’impact de l’inertie des particules sur leur champ de concentration et pour discriminer les stratégies numériques. Une solution analytique est aussi proposée pour ce cas dans la limite des faibles nombres de Stokes. Finalement, la stratégie numérique qui couple le schéma centré d’ordre élevé TTGC et une technique de stabilisation, aussi appelée viscosité artificielle, est celle qui fournit les meilleurs résultats en terme de précision et de robustesse. Les paramètres de viscosité artificielle (c'est-à-dire les senseurs) doivent néanmoins être bien choisis. Ensuite, la question des modèles nécessaires pour d´écrire correctement la dispersion des particules dans une configuration avec un gradient moyen est abordée. Pour ce faire, un des modèles RUM (appelé AXISY-C), proposé par Masi [78] et implémenté dans AVBP par Sierra [120], est validé avec succès dans deux configurations: un jet plan diphasique anisotherme 2D et 3D. Contrairement aux anciens modèles RUM, les principales statistiques de la phase dispersée sont désormais bien prédites au centre et aux bords du jet. Finalement, l’impact de l’inertie thermique des particules sur leur température est étudié. Les résultats montrent un effet important de cette inertie sur les statistiques mettant en évidence la nécessité pour les approches numériques de prendre en compte ce phénomène. Ainsi, l’extension du FEM aux écoulements anisothermes, c’est-à-dire les flux de chaleur RUM (notés RUM HF), est implémentée dans AVBP. L’impact des RUM HF sur les statistiques de température des particules est ensuite évalué sur les configurations des jets 2D et 3D. Les champs Eulériens sont comparés à des solutions Lagrangiennes de référence calculées par B. Leveugle au CORIA et par E. Masi à l’IMFT pour les jets 2D et 3D, respectivement. Les résultats montrent que les RUM HF améliorent la prédiction des fluctuations de température mésoscopique, et dans une moindre mesure la température moyenne des particules en fonction de la configuration. Les statistiques Lagrangiennes sont retrouvées lorsque les RUM HF sont pris en compte alors que les résultats sont dégradés dans le cas contraire. ABSTRACT : This work addresses the Direct Numerical Simulation of non-isothermal turbulent flows laden with solid particles in the dilute regime. The focus is set on the accurate prediction of heat transfer between phases and of particles dispersion. The dispersed phase is described by an Eulerian approach : the Mesoscopic Eulerian Formalism [41, 123], recently extended to non-isothermal flows [78]. The main objective of this work is to assess the ability of this formalism to accurately account for both dynamic and thermal inertia of particles in turbulent sheared flows. The CFD code used in this work is AVBP. The numerical simulation of dilute sprays with an Eulerian approach calls for specific modelling and raises additional numerical issues. First, the numerical methods implemented in AVBP for two-phase flows [69, 103, 109] were tested and revisited. The objective was to propose an accurate and robust numerical strategy that withstands the steep gradients of particle volume fraction due to preferential concentration [132] with a limited numerical diffusion. These numerical strategies have been tested on a series of test cases of increasing complexity and relevant diagnostics were proposed. In particular, the two-dimensional vortex laden with solid particles was suggested as a simple configuration to illustrate the effect of particle inertia on their concentration profile and to test numerical strategies. An analytical solution was also derived in the limit of small inertia. Moreover, dissipations due to numerics and to physical effects were explicitly extracted and quantified. Eventually, the numerical strategy coupling the highorder centered scheme TTGC with a stabilization technique –the so called artificial viscosity– proved to be the most accurate and robust alternative in AVBP if an adequate set-up is used (i.e. sensors). Then, the issue of the accurate prediction of particle dispersion in configurations with a mean shear was adressed. One of the RUM model (denoted AXISY-C), proposed by Masi [78] and implemented by Sierra [120], was successfully validated in a two-dimensional and a three-dimensional non-isothermal jet laden with solid particles. Contrary to the former RUM models [63, 103], the main statistics of the dispersed phase were recovered at both the center and the edges of the jet. Finally, the impact of the thermal inertia of particles on their temperature statistics has been investigated. The results showed a strong dependency of these statistics to thermal inertia, pinpointing the necessity of the numerical approaches to account for this phenomenon. Therefore, the extension of the MEF to non isothermal conditions, i.e. the RUM heat fluxes, has been implemented in AVBP. The impact of the RUM HF terms on the temperature statistics was evaluated in both configurations of 2D and 3D jets. Eulerian solutions were compared with Lagrangian reference computations carried out by B. Leveugle at CORIA and by E. Masi at IMFT for the 2D and 3D jets, respectively. Results showed a strong positive impact of the RUM HF on the fluctuations of mesoscopic temperature, and to a lesser extent on the mean mesoscopic temperature depending of the configuration. Neglecting the RUM HF leads to erroneous results whereas the Lagrangian statistics are recovered when they are accounted for

Direct Numerical Simulation of non-isothermal dilute sprays using the Mesoscopic Eulerian Formalism

This work addresses the Direct Numerical Simulation of non-isothermal turbulent flows laden with solid particles in the dilute regime. The focus is set on the accurate prediction of heat transfer between phases and of particles dispersion. The dispersed phase is described by an Eulerian approach : the Mesoscopic Eulerian Formalism [41, 123], recently extended to non-isothermal flows [78]. The main objective of this work is to assess the ability of this formalism to accurately account for both dynamic and thermal inertia of particles in turbulent sheared flows. The CFD code used in this work is AVBP. The numerical simulation of dilute sprays with an Eulerian approach calls for specific modelling and raises additional numerical issues. First, the numerical methods implemented in AVBP for two-phase flows [69, 103, 109] were tested and revisited. The objective was to propose an accurate and robust numerical strategy that withstands the steep gradients of particle volume fraction due to preferential concentration [132] with a limited numerical diffusion. These numerical strategies have been tested on a series of test cases of increasing complexity and relevant diagnostics were proposed. In particular, the two-dimensional vortex laden with solid particles was suggested as a simple configuration to illustrate the effect of particle inertia on their concentration profile and to test numerical strategies. An analytical solution was also derived in the limit of small inertia. Moreover, dissipations due to numerics and to physical effects were explicitly extracted and quantified. Eventually, the numerical strategy coupling the highorder centered scheme TTGC with a stabilization technique –the so called artificial viscosity– proved to be the most accurate and robust alternative in AVBP if an adequate set-up is used (i.e. sensors). Then, the issue of the accurate prediction of particle dispersion in configurations with a mean shear was adressed. One of the RUM model (denoted AXISY-C), proposed by Masi [78] and implemented by Sierra [120], was successfully validated in a two-dimensional and a three-dimensional non-isothermal jet laden with solid particles. Contrary to the former RUM models [63, 103], the main statistics of the dispersed phase were recovered at both the center and the edges of the jet. Finally, the impact of the thermal inertia of particles on their temperature statistics has been investigated. The results showed a strong dependency of these statistics to thermal inertia, pinpointing the necessity of the numerical approaches to account for this phenomenon. Therefore, the extension of the MEF to non isothermal conditions, i.e. the RUM heat fluxes, has been implemented in AVBP. The impact of the RUM HF terms on the temperature statistics was evaluated in both configurations of 2D and 3D jets. Eulerian solutions were compared with Lagrangian reference computations carried out by B. Leveugle at CORIA and by E. Masi at IMFT for the 2D and 3D jets, respectively. Results showed a strong positive impact of the RUM HF on the fluctuations of mesoscopic temperature, and to a lesser extent on the mean mesoscopic temperature depending of the configuration. Neglecting the RUM HF leads to erroneous results whereas the Lagrangian statistics are recovered when they are accounted for

Quadrature-based moment methods for polydisperse multiphase flow modeling

Polydisperse multiphase flows arise in many applications, and thus there has been considerable interest in the development of numerical methods to find solutions to the kinetic equations used to model such flows. However, the direct numerical solution of the kinetic equations is intractable for most applications due to the large number of independent variables. A useful alternative is to reformulate the problem in terms of the moments of the number density function (NDF), yet the resulting moment transport equations are not closed for flows away from the equilibrium limit. To attain closure, Quadrature-based moment methods (QBMM) is proposed. QBMM reconstruct NDF from a set of moments, which is the key step and named moment-inversion algorithm, then use NDF to close the moment transport equations. By different function approximation, two types of moment-inversion algorithm can be determined. The first type is approximate NDF by Dirac delta function. Quadrature method of moments (QMOM) has been proposed to handle three-dimensional problem for this type of approximation. However, the positivity of NDF cannot be guaranteed by QMOM. Therefore, a novel moment-inversion algorithm, based on 1-D adaptive quadrature of conditional velocity moments, is introduced and shown to yield NDF which is always promise positivity. This conditional quadrature method of moments (CQMOM) can be used to compute exact N-point quadratures for multi-valued solutions, and provides optimal approximations of continuous distributions. In order to control numerical errors arising in volume averaging and spatial transport, an adaptive 1-D quadrature algorithm is formulated for use with CQMOM. The use of adaptive CQMOM in the context of QBMM for the solution of kinetic equations is illustrated by applying it to problems involving particle trajectory crossing, Riemann problem, and granular flow. The drawback of Dirac delta function approximation has two fold, one is when large numbers of nodes are required to achieve the desired accuracy, the moment-inversion problem can become ill-conditioned. Another is value of NDF cannot be provided in QMOM or CQMOM when it is necessary in some applications. To conquer these disadvantages, a new generation of quadrature algorithm is introduced that uses an explicit form for the distribution function. This extended quadrature method of moments (EQMOM) approximates the distribution function by a sum of classical weight functions, which allow unclosed source terms to be computed with great accuracy by increasing the number of quadrature nodes independent of the number of transported moments. EQMOM is used to solve a population balance equations with evaporation, aggregation and breakage terms and compare the results with analytical solutions. This novel quadrature methods EQMOM is then applied to simulate bubbly flow. Bubble-column reactors are widely used in the biological, chemical and petrochemical industries. The accurate design of these reactors depends largely on the complex fluid dynamics of gas-liquid two-phase flows that still remains inadequately understood. Modeling of the fluid dynamics of gas-liquid bubble columns is therefore a challenging task. The Euler-Euler method is widely used in industry to simulate bubble columns. However, accurately predicting polydisperse bubbly flow is a nontrivial task due to the complexity of the bubble number density function, which can involve up to four internal coordinates including size and velocity. To describe polydisperse bubbles, a joint velocity-mass NDF for bubbles is adopted. QBMM is applied to solve the kinetic equation of the joint velocity-mass NDF using EQMOM. It is coupled with an incompressible Navier-Stokes solver for the liquid phase. In this model, transport equations for the joint velocity-mass moments are derived from a kinetic equation for the NDF and closure is attained using a monokinetic NDF, which is valid in the limit of small bubble Stokes number. The pure moments of mass are used to reconstruct mass NDF with EQMOM, while the joint moments determine the conditional velocity. Forces including buoyancy, drag, virtual-mass and lift are accounted for. The injection with a narrow bubble size distribution cases are used to validate the model with experimental data from the literature. Other cases with a continuous size distribution injection show the ability of new model to handle polydisperse bubbles. Results demonstrate that the onset of segregation is sensitive to the bubble size distribution and, thus, an accurate solution for the size-dependent fluxes is required when simulating polydisperse bubbly flows

Simulation aux grandes échelles d'écoulements diphasiques turbulents à phase liquide dispersée

Les écoulements diphasiques turbulents sont présents dans de nombreux systèmes industriels (moteur à piston, turbines à gaz, moteurs fusée...). La compréhension fine de telles configurations s'avèrent de nos jours nécessaire pour limiter notamment les émissions de polluants et de gaz à effet de serre, et la consommation des énergies fossiles. Nous nous intéressons ici à la simulation aux grandes échelles des écoulements diphasiques turbulents, permettant de capturer une large partie du spectre de la turbulence, et ainsi être capable de prédire des phénomènes instables ou transitoires. La phase dispersée est ici modélisée par une approche eulérienne, en raison de ses avantages dans le contexte du calcul haute performance. Le travail de cette thèse a consisté à étendre le formalisme eulérien existant dans le code AVBP à la simulation de sprays polydisperses dans des écoulements turbulents. Pour cela, le Formalisme Eulérien Mésoscopique (FEM) a été couplé à une approche Multi-fluide. Cette nouvelle approche, intitulée Formalisme Eulérien Mésoscopique Multi-fluide (FEMM), a été évaluée sur des cas simples canoniques, permettant de bien caractériser le comportement autant en terme de dynamique turbulente que d'effets polydisperses. Les stratégies numériques disponibles dans le code de calcul AVBP sont aussi analysées, afin d'en cerner les limites pour la simulation eulérienne d'une phase liquide. Ce nouveau formalisme est finalement appliqué à la configuration aéronautique MERCATO, pour laquelle on dispose de résultats numériques obtenus avec d'autres approches (FEM et approche lagrangienne), et de résultats expérimentaux. Un accord satisfaisant avec l'expérience est montré pour toutes les approches, même si le FEM, monodisperse, obtient de moins bon résultats en terme de fluctuations. D'autres résultats expérimentaux s'avèrent nécessaires pour évaluer les approches et déterminer quelle est la plus prédictive pour cette configuration, notamment concernant la fraction massique de kerosene, autant en phase liquide qu'en phase gazeuse. ABSTRACT : Turbulent two-phase flows are encountered in several industrial devices (piston engine, gas turbine, rocket engine...). A fine understanding of such configurations is mandatory to face problems of pollutant emissions, greenhouse gas, and fossil fuel rarefaction. The Large Eddy Simulation seems to be a good candidate. This kind of simulation captures a wide part of turbulence spectrum, and thus allows to predict instabilities and transient phenomena. The dispersed phase is simulated using an Eulerian approach, which seems to be more suitable than lagrangian methods for High Performance Computing. The present work consists in the extension to polydisperse flows of the existing eulerian formalism in the AVBP code. The Mesoscopic Eulerian Formalism (MEF) is coupled with the Multifluid approach. This new formalism, called Multifluid Mesoscopic Eulerian Formalism, is evaluated on simple test cases, showing the ability of such approach to capture turbulent and polydisperse effects. Numerical strategies available in AVBP are also evaluated, in order to emphasize on their limiting aspects for the eulerian simulation of a dispersed phase. The new formalism is finally applied to the simulation of the aeronautical configuration called MERCATO. Several experimental results are available, as well as numerical results using FEM and lagrangian approach. Results show a good agreement between experiments and numerical results, even if FEM results are worse concerning the fluctuations. New experimental results are necessary to determine which is the best approach, especially in terms of liquid and gas kerosene mass fraction

Modélisation et étude de l’évaporation et de la combustion de gouttes dans les moteurs à propergol solide par une approche eulérienne Multi-Fluide

The addition of a significant mass fraction of aluminum particle in the propellant of Solid Rocket Motors improves performance through an increase of the temperature in the combustion chamber. The distributed combustion of aluminum droplets in a portion of the chamber yields a massive amount of disperse aluminum oxide residues with a large size spectrum, called a polydisperse spray, in the entire volume. The spray can have a significant impact on the motor behavior and in particular on the onset/damping of instability. When dealing with aeroacoustical and thermoacoustical instabilities, the faithful prediction of the interactions between the gaseous phase and the spray is a determining step for understanding the physical mechanisms and for future solid rocket motor optimization. In such a harsh environment, experimental measurements have a hard time providing detailed explanation of the physical mechanisms and one has to resort to numerical simulation. For such a purpose, the distributed combustion zone and thermal profile therein, the heat generated by the combustion of the dispersed droplets and the large size distribution of the aluminum oxide residues and its coupling with he gaseous phase hydrodynamic and acoustic fields have to be accurately reproduced through a proper level of modeling and a high fidelity simulation including a precise resolution of size polydispersity, which is a key parameter.In this contribution, we choose a kinetic approach for the description of polydisperse sprays. The Williams-Boltzmann Equation is used to model the disperse phase and we derive a fully Eulerian approach through moment methods. The Multi-Fluid (MF) methods naturally treat droplet size evolution through phenomena such as evaporation and coalescence. These methods rely on the conservation of size moments on fixed intervals called sections and yield systems of conservation laws for a set of "fluids" of droplet of various sizes, which is strongly coupled with the gas phase via source terms. We derive a new optimal and flexible Two Size Moment MF method based on a family of polynomial reconstruction functions to describe the size distribution in the sections, which is second order accurate and particularly efficient at describing accurately the evolution of the size distribution with a moderate number of sections. An original work is also conducted in order to extend this approach to two-component droplets. For size moment MF methods, realizability of the moments is a crucial issue. Thus, we have developed innovative schemes for integrating source terms in moment conservation equations describing transport in phase space. This method enables the use of more representative aluminum droplet combustion models, and leads to more advanced studies of the distributed combustion zone. Moreover, for unsteady two-phase flow simulations, we have developed a robust and accurate coupling strategy between phases that are modeled by a fully Eulerian approach based on operator splitting in order to treat such spatial and temporal very multi-scale problems with reasonable computational time. All the proposed developments have been carried out following two criteria : 1- an attractive cost/accuracy ratio for industrial simulations in the context of high fidelity simulations 2- a preservation of industrial code legacy. Verification of the models and methods have been conducted first using an in-house reseach code and then in the context of a two-phase acoustic study thus emphasizing the relevance of the splitting technique to capture accurately spray-acoustic interactions.En propulsion solide, l'ajout de particules d'aluminium dans le propergol améliore de façon significative les performances du moteur grâce à une augmentation sensible de la température de chambre. La présence de gouttes d'aluminium et de résidus d'alumine de différentes tailles et en quantité importante a un impact notoire sur le fonctionnement du moteur. Dans cette optique, nous souhaitons obtenir une meilleure prévision de la stabilité de fonctionnement en cas de déclenchement d'instabilités d'origine aéroacoustique ou thermoacoustique. Nous visons des calculs plus précis de l'étendue de la zone de combustion, de la chaleur dégagée par la combustion distribuée des gouttes et de la distribution en taille des résidus. Nos efforts ont porté sur la modélisation des échanges entre la phase gazeuse et cette phase dispersée composée de gouttes de nature et de taille très diverses. Le paramètre taille pilotant la dynamique du spray et le couplage avec le gaz, le suivi précis des changements de taille est un enjeu majeur.Dans cette contribution, nous avons choisi une approche cinétique pour la description des sprays polydisperses. L'équation cinétique de Williams-Boltzmann utilisée pour suivre l'évolution des propriétés du spray est résolue par une approche eulérienne. Les méthodes Multi-Fluide (MF) traitent naturellement les changements de taille tels que l'évaporation et la coalescence. Ces méthodes reposent sur une intégration continue de la variable taille sur des intervalles fixes appelés sections sur lesquels nous pouvons dériver des systèmes d'équations de conservation. Chaque système est vu comme un fluide qui est en couplage fort avec la phase gazeuse via des termes sources.Nous avons travaillé sur une méthode MF à deux moments en taille basée sur une famille de fonctions de forme polynomiale pour reconstruire la distribution en taille au sein des sections. Cette approche d'ordre deux en temps et en espace s'avère performante car elle décrit avec précision l'évolution de la distribution avec un nombre modéré de sections. Un travail original a été mené afin d'étendre l'approche MF à des gouttes bicomposants. Cette méthode ouvre la voie à des modèles de combustion des gouttes d'aluminium plus représentatifs. Dans le contexte des simulations instationnaires, nous avons porté une attention particulière à l'emploi d'une stratégie numérique robuste et précise pour le couplage entre les phases modélisées par une approche Euler-Euler. Nous montrons qu'une méthode de splitting séparant le traitement du transport des phases gazeuse/dispersée de celui des termes sources est particulièrement adaptée pour la résolution d'un problème multi-échelle spatial et temporel. Dans la mesure où les conditions de réalisabilité sur les moments en taille des méthodes MF ne sont pas garanties avec des méthodes d'intégration traditionnelles, nous avons développé des schémas innovants pour l'intégration des termes sources. Les travaux proposés dans cette contribution répond à deux exigences : 1- un ratio coût/précision attractif pour des simulations industrielles 2- une facilité d'implémentation des méthodes et une modularité assurant la pérennisation des codes industriels. Ces développements ont d'abord été vérifiés à l'aide d'un code ad hoc ; des cas test d'étude d'acoustique diphasique linéaire ont notamment souligné la pertinence de la technique de splitting pour restituer avec précision les interactions spray-acoustique. Les nouvelles méthodes ont ensuite été implémentées et validées au sein du code multi-physique CEDRE développé à l'ONERA. Des calculs de propulsion solide sur des configurations moteur réalistes ont finalement mis en évidence le niveau de maturité atteint par les méthodes eulériennes pour décrire avec fidélité la dynamique des sprays polydisperses. Les résultats de ces simulations ont mis en avant la sensibilité des niveaux d'instabilités en fonction de la distribution en taille des gouttes d'aluminium et des résidus

Couplage entre modèles diphasiques à « phases séparées » et à « phase dispersée » pour la simulation de l’atomisation primaire en combustion cryotechnique

Two-phase flows play a significant role for the proper functioning of cryogenic liquid-propellant rocketengines, such as those that equip the launchers of the Ariane family. Since the experimental investigationof such propulsion devices is complex and expensive, developing numerical tools able to accuratelysimulate their functioning, is a crucial but nonetheless ambitious objective. The major difficulty is due tothe multiscale nature of the problem, as a result of which there is currently no numerical approach ableto perfectly describe all the liquid scales on its own. Based on this observation the work presented in thisthesis aims at setting up a coupling strategy between models well-adapted to each two-phase flowtopology, in the framework of the ONERA’s multiphysics CEDRE software. The approach adoptedprecisely consists in coupling a 4-equation diffuse interface model for the separated phases and aeulerian kinetic model for the dispersed phase, thus making it possible to describe primary atomization.Besides, the harsh conditions within cryogenic rocket engines, where large temperature, velocity anddensity gradients are encountered, severely challenge the robustness of numerical methods. A newmultislope MUSCL method for general unstructured meshes is thus developed in order to improve therobustness and accuracy of space discretization schemes. The whole coupling strategy is finally appliedto the numerical simulation of the ONERA’s Mascotte test bench for cryogenic combustion research.Les écoulements diphasiques jouent un rôle prépondérant dans les moteurs-fusées à ergols liquides cryogéniques, équipant par exemple les lanceurs de la famille Ariane. L'étude expérimentale de tels engins propulsifs étant complexe et onéreuse, disposer d'outils numériques à même de simuler fidèlement leur fonctionnement se révèle être un objectif aussi important qu'ambitieux. La difficulté majeure réside dans le caractère fortement multi-échelles du problème, si bien qu’aucune approche numérique existante n'est capable à elle seule de décrire parfaitement l'ensemble des échelles liquides. Partant de ce constat, les travaux présentés dans cette thèse visent à mettre en place une stratégie de couplage entre des modèles bien adaptés aux différentes topologies d'écoulement diphasique, et ce dans le cadre de la plateforme logicielle multi-physique CEDRE développée par l'ONERA. La démarche adoptée consiste précisément à coupler un modèle à interface diffuse de type ``4 équations'' pour les zones à phases séparées, et un modèle cinétique eulérien pour la phase dispersée, rendant ainsi possible la description de l’atomisation primaire. Par ailleurs, les conditions sévères qui règnent dans les moteurs cryotechniques, où de forts gradients de température, vitesse et densité sont rencontrés, mettent à l'épreuve la robustesse des méthodes numériques. Une nouvelle méthode MUSCL multipente pour maillages non structurés généraux a ainsi été développée, permettant d’améliorer la robustesse et la précision des schémas de discrétisation spatiale. L’ensemble de la stratégie de couplage est finalement appliquée à la simulation du banc Mascotte de l'ONERA pour la combustion cryotechnique

Generalized averaged Gaussian quadrature and applications

A simple numerical method for constructing the optimal generalized averaged Gaussian quadrature formulas will be presented. These formulas exist in many cases in which real positive GaussKronrod formulas do not exist, and can be used as an adequate alternative in order to estimate the error of a Gaussian rule. We also investigate the conditions under which the optimal averaged Gaussian quadrature formulas and their truncated variants are internal