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

Théorèmes limites pour les sommes de Birkhoff de fonctions d'intégrale nulle en théorie ergodique en mesure infinie

Ce travail est consacré à certaines classes de systèmes dynamiques ergodiques, munis d'une mesure invariante infinie, telles que des applications de l'intervalle avec un point fixe neutre ou des marches aléatoires. Le comportement asymptotique des sommes de Birkhoff d'observables d'intégrale non nulle est assez bien connu, pour peu que le système ait une certaine forme d'hyperbolicité. Une situation particulièrement intéressante est celle des tours au-dessus d'une application Gibbs-Markov. Nous cherchons dans ce contexte à étudier le cas d'observables d'intégrale nulle. Nous obtenons ainsi une forme de théorème central limite pour des systèmes dynamiques munis d'une mesure infinie. Après avoir introduit l'ensemble des notions nécessaires, nous adaptons des résultats de E. Csáki et A. Földes sur les marches aléatoires au cas des applications Gibbs-Markov. Les théorèmes d'indépendance asymptotique qui en découlent forment le cœur de cette thèse, et permettent de démontrer un théorème central limite généralisé. Quelques variations sur l'énoncé de ce théorème sont obtenues. Ensuite, nous abordons les processus en temps continu, tels que des semi-flots et des flots. Un premier travail consiste à étudier les propriété en temps grand du temps de premier retour et du temps local pour des extensions de systèmes dynamiques, ce qui se fait par des méthodes spectrales. Enfin, par réductions successives, nous pouvons obtenir une version du théorème central limite pour des flots périodiques, et en particulier le flot géodésique sur le fibré tangent unitaire de certaines variétés périodiques hyperboliques.This work is focused on some classes of ergodic dynamical systems endowed with an infinite invariant measure, such as transformations of the interval with a neutral fixed point or random walks. The asymptotic behavior of the Birkhoff sums of observables with a non-zero integral is well known, as long as the system shows some kind of hyperbolicity. The towers over a Gibbs-Markov map are especially interesting. In this context, we aim to study the case of observables whose integral is zero. We get the equivalent of a central limit theorem for some dynamical systems endowed with an infinite measure. After we introduce the necessary definitions, we adapt some results by E. Csáki and A. Földes on random walks to the case of Gibbs-Markov maps. We derive a theorem on the asymptotic independence of Birhoff sums, which is the core of this thesis, and from this point we work out a generalised central limit theorem. We also prove a few variations on this generalised central limit theorem. Then, we study dynamical systems in continuous time, such as semi-flows and flows. We first work on the asymptotic properties of the first return time and the local time for extensions of dynamical systems; this is done by spectral methods. Finally, step by step, we extend our generalised central limit theorem to cover some periodic flows, and in particular the geodesic flow on the unitary tangent bundle of some hyperbolic periodic manifolds.RENNES1-Bibl. électronique (352382106) / SudocSudocFranceF