No description available.

Neste trabalho, buscou-se investigar a aplicação de uma modelagem RANS, com a utilização do modelo de turbulência das Tensões de Reynolds (RSTM), na modelagem do escoamento de gás no interior de ciclones, especialmente na simulação do movimento periódico do vórtice, conhecido como Precessing Vortex Core (PVC), com apoio da Fluidodinâmica Computacional. Utilizou-se o pacote de CFD de código aberto OpenFOAM, e as simulações foram realizadas em regime transiente. Dois modelos RSTM foram testados, o modelo LRR com os valores padrões para as constantes e uma modificação das constantes do termo de redistribuição. Compararam-se os resultados obtidos com dados da literatura e verificou-se que ambos os modelos representaram o campo médio de velocidade. No entanto, somente o modelo LRR conseguiu reproduzir o campo de flutuações de velocidade. Detectou-se o movimento de precessão do vórtice, que pode ser visualizado através de imagens sequenciais do escoamento, e uma frequência característica do movimento pode ser calculada através do Espectro de potência. / The flow field of a gas cyclone was simulated with a RANS model using Computational Fluid Dynamics (CFD) and the Precession of the Vortex Core (PVC) was investigated. The CFD simulations were carried out using the open source code OpenFOAM and the simulations were performed in a transient regime. Two Reynolds Stress Turbulence Models (RSTM) were investigated - the LRR model with standard constants and another with modified constants in the pressure-strain term. The results obtained by the models were compared with values of mean velocity taken from the literature and predictions were satisfactory. However, only the LRR model with standard constants predicted with precision the fluctuation velocity field. The precessing vortex was detected and illustrated by the flow field visualization and a characteristic frequency was calculated by Power Spectrum Density.

Neste trabalho, buscou-se investigar a aplicação de uma modelagem RANS, com a utilização do modelo de turbulência das Tensões de Reynolds (RSTM), na modelagem do escoamento de gás no interior de ciclones, especialmente na simulação do movimento periódico do vórtice, conhecido como Precessing Vortex Core (PVC), com apoio da Fluidodinâmica Computacional. Utilizou-se o pacote de CFD de código aberto OpenFOAM, e as simulações foram realizadas em regime transiente. Dois modelos RSTM foram testados, o modelo LRR com os valores padrões para as constantes e uma modificação das constantes do termo de redistribuição. Compararam-se os resultados obtidos com dados da literatura e verificou-se que ambos os modelos representaram o campo médio de velocidade. No entanto, somente o modelo LRR conseguiu reproduzir o campo de flutuações de velocidade. Detectou-se o movimento de precessão do vórtice, que pode ser visualizado através de imagens sequenciais do escoamento, e uma frequência característica do movimento pode ser calculada através do Espectro de potência. / The flow field of a gas cyclone was simulated with a RANS model using Computational Fluid Dynamics (CFD) and the Precession of the Vortex Core (PVC) was investigated. The CFD simulations were carried out using the open source code OpenFOAM and the simulations were performed in a transient regime. Two Reynolds Stress Turbulence Models (RSTM) were investigated - the LRR model with standard constants and another with modified constants in the pressure-strain term. The results obtained by the models were compared with values of mean velocity taken from the literature and predictions were satisfactory. However, only the LRR model with standard constants predicted with precision the fluctuation velocity field. The precessing vortex was detected and illustrated by the flow field visualization and a characteristic frequency was calculated by Power Spectrum Density.

Low-emissions combustor design is crucially important to gas turbine engine manufacturers. Unfortunately, many designs are susceptible to unsteady oscillations that can result in structural fatigue and increased noise. Computational approaches that resolve flow unsteadiness, for example Large Eddy Simulation (LES), are being explored as one avenue to help understand such phenomena. However, in order to quantifY the accuracy of LES predictions, benchmark validation data in suitably chosen test cases are required. Comprehensive experimental data covering both time-averaged and timeresolved features are currently scarce. It was the aim of this thesis, therefore, to provide such data .in a configuration representing the near-field of a typical gas turbine fuel injector. It was decided to focus on the fuel injector since many unsteady events are believed to originate because of the transient interactions between the fuel injector flow and the main combustor flow. A radial fed two-stream fuel injector, based on a preexisting industrial gas-turbine Turbomeca design was used, since this geometry was known to be susceptible to unsteadiness. The fuel injector was investigated under isothermal conditions to place emphasis on the fluid mechanical behaviour of the fuel injector, including detailed capture of any unsteady phenomena present. Light Sheet Imaging (LSI) systems were used as the primary experimental technique to provide high quality spatially and temporally resolved instantaneous velocity and scalar field information in 2D planes (using ParticieImage Velocimetry (PIV) and Planar LaserInduced Fluorescence (PUF) techniques). Several methods were employed to extract information quantifYing the flow unsteadiness and improve visualisation of timedependent large-scale turbulent structures. Proper Orthogonal Decomposition (POD) analysis enabled clear identification of the dominant modes of energy containing structures. The results indicated that periodic high-energy containing vortex structures occurred in the swirl stream shear layer, emerging from the fuel injector. These formed a two-strong two-weak rotating vortex pattern which propagated down the main duct flow path. The formation of these vortices was found to be a function of the swirl number and originated due to an interaction between the forward moving swirl flow and the furthest upstream penetration point ofthe recirculation zone present in the main duct flow. Dependent on the magnitude of the swirl number (influencing the swirl stream cone angle) and the geometry of the fuel injector, the vortex formation point was sometimes found inside the fuel injector itself. If the vortices originated inside the fuel injector they appeared much more coherent in space and time and of higher energy. A second unsteady high energy containing phenomenon was also identified, namely a Precessing Vortex Core (PVC), which was damped out if the fuel injector contained a central jet. The dynamics of the PVC interacted with the dynamics of the swirl stream shear layer vortices to reduce there strength. Transient scalar measurements indicated that there was a clear connection between the unsteady vortex pattern and the rate of mixing, resulting in bursts of high heat release and is therefore identified as one source of combustor oscillations. Future fuel injector designs need to pay close attention to these unsteady features in selecting swirl number and internal geometry parameters.

Current lean-premixed fuel injector designs have shown great potential in terms of reducing emissions of pollutants, but such designs are susceptible to combustion instabilities in which aerodynamic instability plays a major role and also has an effect on mixing of air and fuel. In comparison to prototype testing with combustors running in operating conditions, computational approaches such as Large Eddy Simulations (LES) offer a much more cost-effective alternative in the design stage. However, computational models employed by LES require validation by experimental data. This is one of the main motivations behind the present experimental study. Combined particle image velocimetry (PIV) and planar laser induced fluorescence (PLIF) instrumentation allowed simultaneous measurements of velocity vector and a conserved scalar introduced into the fuel stream. The results show that the inner swirl shear layer features two pairs of vortices, which draw high concentration fuel mixture from the central jet into the swirl stream and causes it to rotate in their wakes. Such periodic entrainment also occurs with the characteristic frequencies of the vortices. This has clear implications for temporal variations in fuel/air ratio in a combusting flow; these bursts of mixing, and hence heat release, could be a possible cause of mixing-induced pressure oscillation in combusting tests. For the first time in such a flow, all 3 components of the turbulent scalar flux were available for validation of LES-based predictions. A careful assessment of experimental errors, particularly the error associated with spatial filtering, was carried out. Comparison of LES predictions with experimental data showed very good agreement for both 1st and 2nd moment statistics, as well as spectra and scalar pdfs. It is particularly noteworthy that comparison between LES computed and measured scalar fluxes was very good; this represents successful validation of the simple (constant Schmidt number) SGS model used for this complex and practically important fuel injector flow. In addition to providing benchmark data for the validation of LES predictions, a new experimental technique has been developed that is capable of providing spatially resolved residence time data. Residence times of combustors have commonly been used to help understand NOx emissions and can also contribute to combustion instabilities. Both the time mean velocity and turbulence fields are important to the residence time, but determining the residence time via analysis of a measured velocity field is difficult due to the inherent unsteadiness and the three dimensional nature of a high-Re swirling flow. A more direct approach to measure residence time is reported here that examines the dynamic response of fuel concentration to a sudden cutoff in the fuel injection. Residence time measurement was mainly taken using a time-resolved PLIF technique, but a second camera for PIV was added to check that the step change does not alter the velocity field and the spectral content of the coherent structures. Characteristic timescales evaluated from the measurements are referred to as convection and half-life times: The former describes the time delay from a fuel injector exit reference point to a downstream point of interest, and the latter describes the rate of decay once the effect of the reduced scalar concentration at the injection source has been transported to the point of interest. Residence time is often defined as the time taken for a conserved scalar to reduce to half its initial value after injection is stopped: this is equivalent to the sum of the convection time and the half-life values. The technique was applied to a high-swirl fuel injector typical of that found in combustor applications. Two test cases have been studied: with central jet (with-jet) and without central jet (no-jet). It was found that the relatively unstable central recirculation zone of the no-jet case resulted in increased transport of fuel into the central region that is dominated by a precessing vortex core, where long half-life times are also found. Based on this, it was inferred that the no-jet case may be more prone to NOx production. The technique is described here for a single-phase isothermal flow field, but with consideration, it could be extended to studying reacting flows to provide more insight into important mixing phenomena and relevant timescales.

No description available.

This thesis aims at improving our knowledge on swirl combustors. The work presented here is based on Large Eddy Simulations (LES) coupled to an advanced combustion model: the Conditional Moment Closure (CMC). Numerical predictions have been systematically compared and validated with detailed experimental datasets. In order to analyze further the physics underlying the large numerical datasets, Proper Orthogonal Decomposition (POD) has also been used throughout the thesis. Various aspects of the aerodynamics of swirling flames are investigated, such as precession or vortex formation caused by flow oscillations, as well as various combustion aspects such as localized extinctions and flame lift-off. All the above affect flame stabilization in different ways and are explored through focused simulations. The first study investigates isothermal air flows behind an enclosed bluff body, with the incoming flow being pulsated. These flows have strong similarities to flows found in combustors experiencing self-excited oscillations and can therefore be considered as canonical problems. At high enough forcing frequencies, double ring vortices are shed from the air pipe exit. Various harmonics of the pulsating frequency are observed in the spectra and their relation with the vortex shedding is investigated through POD. The second study explores the structure of the Delft III piloted turbulent non-premixed flame. The simple configuration allows to analyze further key combustion aspects of combustors, with further insights provided on the dynamics of localized extinctions and re-ignition, as well as the pollutants emissions. The third study presents a comprehensive analysis of the aerodynamics of swirl flows based on the TECFLAM confined non-premixed S09c configuration. A periodic component inside the air inlet pipe and around the central bluff body is observed, for both the inert and reactive flows. POD shows that these flow oscillations are due to single and double helical vortices, similar to Precessing Vortex Cores (PVC), that develop inside the air inlet pipe and whose axes rotate around the burner. The combustion process is found to affect the swirl flow aerodynamics. Finally, the fourth study investigates the TECFLAM configuration again, but here attention is given to the flame lift-off evident in experiments and reproduced by the LES-CMC formulation. The stabilization process and the pollutants emission of the flame are investigated in detail.

[ES] El principal desafío en los motores turbina de gas empleados en aviación reside en aumentar la eficiencia del ciclo termodinámico manteniendo las emisiones contaminantes por debajo de las rigurosas restricciones. Ésto ha conllevado la necesidad de diseñar nuevas estrategias de inyección/combustión que operan en puntos de operación peligrosos por su cercanía al límite inferior de apagado de llama. En este contexto, el concepto Lean Direct Injection (LDI) ha emergido como una tecnología prometedora a la hora de reducir los óxidos de nitrógeno (NOx) emitidos por las plantas propulsoras de los aviones de nueva generación. En este contexto, la presente tesis tiene como objetivos contribuir al conocimiento de los mecanismos físicos que rigen el comportamiento de un quemador LDI y proporcionar herramientas de análisis para una profunda caracterización de las complejas estructuras de flujo de turbulento generadas en el interior de la cámara de combustión. Para ello, se ha desarrollado una metodología numérica basada en CFD capaz de modelar el flujo bifásico no reactivo en el interior de un quemador LDI académico mediante enfoques de turbulencia U-RANS y LES en un marco Euleriano-Lagrangiano. La resolución numérica de este problema multi-escala se aborda mediante la descripción completa del flujo a lo largo de todos los elementos que constituyen la maqueta experimental, incluyendo su paso por el swirler y entrada a la cámara de combustión. Ésto se lleva a cabo través de dos códigos CFD que involucran dos estrategias de mallado diferentes: una basada en algoritmos de generación y refinamiento automático de la malla (AMR) a través de CONVERGE y otra técnica de mallado estático más tradicional mediante OpenFOAM. Por un lado, se ha definido una metodología para obtener una estrategia de mallado óptima mediante el uso del AMR y se han explotado sus beneficios frente a los enfoques tradicionales de malla estática. De esta forma, se ha demostrado que la aplicabilidad de las herramientas de control de malla disponibles en CONVERGE como el refinamiento fijo (fixed embedding) y el AMR son una opción muy interesante para afrontar este tipo de problemas multi-escala. Los resultados destacan una optimización del uso de los recursos computacionales y una mayor precisión en las simulaciones realizadas con la metodología presentada. Por otro lado, el uso de herramientas CFD se ha combinado con la aplicación de técnicas de descomposición modal avanzadas (Proper Orthogonal Decomposition and Dynamic Mode Decomposition). La identificación numérica de los principales modos acústicos en la cámara de combustión ha demostrado el potencial de estas herramientas al permitir caracterizar las estructuras de flujo coherentes generadas como consecuencia de la rotura de los vórtices (VBB) y de los chorros fuertemente torbellinados presentes en el quemador LDI. Además, la implementación de estos procedimientos matemáticos ha permitido tanto recuperar información sobre las características de la dinámica de flujo como proporcionar un enfoque sistemático para identificar los principales mecanismos que sustentan las inestabilidades en la cámara de combustión. Finalmente, la metodología validada ha sido explotada a través de un Diseño de Experimentos (DoE) para cuantificar la influencia de los factores críticos de diseño en el flujo no reactivo. De esta manera, se ha evaluado la contribución individual de algunos parámetros funcionales (el número de palas del swirler, el ángulo de dichas palas, el ancho de la cámara de combustión y la posición axial del orificio del inyector) en los patrones del campo fluido, la distribución del tamaño de gotas del combustible líquido y la aparición de inestabilidades en la cámara de combustión a través de una matriz ortogonal L9 de Taguchi. Este estudio estadístico supone un punto de partida para posteriores estudios de inyección, atomización y combus / [CA] El principal desafiament als motors turbina de gas utilitzats a la aviació resideix en augmentar l'eficiència del cicle termodinàmic mantenint les emissions contaminants per davall de les rigoroses restriccions. Aquest fet comporta la necessitat de dissenyar noves estratègies d'injecció/combustió que radiquen en punts d'operació perillosos per la seva aproximació al límit inferior d'apagat de flama. En aquest context, el concepte Lean Direct Injection (LDI) sorgeix com a eina innovadora a l'hora de reduir els òxids de nitrogen (NOx) emesos per les plantes propulsores dels avions de nova generació. Sota aquest context, aquesta tesis té com a objectius contribuir al coneixement dels mecanismes físics que regeixen el comportament d'un cremador LDI i proporcionar ferramentes d'anàlisi per a una profunda caracterització de les complexes estructures de flux turbulent generades a l'interior de la càmera de combustió. Per tal de dur-ho a terme s'ha desenvolupat una metodología numèrica basada en CFD capaç de modelar el flux bifàsic no reactiu a l'interior d'un cremador LDI acadèmic mitjançant els enfocaments de turbulència U-RANS i LES en un marc Eulerià-Lagrangià. La resolució numèrica d'aquest problema multiescala s'aborda mitjançant la resolució completa del flux al llarg de tots els elements que constitueixen la maqueta experimental, incloent el seu pas pel swirler i l'entrada a la càmera de combustió. Açò es duu a terme a través de dos codis CFD que involucren estratègies de mallat diferents: una basada en la generación automàtica de la malla i en l'algoritme de refinament adaptatiu (AMR) amb CONVERGE i l'altra que es basa en una tècnica de mallat estàtic més tradicional amb OpenFOAM. D'una banda, s'ha definit una metodologia per tal d'obtindre una estrategia de mallat òptima mitjançant l'ús de l'AMR i s'han explotat els seus beneficis front als enfocaments tradicionals de malla estàtica. D'aquesta forma, s'ha demostrat que l'aplicabilitat de les ferramente de control de malla disponibles en CONVERGE com el refinament fixe (fixed embedding) i l'AMR són una opció molt interessant per tal d'afrontar aquest tipus de problemes multiescala. Els resultats destaquen una optimització de l'ús dels recursos computacionals i una major precisió en les simulacions realitzades amb la metodologia presentada. D'altra banda, l'ús d'eines CFD s'ha combinat amb l'aplicació de tècniques de descomposició modal avançades (Proper Orthogonal Decomposition and Dynamic Mode Decomposition). La identificació numèrica dels principals modes acústics a la càmera de combustió ha demostrat el potencial d'aquestes ferramentes al permetre caracteritzar les estructures de flux coherents generades com a conseqüència del trencament dels vòrtex (VBB) i dels raigs fortament arremolinats presents al cremador LDI. A més, la implantació d'estos procediments matemàtics ha permès recuperar informació sobre les característiques de la dinàmica del flux i proporcionar un enfocament sistemàtic per tal d'identificar els principals mecanismes que sustenten les inestabilitats a la càmera de combustió. Finalment, la metodologia validada ha sigut explotada a traves d'un Diseny d'Experiments (DoE) per tal de quantificar la influència dels factors crítics de disseny en el flux no reactiu. D'aquesta manera, s'ha avaluat la contribución individual d'alguns paràmetres funcionals (el nombre de pales del swirler, l'angle de les pales, l'amplada de la càmera de combustió i la posició axial de l'orifici de l'injector) en els patrons del camp fluid, la distribució de la mida de gotes del combustible líquid i l'aparició d'inestabilitats en la càmera de combustió mitjançant una matriu ortogonal L9 de Taguchi. Aquest estudi estadístic és un bon punt de partida per a futurs estudis de injecció, atomització i combustió en cremadors LDI. / [EN] Aeronautical gas turbine engines present the main challenge of increasing the efficiency of the cycle while keeping the pollutant emissions below stringent restrictions. This has led to the design of new injection-combustion strategies working on more risky and problematic operating points such as those close to the lean extinction limit. In this context, the Lean Direct Injection (LDI) concept has emerged as a promising technology to reduce oxides of nitrogen (NOx) for next-generation aircraft power plants In this context, this thesis aims at contributing to the knowledge of the governing physical mechanisms within an LDI burner and to provide analysis tools for a deep characterisation of such complex flows. In order to do so, a numerical CFD methodology capable of reliably modelling the 2-phase nonreacting flow in an academic LDI burner has been developed in an Eulerian-Lagrangian framework, using the U-RANS and LES turbulence approaches. The LDI combustor taken as a reference to carry out the investigation is the laboratory-scale swirled-stabilised CORIA Spray Burner. The multi-scale problem is addressed by solving the complete inlet flow path through the swirl vanes and the combustor through two different CFD codes involving two different meshing strategies: an automatic mesh generation with adaptive mesh refinement (AMR) algorithm through CONVERGE and a more traditional static meshing technique in OpenFOAM. On the one hand, a methodology to obtain an optimal mesh strategy using AMR has been defined, and its benefits against traditional fixed mesh approaches have been exploited. In this way, the applicability of grid control tools available in CONVERGE such as fixed embedding and AMR has been demonstrated to be an interesting option to face this type of multi-scale problem. The results highlight an optimisation of the use of the computational resources and better accuracy in the simulations carried out with the presented methodology. On the other hand, the use of CFD tools has been combined with the application of systematic advanced modal decomposition techniques (i.e., Proper Orthogonal Decomposition and Dynamic Mode Decomposition). The numerical identification of the main acoustic modes in the chamber have proved their potential when studying the characteristics of the most powerful coherent flow structures of strongly swirled jets in a LDI burner undergoing vortex breakdown (VBB). Besides, the implementation of these mathematical procedures has allowed both retrieving information about the flow dynamics features and providing a systematic approach to identify the main mechanisms that sustain instabilities in the combustor. Last, this analysis has also allowed identifying some key features of swirl spray systems such as the complex pulsating, intermittent and cyclical spatial patterns related to the Precessing Vortex Core (PVC). Finally, the validated methodology is exploited through a Design of Experiments (DoE) to quantify the influence of critical design factors on the non-reacting flow. In this way, the individual contribution of some functional parameters (namely the number of swirler vanes, the swirler vane angle, the combustion chamber width and the axial position of the nozzle tip) into both the flow field pattern, the spray size distribution and the occurrence of instabilities in the combustion chamber are evaluated throughout a Taguchi's orthogonal array L9. Such a statistical study has supposed a good starting point for subsequent studies of injection, atomisation and combustion on LDI burners. / Belmar Gil, M. (2020). Computational study on the non-reacting flow in Lean Direct Injection gas turbine combustors through Eulerian-Lagrangian Large-Eddy Simulations [Tesis doctoral]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/159882