No description available.

El cultivo protegido bajo invernadero tiene altos costes energéticos derivados de la ventilación forzada. En cambio, la ventilación natural puede ser una solución barata que reduzca el consumo energético. No obstante, un diseño que permita la ventilación natural del invernadero supone un reto debido a la complejidad del fenómeno físico. Frente a otros métodos de diseño, como los métodos experimentales o analíticos, en los últimos tiempos la dinámica de fluidos computacional (CFD por sus siglas en inglés) se ha convertido en la herramienta más utilizada para estudiar este tipo de fenómenos, gracias a su relativo bajo coste y a la rapidez en la obtención de resultados. No obstante, los modelos CFD deben refrendarse mediante validaciones realizadas a través de datos experimentales. Un análisis bibliográfico detallado del uso de la CFD aplicada a invernaderos muestra, en general, que los trabajos carecen de suficientes datos experimentales, seguramente debido al alto coste de los sensores para adquirirlos y la dificultad que conlleva el trabajo en condiciones de campo. Además, se observa que las simulaciones CFD en invernaderos no se realizan con un procedimiento sistematizado. Por ello, la presente tesis, por un lado describe un dispositivo y unos métodos sencillos y baratos para obtener datos atmosféricos, y por otro, propone una visión crítica sobre la investigación realizada hasta el momento, con el fin de sistematizar la manera de generar modelos CFD aplicados a la ventilación natural de invernaderos. Finalmente, la tesis se complementa con un ejemplo sobre un caso práctico. Para ello, en primer lugar, se realizó una revisión bibliográfica de las diferentes guías de buenas prácticas en diferentes campos de la tecnología principalmente en edificación, para sistematizar y adaptar las recomendaciones para generar modelos CFD en invernaderos. En segundo lugar, se desarrolló un sistema de adquisición de datos sencillo, consistente en una red de sensores, que permite medir simultáneamente la velocidad y dirección de viento en 20 puntos. Este sistema de sensores fue calibrado y probado en campo satisfactoriamente obteniendo una precisión similar a los anemómetros comerciales con un precio 30 veces superior. En tercer lugar, se generaron 24 modelos del flujo de aire alrededor del invernadero, resultantes de la combinación de cuatro modelos de turbulencia (k- ¿; RNG k-¿; k-¿ y RSM); dos esquemas de discretización (primer orden y segundo orden) y tres velocidades de viento exterior (3; 3,5 y 4 m/s) con el fin de analizar sus diferencias y demostrar sus ventajas e inconvenientes. Por este motivo, en cuarto lugar, se comprobó su capacidad de ajustarse a los datos de campo, validando los modelos con un análisis de regresión lineal sobre los datos experimentales. Con este estudió se reveló que los modelos SST k-¿ y el RSM (segundo orden) son los que mejor representan el flujo de ventilación y se demostró que el modelo k-¿ estándar (primer orden), el más utilizado en la bibliografía, no sólo ofrece resultados diferentes al resto de modelos, sino que su rendimiento es pobre para predecir el flujo de ventilación. / Granell Ruiz, R. (2014). Análisis del flujo ambiental y propuesta metodológica para simulaciones CFD aplicadas a la ventilación natural de invernaderos [Tesis doctoral]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/37194

In this dissertation, various VAT parameters are investigated to determine the effect of the overall efficiency of the turbine at a high Reynolds number. To increase the efficiency of the vertical axis turbines, 2D CFD simulations are completed in an effort to better understand the physics behind the operation of these turbines. Specifically, the effect of advance ratio, solidity, and wake interactions were investigated. Simulations were completed in OpenFOAM using the k-ω SST turbulence model at a nominal Reynolds number of 500,000 using a NACA 0015 airfoil. To simulate the motion of the turbine, Arbitrary Mesh Interfacing (AMI) was used. For all of the parameters tested, it was found that the geometric effective angle of attack seen by the turbine blades had a significant impact on the power extracted from the flow. The range of effective angles of attack was found to decrease as the advance ratio increased. In spite of this, a severe loss in the power coefficient occurred at an advance ratio of 2.5 during which the blade experienced dynamic stall. This effect was also seen when the number of turbine blades was changed to four, at a solidity of 1.08. This negative impact on performance was found to be due to the increase in the drag component of the tangential force when dynamic stall occurs. Results indicate that wake interactions between subsequent blades have a large impact on performance especially when the wake interaction alters the flow direction sufficiently to create conditions for dynamic stall. To improve the performance of the VAT in the presence of dynamic stall, calculations were completed of a static twisted blade profile using GenIDLEST and OpenFOAM. There was found to be no improvement in the lift coefficient when comparing the twisted blade profile with a 2D blade at the same median angle of attack as the twisted blade. To further see the effects of the twisted blade, an effective VAT pitching motion was given to the blade and again compared to a 2D blade with the same motion. In this case there was significant improvement seen in the performance of the twisted blade. / Master of Science

The urgent need to reduce green house gas emissions from the transport sector, particularly from heavy-duty trucks, has underscored the importance of developing more efficient internal combustion engines. Using computational fluid dynamics (CFD), this work investigated the impact of valve seat geometry on in-cylinder swirl, addressing a gap in research. Additionally, the suitability of steady-state simulations for providing valid qualitative data on port flow was assessed. To answer both research questions, two approaches were followed: steady-state port flow RANS simulations, and transient RANS simulations in a running engine setup. The results from the steady-state simulations highlighted the limitations of this approach to qualitatively predict swirl, as this quantity is highly dependent on the mesh. Despite these limitations, the steady-state simulations were still able to capture the trade-off between swirl and discharge coefficient, outlined in the literature. Transient simulations revealed that in-cylinder swirl is affected by the geometry of the valve seats. It was found that valve seats that direct the flow towards the liner, while avoiding strong flow separation tend to promote higher swirl, whereas valve seats that induce strong flow separation lead to lower swirl ratios. Despite the trade-off between swirl and volumetric efficiency, the volumetric efficiency losses were found to be practically negligible. The study emphasizes the need for a more comprehensive set of simulations, including more valve lifts and pressure ratios. Given the unsuitability of the steady-state simulations to predict swirl trends, future investigations should focus on replacing this approach by transient simulations with steady-state geometry and boundary conditions, properly addressing flow time-dependency at relatively low computational cost, and facilitating validation with experimental data.

Les objectifs de cette thèse sont d'étudier les capacité prédictive des méthodes statistiques URANS et hybrides RANS-LES à modéliser des écoulements complexes à haut nombre de Reynolds et de réaliser l'analyse physique de la turbulence et des structures cohérentes en proche paroi. Ces travaux traitent de configurations étudiées dans le cadre des projets européens ATAAC (Advanced Turbulent Simulation for Aerodynamics Application Challenges) et TFAST (Transition Location Effect on Shock Wave Boundary Layer Interaction). Premièrement, l'écoulement décollé autour d’une configuration de cylindre en tandem, positionnés l'un derrière l’autre, est étudiée à un nombre de Reynolds de 166000. Un cas statique, correspondant schématique aux support de train d'atterrissage, est d’abord considéré. L’interaction fluide-structure est ensuite étudiée dans le cas dynamique, dans lequel le cylindre aval possède un degré de liberté en translation dans la direction perpendiculaire à l'écoulement. Une étude paramétrique est menée afin d'identifier les différents régimes d'interaction en fonction des paramètres structuraux. Dans un deuxième temps, la physique du tremblement transsonique est étudiée au moyen d’une analyse temps-fréquence et d’une décomposition orthogonale en modes propres (POD), dans l’intervalle de nombre de Mach 0.70–0.75. Les interactions entre le choc principal, la couche limite décollée par intermittence et les tourbillons se développant dans le sillage, sont analysées. Un forçage stochastique, basée sur une réinjection de turbulence synthétique dans les équations de transport de l’énergie cinétique et du taux de dissipation générée à partir de la reconstruction POD, a été introduit dans l’approche OES (organised-eddy simulation). Cette méthode introduit une modélisation de la turbulence “upscale" agissant comme un mécanisme de blocage par tourbillons capable de prendre en compte les interfaces turbulent/non-turbulent et de couches de cisaillement autour des géométries. Cette méthode améliore grandement la prédiction des forces aérodynamiques et ouvre de nouvelles perspectives quant aux approches de type moyennes d’ensemble pour modéliser les processus cohérents et aléatoires à haut nombre de Reynolds. Enfin, l'interaction onde de choc/couche limite (SWBLI) est traitée, dans le cas d’un choc oblique à nombre de Mach 1.7, contribuant aux études de "design d'ailes laminaires" au niveau européen. Les performances des modèles URANS et hybrides RANS-LES ont été analysées en comparant, avec les résultats expérimentaux, les valeurs intégrales de la couche limite (épaisseurs de déplacement et de quantité de mouvement) et les valeurs à la paroi (coefficient de frottement). Les effets de la transition dans la couche limite sur l’interaction choc/couche limite sont caractérisés. / This thesis aims at analysing the predictive capabilities of statistical URANS and hybrid RANS-LES methods to model complex flows at high Reynolds numbers and carrying out a physical analysis of the near-region turbulence and coherent structures. This study handles configurations included in the European research programmes ATAAC (Advanced Turbulent Simulation for Aerodynamics Application Challenges) and TFAST (Transition Location Effect on Shock Wave Boundary Layer Interaction). First, the detached flow in a configuration of a tandem of cylinders, positionned behind one another, is investigated at Reynolds number 166000. A static case, corresponding to the layout of the support of a landing gear, is initially considered. The fluid-structure interaction is then studied in a dynamic case where the downstream cylinder, situated in the wake of the upstream one, is given one degree of freedom in translation in the crosswise direction. A parametric study of the structural parameters is carried out to identify the various regimes of interaction. Secondly, the physics of the transonic buffet is studied by means of time-frequency analysis and proper orthogonal decomposition (POD), in the Mach number range 0.70–0.75. The interactions between the main shock wave, the alternately detached boundary layer and the vortices developing in the wake are analysed. A stochastic forcing, based on reinjection of synthetic turbulence in the transport equations of kinetic energy and dissipation rate by using POD reconstruction, has been introduced in the so-called organised-eddy simulation (OES) approach. This method introduces an upscale turbulence modelling, acting as an eddy-blocking mechanism able to capture thin shear-layer and turbulent/non-turbulent interfaces around the body. This method highly improves the aerodynamic forces prediction and opens new ensemble-averaged approaches able to model the coherent and random processes at high Reynolds number. Finally, the shock-wave/boundary-layer interaction (SWBLI) is investigated in the case of an oblique shock wave at Mach number 1.7 in order to contribute to the so-called "laminar wing design" studies at European level. The performance of statistical URANS and hybrid RANS-LES models is analysed with comparison, with experimental results, of integral boundary-layer values (displacement and momentum thicknesses) and wall quantities (friction coefficient). The influence of a transitional boundary layer on the SWBLI is featured.

Ces travaux de recherche ont pour but d’évaluer la méthode dite de la "Simulation auxEchelles Adaptées" (SAS pour Scale-Adaptive Simulation). Cette approche coïncide avec uneapproche RANS classique dans les zones pariétales attachées et adapte le niveau de viscositéturbulente dans les zones décollées pour y permettre une résolution partielle des structures turbulentes.Dans une première partie, une analyse théorique du modèle SAS original a été menéeet a permis de développer une correction visant à favoriser l’adaptation du niveau de viscositéturbulente dans les zones sièges d’instabilités de type Kelvin-Helmholtz. Le modèle ainsi corrigéest nommé SAS-αL. Les modèles SAS et SAS-αL ont été implantés dans le code de calculNavier-Stokes elsA de l’ONERA. À l’issue de cette étape, trois cas académiques d’écoulementsturbulents instationnaires, cylindre à grand nombre de Reynolds, marche descendante et cavitétranssonique, ont été simulés grâce aux trois modèles de turbulence SST, SAS et SAS-αL. Outreune comparaison aux bases de données expérimentales disponibles, une attention particulièrea été portée à l’influence de paramètres numériques tels que des schémas numériques d’ordreélevé. Enfin, afin d’étudier la viabilité de l’approche SAS dans un contexte industriel, les troismodèles de turbulence ont été testés sur une configuration issue de l’industrie aéronautique etcorrespondant à la sortie d’air chaud d’un système de dégivrage des nacelles d’avion. La comparaisondes prévisions obtenues avec les modèles SST, SAS et SAS-αL aux données expérimentalesobtenues à l’ONERA a permis de montrer un gain de précision grâce à l’emploi de l’approcheSAS et ce pour un coût de calcul compatible avec un cycle de conception industrielle. / This research work is meant to assess an upgraded URANS approach, namely the Scale-Adaptive Simulation (SAS). This method is similar to a conventional RANS approach (namelythe SSTmodel) in attached areas and is able to adapt the eddy-viscosity level in detached areas toensure the resolution, at least partially, of the turbulent structures. In a first part of this researchwork, an improvement of the SAS approach is suggestedto allowa better sensitivity of themodelto instabilities such as Kelvin-Helmholtz ones. This "improved" model is referred to as SAS-αLmodel. Both SAS and SAS-αL models were implemented in the ONERA Navier-Stokes solverelsA and both of themaswell as the SSTmodelwere tested on academic test cases : a cylinder in acrossflowat a high Reynolds number, a backward-facing step flowcorresponding to theDriver&Seegmiller experiment and the transonic flow over the M219 cavity experimentally investigatedby de Henshaw. The influence of the numerical parameters was deeply investigated and particularattention was paid to the high-order space-discretization schemes effects. The reliabilityof the SAS approach in an industrial framework was assessed on an aeronautic configurationnamely a nacelle de-icing device. Comparisons between the threemodels (SST, SAS and SAS-αL)and an experimental database available at ONERA - The French Aerospace Lab have shown thebetter accuracy of the SAS approach as well as the high potential of the SAS-αL model.

Cette thèse présente la construction de schémas distribuant le résidu (RD) d'ordre très élevés, pour la discrétisation d'équations d'advection-diffusion multidimensionnelles et stationnaires sur maillages non structurés. Des schémas linéaires ainsi que des schémas non linéaires sont considérés. Une approximation de la solution polynomiale par morceaux et continue sur chaque élément est adoptée, de plus une procédure de reconstruction du gradient que celle de la solution numérique est utilisée afin d'avoir une représentation continue de la solution numérique et de son gradient. Il est montré que le gradient doit être reconstruit avec la même précision de la solution, sans quoi la précision formel du schéma numérique est perdue dans les cas où les effets de diffusion prévalent sur les effets d'advection, et aussi quand l'advection et la diffusion sont également importants. Ensuite, la méthode est étendue à des systèmes d'équations, en particulier aux équations de Navier-Stokes et aux équations RANS. La précision, l'efficacité et la robustesse du solveur RD implicite sont démontrées sur plusieurs cas tests. / The construction of compact high-order Residual Distribution schemes for the discretizationof steady multidimensional advection-diffusion problems on unstructuredgrids is presented. Linear and non-linear scheme are considered. A piecewise continuouspolynomial approximation of the solution is adopted and a gradient reconstructionprocedure is used in order to have a continuous representation of both thenumerical solution and its gradient. It is shown that the gradient must be reconstructedwith the same accuracy of the solution, otherwise the formal accuracy ofthe numerical scheme is lost in applications in which diffusive effects prevail overthe advective ones, and when advection and diffusion are equally important. Thenthe method is extended to systems of equations, with particular emphasis on theNavier-Stokes and RANS equations. The accuracy, efficiency, and robustness of theimplicit RD solver is demonstrated using a variety of challenging aerodynamic testproblems.

[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

A three dimensional computational model of a full scale axial flow tidal turbine has been used to investigate the effects of a range of realistic environmental conditions on turbine performance. The model, which is based on the Reynolds averaged Navier-Stokes equations, has been developed using the commercial flow solver ANSYS Fluent. A 1:30 scale tidal turbine is simulated in an open channel for comparison to existing experimental data. The rotor blades are directly resolved using a body-fitted, unstructured computational grid. Rotor motion is enabled through a sliding mesh interface between the rotor and the channel boundaries. Reasonably good agreement in thrust and power is observed. The computed performance curves are offset from the measured performance curves by a small increment in rotor speed. Subsequently, a full scale axial flow turbine is modelled in a variety of conditions representative of tidal channel flows. A parametric study is carried out to investigate the effects of flow shear, confinement and alignment on turbine performance, structural loading, and wake recovery. Mean power and thrust are found to be higher in sheared flow, relative to uniform flow of equivalent volumetric flow rate. Large fluctuations in blade thrust and torque occur in sheared flow as the blade passes through the high velocity freestream flow in the upper portion of the profile and the lower velocity flow near the channel bed. A stronger shear layer is formed around the upper portion of the wake in sheared flow, leading to enhanced wake mixing. Mean power and thrust are reduced when the turbine is simulated at a lower position in a sheared velocity profile. However, fluctuations in blade loading are increased due to the higher velocity gradient. The opposite effects are observed when the turbine operates at greater heights in sheared flow. Flow misalignment has a negative impact on mean rotor thrust and power, as well as on unsteady blade loading. Although the range of unsteady loading is not increased significantly, additional perturbations are introduced due to interactions between the blade and the nacelle. A deforming surface is introduced using the volume-of-fluid method. Linear wave theory is combined with the existing free surface model to develop an unsteady inflow boundary condition prescribing combined sheared flow and free surface waves. The relative effects of the sheared profile and wave-induced velocities on turbine loading are identified through frequency analysis. Rotor and blade load fluctuations are found to increase with wave height and wave length. In a separate study, the performance of bi-directional ducted tidal turbines is investigated through a parametric study of a range of duct profiles. A two dimensional axi-symmetric computational model is developed to compare the ducted geometries with an unducted device under consistent blockage conditions. The best-performing ducted device achieves a peak power coefficient of approximately 45% of that of the unducted device. Comparisons of streamtube area, velocity and pressure for the flow through the ducted device shows that the duct limits the pressure drop across the rotor and the mass flow through the rotor, resulting in lower device power.

Le prototypage virtuel des circuits pneumatiques de puissance, par exemple les circuits de freinage des véhicules industriels, constitue un enjeu important en raison de la complexité des écoulements en régime transsonique et des couplages entre les échelles locales et macroscopiques. Ces problèmes sont rencontrés lors de la conception, de la synthèse des commandes et de l'analyse des performances statiques et dynamiques de ces circuits et l'analyse. La mise au point des modèles numériques de ces systèmes induit des coûts et des temps importants par rapport à d'autres systèmes. La démarche proposée dans cette thèse repose sur la construction numérique de bases de données permettant de caractériser le comportement local et macroscopique d'un composant de circuit en fonction de la variation de certains paramètres physiques ou géométriques par rapport à un point de fonctionnement de référence. Les bases de données résultent de l'extrapolation de la solution des équations de Navier Stokes moyennées (RANS) pour le point de référence considéré obtenu à l'aide d’un logiciel de paramétrisation en mécanique des fluides (Turb’Opty). La contribution de cette thèse repose pour l'essentiel dans un travail d'analyse des solutions issues de la paramétrisation dans deux contextes différents: la tuyère De Laval et un élément "coude", des composants élémentaires de circuit. Nous avons montré que ces exemples "simples" conduisent déjà à des difficultés importantes en termes de paramétrisation du problème et du calcul des dérivées des champs aérodynamiques en raison de la taille du problème. Pour pallier cette difficulté, nous avons proposé de déraffiner le maillage et nous avons alors montré que cette démarche conduit parfois à déplacer ou à atténuer certains phénomènes (chocs). La deuxième contribution de ce travail repose sur l'évaluation de la qualité des solutions extrapolées, de leur domaine de validité et la construction des liens entre grandeurs locales et macroscopiques. Nous avons enfin proposé une démarche permettant de reconstruire la caractéristique en débit d'un composant à partir de la détermination de la solution extrapolée pour un nombre limité de points de référence. / Virtual prototyping pneumatic circuits for power transmission, for example braking circuits of trucks, is still a difficulty because of the complexity of the flow behavior in transonic conditions and of the coupling between local and macroscopic scales. These problems are met during system design, control synthesis and for static and dynamic performance analysis. Tuning accurate numerical models requires important costs and time when compared to other systems. The methodology proposed in this PhD thesis relies on numerically determining a data base that characterizes the local and macroscopic behavior of a circuit component according the variation from a reference point of some physical or geometrical parameters. The data bases are obtained from the extrapolation of the Mean Navier Stokes solution (RANS) for a given reference point with the help of a parametrization software dedicated to fluid mechanics (Turb'Flow). The main contribution of this thesis relies io the analysis of the solution obtained from the parametrization in two different cases: the De Laval nozzle and un "elbow" connecting element, which are elementary component in a circuit. We have shown that these two "simple" cases lead already to important difficulties in term of problem parametrization and calculation of the derivatives of the aerodynamic fields because of the problem dimension. In order to tackle this, we proposed to reduce the spatial discretization (mesh derefining) and we showed that this approach could sometimes lead to damp or move some phenomena (shocks). The second contribution of this work relies on evaluating the quality of the extrapolated solution and their validity domain, and on building links between local and macroscopic behavior. Finally, we proposed a method that allows the mass flow rate characteristic of a component to be determined from the calculation of the extrapolated solution issued from a limited number of reference points.