Lors d'un hypothétique accident grave de réacteur à eau sous pression, un mélange de matériaux fondus, appelé corium, issu de la fusion du cœur peut se relocaliser dans le puits de cuve constitué par un radier en béton. Les codes d'évaluation réacteur pour simuler la phénoménologie de l'interaction corium-béton sont basés sur une description à grande échelle des échanges qui soulève de nombreuses questions, tant sur la prise en compte des phénomènes multi-échelles mis en jeu que sur la structure adoptée de la couche limite au voisinage du front d'ablation. Dans ce contexte, l'objectif principal de ce travail consiste à aborder le problème de la structure de la couche limite par simulation numérique directe. Ce travail s'inscrit dans le cadre plus général d'une description et d'une modélisation multi-échelle des échanges, c'est-à-dire de l'échelle locale associée au voisinage du front d'ablation jusqu'à l'échelle du code d'évaluation réacteur. Une telle description multi-échelle des échanges soulève le problème de la description locale de l'écoulement multiphasique multiconstituant mais aussi le problème du changement d'échelle et en particulier le passage de l'échelle locale à l'échelle de description supérieure dite macroscopique associée aux mouvements convectifs dans le bain de corium. Parmi les difficultés associées au changement d'échelle, nous nous intéressons à la problématique de la construction de conditions aux limites effectives ou lois de parois pour les modèles macroscopiques. Devant la complexité du problème multiphasique multiconstituant posé au voisinage du front, cette contribution a été abordée sur un problème modèle. Des conditions aux limites dites effectives ont été construites dans le cadre d'une méthode de décomposition de domaine puis testées pour un problème d'écoulement laminaire de convection naturelle sur parois rugueuses. Mˆeme si le problème traité reste encore éloigné des applications visées, cette contribution offre de nombreuses perspectives et constitue une première étape d'une modélisation multiéchelle des échanges pour la problématique de l'interaction corium-béton. Dans le cas plus complexe des écoulements multiphasiques multiconstituants et devant les difficultés expérimentales associées, le développement de lois de parois pour les outils existants aux échelles de description supérieures nécessite, au préalable, de disposer d'un outil de simulation numérique directe de l'écoulement au voisinage du front d'ablation. L'outil développé dans ce travail correspond à un modèle de Cahn-Hilliard/Navier-Stokes pour un mélange diphasique (liquide-gaz) compositionnel (corium-béton fondu) s'appuyant sur une description du système selon trois paramètres d'ordre associés respectivement aux fractions volumiques du gaz et aux deux espèces miscibles de la phase liquide ainsi que sur une décomposition de l'énergie libre selon une contribution diphasique et compositionnelle. Les équations de transport sont dérivées dans le cadre de la thermodynamique des processus irréversibles et résolues sur la base d'une application éléments finis de la plate-forme PELICANS. Plusieurs expériences numériques illustrent la validité et les potentialités d'application de cet outil sur des problèmes diphasiques et/ou compositionnels. Enfin, à partir de l'outil développé, nous abordons par simulation numérique directe une étude de la structure de la couche limite au voisinage du front d'ablation pour des bétons siliceux et silico-calcaire. / In the late phases of some scenario of hypothetical severe accident in Pressurized Water Reactors, a molten mixture of core and vessel structures, called corium, comes to interact with the concrete basemat. The safety numerical tools are lumped parameter codes. They are based on a large averaged description of heat and mass transfers which raises some uncertainties about the multi-scale description of the exchanges but also about the adopted boundary layer structure in the vicinity of the ablation front. In this context, the aim of this work is to tackle the problem of the boundary layer structure by means of direct numerical simulation. This work joins within the more general framework of a multi-scale description and a multi-scale modeling, namely from the local scale associated with the vicinity of the ablation front to the scale associated with the lumped parameter codes. Such a multi-scale description raises not only the problem of the local description of the multiphase multicomponent flow but also the problem of the upscaling between the local- and the macro-scale which is associated with the convective structures within the pool of corium. Here, we are particularly interested in the building of effective boundary conditions or wall laws for macro-scale models. The difficulty of the multiphase multicomponent problem at the local scale leads us to consider a relatively simplified problem. Effective boundary conditions are built in the frame of a domain decomposition method and numerical experiments are performed for a natural convection problem in a stamp shaped cavity to assess the validity of the proposed wall laws. Even if the treated problem is still far from the target applications, this contribution can be viewed as a first step of a multi-scale modeling of the exchanges for the molten core concrete issue. In the more complicated case of multiphase multicomponent flows, it is necessary to have a direct numerical simulation tool of the flow at the local scale to build wall laws for macro-scale models. Here, the developed tool corresponds to a Cahn-Hilliard/Navier-Stokes model for a two-phase compositional system. It relies on a description of the system by three volume fractions and on a free energy composed by a two-phase part and a compositional part. The governing equations are derived in the frame of the thermodynamic of irreversible processes. They are solved on the basis of a finite element application of the object-oriented software component library PELICANS. Several numerical experiments illustrate the validity and the potentialities of application of this tool on two-phase compositional problems. Finally, using the developed tool, we tackle by means of direct numerical simulation the problem boundary layer structure in the vicinity of the ablation front for limestone-sand and siliceous concretes.

La thèse porte sur la simulation de la thermohydraulique et des transferts thermiques dans un faisceau d'aiguilles d'assemblage combustible de réacteur à neutrons rapides à caloporteur sodium.Des premiers calculs ont été réalisés par une approche moyennée de type RANS à l'aide du code industriel STAR-CCM+. De cette modélisation, il ressort une meilleure compréhension des transferts de chaleur opérés entre les aiguilles et le sodium. Les principales grandeurs macroscopiques de l'écoulement sont en accord avec les corrélations. Cependant, afin d'obtenir une description détaillée des fluctuations de température au niveau des fils espaceur, une approche plus détaillée de type LES et DNS est apparue indispensable. Pour la partie LES, le code TRIO_U a été utilisé. Concernant la partie DNS, un code de recherche a été utilisé. Ces approches requièrent des temps de calculs considérables qui ont nécessité des géométries représentatives mais simplifiées.L'approche DNS permet d'étudier l'écoulement à bas nombre de Prandtl, qui induit un comportement très différent du champ thermique relativement au champ hydraulique. Le calcul LES de l'assemblage montre que la présence du fil espaceur génère l'apparition de points chauds locaux (~20°C) en aval de celui-ci par rapport à l'écoulement sodium, au niveau de son contact avec l'aiguille. Les fluctuations de température au niveau des fils espaceur sont faibles (~1°C-2°C). En régime nominal, l'analyse spectrale montre l'absence de grande amplitude d'oscillations de température à basse fréquence (2-10 Hz); les conséquences sur la tenue mécanique des structures devront être analysées. / The thesis focuses on the numerical simulation of sodium flow in wire wrapped sub-assembly of Sodium-cooled Fast Reactor (SFR).First calculations were carried out by a time averaging approach called RANS (Reynolds- Averaged Navier-Stokes equations) using industrial code STAR-CCM+. This study gives a clear understanding of heat transfer between the fuel pin and sodium. The main variables of the macroscopic flow are in agreement with correlations used hitherto. However, to obtain a detailed description of temperature fluctuations around the spacer wire, more accurate approaches like LES (Large Eddy Simulation) and DNS (Direct Numerical Simulation) are clearly needed. For LES approach, the code TRIO_U was used and for the DNS approach, a research code was used. These approaches require a considerable long calculation time which leads to the need of representative but simplified geometry.The DNS approach enables us to study the thermal hydraulics of sodium that has very low Prandtl number inducing a very different behavior of thermal field in comparison to the hydraulic field. The LES approach is used to study the local region of sub-assembly. This study shows that spacer wire generates the local hot spots (~20°C) on the wake side of spacer wire with respect to the sodium flow at the region of contact with the fuel pin. Temperature fluctuations around the spacer wire are low (~1-2°C). Under nominal operation, the spectral analysis shows the absence of any dominant peak for temperature oscillations at low frequency (2-10Hz). The obtained spectra of temperature oscillations can be used as an input for further mechanical studies to determine its impact on the solid structures.

Turbulent phenomena are found in both natural (e.g. the Earth's oceans, the Sun's corona) and artificial (e.g. flows through pipes, the plasma in a tokamak device) settings; evidence suggests that turbulence is usually the normal behaviour in most cases. Turbulence has been studied extensively for more than a century, but a complete and consistent theoretical description of it has not yet been proposed. It is in this context that the motion of particles under the influence of turbulent fields is studied in this work, with direct numerical simulations. The thesis is structured in three main parts. The first part describes the tools that are used. Methods of integrating particle trajectories are presented, together with a discussion of the properties that these methods should have. The simulation of magnetohydrodynamic (MHD) turbulence is discussed, while also introducing fundamental concepts of fluid turbulence. Particle trajectory integration requires information that is not readily available from simulations of turbulent flows, so the interpolation methods needed to adapt the fluid simulation results are constructed as well. The second part is dedicated to the study of two MHD problems. Simulations of Kolmogorov flow in incompressible MHD are presented and discussed, and also simulations of the dynamo effect in compressible MHD. These two scenarios are chosen because large scale structures are formed spontaneously by the turbulent flow, and there is an interest in studying particle transport in the presence of structures. Studies of particle transport are discussed in the third part. The properties of the overall approach are first analyzed in detail, for stationary predefined fields. Focus is placed on the qualitative properties of the different methods presented. Charged article transport in frozen turbulent fields is then studied. Results concerning transport of particles in fully developed, time-evolving, turbulent fields are presented in the final chapter.<p><p><p>\ / Doctorat en Sciences / info:eu-repo/semantics/nonPublished

[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

Turbulent boundary layers approximating those found on the NASA Orion Multi-Purpose Crew Vehicle (MPCV) thermal protection system during atmospheric reentry from the International Space Station have been studied by direct numerical simulation, with the ultimate goal of reducing aerothermodynamic heating prediction uncertainty. Simulations were performed using a new, well-verified, openly available Fourier/B-spline pseudospectral code called Suzerain equipped with a ``slow growth'' spatiotemporal homogenization approximation recently developed by Topalian et al. A first study aimed to reduce turbulence-driven heating prediction uncertainty by providing high-quality data suitable for calibrating Reynolds-averaged Navier--Stokes turbulence models to address the atypical boundary layer characteristics found in such reentry problems. The two data sets generated were Ma[approximate symbol] 0.9 and 1.15 homogenized boundary layers possessing Re[subscript theta, approximate symbol] 382 and 531, respectively. Edge-to-wall temperature ratios, T[subscript e]/T[subscript w], were close to 4.15 and wall blowing velocities, v[subscript w, superscript plus symbol]= v[subscript w]/u[subscript tau], were about 8 x 10-3 . The favorable pressure gradients had Pohlhausen parameters between 25 and 42. Skin frictions coefficients around 6 x10-3 and Nusselt numbers under 22 were observed. Near-wall vorticity fluctuations show qualitatively different profiles than observed by Spalart (J. Fluid Mech. 187 (1988)) or Guarini et al. (J. Fluid Mech. 414 (2000)). Small or negative displacement effects are evident. Uncertainty estimates and Favre-averaged equation budgets are provided. A second study aimed to reduce transition-driven uncertainty by determining where on the thermal protection system surface the boundary layer could sustain turbulence. Local boundary layer conditions were extracted from a laminar flow solution over the MPCV which included the bow shock, aerothermochemistry, heat shield surface curvature, and ablation. That information, as a function of leeward distance from the stagnation point, was approximated by Re[subscript theta], Ma[subscript e], [mathematical equation], v[subscript w, superscript plus sign], and T[subscript e]/T[subscript w] along with perfect gas assumptions. Homogenized turbulent boundary layers were initialized at those local conditions and evolved until either stationarity, implying the conditions could sustain turbulence, or relaminarization, implying the conditions could not. Fully turbulent fields relaminarized subject to conditions 4.134 m and 3.199 m leeward of the stagnation point. However, different initial conditions produced long-lived fluctuations at leeward position 2.299 m. Locations more than 1.389 m leeward of the stagnation point are predicted to sustain turbulence in this scenario. / text

This thesis details the development, verification and validation of an unsteady unstructured high order (≥ 3) h/p Discontinuous Galerkin - Fourier solver for the incompressible Navier-Stokes equations on static and rotating meshes in two and three dimensions. This general purpose solver is used to provide insight into cross-flow (wind or tidal) turbine physical phenomena. Simulation of this type of turbine for renewable energy generation needs to account for the rotational motion of the blades with respect to the fixed environment. This rotational motion implies azimuthal changes in blade aero/hydro-dynamics that result in complex flow phenomena such as stalled flows, vortex shedding and blade-vortex interactions. Simulation of these flow features necessitates the use of a high order code exhibiting low numerical errors. This thesis presents the development of such a high order solver, which has been conceived and implemented from scratch by the author during his doctoral work. To account for the relative mesh motion, the incompressible Navier-Stokes equations are written in arbitrary Lagrangian-Eulerian form and a non-conformal Discontinuous Galerkin (DG) formulation (i.e. Symmetric Interior Penalty Galerkin) is used for spatial discretisation. The DG method, together with a novel sliding mesh technique, allows direct linking of rotating and static meshes through the numerical fluxes. This technique shows spectral accuracy and no degradation of temporal convergence rates if rotational motion is applied to a region of the mesh. In addition, analytical mappings are introduced to account for curved external boundaries representing circular shapes and NACA foils. To simulate 3D flows, the 2D DG solver is parallelised and extended using Fourier series. This extension allows for laminar and turbulent regimes to be simulated through Direct Numerical Simulation and Large Eddy Simulation (LES) type approaches. Two LES methodologies are proposed. Various 2D and 3D cases are presented for laminar and turbulent regimes. Among others, solutions for: Stokes flows, the Taylor vortex problem, flows around square and circular cylinders, flows around static and rotating NACA foils and flows through rotating cross-flow turbines, are presented.

[ES] El proyecto ESTiMatE nace con el objetivo de mejorar la precisión de las predicciones de hollín para la reducción de emisiones y una mayor sostenibilidad ambiental en la aviación. Para ello es crucial comprender a fondo fenómenos que influyen en la formación de hollín, incluyendo atomización del combustible, mezcla, combustión y formación de emisiones. La presente tesis, como parte del proyecto, se centra en el primer paso para caracterizar las spray flames, el proceso de atomización del líquido. Existen modelos empíricos para ello, pero el objetivo de este trabajo es desarrollar un nuevo modelo fenomenológico a través de simulaciones de alta fidelidad. Se propone acoplar este modelo a los códigos PRECISE-UNS y Alya para inyectar gotas en simulaciones Lagrangianas, reproduciendo de forma efectiva las características clave de atomización primaria en la formación del spray. Realizar estudios sobre la configuración anular de atomizadores prefilming airblast es un gran desafío, por lo que los esfuerzos se han enfocado en la configuración plana como alternativa que proporciona información valiosa. El KIT ha hecho contribuciones notables en este sentido al desarrollar un banco de pruebas de configuración plana, generando una base de datos extensa con resultados que abarcan diversos fluidos de trabajo y condiciones de operación para validar los cálculos. Esta tesis investiga la atomización primaria en estos dispositivos con un enfoque computacional a través de simulaciones eDNS con el código PARIS Simulator. Se utilizan simulaciones VOF-DNS del borde del atomizador para obtener información, y simulaciones LES monofásicas precursoras que permiten tener en cuenta la turbulencia del gas en la entrada. No obstante, se han propuesto dos metodologías derivadas a partir de eDNS, extendiéndola con simulaciones LES bifásicas adicionales. La primera de ellas impone un espesor y velocidad constantes para la película líquida en la entrada, pero sus valores medios se obtienen procesando las nuevas simulaciones LES bifásicas. Esta metodología ha sido validada para una condición de operación de referencia, representativa del reencendido en altitud. La segunda, introduce un espesor y velocidad variables de la película líquida, constituyendo una condición de contorno más realista. Se han realizado dos estudios diferentes. Por un lado, utilizando la primera metodología, se realizó un estudio paramétrico para investigar el impacto de la velocidad media del gas y las propiedades del fluido en la formación de estructuras líquidas. Al variar estos parámetros se observan diferencias en los mecanismos de ruptura basadas en M como en la literatura. La base para un modelo fenomenológico se ha establecido utilizando estos resultados y puede extenderse con más condiciones de operación siguiendo la metodología descrita. Por otra parte, se ha realizado una comparación entre ambas metodologías para considerar la influencia de la historia de la película de combustible en la atomización primaria. En este estudio, al mejorar las condiciones de contorno a la entrada, se ha reproducido con éxito la fenomenología presente en los experimentos, al contrario que con la eDNS estándar. Se ha propuesto una nueva técnica de postprocesado para caracterizar las estructuras líquidas. Los métodos habituales detectan gotas en un dominio 3D, pero solamente detectan ligamentos en una proyección 2D utilizando optimización de funciones algebraicas. Sin embargo, este nuevo método es capaz de encontrar gotas y ligamentos en 3D mediante OpenCV. También ha sido validado con resultados de la literatura. Con todo, la tesis ha sentado las bases para el desarrollo de modelos específicos de atomización airblast mediante simulaciones de alta fidelidad, contribuyendo a una mejor predicción de los fenómenos de inyección, atomización, evaporación, mezcla y combustión de aeromotores que permita realizar un diseño más eficiente de éstos hacia una movilidad aérea sostenible. / [CA] El projecte ESTiMatE naix amb l'objectiu de millorar la precisió de les prediccions de sutge per a la reducció d'emissions i una major sostenibilitat ambiental en l'aviació. Per això, és crucial comprendre a fons fenòmens que influeixen en la formació de sutge, incloent atomització del combustible, mescla, combustió i formació d'emissions. La present tesi, com a part del projecte, se centra en el primer pas per caracteritzar les spray flames, el procés d'atomització del líquid. Existeixen alguns models empírics per això, però l'objectiu d'aquest treball és desenvolupar un nou model fenomenològic a través de simulacions d'alta fildelitat. Es proposa acoblar aquest model als codis PRECISE-UNS i Alya per injectar gotes en simulacions Lagrangianes, reproduint de manera efectiva les característiques clau d'atomització primària en la formació de l'spray. Realitzar estudis sobre la configuració anul·lar d'atomitzadors prefilming airblast és un gran desafiament, per la qual cosa els investigadors han enfocat els seus esforços en la configuració plana com alternativa que proporciona informació valuosa. El KIT ha fet contribucions notables en aquest sentit en desenvolupar un banc de proves de configuració plana, generant una base de dades extensa amb resultats que abasten diversos fluids de treball i condicions d'operació per a validar els càlculs. Esta tesi investiga l'atomització primària en aquests dispositius amb un enfocament computacional a través de simulacions eDNS emprant el codi PARIS Simulator. S'utilitzen simulacions VOF-DNS de la vora de l'atomitzador per a obtindre informació, i simulacions LES monofàsiques precursores que permeten tindre en compte la turbulència del gas a l'entrada. No obstant això, s'han proposat dues noves metodologies derivades a partir de eDNS, estenent-la amb simulacions LES bifàsiques addicionals. La primera d'elles imposa una grossària i velocitat constants per a la pel·lícula líquida en l'entrada, però els seus valors mitjans s'obtenen processant les noves simulacions LES bifàsiques. Aquesta metodologia ha sigut validada per a una condició d'operació de referència, representativa de la reencesa en altitud. La segona, introduïx una grossària i velocitat variables de la pel·lícula líquida, constituint una condició de contorn més realista. S'han realitzat dos estudis diferents. D'una banda, utilitzant la primera metodologia, es va realitzar un estudi paramètric per a investigar l'impacte de la velocitat mitjana del gas i les propietats del fluid en la formació d'estructures líquides. En variar estos paràmetres s'observen diferènciesen els mecanismes de ruptura basades en M com a la literatura. La base per a un model fenomenològic s'ha establit utilitzant aquests resultats i pot estendre's amb més condicions d'operació seguint la metodologia descrita. D'altra banda, s'ha realitzat en paral·lel una comparació entre totes dues metodologies per a considerar la influència de la història de la pel·lícula de combustible en l'atomització primària. En aquest estudi, en millorar les condicionsde contorn a l'entrada, s'ha reproduït amb èxit la fenomenologia present en els experiments, al contrari que amb la eDNS estàndard. S'ha proposat una nova tècnica de postprocessament per a caracteritzar les estructures líquides. Els mètodes habituals detecten gotes en un domini 3D, però solament detecten lligaments en una projecció 2D utilitzant optimització de funcions algebraiques. No obstant això, aquest nou mètode és capaç de trobar gotes i lligaments en 3D mitjançant la llibrería OpenCV. També ha sigut validat amb resultats de la literatura. Amb tot, la tesi ha establit les bases per al desenvolupament de models específics d'atomització airblast mitjançant simulacions d'alta fidelitat, contribuïnt a una millor predicció dels fenòmens d'injecció, atomització, evaporació, mescla i combustió dels aeromotors que permeta realitzar un disseny més eficient dels mateixos cap a una mobilitat aèria sostenible. / [EN] The project ESTiMatE emerges aiming to enhance the precision of soot predictions leading to reduced emissions and improved environmental sustainability in aviation. To achieve this, it is crucial to thoroughly understand various phenomena that influence soot formation, including fuel atomization, mixing, combustion, and subsequent emissions formation. The present thesis, as part of the project, focuses in the first step to characterize spray flames, the liquid atomization process. There exist some empirical models to characterize atomization, but the objective of this specific work package is to develop a new phenomenological model through high-fidelity simulations instead of empirical results. This model is proposed to be coupled to the PRECISE-UNS and Alya codes to inject droplets in Lagrangian simulations, effectively reproducing the key features of the primary atomization process on the spray formation. Conducting fundamental studies on the annular configuration of prefilming airblast atomizers is challenging, so researchers have directed their attention towards planar configurations as an alternative that can provide valuable insights. The KIT has made notable contributions in this regard by developing a dedicated planar test rig, generating a substantial database of results encompassing various working fluids and operating conditions for validation. This thesis investigates primary atomization in these devices with a computational approach using eDNS within the PARIS Simulator code. VOF-DNS simulations of the atomizing Edge are performed to obtain useful data, and precursor one-phase LES simulations in order to account for gas inflow turbulence at the DNS inlet. Nevertheless, two methodologies have been derived from eDNS, extending it with additional precursor two-phase LES simulations. The first one, imposes constant liquid film thickness and velocity at the DNS inlet, but their mean values are obtained processing the new two-phase LES simulations. It has been validated for a reference operating condition, representative of altitude relight. The second one, introduces variable liquid film thickness and velocity instead, constituting a more realistic boundary condition for the DNS. Two different studies have been conducted in this regard. On the one hand, using the former methodology, a parametric study was performed to investigate the impact of mean gas velocity and fluid properties on the formation of liquid structures. By varying these parameters differences on the breakup mechanisms are observed based on M as found in the literature. The basis for a phenomenological model has been established using these results, and it can be extended with more operating conditions following the described methodology. On the other hand, a comparison between both methodologies has been carried out in parallel, in order to account for the influence of film history on the primary atomization process. In this study, by improving the inlet boundary conditions, the phenomenology observed in the experiments has been succesfully reproduced, unlike with the standard eDNS. A novel post-processing technique for characterizing liquid structures has been proposed. Usual post-processing methods detect droplets in a 3D domain, but only detect ligaments in a top view 2D projection of it using algebraic functions optimization. However, this new method is able to find both droplets and ligaments as well as its 3D properties using the Open Computer Vision Library (OpenCV) instead. It has been also validated with the results in the literature. Overall, the thesis has laid the foundations for the development of specific airblast atomization models through the use of high-fidelity simulations,contributing to a better prediction of the injection, atomization, evaporation, mixing and combustion phenomena in aero engines, enabling a more efficient design towards sustainable air mobility. / Moreno Montagud, C. (2024). Computational Modelling of the Primary Atomization in Aero Engine Airblast Atomizers through a DNS Approach [Tesis doctoral]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/211359

粉体の微粒子化に伴い，固液分離操作の中でも高精度かつ多量処理が可能なろ過の重要性はますます高まっている．しかしながら，微粒子スラリーのろ過抵抗は非常に大きく，所望の処理能力を得られないことが多い．本論文では，ろ過操作設計におけるろ材抵抗，スラリー凝集・分散状態，ケークろ過抵抗の予測に関して数値シミュレーションを利用した検討を行うことで，微粒子スラリーのろ過抵抗の低減に有用な知見を得ることができた． / 著者名の「吉」は「土」の下に「口」の置き換え / 博士(工学) / Doctor of Philosophy in Engineering / 同志社大学 / Doshisha University