La simulation aux grandes échelles est devenue un outil d’analyse incontournable pour l’étude des écoulements turbulents dans des géométries complexes. Cependant, à cause de l’augmentation constante des ressources de calcul, le traitement des grandes quantités de données générées par les simulations hautement résolues est devenu un véritable défi qu’il n’est plus possible de relever avec des outils traditionnels. En mécanique des fluides numérique, cette problématique émergente soulève les mêmes questions que celles communément rencontrées en informatique avec des données massives. A ce sujet, certaines méthodes ont déjà été développées telles que le partitionnement et l’ordonnancement des données ou bien encore le traitement en parallèle mais restent insuffisantes pour les simulations numériques modernes. Ainsi, l’objectif de cette thèse est de proposer de nouveaux formalismes permettant de contourner le problème de volume de données en vue des futurs calculs exaflopiques que l’informatique devrait atteindre en 2020. A cette fin, une méthode massivement parallèle de co-traitement, adaptée au formalisme non-structuré, a été développée afin d’extraire les grandes structures des écoulements turbulents. Son principe consiste à introduire une série de grilles de plus en plus grossières réduisant ainsi la quantité de données à traiter tout en gardant intactes les structures cohérentes d’intérêt. Les données sont transférées d’une grille à une autre grâce à l’utilisation de filtres et de méthodes d’interpolation d’ordre élevé. L’efficacité de cette méthodologie a pu être démontrée en appliquant des techniques de décomposition modale lors de la simulation 3D d’une pale de turbine turbulente sur une grille de plusieurs milliards d’éléments. En outre, cette capacité à pouvoir gérer plusieurs niveaux de grilles au sein d’une simulation a été utilisée par la suite pour la mise en place de calculs basés sur une stratégie multi-niveaux. L’objectif de cette méthode est d’évaluer au cours du calcul les erreurs numériques et celles liées à la modélisation en simulant simultanément la même configuration pour deux résolutions différentes. Cette estimation de l’erreur est précieuse car elle permet de générer des grilles optimisées à travers la construction d’une mesure objective de la qualité des grilles. Ainsi, cette méthodologie de multi-résolution tente de limiter le coût de calcul de la simulation en minimisant les erreurs de modélisation en sous-maille, et a été appliquée avec succès à la simulation d’un écoulement turbulent autour d’un cylindre. / Large-Eddy Simulation (LES) has become a major tool for the analysis of highly turbulent flows in complex geometries. However, due to the steadily increase of computational resources, the amount of data generated by well-resolved numerical simulations is such that it has become very challenging to manage them with traditional data processing tools. In Computational Fluid Dynamics (CFD), this emerging problematic leads to the same "Big Data" challenges as in the computer science field. Some techniques have already been developed such as data partitioning and ordering or parallel processing but still remain insufficient for modern numerical simulations. Hence, the objective of this work is to propose new processing formalisms to circumvent the data volume issue for the future 2020 exa-scale computing objectives. To this aim, a massively parallel co-processing method, suited for complex geometries, was developed in order to extract large-scale features in turbulent flows. The principle of the method is to introduce a series of coarser nested grids to reduce the amount of data while keeping the large scales of interest. Data is transferred from one grid level to another using high-order filters and accurate interpolation techniques. This method enabled to apply modal decomposition techniques to a billion-cell LES of a 3D turbulent turbine blade, thus demonstrating its effectiveness. The capability of performing calculations on several embedded grid levels was then used to devise the multi-resolution LES (MR-LES). The aim of the method is to evaluate the modeling and numerical errors during an LES by conducting the same simulation on two different mesh resolutions, simultaneously. This error estimation is highly valuable as it allows to generate optimal grids through the building of an objective grid quality measure. MR-LES intents to limit the computational cost of the simulation while minimizing the sub-grid scale modeling errors. This novel framework was applied successfully to the simulation of a turbulent flow around a 3D cylinder.

No description available.

No description available.

<p>Combustion instability, a complex phenomenon observed in combustion chambers is due to the coupling between heat release and other unsteady flow processes. Combustion instability has long been a topic of interest to rocket scientists and has been extensively investigated experimentally and computationally. However, to date, there is no computational tool that can accurately predict the combustion instabilities in full-size combustors because of the amount of computational power required to perform a high-fidelity simulation of a multi-element chamber. Hence, the focus is shifted to reduced fidelity computational tools which may accurately predict the instability by using the information available from the high-fidelity simulations or experiments of single or few-element combustors. One way of developing reduced fidelity computational tools involves using a reduced fidelity solver together with the flame transfer functions that carry important information about the flame behavior from a high-fidelity simulation or experiment to a reduced fidelity simulation.</p> <p> </p> <p>To date, research has been focused mainly on premixed flames and using acoustic solvers together with the global flame transfer functions that were obtained by integrating over a region. However, in the case of rockets, the flame is non-premixed and distributed in space and time. Further, the mixing of propellants is impacted by the level of flow fluctuations and can lead to non-uniform mean properties and hence, there is a need for reduced fidelity solver that can capture the gas dynamics, nonlinearities and steep-fronted waves accurately. Nonlinear Euler equations have all the required capabilities and are at the bottom of the list in terms of the computational cost among the solvers that can solve for mean flow and allow multi-dimensional modeling of combustion instabilities. Hence, in the current work, nonlinear Euler solver together with the spatially distributed local flame transfer functions that capture the coupling between flame, acoustics, and hydrodynamics is explored.</p> <p> </p> <p>In this thesis, the approach to extract flame transfer functions from high-fidelity simulations and their integration with nonlinear Euler solver is presented. The dynamic mode decomposition (DMD) was used to extract spatially distributed flame transfer function (FTF) from high fidelity simulation of a single element non-premixed flame. Once extracted, the FTF was integrated with nonlinear Euler equations as a fluctuating source term of the energy equation. The time-averaged species destruction rates from the high-fidelity simulation were used as the mean source terms of the species equations. Following a variable gain approach, the local species destruction rates were modified to account for local cell constituents and maintain correct mean conditions at every time step of the nonlinear Euler simulation. The proposed reduced fidelity model was verified using a Rijke tube test case and to further assess the capabilities of the proposed model it was applied to a single element model rocket combustor, the Continuously Variable Resonance Combustor (CVRC), that exhibited self-excited combustion instabilities that are on the order of 10% of the mean pressure. The results showed that the proposed model could reproduce the unsteady behavior of the CVRC predicted by the high-fidelity simulation reasonably well. The effects of control parameters such as the number of modes included in the FTF, the number of sampling points used in the Fourier transform of the unsteady heat release, and mesh size are also studied. The reduced fidelity model could reproduce the limit cycle amplitude within a few percent of the mean pressure. The successful constraints on the model include good spatial resolution and FTF with all modes up to at least one dominant frequency higher than the frequencies of interest. Furthermore, the reduced fidelity model reproduced consistent mode shapes and linear growth rates that reasonably matched the experimental observations, although the apparent ability to match growth rates needs to be better understood. However, the presence of significant heat release near a pressure node of a higher harmonic mode was found to be an issue. This issue was rectified by expanding the pressure node of the higher frequency mode. Analysis of two-dimensional effects and coupling between the local pressure and heat release fluctuations showed that it may be necessary to use two dimensional spatially distributed local FTFs for accurate prediction of combustion instabilities in high energy devices such as rocket combustors. Hybrid RANS/LES-FTF simulation of the CVRC revealed that it might be necessary to use Flame Describing Function (FDF) to capture the growth of pressure fluctuations to limit cycle when Navier-Stokes solver is used.</p> <p> </p> <p>The main objectives of this thesis are:</p> <p>1. Extraction of spatially distributed local flame transfer function from the high fidelity simulation using dynamic mode decomposition and its integration with nonlinear Euler solver</p> <p>2. Verification of the proposed approach and its application to the Continuously Variable Resonance Combustor (CVRC).</p> <p>3. Sensitivity analysis of the reduced fidelity model to control parameters such as the number of modes included in the FTF, the number of sampling points used in the Fourier transform of the unsteady heat release, and mesh size.</p> <p> </p> <p>The goal of this thesis is to contribute towards a reduced fidelity computational tool which can accurately predict the combustion instabilities in practical systems using flame transfer functions, by providing a path way for reduced fidelity multi-element simulation, and by defining the limitations associated with using flame transfer functions and nonlinear Euler equations for non-premixed flames.</p> <p> </p><br>

Le contrôle d'écoulements turbulents est un enjeu majeur en aérodynamique. Cependant, la présence d'un grand nombre de degrés de libertés et d'une dynamique complexe rend délicat la modélisation dynamique de ces écoulements qui est pourtant nécessaire à la conception d'un contrôle efficace. Au cours de cette thèse, différentes directions ont été suivies afin de développer des modèles réduits dans des configurations réalistes d'écoulements et d'utiliser ces modèles pour le contrôle.Premièrement, la décomposition en modes dynamiques (DMD), et certaines de ses variantes, ont été exploitées en tant que base réduite afin d'extraire au mieux le comportement dynamique de l'écoulement. Par la suite, nous nous sommes intéressés à l'assimilation de données 4D-Var qui permet de combiner des informations inhomogènes provenant d'un modèle dynamique, d'observations et de connaissances a priori du système. Nous avons ainsi élaboré des modèles réduits POD et DMD d'un écoulement turbulent autour d'un cylindre à partir de données expérimentales PIV. Finalement, nous avons considéré le contrôle d'écoulement dans un contexte d'interaction fluide/structure. Après avoir montré que les mouvements de solides immergés dans le fluide pouvaient être représentés comme une contrainte supplémentaire dans le modèle réduit, nous avons stabilisé un écoulement de sillage de cylindre par oscillation verticale. / Control of turbulent flows is still today a challenge in aerodynamics. Indeed, the presence of a high number of active degrees of freedom and of a complex dynamics leads to the need of strong modelling efforts for an efficient control design. During this PhD, various directions have been followed in order to develop reduced-order models of flows in realistic situations and to use it for control. First, dynamic mode decomposition (DMD), and some of its variants, have been exploited as reduced basis for extracting at best the dynamical behaviour of the flow. Thereafter, we were interested in 4D-variational data assimilation which combines inhomogeneous informations coming from a dynamical model, observations and an a priori knowledge of the system. POD and DMD reduced-order models of a turbulent cylinder wake flow have been successfully derived using data assimilation of PIV measurements. Finally, we considered flow control in a fluid-structure interaction context. After showing that the immersed body motion can be represented as an additional constraint in the reduced-order model, we stabilized a cylinder wake flow by vertical oscillations.

No description available.

Lors d’un vol hypersonique (Mach 6, 20 km d’altitude) la couche limite se développant sur l’avant-corps d’un véhicule hypersonique est laminaire. Cet état cause un désamorçage du moteur (statoréacteur) assurant la propulsion du véhicule. Pour pallier ce problème, il faut forcer la transition de la couche limite á l’aide d’un dispositif de contrôle dont l’effet est permanent (passif) ou modulable (actif) pendant le vol. Dans ce travail, nous analysons l’efficacité d’un dispositif actif d’injection d’air á la paroi pour forcer la transition de la couche limite sur un avant-corps générique. L’interaction jet d’air/couche limite est simulée numériquement avec une approche aux grandes échelles (LES). Une étude paramétrique sur la pression d’injection permet de quantifier l’efficacité du jet á déstabiliser la couche limite. L’influence des conditions de vol (altitude, Mach) sur la transition est également étudiée. Une analyse des résultats de simulation par Décomposition en Modes Dynamiques (DMD) est menée pour comprendre quels sont les modes dynamiques responsables de la transition et les mécanismes sous-jacents. Des essais dans la soufflerie silencieuse de l’université de Purdue (BAM6QT) ont été effectués pour tester expérimentalement l’efficacité des dispositifs passifs (rugosité isolée en forme de losange) et actifs (mono-injection d’air) pour faire transitionner la couche limite. Une peinture thermo-sensible et des capteurs de pression (PCB, Kulite) ont été utilisés pour déterminer la nature de la couche limite. Les résultats de ce travail montrent qu’une injection sonique suffit pour forcer la couche limite. On observe des essais, que pour une même hauteur de pénétration, les rugosités isolées sont moins efficaces que les jets (mono injection) pour déstabiliser la couche limite. / During a hypersonic flight (Mach 6, 20 km altitude), the boundary layer developing on the forebody of a vehicle is laminar. This state may destabilize the scramjet engine propelling the vehicle. To overcome this problem during the flight, the boundary layer transition has to be forced using a control device whose effect is fixed (passive) or adjustable (active). In this work, we analyze the efficiency of a jet in crossflow in forcing the boundary layer transition on a generic forebody. The flow is computed with a Large Eddy Simulations (LES) approach. A parametric study of the injection pressure allows the efficiency of the jet in tripping the boundary layer to be quantified. The influence of flight conditions (Mach, altitude) on the transition is also studied. Dynamic Mode Decomposition (DMD) is applied to the simulation results to determine the transition leading to dynamic modes and to understand underlying transition mechanisms. Experiments in the Purdue University quiet wind tunnel (BAM6QT) were performed to quantify the efficiency of a passive transition device (diamond roughnesses) and an active transition device (single air jet) in tripping the boundary layer. A thermo-sensitive paint and pressure transducers (Kulite, PCB) were used to determine the state of the boundary layer on the generic forebody. Experimental and numerical results show a sonic injection is sufficient to induce transition. We observe from the experiments that for the same penetration height, a single roughness is less efficient than a single air jet in destabilizing the boundary layer.

<p>Probing the subsurface for evidence related to the degradation of porous mediums and the evolution of damage mechanisms has been a long-standing challenge in geophysics. As such imaging and predicting fracture network development has remained a difficult area for subsurface science for decades despite the seminal and significant works put forward by many researchers. While this has provide great understanding about the behaviours and properties of natural porous media, there is still much that needs to be explored particularly in regard to the mineralogical composition and chemistry of clay-rich rocks. Despite the fact that argillaceous rocks which consist of different types of clays and varied mineral composition are ubiquitous in nature and are often the target of several technologies (e.g. geotechnical engineering, nuclear waste storage and disposal,hydrocarbon exploration and extraction, carbon capture and sequestration, etc.), many studies focus primarily on the bulk properties or the percentage of components in the matrix. For these reason and due to the problems that can be encountered with natural rocks that contain a swelling clay component whether randomly distirbuted or localized in consolidated globs in zones of the matrix, the influence of clay chemistry in relation to fracture development which is not well characterized, especially during desaturation is investigated with analogue rock samples which were systematically fabricated for this purpose.</p> <p><br></p> <p>The research performed in this dissertation investigated, the applicability of  the fabrication protocol for developing synthetic rocks with desirable rock like features and behavior, the impact and relationship between the rock properties, the microstructural composition, water loss, and the macroscopic behavior of the analogue rocks, focusing on the structure and chemistry of the constituent clay materials. Synthetic rocks were fashioned with the necessary geometries, properties, and material compositions. On the macroscopic scale the fracture and drying behavior of the synthetic rocks were examined with 3D X-ray microscopy and further evaluated through the utility of acoustic emission monitoring, water loss monitoring, and unconfined compressive testing. On the finer scale (nano-microscale), the chemical and mechanical properties, and behavior of select clays was explored by exploiting several methods of material characterization which also included cation exchange experiments coupled with inductively coupled plasma – optical emission spectrometry (ICP-OES). </p> <p><br></p> <p>For the finer scale, experiments verified that calcined kaolinite clay had a different mineral structure and negligible to non-existence shrinkage abilities. In contrast, the montmorillonite clays possessed higher and similar moisture contents but, owing to the different principal cations these clays interacted a bit differently in the highly akaline environment, experienced varying degrees of shrinkage, and had observedly minor structural dissimilarities. For the relatively larger scale, the emergence of damage, extent of the damage network, and the patterns of the crack network mainly depended on the microstructural composition of the analogue rocks, particularly it's clay chemistry and/ or distribution. The location of damage depended on the emplacement and percentage of swelling clay in the matrix, and numerical investigations with peridynamics revealed that the observed damage was a consequence of the action of the swelling and non-swelling components of the matrix. Furthermore, if the microstructure consisted of no clay or calcined kaolinite the AE activity was solely attributed to interfacial processes that occurred during fluid front movement. If the microstructure consisted of a particular montmorillonite, the cracks propagated in the direction of the drying front. Conversely, for montmorillonite clay predominated by a different principal cation, the crack network developed and propagated differently during water loss. Additionally, on the laboratory core scale, properties and behavior similar to natural rocks were confirmed and the rock strength, porosity, AE activity, and velocities were primarily affected by the microstructural composition of the analogue rocks. </p> <p><br></p> <p>An added challenge for investigating and monitoring evolving systems and processes, whether on the laboratory or field scale, is the problem of extracting useful information from the physical data that can be used to identify signatures of developing processes, and changes in the properties or the behavior of a system. Here, data driven machine learning modeling and clustering techniques were undertaken to build a mechanistic understanding of the AE activity generated during drying. The intent is for this work to add to the fundamental research aimed at developing methods that will robustly detect and extract signatures related to evolutionary processes or features in the AE signals, and group them according to some degree of similarity. Such research will support reliable interpretations of the physical data for predictions of the behavior of systems, development of engineering controls, and improvement of the understanding of intrinsic dynamics related to complex processes particularly those that occur in clay-rich systems.</p> <p><br></p> <p>Combined chemical and mechanical investigations have great potential for unraveling practical challenges in subsurface science, especially regarding damage processes in clay-rich rock systems, and identifying and interpreting the presence of discontinuities from geophysical data. The present findings are useful for establishing a link between the constituent clay and observed damage, and improving our understanding of the development of damage in clay-bearing systems. These results provide insight on the influence of swelling clay and the chemistry of such clays on the generation of cracks and crack networks in rock like materials which can be useful for the characterization of damage in both laboratory and the field. The work presented here can also be a basis for further experiments that aim to uncover methods and protocols that will help with the indirect characterization of evolutionary processes, damage mechanisms, and damage in clay rich porous media. Additionally support for the use of analogue rocks in experimental rock physics, architected with specific material compositions, pore structures, crack systems, or clay fractions, is provided here. </p>

[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