Les Moteurs à Propergol Solide sont le siège d'instabilités de combustion qui se manifestent par des fluctuations de la pression interne et, en conséquence, de la poussée induite. Les phénomènes oscillatoires observés résultent d'un couplage entre l'acoustique de la chambre de combustion et une ou plusieurs instabilités hydrodynamiques de l'écoulement interne.Dans cette thèse, on s'intéresse à l'étude de l'écoulement induit par injection pariétale en géométrie axisymétrique. Plus particulièrement, on se concentre sur l'interaction qui peut avoir lieu entre cet écoulement et une structure : les Protections Thermiques de Face.Une étude expérimentale reposant sur le montage gaz froid VALDO de l'ONERA, modélisation représentative des MPS P230, a permis d'examiner l'influence d'obstacles, rigides et en élastomère, sur les instationnarités de l'écoulement. Le caractère instable de l'écoulement induit par injection pariétale est alors confirmé et le rôle primordial joué par les PTF sur les fréquences émergeant au sein de l'écoulement a été identifié.L'analyse de stabilité linéaire de l'écoulement, conduite en approche globale, permet l'étude de sa dynamique en tant qu'amplificateur de bruit. Les m¶mécanismes mis en jeu dans l'instabilité de l'écoulement induit par injection pariétale ainsi que son interaction complexe avec la couche de cisaillement issue de l'obstacle sont analysés. La réponse de l'écoulement à un forçage harmonique est alors examinée à la lumière des résultats expérimentaux.Cette étude a été complétée par une approche numérique de l'Interaction Fluide-Structure sur une configuration jugée critique. Une étape préliminaire à la simulation multi-physique est l'élaboration d'un modèle réaliste du comportement des PTF en élastomère. L'identification et la prise en considération des mécanismes propres aux élastomères dans la loi de comportement est ainsi détaillée. Les résultats issus des simulations sont confrontés aux résultats expérimentaux et rendent possible la validation d'un scénario de couplage dédié à l'Interaction Fluide-Structure au sein des MPS. / Fluid Structure Interaction of an inhibitor with the internal flow induced by wall injection was studied in an axisymmetric cold flow apparatus. Experiments were carried out over a wide range of injection velocities in order to underline how the obstacle not only modifies the mean flowfield but also its entire dynamic behaviour.The resulting instability (from the interaction of the unstable shear layer with the Taylor-Culick flow) exhibits a significant shift with respect to the Taylor-Culick instability and therefore emphasizes the strong impact of the inhibitor on hydrodynamics.The mecanisms responsible of such behaviour are studied in this thesis with the help of global linear stability analysis and multi-physics numerical computations.

Les moteurs à propergol solide de certains lanceurs spatiaux peuvent présenter des oscillations de poussée provoquées par des oscillations de la pression interne du moteur. Il est désormais connu qu'une instabilité hydrodynamique propre à l'écoulement et mise en évidence par une approche de stabilité linéaire, appelée Vortex Shedding Parietal, est le phénomène à l'origine de ces oscillations. Cependant plusieurs questions subsistent quant à l'apparition de ces oscillations. En particulier, on constate que ces oscillations ne sont observées que dans la seconde moitié du tir. On suppose que la transition laminaire-turbulent de l'écoulement joue un rôle dans ce phénomène. Par ailleurs, l'utilisation possible de la géométrie étoilée pour les blocs de propergol de géométrie étoilée oblige à reprendre l'étude de stabilité linéaire pour l'étude de ces oscillations. Cette thèse porte ainsi à la fois sur l'étude de la stabilité linéaire de l'écoulement induit par une géométrie étoilée et sur la transition laminaire-turbulent d'une configuration circulaire. Une solution analytique de l'écoulement au sein d'un conduit à paroi débitante de géométrie étoilée est d'abord obtenue. Le profil de vitesse associé peut présenter des points d'inflexion en fonction de la déformation radiale ou du nombre de Reynolds de l'écoulement. Ensuite les approches de stabilité linéaire locale et biglobale sont appliquées à cette solution. Par rapport au cas circulaire, cette étude de stabilité linéaire met en évidence la possible existence d'un ou plusieurs modes amplifiés. Enfin, une simulation numérique LES de la transition laminaire-turbulent de la configuration du montage VALDO est effectuée. / Solid rocket motors may exhibit thrust oscillations induced by pressure oscillations inside the engine. In the scientific literature, the so-called vortex shedding parietal, a hydrodynamic instability specific to the internal flow, has been proven to be the triggering mechanism through linear stability analyses. However, some questions still remain. These oscillations are observed during the second half of Ariane 5 launches only, the laminar-turbulent transition of the intern flow being expected to play a significant role in that behavior. Additionally, non-circular grain shapes used in some engines may impact the stability study. The present PhD work extends the base flow linear stability studies to star-shaped geometries and investigate transition mechanisms in a circular configuration. Firstly, an analytical expression of the base flow in a duct with a star-shaped cross-section is determined. Wall-normal injection is assumed, which is representative of solid rocket engines. The velocity profile may exhibit an inflection point depending on both the Reynolds number and wall radial deformation. Secondly, both local and biglobal stability analyses are applied to this flow. In contrast to the circular case, the flow may exhibit one or several unstable modes. Lastly, a large eddy simulation of the VALDO experimental test bench is performed. The influence of the acoustic impedance of the outflow boundary condition on the amplitude of pressure oscillations is highlighted.

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.

Ce travail s’inscrit dans le cadre de la prévision des oscillations de pression interne des moteurs à propergol solide. Il consiste à étudier la transition laminaire-turbulent de l’écoulement interne,modélisé par celui dans un conduit cylindrique à paroi débitante, et son lien avec l’instabilité naturelle de cet écoulement, le Vortex Shedding Pariétal (VSP). La démarche s’est organisée en trois temps. Des mesures antérieures sur un montage gaz froid, reproduisant l’écoulement modèle,sont analysées afin de mettre en évidence la transition laminaire-turbulent. Cette transition est ensuite imposée dans des simulations URANS afin de permettre l’étude de son influence sur les modes VSP. Enfin, une approche LES est mise en place pour simuler le développement de la transition dans les conditions de l’expérience ; dans ce but, une stratégie de perturbation spatiale de l’écoulement est utilisée. Cette étude met en avant quatre résultats principaux. La transition laminaire-turbulent découle de l’amplification spatiale des modes VSP. La simulation de ce processus met en évidence une forte influence de la perturbation numérique ajoutée à l’écoulement. D’autre part, les simulations URANS montrent que la transition réduit l’amplification des modes VSP et les oscillations de pression interne résultantes. Le rôle de la transition dans l’absence d’oscillations de pression lorsque le domaine a un grand rapport d’aspect, jusqu’alors supposé dans la littérature,est ainsi confirmé. Une particularité importante de cette transition est qu’elle dépend de la position radiale, l’écoulement étant turbulent près de la paroi débitante et laminaire au cœur. / The present work is related to the prediction of oscillations in solid rocket motors inner flow. It consists in a study of the laminar-turbulent transition of the motor’s inner flow, which is represented by a cylindrical injection-driven flow, and the relation between this phenomenon and the natural instability named Parietal Vortex Shedding (PVS). Three aspects have been analyzed.First of all, previous cold-gas experiments – reproducing the injection driven flow – are analyzed in order to highlight the transition laminar-turbulent transition. This transition is then imposedin URANS simulations to enable a study of its influence on the PVS modes. Finally, Large Eddy Simulations are performed to simulate the laminar-turbulent process. A strategy based on spatial steady disturbances is used to ease this process. The mains conclusions of this work are the following ones. The laminar-turbulent transition is a consequence of the spatial amplification of PVS modes. Simulations of this process highlight a strong influence of the injected numerical disturbances. The URANS simulations show that this transition reduces the amplification of PVSmodes, and the resulting pressure oscillations levels. These results confirm the role of the transitionin the absence of pressure oscillations when the motor cavity is long. A distinctive feature ofthis transition is its dependence on the radial position, which leads to the coexistence of a laminar region in the channel core and a turbulent region near the injecting wall at a given axial position.