[pt] Os desafios das economias globais fazem com que as cadeias de suprimento (CS) tenham que aumentar seus processos de colaboração e dependência entre seus nós, gerando um aumento no nível de vulnerabilidade a possíveis impactos e interrupções em suas operações, que por consequência, podem afetar a aplicação da sustentabilidade em seus processos e operações. Com o rápido e extenso desenvolvimento do tema, surgiram pesquisas sobre o crescimento contínuo da complexidade da CS, aumentando a relevância dos estudos sobre suas vulnerabilidades a perturbações. Entretanto, ainda existem lacunas em como definir tais perturbações e como as cadeias afetadas vêm se recuperando delas. Portanto, objetiva-se analisar como os constructos de turbulência, incerteza e resiliência ambiental em cadeias de suprimentos sustentáveis têm evoluído ao longo dos anos. Para chegar a esse objetivo, a presente dissertação foi dividida em duas partes: a primeira teórica e a segunda empírica. Na parte teórica, será realizada uma revisão sistemática da literatura (RSL), tendo como resultado as definições sobre cada constructo e análises bibliométricas da evolução temática do tema, e uma Scoping Review com o objetivo de realizar um comparativo entre as principais escalas já desenvolvidas na área de Supply Chain Management. Na etapa empírica, foi desenvolvida e validada uma escala, através da técnica Fuzzy Delphi, que analise como tais fenômenos definidos anteriormente podem ser medidos, validados e replicados. Esta pesquisa foca em contribuir para a área de Sustainable Supply Chain Management, trazendo uma maior clareza de como os fenômenos de turbulência, incerteza e resiliência podem afetar a trajetória de sustentabilidade em cadeias de suprimento, sendo um tema de alta relevância nos dias atuais. / [en] The challenges of global economies mean that supply chains (CS) have to increase their processes of collaboration and dependence between their nodes, generating an increase in the level of vulnerability to possible impacts and interruptions in their operations, which consequently can affect the application of sustainability in its processes and operations. With the rapid and extensive development of the topic, research has emerged on the continuous growth of SC complexity, increasing the relevance of studies on its vulnerabilities to disturbances. However, there are still gaps in how to define such disturbances and how the affected chains are recovering from them. Therefore, the objective is to analyze how the constructs of turbulence, uncertainty and environmental resilience in sustainable supply chains have evolved over the years. To reach this goal, this dissertation was divided into two parts: the first theoretical and the second empirical. In the theoretical part, a systematic literature review (RSL) was carried out, resulting in the definitions of each construct and bibliometric analyzes of the thematic evolution of the theme, and a Scoping Review with the objective of making a comparison between the main scales already developed in the area of Supply Chain Management. In the empirical stage, a scale was developed and validated, through the Fuzzy Delphi technique, which analyzes how such previously defined phenomena can be measured, validated and replicated. This research focuses on contributing to the area of Sustainable Supply Chain Management, bringing greater clarity to how the phenomena of turbulence, uncertainty and resilience can affect the trajectory of sustainability in supply chains, being a topic of high relevance today.

[pt] Doenças cardiovasculares são responsáveis por um elevado número de óbitos em seres humanos. Muitas dessas patologias são dependentes do ciclo cardíaco e estão localizadas na aorta, a maior e principal artéria do nosso corpo. O conhecimento dos padrões de escoamento e distribuições de tensões nas paredes da aorta podem auxiliar no diagnóstico e prevenção de algumas dessas doenças. Dessa forma, estudou-se numericamente o escoamento do sangue, durante o ciclo cardíaco, em um modelo 3D da aorta de um paciente específico, após a implantação de TAVI (Transcatheter Aortic Valve Implantation). O ciclo cardíaco é formado por dois períodos chamados de sístole e diástole. Durante a sístole, sangue é bombeado do coração para a aorta, apresentado altos valores de vazão, resultando em escoamento turbulento. Por outro lado, na diástole, com o fechamento da válvula aórtica, o sangue escoa com baixas velocidades em regime laminar. Até hoje, cientistas enfrentam um desafio na modelagem da turbulência, pois não existe uma única modelagem que forneça previsibilidade para todas as situações envolvendo o regime turbulento, com esforço computacional razoável. Para seleção do modelo de turbulência mais adequado para análise do escoamento no interior da aorta, na presença da transição de regimes de escoamento durante o ciclo cardíaco, com um custo razoável, selecionou-se a metodologia baseada na Média de Reynolds. Diferentes modelos foram comparados com dados experimentais extraídos do mesmo modelo aórtico em escala real, porém em regime permanente, com vazão correspondente ao pico da sístole. Por fim, avaliou-se o impacto das condições de contorno e dos modelos de turbulência durante o ciclo cardíaco na distribuição e valores de tensões e grandezas turbulentas no endotélio vascular. Mostrou-se que a distribuição espacial das médias temporais de tensão foram qualitativamente e quantitativamente similares, para os dois ciclos cardíacos representativos de diferentes pacientes, porém com pequenas mudanças locais para cada caso. Em termos dos modelos de turbulência, observou-se que o modelo SAS (Scale Adaptive Simulation) foi capaz de representar a relaminarização do escoamento sanguíneo no período diastólico. / [en] Cardiovascular diseases are responsible for a high number of deaths in humans. Many of these pathologies are dependent on the cardiac cycle and are located in the aorta, the largest and main artery in our body. Knowledge of flow patterns and stress distributions in the walls of the aorta can help in the diagnosis and prevention of some of these diseases. Thus, the flow of blood during the cardiac cycle was numerically studied in a 3D model of the aorta of a specific patient, after TAVI (Transcatheter Aortic Valve Implantation) implantation. The cardiac cycle consists of two periods called systole and diastole. During the systole, blood is pumped from the heart to the aorta, presenting high flow rates, resulting in a turbulent flow. On the other hand, in diastole, with the closure of the aortic valve, the blood flows with low velocities in laminar regime. Until today, scientists face a challenge in turbulence modeling, as there is no single model that provides predictability for all situations involving the turbulent regime, with reasonable computational effort. In order to select the most suitable turbulence model for the analysis of the flow inside the aorta, in the presence of the transition of flow regimes during the cardiac cycle, with a reasonable cost, the methodology based on the Reynolds Average was selected. Different models were compared with experimental data extracted from the same real-scale aortic model, but a in steady state, with flow corresponding to the systolic peak. Finally, the impact of boundary conditions and turbulence models during the cardiac cycle on the distribution and values of stresses and turbulent quantities in the vascular endothelium were evaluated. It was shown that the spatial distribution of the temporal averages of tension was qualitatively and quantitatively similar, for the two cardiac cycles representative of different patients, but with small local changes for each case. In terms of turbulence models, it was observed that the SAS (Scale Adaptive Simulation) model was able to represent the relaminarization of blood flow in the diastolic period.

We study the formation and spreading of a turbulent jet of bubbles in microgravity. This has been analyzed from the recordings obtained in previous experimental campaigns of microgravity. Results have been compared with a simplified model of passive bubbles, in which bubbles are advected by the mean flow and dispersed due to the local degree of turbulence at each point of the jet. Thanks to the expertise obtained with this part of the thesis, we have designed and built a new experiment that has been used 36 times in the 4.7 s drop tower of ZARM (“Center of Applied Space Technology and Microgravity”) in Bremen. With this experiment we have obtained, for the first time in history, a monodisperse suspension of bubbles, within a turbulent flow, in microgravity. From the resulting measures we have characterized the relaxation time of pseudo-turbulence (previously generated due to the effect of buoyancy forces upon the injected bubbles in normal gravity conditions). We have also studied the interaction between bubbles and the turbulent medium. Results have been compared with Lattice-Boltzmann simulations of the flow. On the other hand, we have also studied the impact of residual gravitational vibrations (known as g-jitters) upon the quality of semiconductors solidified in microgravity. The quality of the resulting crystals has been studied from the analysis of the inhomogeneities in their dopant concentration. This study has been based entirely on simulations, but g-jitters have been modeled from acceleration signals measured in real space missions. / En la present tesi s’estudia, en primer lloc, la formació d’un doll turbulent de bombolles en condicions de microgravetat. Aquest ha sigut analitzat a partir del tractament de les gravacions obtingudes per altres investigadors en experiments de microgravetat. Els resultats s’han comparat amb un model simplificat de bombolles passives, en el que aquestes són arrossegades pel flux mitjà i, simultàniament, són dispersades degut al grau local de turbulència a cada punt. Gràcies a la experiència obtinguda en aquest anàlisi, s’ha dissenyat un nou experiment que ha sigut utilitzat en 36 llançaments de la torre de caiguda de 4.7 segons del ZARM (“Centre de Tecnologia Espacial Aplicada i Microgravetat”) a Bremen. Amb aquest experiment s’ha aconseguit, per primera vegada a la història, una suspensió monodispersa de bombolles, en el sí d’un flux turbulent, en condicions de microgravetat. A partir dels resultats obtinguts, s’ha caracteritzat per primera vegada el temps de relaxació de la pseudo-turulència (generada prèviament degut a l’efecte de les forces de flotació sobre les bombolles injectades en gravetat normal). També s’ha estudiat l’efecte causat per les bombolles en el medi turbulent. Els resultats han sigut comparats amb simulacions realitzades mitjançant el model de Lattice-Boltzmann. Per altra banda, s’ha estudiat també l’efecte que tenen les vibracions gravitatòries residuals sobre la qualitat de semiconductors solidificats en microgravetat. S’ha analitzat la qualitat dels cristalls resultants a partir de l’estudi de les inhomogeneïtats en la concentració de dopant. Aquest estudi ha sigut realitzat íntegrament a base de simulacions, però s’han establert els paràmetres dominants del soroll gravitatori a partir de valors mesurats en missions espacials reals.

Coastal zones are governed by physical forces originating from tidal currents, waves, winds and night convection, amongst others, and are characterized by the presence of canopy meadows. This thesis studies the hydrodynamic in a fluid dominated by: nearly isotropic turbulence, progressive waves and breaking waves in different canopy models. Under nearly isotropic turbulence sheltering is enhanced by a reduction in the plant-to-plant distance. Under progressive waves sheltering is associated with the reduction of wave velocity at the top of submerged rigid canopy. Sheltering observed in the submerged flexible model is caused by blade movement which absorbes the energy. Emergent rigid vegetation shows sheltering. On the other hand, for some specific progressive wave conditions and plant densities and under a fluid dominated by breaking waves, turbulence increases within the meadow if Reynolds number, based on wave velocity, is larger than 300 / Les zones costeres estan governades per forces físiques originades per la marea, l’onatge, i la convecció nocturna, entre d’altres, i caracteritzades per la presència de plantes aquàtiques. Aquesta tesi té per objectiu estudiar la hidrodinàmica d’un fluid dominat per: turbulència quasi-isotròpica, onades progressives i onades trencants en diferents models de vegetació. En un fluid dominat per turbulència isotròpica, la reducció de la turbulència augmenta al reduir la distància entre plantes. En un fluid dominat per onades progressives, la disminució de la turbulència està associada a la reducció de la velocitat d’onada al cim de la vegetació. Aquesta reducció s’observa en plantes submergides flexibles degut a que el moviment de les plantes absorbeix l’energia. La vegetació emergent també redueix la turbulència. Per condicions específiques d’onades progressives i densitat de plantes i sota un fluid dominat per onades trencants, la turbulència incrementa si el número de Reynolds és més gran que 300

A T E N C I Ó ! ! !Per consultar aquesta tesi cal navegar partint del fitxer 0INICI.pdfRESUMENLas variables físicas que describen la mezcla, son perfectamente comparables para experimentos de laboratorio estudiados por un nutrido grupo de investigadores en la dinámica de fluidos medioambientales. En este trabajo se relacionan variables de tipo físicas y biológicas y se destaca el uso de las primeras como ayuda en la interpretación de las variables biológicas, demostrándose su utilidad. Se utilizan además distintos experimentos de laboratorio enfocados al estudio de la dinámica de la mezcla turbulenta en un flujo estratificado inducida por cizalla y a las escalas de las estructuras inducidas por dinámica de la pluma en el delta. Este objetivo se lleva a cabo describiendo las características de un ambientes estuarino y de la pluma por medio de indicadores como el número de Richardson y de Reynolds. Para ello se emplea una metodología avanzada de toma de muestras, técnicas experimentales de laboratorio como procesamiento de imágenes y técnicas numéricas para la modelización de la pluma en la desembocadura del estuario.El estuario bajo estudio se encuentra presente en el delta del Ebro, donde existe la mayor parte de año una cuña estable que es controlada principalmente por los caudales locales del río. Con la finalidad de observar una cuña estable, se realizo una campaña en el mes de julio de 1997. Adicionalmente para la obtención de variable a partir de las cuales puede hacerse una descripción de la dinámica de la mezcla y la obtención de variables biológicas, se desarrollaron una serie de campañas donde se buscaron diferentes condiciones del flujo del río que determina la mezcla en el estuario; estas campañas se realizaron las fechas del 5 de abril de 1998, 12 de julio, 5 de octubre de 1999 y 5 de febrero de 2000. En este trabajo junto al trabajo de campo, se realizan experimentos de laboratorio y simulaciones numéricas; se establece además un criterio cuantitativo del inicio y final de una capa de una estratificación en la columna de agua, se establecen zonas características de la mezcla en el estuario, la distancia de la cabeza de la pluma en cualquier punto del estuario incluyendo la boca a partir de la medida de la densidad mínima en la columna de agua, que la medida de los caudales locales maneja dinámica de la mezcla tanto en el cuerpo del estuario como en la pluma, además de medidas de eficiencia de mezcla, relaciones de la abundancia de organismos típicos de agua marina en el estuario de acuerdo a uno de los descriptores de la mezcla como el número de Reynolds, se describe la dinámica general de la circulación superficial inducida por la pluma en la zona costera adyacente del delta del Ebro por medio del seguimiento de partículas trazadoras y de la circulación en la columna de agua por seguimiento de colorante, por último se comparan las estructuras y escalas dominantes observadas experimentalmente con estructuras observadas por medio de percepción remota y simulaciones numéricas obteniéndose una buena correlación con los tendencias las medidas naturales y experimentales de laboratorio y numéricas. / Physical variables that describe mixing are compared with experimental laboratory results from a large number of studies of mixing dynamics in environmental fluids to find habitat regimes for primary production. The work relates physical and biological variables, and emphasises the utility of laboratory studies. Several laboratory experiments that focused on turbulent mixing dynamics in stratified shear flows are used to describe (a) mixing in the estuary and (b) induced circulation in the river plume. Mixing descriptors as entrainment, Richardson number and Reynolds number and field data were employed, and advanced techniques of laboratory simulations, image processing and numerical modelling were used to match (a) to (b). Four kinds of experiments were used to describe the dynamics in the whole estuary. 1) Mixing turbulence across a density interface generated by an oscillating grid inside a mixing-box. 2) The horizontal advance of a turbulent front in a stratified system with a lateral current, inside a 1 m x 1 m square box. 3) Induced circulation in the delta del Ebro slope and shelf that were performed with an experimental model in a 2 m x 4 m rectangular tank on a 5-m diameter turntable. 4) Dispersion simulations in the river plume with the OCK3D code. Experiments 1 and 2 were developed in the UPC applied physic laboratories in Barcelona, experiments 3 were performed in the SINTEF laboratories in Trondheim and numerical experiments 4 were realized in the LSEET laboratories in Toulon. The estuary under study is the delta of the del Ebro, where a stable saline wedge is present the greater part of the year, driven by the river flow. A field campaign was carried out to observe a complete large stable saline wedge during the month of July of 1997. Additionally, to obtain the physical descriptors, cruises were made on 5 April, 1998, 12 July and 5 October, 1999 and 5 February, 2000 to observe several flow conditions in the river estuary. Particle tracking and optical measurements in the water column were used to study the induced length scale dynamics in the shelf and slope of the delta del Ebro. The dominant structures and scales from SAR images and numerical and laboratory simulations were compared with the field data.With the help of laboratory and field experiments and numerical simulations, qualitative criteria were suggested to determine the extent of the interface ending in a stratified water column, the different mixing characteristic zones within the estuary, and the distance of the estuary head, since elsewhere points in the estuary (including the river mouth) are obtained by means of the minimum density in the water column. Both estuary mixing and plume dynamics are driven by the local flow. In the present analysis of mixing efficiency, the phytoplankton abundance is in accord with the Reynolds number as a mixing descriptor.

[ES] El principal objetivo de esta tesis es el estudio de flujos térmicos turbulentos en canales para obtener un mayor conocimiento sobre el fenómeno de la turbulencia. Para ello, se ha realizado un estudio desde el punto de vista de la mecánica de fluidos computacional, en concreto, se ha utilizado la técnica de las simulaciones numéricas directas (DNS de sus siglas en inglés). La idea principal de las simulaciones realizadas ha sido ampliar el estado del arte actual, en lo referente a los dos parámetros principales que caracterizan el flujo: el número de Reynolds de fricción, Reτ, y el número de Prandtl, Pr. Dos configuraciones del flujo han sido utilizadas: flujo de Poiseuille y flujo de Couette, siendo la primera el principal foco del estudio. En cuanto al campo de temperaturas, se ha utilizado una condición de contorno mixta y se ha considerado como un escalar pasivo. Así pues, los números de Reynolds de fricción simulados para un flujo de Poisuille han sido Reτ = 500, 1000 y 2000, para números de Prandtl que varían desde 0.007 (metales fundidos) hasta 10 (agua), pasando por 0.71 que es el valor más utilizado por ser éste el número de Prandtl del aire. Además, se ha realizado una simulación con Reτ = 5000 y Pr = 0.71, la cual es la DNS térmica con el número de Reynolds de fricción más alto hasta la fecha. Destacar que para los números de Prandtl más altos, se ha observado que el valor máximo de la varianza de la temperatura es constante. Esto tiene un importante beneficio en el escalado cerca de la pared de los términos de disipación y difusión viscosa del balance de energía de θ′'+. Por último en lo referente a simulaciones de flujos de Poiseuille, se ha estudiado el caso isotérmico con Reτ = 10000, la cual es la mayor DNS de un canal turbulento, obteniendo por primera vez en una DNS una capa logarítmica perfectamente desarrollada en el campo de velocidades. Un estudio teórico, basado en las simetrías de Lie, ha sido llevado a cabo en paralelo a las simulaciones. El principal objetivo ha sido la generación de leyes de escaldo, basadas en primeros principios, del campo de velocidades, temperatura y momentos de altos órdenes de ambos campos. El resultado es que para números de Reynolds y Péclet suficientemente altos, dichos campos escalan como leyes de defecto de funciones de potencia de la distancia de la pared en el centro del canal. De la misma forma, se ha obtenido un escalado de la velocidad en la capa logarítmica para el caso de Reτ = 10000, obteniendo la clásica función logarítmica para la velocidad media y una función potencial para los momentos de órdenes superiores. Las leyes de escalado han sido validadas con los datos obtenidos en las DNS, obteniendo una precisión excelente. Por último, se han realizado una simulaciones de flujo de Couette con el número de Prandtl de aire, Pr = 0.71, y números de Reynolds de fricción de valores Reτ = 180, 250 y 500. El principal objetivo era el estudio de las estructuras coherentes que se forman en estos flujos de Couette. En concreto, se ha visto que las intensidades turbulentas dependen del número y tamaño de las estructuras. Es por ello que se necesita como mínimo una anchura del dominio computacional de 6πh para que las estadísticas sean independientes. Una última serie de simulaciones ha sido llevada a cabo considerando flujo estratificado. El objetivo era estudiar si las estructuras de Couette persisten en este tipo de flujos. Para un Reτ = 500 y el número de Prandtl de aire, Pr = 0.71, se han tomado valores del número de Richardson de fricción, Riτ = 0.5, 1.65 y 2.90. Para los dos casos con el número de Richardson de fricción más alto, las estructuras del flujo de Couette se debilitan hasta el punto de ser casi inexistentes. Las principales estadísticas de las simulaciones se encuentran disponibles en la base de datos del grupo de investigación, la cual está abierta a la comunidad científica y se puede acceder desde el siguiente enlace http://personales.upv.es/serhocal/ / [CA] El principal objectiu d'aquesta tesis és l'estudi de fluids tèrmics turbulents en canals per obtenir un major coneixement sobre el fenomen de la turbulència. Per a això, s'ha realitzat un estudi des de el punt de vista de la mecànica de fluids computacional, més concretament, s'ha utilitzat la tècnica de les simulacions numèriques directes (DNS de les seues sigles en anglès). La idea principal de les simulacions realitzades ha sigut ampliar l'estat de l'art actual, en lo referent al dos paràmetres principals que caracteritzen un flux: el número de Reynolds de fricció, Reτ , i el número de Prandtl, Pr. Dos configuracions del fluid han sigut utilitzades: flux de Poiseuille i flux de Couette, sent la primera el principal focus de l'estudi. En quant al camp de temperatures, s'ha utilitzat una condició de contorn mixta i s'ha considerat com a un escalar passiu. Així doncs, els primers números de Reynolds de fricció simulats per a un flux de Poisuille han sigut Reτ = 500, 1000 i 2000, per a números de Prandtl que varien des de 0.007 (metalls fosos) fins 10 (aigua), passant per 0.71 que és el valor més utilitzat per se aquest el número de Prandtl de l'aire. A més, s'ha realitzat una simulació con Reτ = 5000 i Pr = 0.71, la qual és la DNS tèrmica con el número de Reynolds de fricció més alt fins avui. Destacar que per a números de Prandtl més alts, s'ha observat que el valor màxim de la variància de la temperatura és constant. Això té un important benefici en l'escalat cerca de la paret dels terminis de dissipació i difusió viscosa del balanç d'energia de θ′'+. Per últim, en lo referent a simulacions de fluxos de Poiseuille, s'ha estudiat el cas isotèrmic amb Reτ = 10000, el qual és el major DNS de un canal turbulent, obtenint per primera vegada en una DNS una capa logarítmica perfectament desenvolupada en el camp de les velocitats. Un estudi teòric , basat en les simetries de Lie, ha sigut portat a terme en paral·lel a les simulacions. El principal objectiu ha sigut la generació de lleis d'escalat, basades en primers principis, del camp de velocitats, temperatura i els moments d'altos ordres de ambdós camps. El resultat és que per a números de Reynolds i Péclet suficientment alts, aquests camps escalen com lleis de defecte de funcions de potència de la distància de la paret en el centre del canal. De la mateixa manera, s'ha obtingut un escalat de la velocitat en la capa logarítmica per al caso de Reτ = 10000, obtenint la clàssica funció logarítmica per a la velocitat mitjana i una funció potencial per als moments d'ordres superiors. Les lleis d'escalat han sigut validades amb les dades obtingudes en les DNS, obtenint una precisió excel·lent. Per últim, s'ha realitzat una simulació de fluxos de Couette amb el número de Prandtl d l'aire, Pr = 0.71, i números de Reynolds de fricció de valors Reτ = 180, 250 i 500. El principal objectiu era l'estudi de les estructures coherents que es formen en aquests fluxos de Couette. Concretament, s'ha vist que les intensitats turbulentes depenen del número i mesura de les estructures. Es per això que es necessita com a mínim una amplada del domini computacional de 6π h per a que les estadístiques siguin independents. Una última sèrie de simulacions ha sigut feta considerant el flux estratificat. L'objectiu era estudiar si les estructures de Couette persisteixen en aquest tipus de fluxos. Per a un Reτ = 500 i el número de Prandtl d'aire, Pr = 0.71, s'han agafat valors del número de Richardson de fricció, Riτ = 0.5, 1.65 i 2.90. Per als dos casos amb el número de Richardson de fricció més alt, les estructures de flux de Couette es debiliten fins al punt de ser casi inexistents. Les principals estadístiques de les simulacions es troben disponibles en les bases de dades del grup d'investigació, el qual està obert a la comunitat científica i es pot accedir des de el següent enllaç http://personales.upv.es/serhocal/ / [EN] The main objective of this thesis is the study of thermal turbulent channel flows to obtain a greater knowledge about the phenomenon of turbulence. For this, a study has been carried out from the point of view of computational fluid mechanics, specifically, the technique of direct numerical simulations (DNS) has been used. The main idea of the simulations conducted has been to expand the current state of the art, in relation to the two main parameters that characterize the flow: the friction Reynolds number, Reτ, and the Prandtl number, Pr. Two flow configurations have been used: Poiseuille flow and Couette flow, the former being the main focus of the study. Regarding the temperature field, a mixed boundary condition has been used and it has been considered as a passive scalar. Thus, the simulated friction Reynolds numbers for a Poisuille flow have been Reτ = 500, 1000 and 2000, for Prandtl numbers that vary from 0.007 (molten metals) to 10 (water), passing through 0.71 which is the value more used because this is the Prandtl number of the air. In addition, a simulation has been carried out with Reτ = 5000 and Pr = 0.71, which is the thermal DNS with the highest friction Reynolds number to date. It should be noted that for the highest Prandtl numbers, it has been observed that the maximum value of the variance of the temperature is constant. This has an important benefit in scaling near the wall of the dissipation and viscous diffusion budget terms of the τ'+. Finally, with regard to Poiseuille flow simulations, the isothermal case has been studied with Reτ = 10000, which is the highest DNS of a turbulent channel flow, obtaining for the first time in a DNS a perfectly developed logarithmic layer in the velocity field. A theoretical study, based on Lie symmetries, has been carried out in parallel to the simulations. The main objective has been the generation of scald laws, based on first principles, of the field of velocity, temperature and high order moments of both fields. The result is that for sufficiently high Reynolds and Péclet numbers, these fields scale as defect laws of power functions of the distance to the wall in the center of the channel. In the same way, a scaling of the speed in the logarithmic layer has been obtained for the case of Reτ = 10000, obtaining the classic logarithmic function for the average velocity and a potential function for the moments of higher orders. The scaling laws have been validated with the data obtained in the DNS, obtaining excellent precision. Finally, a set of Couette flow simulations have been carried out with the Prandtl number of air, Pr = 0.71, and Reynolds friction numbers of Reτ = 180, 250 and 500. The main objective was the study of coherent structures that are formed in these Couette flows. Specifically, it has been seen that turbulent intensities depend on the number and size of the structures. For this reason, a minimum width of the computational domain of 6πh is required for the statistics to be independent. A last series of simulations has been carried out considering stratified flow. The objective was to study whether Couette structures persist in this type of flow. Thus, for a Reτ = 500 and the Prandtl number of air, Pr = 0.71, the values of the friction Richardson number have been varied according to, Riτ = 0.5, 1.65 and 2.90, for each simulation. For the two cases with the highest friction Richardson number, the Couette flow structures weaken to the point of being almost non-existent. The main statistics of the simulations are available in the research group's database, which is open to the scientific community and can be accessed from the following link http://personales.upv.es/serhocal/ / Esta tesis ha recibido la ayuda de MINECO/FEDER proyecto ENE2015-71333-R. / Alcántara Ávila, F. (2021). Study of the thermal field of turbulent channel flows via Direct Numerical Simulations [Tesis doctoral]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/180122 / TESIS

[ES] Debido a las regulaciones en materia de emisiones y CO2 la industria automotriz a desarrollado diferentes tecnologías innovadoras. Estas tecnologías incluyen combustibles alternativos y nuevos modos de combustión, entre otros. De aquí surge la necesidad del desarrollo de nuevos métodos para el control de la combustión en estas condiciones mencionadas. Por este motivo, en este trabajo se han desarrollado diferentes modelos e indicadores orientados al diagnóstico y control de la combustión tanto en condiciones normales como anormales. Para los casos de combustión normal, se ha desarrollado un modelo de combustión, cuyo objetivo es estimar la media de la evolución de la fracción de la masa quemada y la presión dentro del cilindro. Se implementó un observador, basado en la señal de knock, con la finalidad de mejorar la estimación en condiciones transitorias y poder aplicar así el modelo a diferentes tipos de combustibles. También se presenta un modelo de variabilidad cíclica, en el cual, a partir del modelo de combustión, se propaga una distribución en dos de los parámetros de dicho modelo. Ambos modelos han sido aplicados para un motor de encendido provocado y un motor de combustión de encendido por chorro turbulento. En los casos de combustión anormal, se ha incluido un análisis de la resonancia dentro de la cámara de combustión, en donde también se desarrollaron dos modelos capaces de estimar la evolución de la resonancia. Estos modelos, tanto para condiciones normales como anormales, se utilizaron para el diagnóstico de la combustión. Por una parte, para la detección de knock, en donde tres estrategias de detección de knock fueron desarrolladas: dos basadas en el sensor de presión en cámara y una en el sensor de knock. Por otra parte, se realizó una aplicación de un modelo de resonancia para la mejora de la estimación de la masa atrapada a partir de la resonancia. Finalmente, para mostrar el potencial de los modelos de diagnóstico, dos aplicaciones a control se desarrollaron: una para el control de knock a través de la actuación de la chispa, y otra para el control de gases residuales, a través de la actuación de la distribución variable, realizando paralelamente una optimización de la combustión a través de la actuación de la chispa. / [CA] Impulsada per les regulacions en matèria d'emissions i CO2 la indústria automotriu a desenvolupat diferents tecnologies inovadore. Aquestes tecnologies inclouen combustibles alternatius i nous modes de combustió, entre altres. D'ací sorgix la necessitat posar en pràctica nous mètodes per al control de la combustió. En aquest context, el present trevall proposa diferents models i indicadors orientats al diagnòstic i control de la combustió tant en condicions normals com anormals. Per als casos de combustió normal, es proposa un model de combustió, l'objectiu del qual és estimar la mitjana de l'evolució de la fracció de la massa cremada i la pressió dins del cilindre. Es va implementar un observador, basat en la senyal de knock, amb la finalitat de millorar l'estimació en condicions transitòries i poder aplicar així el model a diferents tipus de combustibles. També es presenta un model de variabilitat cíclica, en el qual, a partir del model de combustió, es propaga una distribució en dos dels parametres del dit model. Ambdós models han sigut aplicats a un motor d'encesa provocada i un motor de combustió d'encesa per doll turbulent. Als casos de combustió anormal, s'ha inclos un anàlisi de la ressonància dins de la cambra de combustió, on també es van desenvolupar dos models capaços d'estimar l'evolució de la ressonància. Aquests models, tant per a condicions normals com anormals, s'utilitzen per al diagnòstic de la combustió. Per una part, per a la detecció de knock, on tres estratègies de detecció de knock s'han desenvolupat: dues basades en el sensor de pressió a la cambra de combustió i una altra basada en el sensor de knock. Per altra part, es va realitzar una aplicació d'un model de ressonància per a la millora de l'estimació de la massa atrapada a partir de la ressonància. Finalment, per a mostrar el potencial dels models de diagnòstic, dos aplicacions de control es van desenvolupar: una per al control de knock a través de l'actuació de l'espurna, i una altra per al control de gasos residuals, a través de l'actuació de la distribució variable, realitzant paral·lelament una optimització de la combustió a través de l'actuació de l'espurna. / [EN] The need to satisfy emissions and CO2 regulations is pushing the automotive industry to develop different innovative technologies. These technologies include alternative fuels and new modes of combustion, among others. Therefore, the need for the development of new methods for combustion control in these mentioned conditions arises. For this reason, in this work different models and indicators have been developed aimed at the diagnosis and control of combustion in both normal and abnormal conditions. For normal combustion cases, a combustion model has been developed, the objective of this model is to estimate the mean of evolution of the mass fraction burned and the in-cylinder pressure. An observer had been implemented, based on knock sensor signal, in order to improve the estimation in transient conditions and also to be able to make use of the model with different fuels. A cyclic variability model is also presented, where from the combustion model, a probability distribution is propagated over two of the parameters of such model. Both models had been applied for a spark ignition engine and a turbulent jet ignition combustion engine. For the abnormal combustion cases, an analysis of the resonance within the combustion chamber had been included, where two models capable of estimating the evolution of the resonance were also developed. These models, for both normal and abnormal conditions, were used for the diagnosis of combustion: on the one hand, for knock recognition, where three knock detection strategies were developed: two based on the in-cylinder pressure sensor and one on the knock sensor. On the other hand, an application of a resonance model was carried out in order to improve the estimation of the trapped mass from the resonance excitation. Finally, to show the potential of such models and applications, two control strategies were developed: one for the control of knock through the actuation of the spark advance, and a second for the control of residual gases, through the actuation of the variable valve timing, while optimizing the combustion through the actuation of the spark advance. / El trabajo desarrollado en esta tesis ha sido posible gracias a la financiación de la Generalitat Valenciana y el fondo social europeo a través de la beca 132 GRISO- LIAP/2018/132 y BEFPI/2021/042. / Jiménez, IA. (2022). Development of Combustion Indicators for Control of Multi-Fuel Engines Based on New Combustion Concepts [Tesis doctoral]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/181561 / TESIS

[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

<p>This is a thesis work proposed by Sweco System in order to carry out a study related to the heating system of a circular fuel oil storage tank or cistern. The study tank is a 23m diameter and 18m height with a storage capacity of around 7500m3 of Eo5 heavy fuel oil. The content ought to be at a minimum storage temperature of 50ºC so that the fuel oil is fluid enough and operation labors can be adequately performed. In fact, these types of heavy fuel oils have fairly high viscosities at lower temperatures and the heating and pumping system can be compromised at temperatures below the pour point. For this purpose a heating system is installed to maintain the fluid warm. So far the system was operated by an oil burner but there are plans to its replacement by a District Heating-heat exchanger combo. Thereby, tank heating needs, flow and thermal patterns and heat transfer within it are principally studied.</p><p> </p><p>Tank boundaries are studied and their thermal resistances are calculated in order to dimension heat supply capacity. The study implies Finite Elements (Comsol Multiphysics) and Finite Volume (Fluent) analysis to work out some stationary heat transfer by conduction cases on some parts and thermal bridges present on these boundaries. Afterwards both cooling and heating processes of the fuel oil are studied using several strategies: basic models and Computational Fluid Dynamics (CFD). CFD work with Fluent is focused on optimizing inlet and outlet topologies. Understanding the cooling process is sought as well; Fluent CFD transient models are simulated in this way as well. Additionally the effect of filling levels is taken into account leading to a multiphase (fuel oil and air) flow cases where especially heating coupling of both phases is analyzed.</p><p> </p><p>Results show that maximum heat supply needs are around 80kW when the tank temperature is around 60ºC and 70kW when it is around 50ºC. Expectedly the main characteristic of the flow turns out to be the buoyancy driven convective pattern. K-ε turbulence viscous models are applied to both heating and cooling processes showing thermal stratification, especially at the bottom of the tank. Hotter fluid above follows very complex flow patterns. During the heating processes models used predict fairly well mixed and homogenous temperature distribution regardless small stratification at the bottom of the tank. In this way no concrete inlet-outlet configuration shows clear advantages over the rest. Due to the insulation of the tank, low thermal conductivity of the fluid and vast amount of mass present in the tank, the cooling process is slow (fluid average temperature drops around 5.7 ºC from 60ºC in 15 days when the tank is full and ambient temperature is considered to be at -20ºC) and lies somewhere in the middle between the solid rigid and perfect mixture cooling processes. However, due to stratification some parts of the fluid reach minimum admissible temperatures much faster than average temperature does. On the other hand, as expected, air phase acts as an additional thermal resistance; anyhow the cooling process is still faster for lower filling levels than the full one.</p> / <p>El presente proyecto fue propuesto por Sweco Systems para llevar a cabo un estudio relacionado con el sistema de calefacción de una cisterna o tanque de almacenamiento de fuel oil circular. Dicho tanque tiene 23 m de diámetro  y 18 m de altura con una capacidad de almacenamiento de alrededor de 7500 m³ de Eo5 fuel oil pesado. El contenido mantenerse a una temperatura mínima de 50 ºC de manera que el fuel oil es suficientemente fluido para que las labores de operación puedan ser ejecutadas adecuadamente. De hecho, estos tipos de fuel oil pesado tienen altas viscosidades a bajas temperaturas y, por tanto, tanto los sistemas de calefacción y como el de bombeo pueden verse comprometidosr a temperaturas por debajo del pour point. Con este fin un sistema de calefacción es instalado para mantener el fluido suficientemente caliente. Hasta el momento, el sistema era operado por un quemador de fuel, sin embargo, hay planes que éste sea sustituido por un combo intercambiador de calor-District Heating. Por lo tanto, principalmente son estudiadas las necesidades de calefacción así como los flujos térmicos y fluidos.</p><p>Se estudian las fronteras del tanque, y sus respectivas resistencias térmicas son calculadas con el fin de dimensionar la capacidad necesaria de suministro de calor. El estudio implica Elementos Finitos (Comsol Multiphysics) y Volúmenes Finitos (Fluent) para elaborar análisis estacionarios de transferencia de calor por conducción en algunos casos. Existen puentes térmicos en las paredes y su importancia es también anallizada. Posteriormente se estudian tanto los procesos de calentamiento y enfriamiento del fuel oil utilizando diversas estrategias: modelos básicos y Dinámica de Fluidos Computacional (CFD). El trabajo con CFD se centra en la optimización de topologías de entradas y salidas del sistema. También es solicitado entender el proceso de enfriamiento; En este sentido, se simulan modelos CFD transitorios de Fluent. Además, el efecto de los niveles de llenado se tiene en cuenta dando lugar a estudios de flujo multifase (fuel oil y aire), haciendo hincapié en el análisis de acoplamiento de transferencia de calor entre las dos fases.</p><p>Los resultados muestran que las necesidades de calefacción máximas son de alrededor de 80kW cuando la temperatura del tanque es de alrededor de 60 º C y 70kW cuando está alrededor de 50 ºC. Como era de esperar, la principal característica de este tipo de flujos es la convección natural resultante de las fuerzas de flotabilidad. Se aplican modelos turbulentos k-ε a los procesos de calentamiento y enfriamiento, mostrando estratificación térmica, sobre todo en la parte inferior de la cisterna. El líquido más caliente que se sitúa encima muestra complejos patrones de flujo. Durante los procesos de calentamiento, los modelos utilizados predicen un buen mezclado y distribución homogénea de la temperatura independientemente de esta pequeña estratificación en la parte inferior de la cisterna. De esta manera, ninguna concreta configuración de entradas-salidas simuladas muestra claras ventajas sobre el resto. Debido al aislamiento de la cisterna, la baja conductividad térmica del fluido y la gran cantidad de masa presente en el tanque el proceso de enfriamiento es lento (la temperatura media del fluido desciende 5.7 º C desde 60 º C en 15 días cuando el tanque está lleno y la temperatura ambiente es de -20 º C) y se encuentra en algún lugar en medio de los procesos de enfriamiento del sólido rígido y perfecta mezcla. Sin embargo, debido a la estratificación, algunas partes el líquido alcanzan la temperatura mínima admisible mucho más rápido que la media de temperatura. Por otra parte, como se esperaba, la fase de aire actúa como una resistencia térmica adicional, de todos modos, el proceso de enfriamiento es aún más rápido para niveles de llenado más bajos que el lleno.</p>

This is a thesis work proposed by Sweco System in order to carry out a study related to the heating system of a circular fuel oil storage tank or cistern. The study tank is a 23m diameter and 18m height with a storage capacity of around 7500m3 of Eo5 heavy fuel oil. The content ought to be at a minimum storage temperature of 50ºC so that the fuel oil is fluid enough and operation labors can be adequately performed. In fact, these types of heavy fuel oils have fairly high viscosities at lower temperatures and the heating and pumping system can be compromised at temperatures below the pour point. For this purpose a heating system is installed to maintain the fluid warm. So far the system was operated by an oil burner but there are plans to its replacement by a District Heating-heat exchanger combo. Thereby, tank heating needs, flow and thermal patterns and heat transfer within it are principally studied.   Tank boundaries are studied and their thermal resistances are calculated in order to dimension heat supply capacity. The study implies Finite Elements (Comsol Multiphysics) and Finite Volume (Fluent) analysis to work out some stationary heat transfer by conduction cases on some parts and thermal bridges present on these boundaries. Afterwards both cooling and heating processes of the fuel oil are studied using several strategies: basic models and Computational Fluid Dynamics (CFD). CFD work with Fluent is focused on optimizing inlet and outlet topologies. Understanding the cooling process is sought as well; Fluent CFD transient models are simulated in this way as well. Additionally the effect of filling levels is taken into account leading to a multiphase (fuel oil and air) flow cases where especially heating coupling of both phases is analyzed.   Results show that maximum heat supply needs are around 80kW when the tank temperature is around 60ºC and 70kW when it is around 50ºC. Expectedly the main characteristic of the flow turns out to be the buoyancy driven convective pattern. K-ε turbulence viscous models are applied to both heating and cooling processes showing thermal stratification, especially at the bottom of the tank. Hotter fluid above follows very complex flow patterns. During the heating processes models used predict fairly well mixed and homogenous temperature distribution regardless small stratification at the bottom of the tank. In this way no concrete inlet-outlet configuration shows clear advantages over the rest. Due to the insulation of the tank, low thermal conductivity of the fluid and vast amount of mass present in the tank, the cooling process is slow (fluid average temperature drops around 5.7 ºC from 60ºC in 15 days when the tank is full and ambient temperature is considered to be at -20ºC) and lies somewhere in the middle between the solid rigid and perfect mixture cooling processes. However, due to stratification some parts of the fluid reach minimum admissible temperatures much faster than average temperature does. On the other hand, as expected, air phase acts as an additional thermal resistance; anyhow the cooling process is still faster for lower filling levels than the full one. / El presente proyecto fue propuesto por Sweco Systems para llevar a cabo un estudio relacionado con el sistema de calefacción de una cisterna o tanque de almacenamiento de fuel oil circular. Dicho tanque tiene 23 m de diámetro  y 18 m de altura con una capacidad de almacenamiento de alrededor de 7500 m3 de Eo5 fuel oil pesado. El contenido mantenerse a una temperatura mínima de 50 ºC de manera que el fuel oil es suficientemente fluido para que las labores de operación puedan ser ejecutadas adecuadamente. De hecho, estos tipos de fuel oil pesado tienen altas viscosidades a bajas temperaturas y, por tanto, tanto los sistemas de calefacción y como el de bombeo pueden verse comprometidosr a temperaturas por debajo del pour point. Con este fin un sistema de calefacción es instalado para mantener el fluido suficientemente caliente. Hasta el momento, el sistema era operado por un quemador de fuel, sin embargo, hay planes que éste sea sustituido por un combo intercambiador de calor-District Heating. Por lo tanto, principalmente son estudiadas las necesidades de calefacción así como los flujos térmicos y fluidos. Se estudian las fronteras del tanque, y sus respectivas resistencias térmicas son calculadas con el fin de dimensionar la capacidad necesaria de suministro de calor. El estudio implica Elementos Finitos (Comsol Multiphysics) y Volúmenes Finitos (Fluent) para elaborar análisis estacionarios de transferencia de calor por conducción en algunos casos. Existen puentes térmicos en las paredes y su importancia es también anallizada. Posteriormente se estudian tanto los procesos de calentamiento y enfriamiento del fuel oil utilizando diversas estrategias: modelos básicos y Dinámica de Fluidos Computacional (CFD). El trabajo con CFD se centra en la optimización de topologías de entradas y salidas del sistema. También es solicitado entender el proceso de enfriamiento; En este sentido, se simulan modelos CFD transitorios de Fluent. Además, el efecto de los niveles de llenado se tiene en cuenta dando lugar a estudios de flujo multifase (fuel oil y aire), haciendo hincapié en el análisis de acoplamiento de transferencia de calor entre las dos fases. Los resultados muestran que las necesidades de calefacción máximas son de alrededor de 80kW cuando la temperatura del tanque es de alrededor de 60 º C y 70kW cuando está alrededor de 50 ºC. Como era de esperar, la principal característica de este tipo de flujos es la convección natural resultante de las fuerzas de flotabilidad. Se aplican modelos turbulentos k-ε a los procesos de calentamiento y enfriamiento, mostrando estratificación térmica, sobre todo en la parte inferior de la cisterna. El líquido más caliente que se sitúa encima muestra complejos patrones de flujo. Durante los procesos de calentamiento, los modelos utilizados predicen un buen mezclado y distribución homogénea de la temperatura independientemente de esta pequeña estratificación en la parte inferior de la cisterna. De esta manera, ninguna concreta configuración de entradas-salidas simuladas muestra claras ventajas sobre el resto. Debido al aislamiento de la cisterna, la baja conductividad térmica del fluido y la gran cantidad de masa presente en el tanque el proceso de enfriamiento es lento (la temperatura media del fluido desciende 5.7 º C desde 60 º C en 15 días cuando el tanque está lleno y la temperatura ambiente es de -20 º C) y se encuentra en algún lugar en medio de los procesos de enfriamiento del sólido rígido y perfecta mezcla. Sin embargo, debido a la estratificación, algunas partes el líquido alcanzan la temperatura mínima admisible mucho más rápido que la media de temperatura. Por otra parte, como se esperaba, la fase de aire actúa como una resistencia térmica adicional, de todos modos, el proceso de enfriamiento es aún más rápido para niveles de llenado más bajos que el lleno.