Des crises vibratoires ont été constatées dans plusieurs centrales thermiques d’EDF opérant avec du fioul lourd, certaines ayant entraîné l’arrêt du foyer. Ce travail traite des instabilités de combustion pouvant se déclencher dans ce type de système où le combustible liquide est injecté avec de la vapeur d’eau et où l’écoulement d’air est mis en rotation. Ces phénomènes vibratoires résultent d’un couplage résonant entre la dynamique de la combustion et l’acoustique du foyer. La réponse acoustique des flammes diphasiques non-prémélangées swirlées reste largement méconnue et est difficilement analysable sur le foyer réel. L’objectif de ce travail est donc d’étudier la stabilité des chaudières EDF à partir de l’analyse de la réponse d’une flamme diphasique non-prémélangée swirlée issue d’un injecteur générique et soumise à des perturbations de la vitesse acoustique. Cette réponse est déterminée sur un dispositif original (DIFAV) équipé d’un swirler et d’un injecteur bi-fluides fonctionnant à la vapeur d’eau et au dodécane. Ce système est constitué des principaux éléments des brûleurs utilisés sur les centrales thermiques EDF à une échelle 1/7000. Le dispositif est conçu pour facilement modifier la géométrie de la tête d’injection, les conditions d’injection de combustible et de vapeur et ainsi contrôler le spray généré. Des visualisations à la sortie d’une buse d’injection montrent l’influence de la topologie de l’écoulement diphasique dans l’injecteur sur la taille des gouttes mesurées dans le spray. Des mesures de taille et de vitesse des gouttes lorsque le rapport des débits de vapeur et de combustible (GLR) est modifié sont réalisées. Ces données comparées à des modèles ont permis d’estimer l’évolution de la taille des gouttes générées par l’injecteur qui équipe les centrales thermiques lorsque le GLR varie. Une analyse modale du foyer DIFAV et d’un modèle simplifié de la chaudière réelle est ensuite menée. Les fréquences propres et les taux d’amortissement du foyer DIFAV sont déterminés expérimentalement en soumettant le système à une modulation acoustique externe. Un modèle acoustique simplifié composé de trois cavités couplées représentatif du brûleur DIFAV est également développé. Des simulations acoustiques réalisées avec COMSOL Multiphysics sur une coupe transverse d’une chaudière générique représentative de la chaudière industrielle permettent d’identifier trois modes à basses fréquences établis entre les plenums et la chambre de combustion qui sont susceptibles d’être instables. La sensibilité de ces modes à la géométrie du foyer et aux conditions limites est étudiée. La réponse de la flamme générique lorsqu’elle est soumise à des modulations acoustiques de l’écoulement d’air en amont du brûleur est ensuite mesurée sur le banc DIFAV pour différents niveaux d’excitation et deux topologies de flamme lorsque les conditions d’injection sont modifiées. Les mécanismes qui pilotent l’évolution du gain de l’une des fonctions de transfert généralisées (FDF) de la flamme sont étudiés à l’aide de visualisations en moyenne de phase de l’écoulement et de mesures des vitesses axiale et azimutale de l’écoulement d’air au cours d’un cycle de modulation. Une forte sensibilité de la phase de la FDF à l’amplitude des perturbations acoustiques est observée. Un adimensionnement par le nombre de Strouhal basé sur la vitesse débitante et la longueur efficace de la flamme est proposé pour transposer ces FDFs sur le brûleur réel. Une analyse de stabilité du foyer DIFAV est réalisée en intégrant les FDF au modèle acoustique afin de déterminer les cycles limites des oscillations lorsque la longueur de la chambre de combustion varie. Ces calculs sont comparés aux fréquences des instabilités auto-entretenues mesurées aux cycles limites dans le foyer DIFAV. [...] / Vibratory crises have been observed in EDF thermal power plants operating with heavy fuel oil. Such instabilities may lead to shutdown and damage the boiler. This work deals with combustion instabilities that can take place in boilers equipped with steam-assisted atomizers and where the airflow is swirled. These vibratory phenomena result from a resonant coupling between the combustion dynamics and the boiler acoustics. Analyses of combustion dynamics of non-premixed swirling spray flames remain rare and are difficult to realize on the real system. The objective of this work is to analyze the stability of EDF boilers using the response of generic non-premixed swirling spray flames submitted to acoustic velocity disturbances. This response is determined on an original device (DIFAV) equipped with a swirling vane and a twin-fluid atomizer operated with steam and dodecane. This burner is equipped with the main elements of those used in the thermal power plant, but has a reduced scale of 1/7000. The influence of the injector geometry and of the operating conditions on the spray generated by the injector can be studied. Spray visualizations at the outlet of the injector reveal the relationship between the topology of the two-phase flow in the injector and the measured droplet size. Measurements of the droplet diameter and velocity as a function of the gas-to-liquid ratio (GLR) have been performed at the outlet of the injector. These data have been compared to models and were used to estimate the evolution of the droplets diameter as a function of the GLR generated by the industrial injector. A modal analysis of the DIFAV combustor is then carried out and a simplified acoustic model made of three coupled cavities is developed. The natural frequencies and damping rates of the DIFAV combustor are determined experimentally when it is submitted to acoustic modulation. Acoustic simulations are performed with COMSOL Multiphysics on a simplified geometrical model of the industrial boiler. Three low frequency modes established between the plenums and the combustion chamber have been identified and may be unstable. Their sensitivity to modifications of the boiler geometry and boundary conditions are studied. Flame responses subjected to acoustic modulations of the airflow rate are then measured on the DIFAV combustor for several amplitudes and two flames topologies obtained at globally lean condition. Phase-conditioned flame visualizations and measurements of swirl number fluctuations during an acoustic forcing cycle are conducted to explain the mechanisms that control the evolution of gain of the Flame Describing Function (FDF). A high sensitivity of the phase of the FDF to the amplitude of the acoustic disturbance is observed. The Strouhal number based on the airflow velocity and the effective length of the flame is used to transpose these FDF on the industrial burner. FDF are integrated in the acoustic model of the DIFAV setup to carry out a stability analysis and predict the limit cycle oscillations as a function of the combustion chamber length. These calculations are compared to frequencies of self-sustained instability measured at the limit cycles in the DIFAV combustor. A reasonable agreement is obtained showing the validity of the stability analysis for the non-premixed two-phase flames investigated based on the knowledge of their FDF. Finally, a stability analysis of the EDF boiler is conducted with the COMSOL Multiphysics model by including the acoustic flame response of the industrial burner in the simulation. This FDF is deducted from the dimensionless FDF measured on the generic burner. The Rayleigh criterion is used to analyze the stability of the combustor as a function of the flame length for different boundary conditions. Indications are given to improve the stability of the EDF boiler.

Pour des raisons évidentes d’opération et de sécurité, l’allumage est un problème essentiel dans les moteurs aéronautiques. La conception d’une chambre de combustion de turbine à gaz intègre de multiples objectifs contradictoires, l’un d’entre eux étant un allumage ou ré-allumage efficace des brûleurs. Parmi les paramètres dont disposent les ingénieurs dans la phase d’optimisation du design, le nombre de systèmes d’injection de carburant et leur espacement sont des points cruciaux qui doivent être fixés dès le début. En effet, de tels choix ont non seulement un impact sur le coût de fabrication et la taille de la chambre mais ils affectent aussi l’efficacité d’un moteur ainsi que ses caractéristiques d’allumage. Afin d’améliorer les connaissances relatives au processus l’allumage dans des moteurs réels, la recherche actuelle combine des expériences fondamentales de plus en plus complexes et des simulations numériques de haute fidélité. Ces actions se concentrent d’une part sur les premiers instants où le noyau de flamme apparaît et d’autre part sur la phase de propagation entre les différents brûleurs. Ces deux phases sont capitales mais restent difficiles à étudier simultanément. Le premier objectif de cette thèse vise à évaluer les modèles SGE sur un seul brûleur expérimental situé au CORIA (France) pour mettre en place une méthodologie fiable afin de réaliser numériquement une séquence d’allumage dans des conditions d’opération réelles et équivalentes aux premiers instants. Une telle étude met en jeu plusieurs phénomènes tels que les écoulement swirlés, l’allumage, l’extinction, la propagation de flamme et les interactions flamme/turbulence. Tous ces processus et mécanismes interagissent et augmentent de façon significative le niveau de difficulté, notamment pour modéliser la combustion turbulente d’un tel allumage. Ces modèles requièrent donc d’être évalués précisément. Ensuite, ce travail examine par la simulation numérique la phase de propagation en utilisant les expériences réalisées sur une chambre composée de plusieurs injecteurs. La comparaison des séquences d’allumage obtenues numériquement avec celles des données expérimentales montre que la SGE reproduit les bonnes tendances et s’avère prédictive. D’un point de vue global, les caractéristiques de propagation du front de flamme en direction des injecteurs voisins sont bien capturées par le numériquemontrant desmodes de propagation identiques à ceux obtenus expérimentalement (radial ou axial) et des temps d’allumage similaires. Pour finir, l’analyse détaillée de ces données numériques a permis d’identifier les mécanismes principaux qui sont à l’origine des différents modes de propagation. / For evident operational and safety reasons, ignition is a key feature of aeronautical gas turbine applications. In fact the design of a gas turbine combustion chamber imposes multiple contradicting objectives one of them being efficient ignition or re-ignition. Among all the parameters available to the engineers, the number of fuel injection systems and their spacing are crucial elements, that must be fixed early on in the design phase. Such choices however not only impact the manufacturing cost and size of the combustor but they also affect the operability of the engine as well as its ignition. To improve knowledge of the ignition process occurring in real engines, current research combines fundamental and increasingly complex experiments complemented by high fidelity numerical simulations. These actions focus on the one hand on the initial instants where the first flame kernel appears as well as the follow-on instants corresponding to the light-around phase or burner to burner flame propagation phase. Both phases are clearly important but are difficult to study simultaneously. The first purpose of this thesis aims at assessing LES models on a single experimental burner located at CORIA (France) to provide a reliable numerical methodology to achieve an ignition sequence in real engines. Indeed, various phenomena are involved in such numerical studies dedicated to real aeronautical combustion chambers and all need to be reproduced by numerics: swirling flows, ignition, quenching, flame propagation, flame/turbulence interactions. All of these processes interact and clearly raise the level of difficulty notably in terms of turbulent combustion modeling of an ignition transient. Having assessed the method on a single burner configuration, the work then investigates the second phase, using a multi-injector experiment simulated by LES to study the flame propagation during ignition. The comparison of numerical fully transient ignition sequences with experimental data shows that LES recovers features found in the experiment. Global events such as the propagation of the flame front to neighboring swirlers are well captured and correct propagation modes (radial or axial) as well as correct overall ignition time delay are obtained. Finally the detailed analysis of LES data allows to identify the driving mechanisms governing each of these propagation modes.