Cette thèse s’inscrit dans le domaine de la simulation numérique appliquée à la propulsion. Le moteur à détonation rotative (RDE) fait partie des candidats susceptibles de remplacer nos actuels moyens de propulsion grâce à l’augmentation du rendement thermodynamique du moteur. Pour conserver l’avantage de la détonation, l’injecteur doit fournir un mélange dont la qualité doit être la meilleure possible tout en limitant les pertes de pression totale. La présente étude porte sur le développement et l’optimisation numérique d’un injecteur adapté au fonctionnement d’un RDE. L’injection d’hydrogène et d’oxygène gazeux en rapport stoechiométrique est considérée pour une utilisation en propulsion fusée. Le premier objectif est de proposer un concept réaliste d’injecteur permettant de maximiser le mélange des ergols. Le second objectif est de réaliser des études du mélange dans la chambre par des simulations LES (Large Eddy Simulation). Le troisième objectif est de simuler la propagation d’une détonation rotative (RD) alimentée par différents injecteurs en régimes prémélangé et séparé. Deux éléments d’injection sont mis en concurrence. Le premier utilise le principe de jets semi-impactants de H2 et de O2. Le deuxième représente une configuration améliorée. Les simulations de RD avec les deux injecteurs donnent des résultats similaires lorsque l’injection est prémélangée. La part du mélange injecté perdu par déflagration est de 30% et la vitesse de propagation de la RD est proche de la vitesse théorique CJ. Pour les injections séparées de H2 et O2, l’injecteur amélioré permet de conserver un bon niveau de mélange dans la chambre, contrairement à l’injecteur à semi-impact qui produit une forte stratification des ergols dans la chambre. En conséquence, la vitesse de propagation de la RD est proche de la vitesse CJ avec l’injecteur amélioré et limitée à 80% de la vitesse CJ avec l’injecteur à semi-impact. / This thesis pertains to the domain of numerical simulation for propulsion applications. The rotating detonation engine (RDE) appears to be a good candidate to replace our current means of propulsion thanks to the increase of the thermodynamic efficiency. To preserve the advantage given by the detonation mode, the injector must provide the best possible mixing of the propellants together with acceptable total pressure losses. This numerical study deals with developing and optimizing an injector adapted to the operation of a RDE. Injection of gaseous H2 and O2 at stoichiometric ratio is considered to be suitable for rocket propulsion application. The first goal is to propose an efficient injector design so that the mixing between the propellants is maximized. The second goal is to perform simulations of the mixing process in the chamber by LES (Large Eddy Simulation) computations. The third goal consists in computing the propagation of a rotating detonation (RD) fed by different injectors in premixed and separate regimes. This study allows the comparison of two injection elements. The first one uses the principle of semi-impinging jets of H2 and O2. The second one represents an improved configuration. RD simulations with both injectors provide similar results when premixed injection is considered. The part of the injected mixture that burns by deflagration is 30% and the detonation velocity remains close the theoretical CJ velocity. In the regime of separate injection of H2 and O2, the improved injector enables to keep a high mixing efficiency in the chamber whereas the semi-impingement injector produces a strong stratification of the propellants in the chamber. As a consequence, the detonation velocity is close to the CJ velocity with the improved injector and limited to 80% of the CJ velocity with the semi-impingement injector.

De nombreuses recherches sont actuellement menées afin de réduire les émissions polluantesdes aéronefs. Le contrôle actif des écoulements aérodynamiques est une piste envisagéepour répondre à ces enjeux. Parmi les technologies de contrôle en développement, lestechnologies plasma offrent plusieurs avantages, dont la compacité, la simplicité de mise enoeuvre et la réactivité. Ce travail de thèse a été consacré à l’étude d’un actionneur plasmade type jet synthétique. Il se présente sous la forme d’une petite cavité insérée en paroiet reliée à l’extérieur par une tuyère. Un arc électrique est généré dans la cavité, ce quientraîne une augmentation de la pression de l’air dans celle-ci. Par conséquent, un jet estproduit à la sortie de la tuyère, et celui-ci peut interagir avec l’écoulement extérieur. A lasuite de cette phase d’éjection, de l’air extérieur est naturellement aspiré par la cavité, cequi permet au processus d’être répété à des fréquences pouvant atteindre plusieurs kilohertz.L’objectif de ce travail de thèse est de construire un modèle numérique capable dereproduire ces phénomènes physiques. Pour cela, un calcul aérodynamique de type LargeEddy Simulation est mis en oeuvre. L’action du plasma d’arc est prise en compte au traversde termes sources dans l’équation de l’énergie. Ces derniers sont notamment calculés grâceà l’hypothèse d’équilibre thermodynamique local dans le plasma. De plus, l’augmentationde la température dans la partie solide de l’actionneur est simulée lorsque celui-ci est opéréà haute-fréquence. Les résultats du modèle numérique sont comparés à des mesures de vitesseeffectuées lorsque l’actionneur fonctionne dans un environnement extérieur au reposet lorsque celui-ci interagit avec une couche limite. / Nowadays, many studies are conducted in order to decrease greenhouse gases and noiseemissions from aircrafts. Active aerodynamic flow control is a way considered to meet thesechallenges. Among developed technologies of control, plasma actuators offer several advantages,including compactness, easy implementation and fast response. This thesis is devotedto the study of a plasma synthetic jet actuator which comes in the form of a small cavityinserted in wall and connected to the environment through a nozzle. An electrical arc isgenerated in the cavity, resulting in an increase of the cavity air pressure. Consequently,a jet is produced at the nozzle exhaust that can interact with the external flow. After thisejection phase, an aspiration phase naturally occurs, allowing the process to be repeatedat a frequency that can reach several kilohertz. This thesis aims at building a numericalmodel able to reproduce these physical phenomena. To do so, a Large Eddy Simulationis performed. The action of the plasma arc is taken into account through source terms inthe energy equation. In particular, these source terms are calculated using the assumptionof local thermodynamic equilibrium in the plasma. In addition, the heating of the actuatormaterials is simulated when it is operated at high frequency. The results from the numericalmodel are compared with velocity measurements, carried out when the actuator works ina quiescent environment and when it interacts with a boundary layer.

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. / 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.

" Le sujet de cette thèse traitera du Liban entre les deux retraits israélien et syrien, 2000 - 2005.En effet, cette récente période s'inscrit comme l'une des plus critiques dans l'histoire contemporaine du Liban. La guerre civile déclenchée en avril 1975 est officiellement terminée en octobre 1990 par une " paix syrienne " imposée aux belligérants Libanais, soutenue par la Ligue Arabe et la communauté internationale. Cette période de paix entre 1990 et 2005, fut d'une part une période de reconstruction et de développement du Liban. Mais d'autre part, ces quinze années témoignent d'un profond malaise et mécontentement inscrits dans une situation globale apparemment paisible mais en fait marquée de fragilités politiques, sociales et communautaires profondes.Deux évènements majeurs - les retraits : israélien en 2000 et syrien en 2005 - viennent bouleverser ce statut quo libanais et déclenchent le processus inévitable d'une nouvelle " démocratisation " du Liban dont les conséquences ne cessent de se ressentir et de rebondir à ce jour. Entre ces deux retraits, un troisième " évènement choc ", le 11 septembre 2001 semble au regard de certains Libanais la perspective d'une nouvelle politique américaine et européenne dans la région.Cette période sera marquée au Liban, par les événements d'août 2000, arrestations arbitraires à l'encontre des opposants au régime prosyrien, la fermeture forcée de la chaîne de télévision libanaise anti syrienne MTV, les discours virulents des prélats maronites de l'église du Liban et du patriarche du Liban contre le régime en place, la nouvelle politique du leader druze Walid Joumblatt, l'éloignement du Premier Ministre libanais Rafic Hariri de la politique syrienne, et le début de la création d'une opposition multiconfessionnelle contre la reconduction du président libanais prosyrien Emile Lahoud, le vote de la résolution 1559 au conseil de sécurité de l'ONU, et finalement le " 11 septembre libanais " ou l'assassinat de Rafic Hariri, en février 2005 et la " révolution du Cèdre " qui en suit.

Celui-ci a pour objectif de valider l'ensemble des modèles physiques utilisés dans un code de simulation numérique pour simuler un écoulement de type annulaire dispersé en conduite rencontré lors d'une procédure de lavage à eau utilisé dans les raffineries. Pour ce faire une banque de données expérimentale est mise en place sur des configurations représentatives de celles utilisées en condition industrielle. La géométrie retenue comporte une zone horizontale d'injection rectiligne avec un injecteur central, suivi d'un coude à 90° situé dans un plan vertical. Différentes conditions expérimentales permettent d'étudier l'influence de la vitesse du gaz, de la condition d'injection du brouillard et de la pression sur les différents processus physiques. Ces résultats comprenant des visualisations du brouillard et du film pariétale, des mesures de taille et de distribution de gouttes,des mesures de débit et d'épaisseur de film, sont analysés pour faire ressortir les principaux mécanismes d'interaction entre le gaz et la phase dispersée, le gaz et le film liquide pariétal et la phase dispersée et le film pariétal. En parallèle, des premières simulations, avec une approche RANS, sont réalisées avec le code CEDRE de l'ONERA et les résultats sont confrontés aux mesures. / This work has been realised within a CIFRE contract with TOTAL. Its aim was to validate all the physical models used in a computation, which simulates an annular dispersed flow through a pipe used in a water washing process in refinery plants. That is why, a whole set of data has been gathered using experimental boundary conditions which are representative to those used in industrial configurations. The geometry is made of a horizontal pipe with a centred nozzle followed by a 90º elbow in the vertical plane. Several experimental boundary conditions enable one to study the influence of the gas velocity, the type of the spray injection and the pressure on the different physical phenomena. These results including spray and liquid film visualisations, droplets distribution and size measurements as well as liquid film thickness and mass flow measurements were analysed in order to extract the main interaction mechanism between the gas and the dispersed phase, the gas and the liquid film, and the dispersed phase and the annular liquid film. Meanwhile, simulations using a RANS approach were realized with the ONERA code named CEDRE and its results were compared to the gathered measurements.

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. / 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.

« Le sujet de cette thèse traitera du Liban entre les deux retraits israélien et syrien, 2000 – 2005.En effet, cette récente période s’inscrit comme l’une des plus critiques dans l’histoire contemporaine du Liban. La guerre civile déclenchée en avril 1975 est officiellement terminée en octobre 1990 par une « paix syrienne » imposée aux belligérants Libanais, soutenue par la Ligue Arabe et la communauté internationale. Cette période de paix entre 1990 et 2005, fut d’une part une période de reconstruction et de développement du Liban. Mais d’autre part, ces quinze années témoignent d’un profond malaise et mécontentement inscrits dans une situation globale apparemment paisible mais en fait marquée de fragilités politiques, sociales et communautaires profondes.Deux évènements majeurs - les retraits : israélien en 2000 et syrien en 2005 - viennent bouleverser ce statut quo libanais et déclenchent le processus inévitable d’une nouvelle « démocratisation » du Liban dont les conséquences ne cessent de se ressentir et de rebondir à ce jour. Entre ces deux retraits, un troisième « évènement choc », le 11 septembre 2001 semble au regard de certains Libanais la perspective d’une nouvelle politique américaine et européenne dans la région.Cette période sera marquée au Liban, par les événements d’août 2000, arrestations arbitraires à l’encontre des opposants au régime prosyrien, la fermeture forcée de la chaîne de télévision libanaise anti syrienne MTV, les discours virulents des prélats maronites de l’église du Liban et du patriarche du Liban contre le régime en place, la nouvelle politique du leader druze Walid Joumblatt, l’éloignement du Premier Ministre libanais Rafic Hariri de la politique syrienne, et le début de la création d’une opposition multiconfessionnelle contre la reconduction du président libanais prosyrien Emile Lahoud, le vote de la résolution 1559 au conseil de sécurité de l’ONU, et finalement le « 11 septembre libanais » ou l’assassinat de Rafic Hariri, en février 2005 et la « révolution du Cèdre » qui en suit. / «The subject of that thesis will deal with Lebanon between both Israeli and Syrian withdrawals, on 2000 - 2005. Indeed, this recent period joins as one of the most critical in the contemporary history of Lebanon. The civil war started in April, 1975 is officially ended in October, 1990 with a " Syrian peace " compulsory for the Lebanese belligerent parties, supported by the Arab League and the international community. This peace period between 1990 and 2005, was on one hand a period of reconstruction and development of the Lebanon. But on the other hand, these fifteen years testify of a profound faintness and a dissatisfaction registered in an apparently peaceful global situation but in fact marked with deep political, social and community fragilities:Two major events - withdrawals: Israeli in 2000 and Syrian in 2005 come to upset this Lebanese status quo and activate the inevitable process of a new "democratization" of Lebanon of which the consequences do not stop feeling the effects and bouncing to this day. Between these two withdrawals, the third " event shock ", September 11th, 2001 seems with regard to certain Lebanese the perspective of a new American and European policy in the region.This period will be marked in Lebanon, by the events of August 2000, arbitrary arrests against the opponents of the proSyrian regime, the forced closure of the Lebanese anti Syrian television channel MTV, the virulent speeches of the Maronite prelates of the church and the patriarch of Lebanon against the regime in position, the new policy of the Druze leader Walid Joumblatt, the estrangement of the Lebanese Prime Minister Rafic Hariri of the Syrian politics, and the beginning of the creation of a multidenominational opposition against the renewal of the proSyrian Lebanese president Emile Lahoud, the vote of the resolution 1559 in the United Nations Security Council, and finally on the " Lebanese September 11th " or the murder of Rafic Hariri, in February, 2005 and the " revolution of the Cedar ".