Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Este trabalho apresenta uma família de novos métodos numéricos euleriano-lagrangeanos localmente conservativos para leis de conservação hiperbólicas escalares. Estes métodos não utilizam soluções analíticas de problemas de Riemann e são bastante precisos na captura de saltos nas soluções. Estes métodos foram introduzidos, implementados computacionalmente e testados para leis de conservação em uma e duas dimensões espaciais. Foram consideradas as equações de Burgers e Buckley-Leverett. Nossos experimentos numéricos indicaram que os métodos são pouco difusivos e que as soluções não apresentam oscilações espúrias.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Este trabalho apresenta uma família de novos métodos numéricos euleriano-lagrangeanos localmente conservativos para leis de conservação hiperbólicas escalares. Estes métodos não utilizam soluções analíticas de problemas de Riemann e são bastante precisos na captura de saltos nas soluções. Estes métodos foram introduzidos, implementados computacionalmente e testados para leis de conservação em uma e duas dimensões espaciais. Foram consideradas as equações de Burgers e Buckley-Leverett. Nossos experimentos numéricos indicaram que os métodos são pouco difusivos e que as soluções não apresentam oscilações espúrias.

In power plants using large utility coal-fired boilers for generation of electricity the coal is pulverised in coal mills and then it has to be pneumatically transported and distributed to a larger number of burners (e.g. 30-40) circumferentially arranged in several rows around the burning chamber of the boiler. Besides the large pipework flow splitting devices are necessary for distribution of an equal amount of pulverised fuel (PF) to each of the burners. So called trifurcators (without inner fittings or guiding vanes) and ''riffle'' type bifurcators are commonly used to split the gas-coal particle flow into two or three pipes/channels with an equal amount of PF mass flow rate in each outflow cross section of the flow splitting device. These PF flow splitting devices are subject of a number of problems. First of all an uneven distribution of PF over the burners of a large utility boiler leads to operational and maintenance problems, increased level of unburned carbon and higher rates of NOX emissions. Maldistribution of fuel between burners caused by non uniform concentration of the PF (particle roping) in pipe and channel bends prior to flow splitting devices leads to uncontrolled differences in the fuel to air ratio between burners. This results in localised regions in the furnace which are fuel rich, where insufficient air causes incomplete combustion of the fuel. Other regions in the furnace become fuel lean, forming high local concentrations of NOX due to the high local concentrations of O2. Otherwise PF maldistribution can impact on power plant maintenance in terms of uneven wear on PF pipework, flow splitters as well as the effects on boiler panels (PF deposition, corrosion, slagging). In order to address these problems in establishing uniform PF distribution over the outlet cross sections of flow splitting devices in the pipework of coal-fired power plants the present paper deals with numerical prediction and analysis of the complex gas and coal particle (PF) flow through trifurcators and ''riffle'' type bifurcators. The numerical investigation is based on a 3-dimensional Eulerian- Lagrangian approach (MISTRAL/PartFlow-3D) developed by Frank et al. The numerical method is capable to predict isothermal, incompressible, steady gas- particle flows in 3-dimensional, geometrically complex flow geometries using boundary fitted, block-structured, numerical grids. Due to the very high numerical effort of the investigated gas-particle flows the numerical approach has been developed with special emphasis on efficient parallel computing on clusters of workstations or other high performance computing architectures. Besides the aerodynamically interaction between the carrier fluid phase and the PF particles the gas-particle flow is mainly influenced by particle-wall interactions with the outer wall boundaries and the inner fittings and guiding vanes of the investigated flow splitting devices. In order to allow accurate quantitative prediction of the motion of the disperse phase the numerical model requires detailed information about the particle-wall collision process. In commonly used physical models of the particle-wall interaction this is the knowledge or experimental prediction of the restitution coefficients (dynamic friction coefficient, coefficient of restitution) for the used combination of particle and wall material, e.g. PF particles on steel. In the present investigation these parameters of the particle-wall interaction model have been obtained from special experiments in two test facilities. Basic experiments to clarify the details of the particle-wall interaction process were made in a test facility with a spherical disk accelerator. This test facility furthermore provides the opportunity to investigate the bouncing process under normal pressure as well as under vacuum conditions, thus excluding aerodynamically influences on the motion of small particles in the near vicinity of solid wall surfaces (especially under small angles of attack). In this experiments spherical glass beads were used as particle material. In a second test facility we have investigated the real impact of non-spherical pulverised fuel particles on a steel/ceramic target. In this experiments PF particles were accelerated by an injector using inert gas like e.g. CO2 or N2 as the carrier phase in order to avoid dust explosion hazards. The obtained data for the particle-wall collision models were compared to those obtained for glass spheres, where bouncing models are proofed to be valid. Furthermore the second test facility was used to obtain particle erosion rates for PF particles on steel targets as a function of impact angles and velocities. The results of experimental investigations has been incorporated into the numerical model. Hereafter the numerical approach MISTRAL/PartFlow-3D has been applied to the PF flow through a ''riffle'' type bifurcator. Using ICEM/CFD-Hexa as grid generator a numerical mesh with approximately 4 million grid cells has been designed for approximation of the complex geometry of the flow splitting device with all its interior fittings and guiding vanes. Based on a predicted gas flow field a large number of PF particles are tracked throughout the flow geometry of the flow-splitter. Besides mean quantities of the particle flow field like e.g. local particle concentrations, mean particle velocities, distribution of mean particle diameter, etc. it is now possible to obtain information about particle erosion on riffle plates and guiding vanes of the flow splitting device. Furthermore the influence of different roping patterns in front of the flow splitter on the uniformness of PF mass flow rate splitting after the bifurcator has been investigated numerically. Results show the efficient operation of the investigated bifurcator in absence of particle roping, this means under conditions of an uniform PF particle concentration distribution in the inflow cross section of the bifurcator. If particle roping occurs and particle concentration differs over the pipe cross section in front of the bifurcator the equal PF particle mass flow rate splitting can be strongly deteriorated in dependence on the location and intensity of the particle rope or particle concentration irregularities. The presented results show the importance of further development of efficient rope splitting devices for applications in coal-fired power plants. Numerical analysis can be used as an efficient tool for their investigation and further optimisation under various operating and flow conditions.

The transport or the separation of solid particles or droplets suspended in a fluid flow is a common task in mechanical and process engineering. To improve machinery and physical processes (e.g. for coal combustion, reduction of NO_x and soot) an optimization of complex phenomena by simulation applying the fundamental conservation equations is required. Fluid-particle flows are characterized by the ratio of density of the two phases gamma=rho_P/rho_F, by the Stokes number St=tau_P/tau_F and by the loading in terms of void and mass fraction. Those numbers (Stokes number, gamma) define the flow regime and which relevant forces are acting on the particle. Dependent on the geometrical configuration the particle-wall interaction might have a heavy impact on the mean flow structure. The occurrence of particle-particle collisions becomes also more and more important with the increase of the local void fraction of the particulate phase. With increase of the particle loading the interaction with the fluid phase can not been neglected and 2-way or even 4-way coupling between the continous and disperse phases has to be taken into account. For dilute to moderate dense particle flows the Euler-Lagrange method is capable to resolve the main flow mechanism. An accurate computation needs unfortunately a high number of numerical particles (1,...,10^7) to get the reliable statistics for the underlying modelling correlations. Due to the fact that a Lagrangian algorithm cannot be vectorized for complex meshes the only way to finish those simulations in a reasonable time is the parallization applying the message passing paradigma. Frank et al. describes the basic ideas for a parallel Eulererian-Lagrangian solver, which uses multigrid for acceleration of the flow equations. The performance figures are quite good, though only steady problems are tackled. The presented paper is aimed to the numerical prediction of time-dependend fluid-particle flows using the simultanous particle tracking approach based on the Eulerian-Lagrangian and the particle-source-in-cell (PSI-Cell) approach. It is shown in the paper that for the unsteady flow prediction efficiency and load balancing of the parallel numerical simulation is an even more pronounced problem in comparison with the steady flow calculations, because the time steps for the time integration along one particle trajectory are very small per one time step of fluid flow integration and so the floating point workload on a single processor node is usualy rather low. Much time is spent for communication and waiting time of the processors, because for cold flow particle convection not very extensive calculations are necessary. One remedy might be a highspeed switch like Myrinet or Dolphin PCI/SCI (500 MByte/s), which could balance the relative high floating point performance of INTEL PIII processors and the weak capacity of the Fast-Ethernet communication network (100 Mbit/s) of the Chemnitz Linux Cluster (CLIC) used for the presented calculations. Corresponding to the discussed examples calculation times and parallel performance will be presented. Another point is the communication of many small packages, which should be summed up to bigger messages, because each message requires a startup time independently of its size. Summarising the potential of such a parallel algorithm, it will be shown that a Beowulf-type cluster computer is a highly competitve alternative to the classical main frame computer for the investigated Eulerian-Lagrangian simultanous particle tracking approach.

[ES] El proceso de atomización primaria es el mecanismo por el cual una vena líquida se disgrega en un ambiente gaseoso. Este proceso está presente en muchas aplicaciones de ingeniería realizando diferentes tareas. En ocasiones es un paso previo antes de ser quemado, como en la industria energética o de propulsión, donde el objetivo es extraer la energía específica del líquido. En otros sectores, como el revestimiento o la extinción de incendios, el objetivo es maximizar el área cubierta por el chorro. Sin embargo, aunque la atomización es una parte fundamental de varios procesos industriales, está lejos de comprenderse por completo. El proceso de atomización es una mezcla de fenómenos de interacción gas-líquido dentro de un campo turbulento que tiene lugar en el campo cercano, que es la región más densa del chorro. Cuando se trata de arrojar luz sobre el proceso de atomización primaria, el problema principal es la falta de teorías físicas definitivas capaces de vincular los complejos eventos de ruptura con la turbulencia. El principal obstáculo que impide investigar el proceso de atomización primaria es la incapacidad de las técnicas ópticas clásicas para proporcionar información de la región densa del chorro. Solo en los últimos años, las nuevas técnicas basadas en rayos X podrían proporcionar nueva información sobre las características de la atomización cerca de la salida de la tobera. Esto también afecta a los modelos computacionales de atomización primaria que, al no disponer de información experimental sobre la región densa, requieren una calibración precisa de sus constantes para proporcionar resultados fiables en el campo lejano. Esta tesis se centra en mejorar el conocimiento del proceso de atomización primaria, especialmente en cómo las condiciones de inyección afectan el desarrollo del chorro en el campo cercano desde dos puntos de vista diferentes. Por un lado, con un enfoque computacional usando Direct Numerical Simulations y, por otro lado, experimentalmente usando Near-Field Microscopy. El estudio computacional se centra en variar los números de Reynolds y Weber de inyección. Los resultados muestran que aumentar el número de Reynolds mejora la desintegración del líquido, mostrando un aumento de las gotas generadas y una nube de gotas más fina. Sin embargo, la falta de un perfil turbulento de flujo de entrada completamente desarrollado conduce a comportamientos inesperados en la longitud de ruptura de la vena líquida que también aumenta con el número de Reynolds. El número de gotas también aumenta cuando aumenta el número de Weber, pero los tamaños característicos de las gotas siguen siendo los mismos. La longitud de ruptura no varía, lo que sugiere que las variaciones de la tensión superficial afectan la ruptura de las gotas y los ligamentos, pero no la desintegración del núcleo líquido en sí. Con los resultados obtenidos de ambos estudios, se propone un modelo fenomenológico que predice la distribución del tamaño de gota en función de las condiciones de inyección. Además, también se ha estudiado el efecto de usar toberas elípticas. Se ha obtenido que el número de gotas detectadas aumenta en comparación con el chorro redondo manteniendo ángulos de apertura del chorro similares. Sin embargo, cuando se utilizan toberas extremadamente excéntricas, la disminución de la turbulencia del flujo de entrada contrarresta los beneficios de este tipo de inyectores. En cuanto al análisis experimental, usar Near-Field Microscopy permite magnificar la región densa y analizar las características macroscópicas del chorro. Por lo tanto, se varían las presiones de inyección y descarga, centrándose en el ángulo de apertura del chorro. Se observa el aumento esperado en el ángulo al aumentar tanto la presión de inyección como la de descarga. Sin embargo, adicionalmente, se realiza un análisis de las perturbaciones del contorno del chorro, concluyendo que, al aumentar la presión de inyección, y por lo tanto la turbulencia del flujo de / [CA] El procés d'atomització primària és el mecanisme pel qual una vena líquida es disgrega en un ambient gasós. Aquest procés és present en moltes aplicacions d'enginyeria fent diferents tasques. De vegades és un pas previ abans de ser cremat, com ara en la indústria energètica o de propulsió, on l'objectiu és extraure l'energia específica del líquid. En altres sectors, com ara el revestiment o l'extinció d'incendis, l'objectiu és maximitzar l'àrea coberta pel doll. No obstant això, tot i que l'atomització és una part fonamental de diversos processos industrials, està lluny de comprendre's per complet. El procés d'atomització és una barreja de fenòmens d'interacció gas-líquid dins d'un camp turbulent que té lloc en el camp pròxim, que és la regió més densa del doll. Quan es tracta de donar llum sobre el procés d'atomització primària, el problema principal és la falta de teories físiques definitives capaces de vincular els complexos esdeveniments de trencament amb la turbulència. El principal obstacle que impedeix investigar el procés d'atomització primària és la incapacitat de les tècniques òptiques clàssiques per a proporcionar informació de la regió densa del doll. Només en els últims anys, les noves tècniques basades en raigs X podrien proporcionar nova informació sobre les característiques de l'atomització prop de l'eixida de la tovera. Això també afecta els models computacionals d'atomització primària que, en no disposar d'informació experimental sobre la regió densa, requereixen un calibratge precís de les seues constants per a proporcionar resultats fiables en el camp llunyà. Aquesta tesi se centra a millorar el coneixement del procés d'atomització primària, especialment en com les condicions d'injecció afecten el desenvolupament del doll en el camp pròxim des de dos punts de vista diferents. D'una banda, amb un enfocament computacional usant Direct Numerical Simulations i, d'altra banda, experimentalment usant Near-Field Microscopy. L'estudi computacional se centra a variar els nombres de Reynolds i Weber d'injecció. Els resultats mostren que augmentar el nombre de Reynolds millora la desintegració del líquid, tot mostrant un augment de les gotes generades i un núvol de gotes més fi. No obstant això, la falta d'un perfil turbulent de flux d'entrada completament desenvolupat condueix a comportaments inesperats en la longitud de ruptura de la vena líquida que també augmenta amb el nombre de Reynolds. El nombre de gotes també augmenta quan creix el nombre de Weber, però les grandàries característiques de les gotes continuen sent les mateixes. La longitud de ruptura no varia, la qual cosa suggereix que les variacions de la tensió superficial afecten la ruptura de les gotes i els lligaments, però no la desintegració del nucli líquid en ell mateix. Amb els resultats obtinguts de tots dos estudis, es proposa un model fenomenològic que prediu la distribució de la grandària de gota en funció de les condicions d'injecció. A més, també s'ha estudiat l'efecte d'usar toveres el·líptiques. S'ha obtingut que el nombre de gotes detectades augmenta en comparació amb el doll redó tot mantenint angles d'obertura del doll similars. No obstant això, quan s'utilitzen toveres extremadament excèntriques, la disminució de la turbulència del flux d'entrada contraresta els beneficis d'aquesta mena d'injectors. Quant a l'anàlisi experimental, usar Near-Field Microscopy permet magnificar la regió densa i analitzar les característiques macroscòpiques del doll. Per tant, es varien les pressions d'injecció i descàrrega, tot centrant-se en l'angle d'obertura del doll. S'observa l'augment esperat en l'angle en augmentar tant la pressió d'injecció com la de descàrrega. No obstant això, addicionalment, es realitza una anàlisi de les pertorbacions del contorn del doll i es conclou que en augmentar la pressió d'injecció, i per tant la turbulència del flux d'entrada, augmenten les pertorbacions en el contorn del ruixat, especialment a pressions de descàrrega mé / [EN] The primary atomisation process is the mechanism by which a liquid vein breaks into droplets in a gaseous ambient. This process is present in many engineering applications accomplishing different tasks. Sometimes it is a previous step before being burned, as in the energy or propulsion industry, where the objective is to extract the specific energy of the liquid. In other sectors, such as the coating or fire extinction, the objective is to maximise the area covered by the droplet cloud. However, although atomisation is a fundamental part of several industrial processes, it is far from fully understood. The atomisation process is a mixture of gas-liquid interaction phenomena within a turbulent field that takes place in the near-field, which is the denser region of the spray. When trying to shed light on the primary atomisation process, the main issue is the lack of definitive physical theories able to link the complex breakup events and the turbulence. The principal impediment that prevents the investigation from breaking through the atomisation process is the inability of the classic optical techniques to provide information from the dense region of the spray. Only in the last years, newer techniques based on X-Ray could provide new information on spray characteristics near the nozzle outlet. This also affects the computational primary atomisation models that, as there is no available experimental information on the dense region, require an accurate calibration of their constants to provide reliable results on the far-field. This thesis focuses on improving the knowledge of the primary atomisation process, especially on how the injection conditions affect the spray development in the near field from two different standpoints. On the one hand, with a computational approach using Direct Numerical Simulations and on the other hand, experimentally using Near-Field Microscopy. The computational study is focused on varying the inflow Reynolds and Weber numbers. Results show that increasing the Reynolds number improves the liquid disintegration, exhibiting an increase of generated droplets and a finer droplet cloud. However, the lack of a fully developed inflow turbulent profile leads to characteristic behaviours on the breakup length of the spray that also increases with the Reynolds number. The number of droplets increases when the Weber number increases, but the characteristic droplet sizes remain the same. The breakup length does not vary, suggesting that the surface tension variations affect the droplet and ligament breakup but not the core disintegration itself. With the results obtained from both studies, a phenomenological model is proposed to predict the droplet size distribution depending on the injection conditions. Additionally, using elliptical nozzles, the number of detected droplets increases compared with the round spray and maintain similar spray apertures. However, when using extremely eccentric nozzles, the inflow turbulence decrease counteracts the elliptical sprays' benefits. Regarding the experimental analysis, the Near-Field Microscopy magnifies the dense region and analyses the macroscopic features on the spray. So the injection and discharge pressure are varied, and the spotlight is put on the spray angle. The expected increase in the spray angle when increasing both the injection and discharge pressure is observed. Nevertheless, additionally, an analysis of the spray contour perturbations is performed, concluding that increasing the injection pressure, and thus the inflow turbulence, increases the perturbations on the spray contour, especially at lower discharge pressures. / González Montero, LA. (2022). Computational and Experimental Study of the Primary Atomisation Process under Different Injection Conditions [Tesis doctoral]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/190635

[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