Les réglementations toujours plus drastiques sur les émissions de polluants ont conduit au développement de systèmes de combustion opérant en régimes pauvres qui sont malheureusement sujet aux instabilités thermo acoustiques. La capacité de la Simulation aux Grandes Echelles (SGE) à simuler des turbines à gaz industrielles complexes de grande puissance est mise en évidence au cours de ce travail de thèse. Tout d’abord, la SGE est appliquée à un brûleur académique et validée par comparaison à des mesures effectuées à l’Université de Berlin ainsi qu’à des simulations SGE effectuées avec OpenFOAM chez Siemens. Afin de déterminer la stabilité de ce bruleur le couplage entre l’acoustique et la combustion est modélisé par l’approche de type fonction de transfert de flamme (FTF). Suite à ces calcules et l’évaluation de la FTF les fluctuations du nombre de swirl sont identifiées comme un paramètre à même de modifier cette réponse de flamme. Après cette première étape de validation, une turbine à gaz industrielle est simulée en SGE pour deux géométries différentes du brûleur et pour deux points de fonctionnement. La FTF issue de ces calculs est peu influencée par les deux points de fonctionnement. A l’inverse, des légères modifications de la géométrie du swirler modifient les caractéristiques de la FTF montrant que plusieurs mécanismes sont en jeu. Ces mécanismes sont identifiés comme étant la vitesse d’entrée, les fluctuations de swirl et les fluctuations de fraction de mélange. Cette dernière est causée par: 1) la pulsation du débit de carburant injecté et 2) la trajectoire fluctuante des jets de carburant. Bien que le swirler soit conçu pour fournir un mélange le plus homogène possible, d’importantes hétérogénéités de mélange à l’entrée de la chambre de combustion sont présentes. Les perturbations de mélange se combinent avec les fluctuations de vitesse (et donc avec les fluctuations de swirl) aboutissant à des résultats de FTF différents. Un modèle étendu pour la FTF reliant le dégagement de chaleur à la vitesse d’entrée et à la fluctuation de fraction de mélange (modèle MISO) se révèle être une bonne solution pour ces systèmes complexes. Une analyse non linéaire montre en outre que l’amplitude de forçage conduit non seulement à une saturation de la flamme, mais aussi à un changement de la réponse de flamme. La saturation de la flamme n’est vérifiée que pour la FTF globale et le gain augmente localement avec une amplitude croissante. Pour ce système on notera enfin que la flamme linéaire, comme la flamme non linéaire, ne sont pas compactes: certaines zones pourtant situées l’une à coté de l’autre, ont des différences significatives de délai de FTF, montrant que certaines parties de la flamme amortissent l’excitation alors que d’autres l’amplifient. / Modern pollutant regulation have led to a trend towards lean combustion systems which are prone to thermo-acoustic instabilities. The ability of Large Eddy Simulation (LES) to handle complex industrial heavy-duty gas turbines is evidenced during this thesis work. First, LES is applied to an academic single burner in order to validate the modeling against measurements performed at TU Berlin and against OpenFoam LES simulations done at Siemens. The coupling between acoustic and combustion is modeled with the Flame Transfer Function (FTF) approach and swirl number fluctuations are identified changing the FTF amplitude response of the flame. Then, an industrial gas turbine is analyzed for two different burner geometries and operating conditions. The FTF is only slightly influenced for the two operating points but slight modifications of the swirler geometry do modify the characteristics of the FTF showing that a simple model taking only into account the flight time is not appropriate and additional mechanisms are at play. Those mechanisms are identified being the inlet velocity, the swirl and the inlet mixture fraction fluctuations. The latter is caused by two mechanisms: 1) the pulsating injected fuel flow rate and 2) the fluctuating trajectory of the fuel jets. Although the diagonal swirler is designed to provide good mixing, effects of mixing heterogeneities at the combustion chamber inlet occur. Mixture perturbations phase with velocity (and hence with swirl) fluctuations and combine with them to lead to different FTF results. Another FTF approach linking heat release to inlet velocity and mixture fraction fluctuation (MISO model) shows further to be a good solution for complex systems. A nonlinear analysis shows that the forcing amplitude not only leads to a saturation of the flame but also to changes of the delay response. Flame saturation is only true for the global FTF and the gain increases locally with increasing forcing amplitude. Both, the linear and the nonlinear flames, are not compact: flame regions located right next to each other exhibited significant differences in delay meaning that at the same instant certain parts of the flame damp the excitation while others feed it.

L'allumage et le réallumage de haute altitude présentent de grandes difficultés dans le cadre des chambres de combustion aéronautiques. Le succès d'un allumage dépend de multiples facteurs, des caractéristiques de l'allumeur à la taille des gouttes du spray en passant par le niveau de turbulence au point d'allumage. Déterminer la position optimale de l'allumeur ou le potentiel d'allumage d'une source d'énergie donnée à une position donnée sont ainsi des paramètres essentiels lors du design de chambre de combustion. Le but de ces travaux de thèse est d'étudier l'allumage forcé des chambres de combustion aéronautiques. Pour cela, des Simulation numériques aux Grandes Echelles (SGE) d'écoulements diphasiques réactifs sont utilisées et analysées. Afin de les valider, des données expérimentales issues du banc MERCATO installé à l'ONERA Fauga-Mauzac sont utilisées. Cela permet dans un premier temps de valider la méthodologie ainsi que les modèles utilisés pour les SGE diphasiques évaporantes avant leur utilisation dans d'autres conditions d'écoulement. Le cas diphasique réactif statistiquement stationnaire est ensuite comparé aux données disponibles pour évaluer les modèles en condition réactives. Ce cas est étudié plus en détail à travers l'analyse de caractéristiques de la flamme. Celle-ci semble être le théâtre de régimes de combustion très différents. On note aussi que la détermination de la méthode numérique la plus appropriée pour le calcul d'écoulements diphasiques n'est pas évidente. De plus, deux méthodes numériques différentes peuvent donner des résultats en bon accord avec l'expérience et pourtant avoir des modes de combustion différents. Les capacités de la SGE à correctement calculer un écoulement diphasique réactif étant validé, des SGE du phénomène transitoire d'allumage sont effectuées. La sensibilité observée expérimentalement de l'allumage aux conditions initiales, i.e. à l'instant de claquage, est retrouvé par les SGE. L'analyse met en évidence le rôle prépondérant de la dispersion du spray dans le développement initial du noyau de flamme. L'utilisation des SGE pour calculer les séquences d'allumage fournie de nombreuses informations sur le phénomène d'allumage, cependant d'un point de vue industriel, cela ne donne pas de résultat optimal, à moins de ne tester toutes les positions, ce qui rendrait le coût CPU déraisonnable. Des alternatives sont donc nécessaires et font l'objet de la dernière partie de ces travaux. On propose de dériver un critère local d'allumage, donnant la probabilité d'allumage à partir d'un écoulement diphasique (air et carburant) non réactif instationnaire. Ce modèle est basé sur des critères liés aux différentes phases menant à un allumage réussi, de la formation d'un premier noyau à la propagation de la flamme vers l'injecteur. Enfin, des comparaisons avec des données expérimentales sur des chambres aéronautiques sont présentées et sont en bon accord, indiquant que le critère d'allumage proposé, couplé avec une SGE d'écoulement diphasique non réactif, peut être utilisé pour optimiser la puissance et la position du système d'allumage. / Ignition and altitude reignition are critical issues for aeronautical combustion chambers. The success of ignition depends on multiple factors, from the characteristics of the igniter to the spray droplet size or the level of turbulence at the ignition site. Finding the optimal location of the igniter or the potential of ignition success of a given energy source at a given location are therefore parameters of primary importance in the design of combustion chambers. The purpose of this thesis is to study forced ignition of aeronautical combustion chambers. To do so, Large Eddy Simulations (LES) of two-phase reacting flows are performed and analyzed. First, the equations of the Eulerian formalism used to describe the dispersed phase are presented. To validate the successive LES, experimental data from the MERCATO bench installed at ONERA Fauga-Mauzac are used. It allows to validate the two-phase evaporating flow LES methodology and models prior to its use to other flow conditions. The statistically stationary two-phase flow reacting case is then compared to available data to evaluate the model in reacting conditions. This case is more deeply studied through the analysis of the characteristics of the flame. This last one appears to experience very different combustion regimes. It is also seen that the determination of the most appropriate methodology to compute two-phase flow flame is not obvious. Furthermore, two different methodologies may both agree with the data and still have different burning modes. The ability of the LES to correctly compute burning two-phase flow being validated, LES of the transient ignition phenomena are performed. The experimentally observed sensitivity of ignition to initial conditions, i.e. to sparking time, is recovered with LES. The analysis highlights the major role played by the spray dispersion in the development of the initial flame kernel. The use of LES to compute ignition sequences provides a lot of information about the ignition phenomena, however from an industrial point of view, it does not give an optimal result, unless all locations are tested, which brings the CPU cost to unreasonable values. Alternatives are hence needed and are the objective of the last part of this work. It is proposed to derive a local ignition criterion, giving the probability of ignition from the knowledge of the unsteady non-reacting two-phase (air and fuel) flow. This model is based on criteria for the phases of a successful ignition process, from the first kernel formation to the flame propagation towards the injector. Then, comparisons with experimental data on aeronautical chambers are done and show good agreement, indicating that the proposed ignition criterion, coupled to a Large Eddy Simulation of the stationary evaporating two-phase non-reacting flow, can be used to optimize the igniter location and power.

L'amélioration des moteurs à allumage commandé représente un défi de première importance pour les ingénieurs afin de produire plus d'énergie, de consommer moins de matière première et de réduire les émissions polluantes. Les nouvelles technologies apparues ces dernières années amènent les moteurs de plus en plus proches de leurs limites de fonctionnement, favorisant ainsi des phénomènes néfastes qui doivent être contrôlés. Parmi ces phénomènes, les variations cycle-à-cycle (VCC) doivent être minimisées pour garder une performance optimale et éviter une dégradation rapide du moteur. La Simulation aux Grandes Echelles (SGE) est un outil prometteur afin de prédire numériquement les niveaux de variabilités obtenues lors du design d'un moteur (limitant ainsi les coûteuses campagnes de mesures expérimentales). Ce manuscrit s'est attaché à développer une méthodologie numérique pour la prédiction des variabilités cycliques, à simuler un nombre suffisant de cycles pour pouvoir estimer les niveaux de VCC et à valider les résultats obtenus par rapport aux résultats expérimentaux. La SGE semble capter les points de fonctionnements stable et instable étudiés. Les sources qui provoquent ces VCC ont aussi été analysées et une modification du fonctionnement du moteur a été proposée afin de réduire les VCC. / The improvement of the spark ignition engines is a major challenge for engineers in order to produce more energy, to minimize fuel consumption and to reduce the pollutant emissions. The new technologies which appear in the last years bring the engines closer to their stability limit while increasing various unwanted phenomena. Among these phenomena, cycle-to-cycle variation (CCV) need to be minimized in order to keep the performances as high as possible and avoid damages on the engines. Large Eddy Simulation (LES), which is a very promising tool in order to predict the level of CCV of a given engine, has been used in the present document to simulate a mono-cylinder spark ignition engine. The present document presented a numerical methodology for the prediction of CCV, numerous engine cycles were simulated by LES in order to validate the results in comparison to the experimental findings. The LES seems to be able to capture stable and instable (in terms of CCV) operating points of the engine. In addition, the sources of CCV were also analyzed and a modification of the engine has been proposed to reduce CCV.

The current demands on Swedish steel industry to produce low quantity batches of specialized products requires research on steel casting processes. There are several physical processes that need be taken into account for this problem to be viewed in full light such as thermal-processes, solidification and fluid dynamics. This work focuses on the fluid-dynamics part; more specifically, the dependence of flow quality within the caster on nozzle and mould geometry. The simulations are carried out using a scale-resolving method, in specific LES (Large Eddy Simulation) which is coupled with a DPM (Discrete Phase Model) to model Argon behaviour. The results of these simulations are presented and validated against physical experiment and data from industrial trials. Conclusions are drawn regarding optimal nozzle types in respect to different mould geometries. The mould eigenfrequencies are shown to exhibit a connection with the casting velocity. This results in so called sweet spots in casting velocity where flow irregularities due to sloshing is minimal. It is shown that the mountain type nozzle is preferable for smaller geometries whilst comparatively larger geometries benefit from a cup type. / FLOWFLEX CC

Hot wire anemometers have been used in several wind velocity sensors deployed in Mars. They are based in keeping the temperature of a surface at a constant value, above the ambient. This is done by means of a heater controlled with an electronic system. The cooling rate of each point at the sensor surface can be used to calculate the wind velocity and direction. However, due to turbulent fluctuations, the cooling rate is not constant even in the case of constant velocity. Moreover, RANS simulations cannot estimate such fluctuations as they only provide an estimation of the averaged flow field. The goal of this work has been to estimate such fluctuations and the e↵ect they might have on the sensor readings. To do so, the turbulent cooling rate (Nusselt number) of a sensor with a generic shape, under the typical conditions to be find in Mars, has been simulated using high performance LES (Large Eddy Simulation) simulations and compared with RANS and URANS simulations.

Combustion research is at the forefront of development of clean and efficient IC engines, gas turbines, rocket propulsion systems etc. With the advent of faster computers and parallel programming, computational studies of turbulent combustion is increasing rapidly. Many turbulent combustion models have been previously developed based on certain underlying assumptions. One of the major assumptions of the models is the regime it can be used for: either premixed or non-premixed combustion. However in reality, combustion systems are multi-regime in nature, i.e.,\ co-existence of premixed and non-premixed modes. Thus, there is a need for development of multi-regime combustion models which closely follows the physics of combustion phenomena. Much of previous modeling efforts for multi-regime combustion was done using flamelet-type models. As a first, the current study uses the highly robust transported Probability Density Function (PDF) method coupled with Large Eddy Simulation (LES) to develop a multi-regime model. The model performance is tested for Sydney Flame L, a piloted methane-air turbulent flame. The concept of flame index is used to detect the extent of premixed and non-premixed combustion modes. The drawbacks of using the traditional flame index definition in the context of PDF method are identified. Necessary refinements to this definition, which are based on the species gradient magnitudes, are proposed for the multi-regime model development. This results in identifying a new model parameter beta which defines a gradient threshold for the calculation of flame index. A parametric study is done to determine a suitable value for beta, using which the multi-regime model performance is assessed for Flame L by comparing it against the widely used non-premixed PDF model for three mixing models: Modified Curl (MCurl), Interaction by Exchange with Mean (IEM) and Euclidean Minimum Spanning Trees (EMST). The multi-regime model shows a significant improvement in prediction of mean scalar quantities compared to the non-premixed PDF model when MCurl mixing model is used. Similar improvements are observed in the multi-regime model when IEM and EMST mixing models are used. The results show potential foundation for further multi-regime model development using PDF model.

Les réseaux de tuyauterie industriels sont le siège de niveaux importants de vibrations induites par l’écoulement qui peuvent mener à la rupture par fatigue des installations. La présente étude se concentre sur l’analyse et modélisation simplifiée de la source vibratoire associée au passage d’un écoulement liquide turbulent par un coude à 90°. Une approche combinant expériences et simulation a été conduite. Une boucle de circulation d’eau munie d’un coude transparent a été conçue afin de permettre des me-sures de vitesse à l’intérieur du coude. Pour cette finalité, les techniques de Vélocimétrie par Imagerie de Particules (PIV) plane et stéréoscopique ont été utilisées. La pression pariétale et les vibrations du banc d’essais ont été mesurées simultanément. Plusieurs configurations d’écoulement ont été testées afin d’obtenir une riche base de données couplées reliant l’écoulement fluide à l’excitation dynamique des parois et, finalement, à la réponse vibratoire de la structure. En parallèle, l’écoulement instationnaire d’eau dans le coude a été simulé au moyen d’une approche du type Simulation des Grandes Echelles (LES). La simulation fluide a permis d’étudier en détails la topologie de l’écoulement turbulent au passage du coude ainsi que le champ instationnaire de pression fluctuante induit sur la paroi. Finale-ment, un ensemble d’outils statistiques a été appliqué aux données expérimentales et numériques afin de proposer un modèle simplifié des transferts qui relient l’écoulement turbulent à la sollicitation dynamique de la structure contenant le coude. / The flow of dense fluids within thin-walled piping systems may lead to significant levels of Flow-Induced Vibration, mainly in the vicinities of singularities such as obstacles inserted into the flow, sudden changes of cross-sectional area or flow direction. This study focuses on the analysis and reduced-order modelling of the vibrational source associated with the turbulent flow of liquids through a 90° elbow.A mixed experimental-computational approach is undertaken. A closed water loop containing a transparent elbow was designed in order to allow for fluid velocity measurements inside the singularity. To this purpose, planar and stereoscopic Particle Image Velocimetry (PIV) were employed. Wall pressure fluctuations and structural vibrations were measured simultaneously. Several flow configurations were tested in order to obtain a large coupled database linking the flow to the dynamic excitation, and then to the vibration response of the structure.In parallel, the unsteady water flow through the elbow was computed using Large-Eddy Simulation (LES). The fluid simulation allowed for a detailed study of the turbulent flow through the singularity and the unsteady pressure field induced on the piping walls. Finally, a set of statistical tools was applied to both experimental and computational data in order to propose a reduced-order model of the transfer function that links the tur-bulent flow to the dynamic excitation of the elbowed piping structure.

Interactions between flames and thermal phenomena are the guiding thread of this work. Flamesproduce heat indeed, but can also be affected by it. Large Eddy Simulations (LES) are used hereto investigate these interactions, with a focus on two main topics: wall heat transfer andcombustion control. In a first part, wall heat transfer in a rocket engine sub-scale CH4/O2 burner isstudied. In the context of launchers re-usability and cost reduction, which are major challenges,new propellant combinations are considered and wall heat fluxes have to be precisely predicted.The aim of this work is to evaluate LES needs and performances to simulate this kind ofconfiguration and provide a computational methodology permitting to simulate variousconfigurations. Numerical results are compared to experimental data provided by the TechnischeUniversität München (Germany). In a second part, combustion control by means of NanosecondRepetitively Pulsed (NRP) plasma discharges is studied. Modern gas turbine systems use indeedlean combustion with the aim of reducing fuel consumption and pollutant emissions. Lean flamesare however known to be prone to instabilities and combustion control can play a major role in thisdomain. A phenomenological model which considers the plasma discharges as a heat source isdeveloped and applied to a swirl-stabilized CH4/Air premixed lean burner. LES are performed inorder to evaluate the effects of the NRP discharges on the flame. Numerical results are comparedwith experimental observations made at the King Abdulla University of Science and Technology(Saudi Arabia).

Combustion in cryogenic engines is a complex phenomenon, involving either liquid or supercritical fluids at high pressure, strong and fast oxidation chemistry, and high turbulence intensity. Due to extreme operating conditions, a particularly critical issue in rocket engine is wall heat transfer which requires efficient cooling of the combustor walls. The concern goes beyond material resistance: heat fluxes extracted through the chamber walls may be reused to reduce ergol mass or increase the power of the engine. In expander-type engine cycle, this is even more important since the heat extracted by the cooling system is used to drive the turbo-pumps that feed the chamber in fuel and oxidizer. The design of rocket combustors requires therefore an accurate prediction of wall heat flux. To understand and control the physics at play in such combustor, the Large Eddy Simulation (LES) approach is an efficient and reliable numerical tool. In this thesis work, the objective is to predict wall fluxes in a subcritical rocket engine configuration by means of LES. In such condition, ergols may be in their liquid state and it is necessary to model liquid jet atomization, dispersion and evaporation.The physics that have to be treated in such engine are: highly turbulent reactive flow, liquid jet atomization, fast and strong kinetic chemistry and finally important wall heat fluxes. This work first focuses on several modeling aspects that are needed to perform the target simulations. H2/O2 flames are driven by a very fast chemistry, modeled with a reduced mechanism validated on academic configurations for a large range of operating conditions in laminar pre- mixed and non-premixed flames. To form the spray issued from the atomization of liquid oxygen (LOx) an injection model is proposed based on empirical correlations. Finally, a wall law is employed to recover the wall fluxes without resolving directly the boundary layer. It has been specifically developed for important temperature gradients at the wall and validated on turbulent channel configurations by comparison with wall resolved LES. The above models are then applied first to the simulation of the CONFORTH sub-scale thrust chamber. This configuration studied on the MASCOTTE test facility (ONERA) has been measured in terms of wall temperature and heat flux. The LES shows a good agreement compared to experiment, which demonstrates the capability of LES to predict heat fluxes in rocket combustion chambers. Finally, the JAXA experiment conducted at JAXA/Kakuda space center to observe heat transfer enhancement brought by longitudinal ribs along the chamber inner walls is also simulated with the same methodology. Temperature and wall fluxes measured with smooth walls and ribbed walls are well recovered by LES. This confirms that the LES methodology proposed in this work is able to handle wall fluxes in complex geometries for rocket operating conditions.

The marine current turbine (MCT) is an exciting proposition for the extraction of renewable tidal and marine current power. However, the numerical prediction of the performance of the MCT is difficult due to its complex geometry, the surrounding turbulent flow and the free surface. The main purpose of this research is to develop a computational tool for the simulation of a MCT in turbulent flow and in this thesis, the author has modified a 3D Large Eddy Simulation (LES) numerical code to simulate a three blade MCT under a variety of operating conditions based on the Immersed Boundary Method (IBM) and the Conservative Level Set Method (CLS). The interaction between the solid structure and surrounding fluid is modelled by the immersed boundary method, which the author modified to handle the complex geometrical conditions. The conservative free surface (CLS) scheme was implemented in the original Cgles code to capture the free surface effect. A series of simulations of turbulent flow in an open channel with different slope conditions were conducted using the modified free surface code. Supercritical flow with Froude number up to 1.94 was simulated and a decrease of the integral constant in the law of the wall has been noticed which matches well with the experimental data. Further simulations of the marine current turbine in turbulent flow have been carried out for different operating conditions and good match with experimental data was observed for all flow conditions. The effect of waves on the performance of the turbine was also investigated and it has been noticed that this existence will increase the power performance of the turbine due to the increase of free stream velocity.