This work deals with study of swirling flows where the spiral vortical structure appears. The main relation is to flow seen in the draft tube cone of hydraulic turbines operated out of the design point (i.e. best efficiency point). In this cases large coherent vortex structure (vortex rope) appears and consequently high pressure pulsations are propagated to the whole machine system leading to possible restriction of turbine operation. This flow features are consequence of flow instability called vortex breakdown in case of Francis turbine operated at part load (flow rate lower than optimal one). The present study is carried out using simplified device of swirl generator in order to access similar flow conditions as can be found in real hydraulic turbines. Both the dynamic and dissipation effect of spiral vortex breakdown are investigated. The first part of thesis deals with spiral form of vortex breakdown. The experimentally measured velocity profiles (LDA) and wall static pressures are correlated with numerical simulations carried out using open-source CFD package OpenFOAM 2.2.2. The high speed camera recording of cavitating vortex core is used to obtain image ensemble for further post-processing. The dissipation effect of spiral vortex structure is in detail discussed based on computed flow fields. The second part of thesis is dedicated to the application of POD decomposition to the study of spatio-temporal features of spiral vortex dynamics. Firstly the POD is applied to the both the experimentally obtained image ensemble of cavitating vortex and numerically computed static pressure fields. Secondly the comprehensive analysis of spiral vortex mitigation effect by the axial water jet is analyzed. The collaborative study employing the swirl generator apparatus designed by the researchers from Politehnica University of Timisoara in Romania is performed and changes in spatio-temporal vortex dynamic are studied. In this study the numerical data (in a form of three-dimensional pressure and velocity fields) are obtained using commercial CFD software ANSYS Fluent R14.

The natural unsteadiness in the post-breakdown flowfield of a 75° sweep delta wing at 40° angle of attack was studied with dual and single point hot-wire anemometry in the Engineering Science and Mechanics (ESM) Wind Tunnel at a Reynolds number Re = 210,000. Data were taken in five crossflow planes surrounding the wing's trailing edge. Results showed a dominant narrowband Strouhal frequency of St = 1.5 covering approximately 80% of the area with lower-intensity broadband secondary frequencies over 15% of that region. Cross-correlations between a fixed and traversing wire were calculated and phase and coherences mapped to determine the convection speed and trajectory of the helical mode instability. High-speed Particle Image Velocimetry (PIV) was conducted over a 75° sweep delta wing at 40° angle of attack in the ESM Water Tunnel II at Re = 45,000. Data were taken along the axis of the vortex in the breakdown flowfield at a speed of 0.1% of the convective time scale of the flow. Animations of instantaneous streamlines and velocity vectors revealed the impression of a helically spiralling vortex core on the measurement plane. Spectral analysis of the PIV data showed reduced frequencies which confirmed those found with the single-point measurements made in the ESM Wind Tunnel. The effect of four novel control surfaces on the breakdown flowfield of the delta wing was studied with surface pressure measurements along the axis of the vortex in the ESM Wind Tunnel. The apex flap was found effective and delayed vortex breakdown by 8° for a fixed wing. The flowfield was characterized over the delta wing executing a pitch-up maneuver at a reduced frequency of 0.06. Surface pressure measurements were taken in the ESM Wind Tunnel and Laser Doppler Velocimetry (LDV) was employed in the ESM Water Tunnel I as both the unmodified wing and then the wing with an apex flap deployed at an optimal angle <font face="Symbol">b</font> = 15° executed the pitch-up. Both sets of data confirmed the hysteresis of the flowfield. The LDV data, taken in two crossflow planes throughout the maneuver, showed an asymmetric breakdown development. As a practical extension of the study of the breakdown wake flowfield, hot-wire measurements were made over an F/A-18 model to determine the spectral characteristics of the flowfield. Three-dimensional vortex interactions were investigated with helium bubble flow visualization in the VPI Stability Tunnel. / Ph. D.

Vortex breakdown is an important phenomenon observed in swirling flows involving the development of a stagnation point on the axis of the vortex followed by a region of recirculation when the swirl increases beyond a particular level. It has been studied extensively over past 50 years and various theories have been proposed to explain its various aspects. However, a single model explaining all the aspects together is yet to emerge. Numerical simulations of breakdown have been performed using a variety of grid-based as well as vortex methods. Vortex methods are a Lagrangian alternative to grid-based methods wherein the motion of the vorticity is determined by the local fluid velocity convection, with models for viscous effects when considered. The fluid velocity is obtained from the vorticity field. Only the rotational regions of the flow need to be considered leading to significant economy of computational effort for simulations of vorticity dominated flows, such as vortex breakdown. The inviscid vortex filament method has been used to simulate several aspects of the vortex breakdown phenomenon. The vortex filament method however, cannot easily simulate viscous effects. To simulate the viscous effects the viscous vortex particle method needs to be used. This work was intended to be a first step towards this end by initially evaluating the effectiveness of the inviscid version of the vortex particle method in simulating the breakdown phenomenon. The inviscid vortex particle method was found to satisfactorily simulate most qualitative aspects involved in the formation of vortex breakdown such as the retardation of axial velocity along centerline, radial swelling of the vortex core, formation of stagnation points, creation of azimuthal vorticity gradient from axial vorticity gradient and the turning of vortex lines along with the formation of a bubble-like structure with recirculating flow within. The effect of a wall placed adjacent to the vortex core was simulated by using image vortices. The wall was not found to influence the location of breakdown. However, the initiation of the spiral mode was found to occur earlier when a wall was present. For a quantitative assessment, a simulation of the experimental results of Faler and Leibovich (1978) was attempted. The simulation managed to predict the location of the breakdown and the extent of the bubble. The shape and height of the bubble obtained however were not in accord with the experimental observations. A single vortical cell was obtained in the interior of the bubble.

The efficient and enhanced mixing of heat and incoming reactants is achieved in modern gas turbine systems by employing swirling flows. This is realized by a low velocity region (internal recirculation zone -IRZ) zone resulting from vortex breakdown phenomenon. Besides, IRZ acts as effective flame holder/stabilization mode. Double concentric swirling jet is employed in plethora of industrial applications such as heat exchange, spray drying and combustion. As such, understanding essential features of vortex breakdown induced IRZ and its acoustic response in swirling flow/flame is important in thermo-acoustic instability studies. The key results of the present experimental investigation are discussed in four parts. In the first part, primary transition (sub-critical states) from a pre-vortex breakdown (Pre-VB) flow reversal to a fully-developed central toroidal recirculation zone (CTRZ) in a non-reacting, double-concentric swirling jet configuration is discussed when the swirl number is varied in the range 0.592 S 0.801. This transition proceeds with the formation of two intermediate, critical flow regimes. First, a partially-penetrated vortex breakdown bubble (VBB) is formed that indicates the first occurrence of an enclosed structure resulting in an opposed flow stagnation region. Second, a metastable transition structure is formed that marks the collapse of inner mixing vortices. In this study, the time-averaged topological changes in the coherent recirculation structures are discussed based on the non-dimensional modified Rossby number (Rom) which appears to describe the spreading of the zone of swirl influence in different flow regimes. The second part describes a secondary transition from an open-bubble type axisymmetric vortex breakdown (sub-critical states) to partially-open bubble mode (super-critical states) through an intermediate, critical regime of conical sheet formation for flow modes Rom ≤ 1 is discussed when the swirl number (S) is increased beyond 0.801. In the third part, amplitude dependent acoustic response of above mentioned sub and supercritical flow states is discussed. It was observed that the global acoustic response of the sub-critical VB states was fundamentally different from their corresponding super-critical modes. In particular, with a stepwise increase in excitation amplitude till a critical value, the sub-critical VB topology moved downstream and radially outward. Beyond a critical magnitude, the VB bubble transited back upstream and finally underwent radial shrinkage at the threshold excitation amplitude. On the other hand, the topology of the super-critical VB state continuously moved downstream and radially outwards and finally widened/fanned-out at threshold amplitude. In the final part, transition in time-averaged flame global flame structure is reported as a function of geometric swirl number. In particular, with a stepwise increase in swirl intensity, primary transition is depicted as a transformation from zero-swirl straight jet flame to lifted flame with blue base and finally to swirling seated flame. Further, a secondary transition is reported which consists of transformation from swirling seated flame to swirling flame with a conical tailpiece and finally to highly-swirled near blowout ultra-lean flame. For this purpose, CH* chemiluminescence imaging and 2D PIV in meridional planes were employed. Three baseline fuel flow rates through the central fuel injection pipe were considered. For each of the fuel flow cases (Ref), six different co-airflow rate settings (Rea) were employed. The geometric swirl number (SG) was increased in steps from zero till blowout for a particular fuel and co-airflow setting. A regime map (SG vs Rea) depicting different regions of flame stabilization were then drawn for each fuel flow case. The secondary transformation is explained on the basis of physical significance of Rom.

La thèse est une contribution à l’étude des tourbillons longitudinaux qui se développent sur la lunette arrière des véhicules automobiles, dans l’idée de provoquer leur éclatement afin de réduire la traînée aérodynamique. On conçoit tout d’abord un système d’acquisition dénommé A-SPIV, permettant de reconstruire le champ moyen de vitesse 3D d’un sillage turbulent à partir de plans stéréo-PIV acquis par translation du système caméras-laser, sans qu’il soit nécessaire de le recalibrer. En complément, on propose une méthode de reconstruction de la pression moyenne reposant sur les données A-SPIV et sur une mesure de la pression pariétale. L’ensemble forme un nouveau protocole expérimental, validé dans le sillage d’un corps d’Ahmed 25° et dont les résultats à haut Reynolds sont comparés à la littérature existante. L'analyse topologique des tourbillons longitudinaux suggère l’existence d’un éclatement tourbillonnaire spontané dans le sillage proche, malgré l’absence de point de stagnation. On démontre formellement l’existence de cet éclatement par deux critères théoriques qui considèrent ce phénomène, soit comme la conséquence d’une réorganisation de la vorticité, soit comme la résultante d’une accumulation, en un point critique, d’ondes inertielles se propageant le long du tourbillon. Les analyses sous-jacentes sont menées dans un repère cylindrique attaché à l’axe tourbillonnaire et prédisent une même position d’éclatement, en bon accord avec la position singulière issue des mesures A-SPIV. La thèse se conclut par une analyse de stabilité globale de l'écoulement moyen qui suggère un lien possible entre l’éclatement et une instabilité globale de l’écoulement tourbillonnaire. / This thesis is a contribution to the study of the longitudinal vortices developing in the near wake of ground vehicles, with the general purpose of reducing the aerodynamic drag by triggering the vortex breakdown phenomenon. We present a new data acquisition system called A-SPIV, allowing to reconstruct a 3D turbulent time-averaged velocity field from stereo-PIV planes measured by translation of the whole cameras-laser system, with no need to recalibrate. We also propose a method to reconstruct the mean pressure in the bulk from the A-SPIV data and from a dedicated wall static pressure measurement. This new overall experimental protocol is applied to a standard aerodynamic test-case, the 25° Ahmed body, all results being compared and validated at high turbulent Reynolds number against existing data from the literature. A thorough analysis of the longitudinal vortices suggests the occurrence of a spontaneous vortex breakdown in the near-wake, although there exist no stagnation point in the experimental data. Such vortex breakdown is therefore evidenced using two different theoretical criteria considering the phenomenon as the consequence of either a reorganization of the vorticity, or an accumulation of inertial waves propagating along the vortex core. The underlying analyses are carried out in a cylindrical system attached to the vortex axis and predict a single breakdown position, in good agreement with the singular position initially inferred from the A-SPIV data. The thesis ends with a global stability analysis of the turbulent mean flow suggesting a possible connection between the occurrence of vortex breakdown and a global instability of the longitudinal vortex.

[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 / TESIS