Etude numérique de la stabilisation de flamme et des instabilités de combustions dans les systèmes de propulsion

Cette thèse se focalise sur la compréhension et la prédiction des instabilités de combustion dans les systèmes à haute pression. Elle s'oriente autour de la simulation numérique d’un banc d'essai, opéré à l'université de Purdue, comprenant un injecteur caractéristique des moteurs-fusées et dont les propriétés acoustiques peuvent varier à l'aide d’un tube d'injection mobile. Une méthode d'initialisation et d'allumage pour les calculs LES de chambres de combustions terminées par une tuyère est présentée. Un point de fonctionnement instable est choisi pour étudier le mécanisme de l'instabilité. Les simulations sont comparées aux résultats expérimentaux en terme de fréquence et structure du mode instable. La fonction de transfert de flamme est calculée à l'aide du modèle n − τ puis implémentée dans un solveur acoustique (ne résolvant que les perturbations acoustiques à partir de l'équation de Helmholtz en écoulement réactif). Différents modèles d'impédance de tuyère, extraits de la littérature, sont comparés et leurs impacts sur les résultats de stabilité sont analysés. Le théorème d’impédance translatée est implémenté dans le solveur acoustique pour analyser, à faible coût de calcul, l’influence de la variation de la longueur du tube d'injection. Des écarts entre les fréquences prédites et celles trouvées expérimentalement subsistent mais la carte de stabilité de l’expérience est bien reproduite.This thesis focuses on the understanding and the prediction of combustion instability in high-pressure devices. A model rocket combustor, tested experimentally at Purdue University, with continuously variable acoustic properties, thanks to a variable-length injector tube, is simulated. A method to initialize and ignite Large-Eddy-Simulation (LES) calculation of combustion chamber surrounded by nozzle is proposed. An unstable operating point is then chosen to investigate the mechanism of the instability. The simulations are compared to experimental results in terms of frequency and mode structure. The flame transfer function is calculated using the n − τ model to feed an acoustic solver which solves only the acoustic perturbation using a Helmholtz equation in reacting flows. The importance of the modeling of the nozzles impedance is studied through the main theories in the literature. The impedance translation theorem is implemented in the acoustic solver to analyze at low cost the influence of the variation of the injector tube. Despite differences in frequency of the instability, the stability map of the experiment is well reproduced