Dissertations / Theses on the topic 'Computational Fluid Dynamics (CFD) model'

Thesis (MEng)--Stellenbosch University, 2015.<br>ENGLISH ABSTRACT: Modern Computational Fluid Dynamics modelling (CFD) are becoming common design and analysis tools in the engineering field. Nowadays, project designs involve the use of CFD techniques along with physical scale modelling to analyse the complex rapidly varied and turbulent flows which would not be easily analysed by physical modelling. In particular, the consideration and/or use of CFD modelling in the Hydraulic Engineering field remains on the increase. Apart from being used for comparison with other design techniques, CFD may in future become a standalone modelling technique in hydraulic structures design. This research aims to use CFD models to validate the simulation of the flow over two ogee dam spillways which are installed in the Hydraulic Laboratory of Stellenbosch University. To achieve this simulation of the flow which involves an interaction between water and air, the flow behaviour has been mapped by the Volume of Fluid (VOF) and the realisable "𝑘−𝜀" turbulence numerical models. The Volume of Fluid (VOF) and the realisable "𝑘−𝜀" models simulate the free surface of two-phase flow and the flow turbulence, respectively. Firstly, the study embarks with details on the actual design approaches of a typical ogee dam spillway. It subsequently presents the geometry and dimensions of the physical models, the testing procedure and the experimental test results achieved from this modelling exercise. For CFD modelling, a commercially available Computational Fluid Dynamics (CFD) package, Ansys-Fluent, was used. To model the physical model, the use of Reynolds-averaged Navier-Stokes equations in combination with the realisable k-ε eddy-viscosity closure model was adopted. The process of CFD model development and the underlying theory of it are discussed in this thesis. Different test scenarios including steady and fully hydrodynamic states simulation for two and three-dimensional geometries were considered in this simulation to achieve the most accurate results. In order to determine the required mesh size, the mesh sensitivity tests were conducted on the 2 dimensional and 3 dimensional models. Finally, the pressure readings and water levels produced by numerical models are discussed through a validation process by comparing the CFD model results with the results obtained from physical models. The outcome proved that CFD models are able to map the behaviour of both flow phases since they exhibited a close correlation to those achieved in the physical models. Even though some slight differences in values were revealed, the graphical trend remains reasonably similar for all test results.<br>AFRIKAANSE OPSOMMING: Moderne gerekenariseerde vloeidinamika numeriese modelle (CFD) word deesdae dikwels deur ingenieurs gebruik. Projekontwerpe sluit tans die gebruik van CFD tegnieke asook fisiese skaalmodellering in om komplekse, vinnig-veranderede en turbulente vloei te ontleed. Hierdie tipe vloeie is moeilik om met fisiese modellering te ontleed. Die gebruik van CFD numeriese modelle in hidrouliese ingenieurswese is besig om toe te neem, Die bevindinge van CFD modelering word tans vergelyk met die bevindinge van ander ontwerptegnieke, maar in die toekoms mag dit moontlik gebruik word as die enigste modelleringstegniek in hidrouliese struktuurontwerp. Die doel met hierdie navorsing is om CFD modelering te gebruik om die vloei oor twee ogee-vormige afvoergeute wat in die hidrouliese labrotorium van die Universiteit van Stellenbosch ge-installeer is, te ondersoek. Ten einde hierdie vloei, wat die interaksie tussen water en lug insluit, te simuleer, is die vloeigedrag deur ”volume van vloeistof” (VOF) en die "𝑘−𝜀" turbulensie numeriese modules, gemodeleer. Die VOF en "𝑘−𝜀" numeriese modules simuleer onderskeidelik die vry oppervlakte vloei van die twee-fase vloei en turbulente vloei. Die ontwerp van ’n tipiese ”ogee”-tipe dam oorloop word bespreek, gevolg deur die beskrywing van die geometrie van die fisiese modelle, die toetsprosedure en die eksperimentele toetsresultate. Vir die CFD modellering is die CFD pakket, Ansys-Fluent, gebruik. Vir die simulering van die fisiese model is die Reynolds-gemiddeld Navier-Stokes vergelykings tesame met die k-ε eddy-viskositeit geslote module gebruik. Die proses van CFD ontwikkeling en die onderliggende teorie daarvan word bespreek. Verskillende toets-scenario’s wat 2D en 3D simulasies insluit, uitgevoer. Ten einde die toepaslike berekeningsrooster grootte vir die numeriese model te verkry, is sensitiewiteitstoetse uitgevoer op die twee- en drie-dimensionele numeriese modelle. Laastens is die CFD numeries gesimuleerde drukke en die watervlakke met die van die fisiese modelle vergelyk om die akkuraatheid van die CFD resultate te verkry. Die uitkomstes het getoon dat CFD modelle gebruik kan word om die gedrag van albei vloei fases te simuleer aangesien dit goed vergelyk het met die uitkomstes van die fisiese modellering. Daar was wel klein verskille in die druk waardes, maar die tendense in drukverspreiding was ooreenstemmend.

In this thesis a methodology is presented to optimise non–linear mathematical models in numerical engineering applications. The method is based on biological evolution and uses known concepts of genetic algorithms and evolutionary compu- tation. The working principle is explained in detail, the implementation is outlined and alternative approaches are mentioned. The optimisation is then tested on a series of benchmark cases to prove its validity. It is then applied to two different types of problems in computational engineering. The first application is the mathematical modeling of turbulence. An overview of existing turbulence models is followed by a series of tests of different models applied to various types of flows. In this thesis the optimisation method is used to find improved coefficient values for the k–ε, the k–ω-SST and the Spalart–Allmaras models. In a second application optimisation is used to improve the quality of a computational mesh automatically generated by a third party software tool. This generation can be controlled by a set of parameters, which are subject to the optimisation. The results obtained in this work show an improvement when compared to non–optimised results. While computationally expensive, the genetic optimisation method can still be used in engineering applications to tune predefined settings with the aim to produce results of higher quality. The implementation is modular and allows for further extensions and modifications for future applications.

We have studied gas flow and particle deposition in a realistic three-dimensional model of the bronchial tree, extending from the trachea to the segmental bronchi (7th airway generation for the most distal ones) using Computational Fluid Dynamics (CFD). The model is based on the morphometrical data of Horsfield et al. [J. Appl. Physiol. 31: 207-217, 1971] and on bronchoscopic and CT images, which give the spatial 3D-orientation of the curved ducts. It incorporates realistic angles of successive branching planes. Steady inspiratory flow varying between 50cm³/s and 500cm³/s was simulated as well as deposition of spherical aerosol particles (1 to 7&<br>Doctorat en sciences appliquées<br>info:eu-repo/semantics/nonPublished

Heart failure occurs when the heart is not capable to pump blood at a sufficient rate to meet the demands of the body. Depending on the health of the heart, doctors may recommend a heart transplant, but finding a suitable donor is often a long duration process and the patient might be at an advance condition or the patient is not adequate for a heart transplant. In such cases Ventricular assist devices (VAD) are implemented. The purpose of a VAD is to aid the heart to pump the correct amount of blood, by doing so it relives the load that is put on the heart while giving the patient a chance for recovery. This study focuses on observing the hemodynamic effects of implementing a left ventricular assist device (LVAD) along the aortic arch and main arteries. Thrombi creation and transportation is other subject included in the study, due to the fact that thrombi can obstruct blood flow to critical arteries, manly carotid and vertebral. Occlusion of these can lead to a stroke with devastating effects on the neurocognitive functions and even death. A multi-scale CFD analysis a patient specific geometry model is used as well as a lumped system which provides the correct conditions in order to simulate the whole cardiovascular system. The main goal of the study is to understand the difference in flow behavior created by the unsteady pulsatile boundary conditions. The model described in this work has a total cardiac output of 7.0 Liters/ minute, this for a healthy heart. Two cardiac output splits are used to simulate heart failure conditions. The first split consists of 5 Liters/minute flowing through the LVAD cannula and 2 Liters/minute via the aortic root. The second scenario is when heart iv failure is critical, meaning that zero flow is being output by the left ventricle, thus a split of 7 Liter/minute trough the LVAD cannula and 0 Liters/minute traveling through the aortic root. A statistical analysis for the thrombi motion throughout the patient aortic arch was performed in order to quantify the influence that pulsatile flow has on the particles being track. Spherical particles of 2mm, 4mm and 5mm were released and accounted in the statistical analysis for each of the two split configurations. The study focuses on particles that escaped on the outlet boundaries of the upper arteries (Right Carotid, Left Carotid, and Vertebral). Results exhibit the statistical comparison of means for each particle diameter as well as for the overall probability for the steady and unsteady flow condition.<br>M.S.M.E.<br>Masters<br>Mechanical and Aerospace Engineering<br>Engineering and Computer Science<br>Mechanical Engineering; Thermo-Fluids

Although oscillating water column type wave energy devices are nearing the stage of commercial exploitation, there is still much to be learnt about many facets of their hydrodynamic performance. This research uses the commercially available FLUENT computational fluid dynamics flow solver to model a complete OWC system in a two dimensional numerical wave tank. A key feature of the numerical modelling is the focus on the influence of the front wall geometry and in particular the effect of the front wall aperture shape on the hydrodynamic conversion efficiency. In order to validate the numerical modelling, a 1:12.5 scale experimental model has been tested in a wave tank under regular wave conditions. The effects of the front lip shape on the hydrodynamic efficiency are investigated both numerically and experimentally and the results compared. The results obtained show that with careful consideration of key modelling parameters as well as ensuring sufficient data resolution, there is good agreement between the two methods. The results of the testing have also illustrated that simple changes to the front wall aperture shape can provide marked improvements in the efficiency of energy capture for OWC type devices.

A high fidelity computational fluid dynamics approach to perform direct simulations of ship maneuvers is presented in this thesis. The approach uses dynamic overset grids with a hierarchy of bodies to enable arbitrary motions between objects, and overcome the difficulties in simulation of the moving rudder and rotating propeller. To better resolve propeller/rudder interaction a Delayed Detached Eddy Simulation turbulence model based on Menter’s SST is used. The methodology was implemented in the general purpose RANS/DES/DDES research code REX, and is applied to the KRISO Container Ship (KCS) with moving rudder and rotating propeller in deep and shallow water. For the first time, a grid study is conducted for the self-propulsion condition for the propeller RPM, thrust, torque and lateral force, and for the roll and pitch motions, using grids of 8.7 (coarse), 24.6 (medium) and 71.3 (fine) million points. A grid study is also performed for the zigzag maneuver evaluating the maximum and minimum values of propeller thrust, torque and lateral force roll, pitch, yaw, roll rate, yaw rate and drift throughout the maneuver. An extensive comparison between predicted motions and forces of the direct simulations and the experimental data collected by Schiffbau-Versuchsanstalt Potsdam GmbH (SVA) and Flanders Hydraulics Research (FHR) are presented. While the results and comparisons with experimental data show that using direct CFD to compute modified and standard maneuvers with moving rudder and rotating discretized propeller is feasible, computational cost remains an impediment for many practical applications. Coupling a dynamic overset CFD solver with a potential propeller code can dramatically reduce the computational time to perform maneuvering simulations by using one order of magnitude larger time step than direct simulation. This thesis investigates the ability of a coupled CFD/potential propeller code approach to simulate maneuvers in ships, where the rudder is located downstream of the propeller. While the approach has been successfully applied to submarine maneuvers, in which the propeller wake is free of interference, the concept had not been evaluated before for cases where an object (the rudder) is immersed in the wake. The study is performed using the CFD code REX and the propeller code PUF-14. Performance of the coupled REX/PUF-14 approach is first tested studying propeller/rudder interaction, evaluating influence of the propeller/rudder gap size and rudder deflection on propeller performance curves and rudder forces, comparing against DDES simulations with a discretized rotating propeller. A grid study was performed for advance coefficient J=0.6 and a rudder angle δ=20 degrees for a propeller rudder gap of 0.2 times the rudder radius, with the resulting grid uncertainties for propeller thrust and torque coefficients suggesting that the effects of the grid changes are small for the present range of grid sizes. A 15/1 zigzag maneuver for the KCS container ship, in which case the rudder is very close downstream of the propeller, is then analyzed, and compared against discretized propeller simulations and experimental data. Self-propulsion coupled REX/PUF-14 results agree very well with experiments and discretized propeller simulations. Prediction of motions, forces and moments, and mean flow field with the coupled REX/PUF-14 approach are comparable to results obtained with discretized propeller simulations and agree with experiments well, though as implemented the coupled approach is unable to resolve tip vortices and other flow structures that interact with the rudder, potentially affecting prediction of flow separation. It can be concluded that coupled CFD/potential flow propeller approaches are an effective and economical way to perform direct simulation of surface ship maneuvers with CFD.

Thesis (MScEng)--Stellenbosch University, 2014.<br>ENGLISH ABSTRACT: The designs of pump intake structures are generally verified by means of a hydraulic model study to ensure the structure will operate without large air or swirl entraining vortices. This is an expensive and time consuming process. The popularity of Computational Fluid Dynamics (CFD) is ever increasing in the engineering world and the use thereof in the modelling of flows in pump intake structures is just one of the many areas of focus. The use of a validated CFD model in the design phase of pump intakes would be beneficial, and ideally would reduce the frequency of requiring hydraulic model studies. This study investigates the applicability of using a 3 dimensional CFD model as a design tool by validation with a physical model. A physical model of a simple pump intake was available and with a few modifications the flow conditions were significantly improved. The quantitative validation of the CFD model was carried out by means of comparing corresponding point velocity magnitudes. An Acoustic Doppler Velocimeter (ADV) was used in the physical model. The data obtained were successfully filtered to remove noise and other disturbances. The qualitative validation was done by means of photographs and observations. The photographs made use of illuminated tracer particles to identify flow patterns. The CFD model compares well qualitatively, but the velocity magnitudes are not yet sufficiently reliable. It is recommended that the CFD model can be used for qualitative studies, but future research should focus on the accuracy of the CFD model. Using higher resolution velocity measurements in the physical model by means of other types of instruments, a better comparison can be made, as well as enabling validation of the ADV readings.<br>AFRIKAANSE OPSOMMING: Die ontwerpe van die pomp inname strukture word gewoonlik bevestig deur middel van 'n hidrouliese model studie om lug-intrekking te minimeer. Dit is 'n duur en tydrowende proses . Die gebruik van Computational Fluid Dynamics ( CFD ), ‘n tipe numeriese modelle neem toe in die ingenieurswese wêreld. Die gebruik daarvan in die modellering van vloei in pomp inname strukture is net een van die vele areas van toepassing . Die gebruik van 'n betroubare CFD model in die ontwerp fase van pomp innames sal voordelig wees , en ideaal sou die noodaaklikheid van fisiese hidrouliese model studies verminder. Hierdie studie ondersoek die toepaslikheid van die gebruik van 'n 3 -dimensionele model CFD model as 'n ontwerp instrument deur bevestiging met 'n fisiese model . 'n Eenvoudige pomp-inname model was beskikbaar en het met ‘n paar veranderinge die vloei toestande aansienlik verbeter . Die kwantitatiewe bevestiging van die CFD model is uitgevoer deur middel van vergelyking van korresponderende punt snelheid groottes. In die fisiese model is die snelhede met 'n Akoestiese Doppler Velocimeter (ADV) gemeet. Die data wat verkry is suksesvol gefiltreer om geraas en ander versteurings te verwyder . Die kwalitatiewe bevestiging is gedoen deur middel van foto's en waarnemings. In die foto's is gebruik gemaak van verligte spoor deeltjies om vloei patrone te identifiseer. Die CFD model vergelyk kwalitatief goed, maar die snelheid groottes is nog nie voldoende betroubaar nie. Dit word aanbeveel dat die CFD model gebruik kan word vir kwalitatiewe studies, maar toekomstige navorsing behoort te fokus op die verbeteringe van dieakkuraatheid van die CFD model. Met behulp van 'n hoër resolusie snelheid metings van die fisiese model, deur middel van ander tipe instrumente, kan 'n beter vergelyking gemaak word , sowel as die geldigheid van die ADV lesings in staat te stel.

Hydropower plants with spillways exposed to high pressures because of large discharges can result in a dam failure which ultimately results in devastating consequences. Therefore properly designed and constructed spillways are important. In this project there is an interest in examining the pressure distribution on a guided wall, measuring the water level inside the spillway and examine the water flow characteristics. The main purpose was to numerically model the flow behaviors at two upper pool levels +29.82 m and +30.92 m in order to achieve a comparison between a numerical model and a physical model. A previous experiment was carried out at Vattenfall R&amp;D in Älvkarleby where a physical model was built. The main purpose of that study was to find an optimal form of a guide wall for free weir flow. This study served as guidelines for what should be investigated in this project. The methodology of the project began with a development of an initial CAD model created by Vattenfall R&amp;D. The CFD model could be constructed into a volume model containing also an upstream reservoir. Next, a mesh that could numerically calculate the flow behavior was created in order to enable execution of the calculations. Furthermore, a big part of the method was to adjust the settings in Fluent so that the numerical model could recreate the flow behaviors of the physical model tests. The results indicate that the numerical model is an accurate replica of the physical model. The deviation comprising the mass flow rate came as close as 0.6 % from the physical model for the lower upper pool level. Highest calculated pressure was positioned at the bottom of the guided wall for both upper pool levels, which the results conducted from the physical model also indicated. Furthermore, results regarding the water level at the guided wall confirmed same statement Vattenfall's examinations pointed out, that the roof placed at the guided wall is needed in order to avoid major splashing. The present project has shown great results regarding the flow characteristics, pressure distribution, and the water level at the guided wall. Therefore, it is considered trustworthy enough to be used as a tool for Vattenfall R&amp;D in future studies comprising the targeted spillway.

Orientador: Sávio Souza Venâncio Vianna<br>Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Química<br>Made available in DSpace on 2018-08-24T13:25:34Z (GMT). No. of bitstreams: 1 Ferreira_TatieleDalfior_M.pdf: 4562241 bytes, checksum: 69c742236e806f040cfc237f2ba91cf4 (MD5) Previous issue date: 2014<br>Resumo: Este trabalho tem como objetivo desenvolver um modelo matemático capaz de prever o tamanho de nuvem de gás inflamável formada em uma típica plataforma de petróleo considerando condições reais de ventilação e de operação de uma planta de processo. Para tanto, foi realizado um estudo de dispersão de gás inflamável (gás natural) na plataforma em questão utilizando Fluidodinâmica Computacional (CFD). Os resultados deste estudo de dispersão serviram como base para a construção do modelo matemático utilizando Metodologia de Superfície de Resposta. Tal modelo permite o cálculo do tamanho de nuvem de gás inflamável no ambiente estudado usando duas variáveis principais: a taxa não-dimensional de vazamento (que contabiliza a relação entre a taxa de vazamento de gás e a taxa de ventilação na plataforma) e a direção adimensional de vazamento (que computa a relação entre as direções de vazamento de gás e do vento). O modelo desenvolvido mostrou-se eficaz, pois foi capaz de prever com considerável grau de confiabilidade os tamanhos de nuvem de gás inflamável quando comparados aos valores fornecidos por simulações com CFD<br>Abstract: This work proposes the development of a mathematical correlation for prediction of flammable gas cloud size in a typical offshore module. Real conditions regarding the ventilation and process plant operation were considered. A dispersion study of natural gas release in the module was conducted using Computational Fluid Dynamics (CFD) and the state of art as far as the gas dispersion modelling is concerned. A mathematical model was built based on the numerical results and Response Surface Methodology (RSM). The approach comprises into a single mathematical model the most relevant independent variables. The response surface curves calculate the flammable gas cloud volume as a function of the non-dimensional leak rate (that concerns the ventilation and the gas release rate) and the non-dimensional leak direction (which comprises the wind direction and the leak direction). The developed model had proved to be effective. It was able to predict flammable gas volume and good agreement with CFD results was observed<br>Mestrado<br>Sistemas de Processos Quimicos e Informatica<br>Mestra em Engenharia Química

The production and transportation of hydrocarbon fluids in multiphase pipelines could be severely hindered by particulate solids deposit such as produced sand particles which accompany hydrocarbon production. Knowledge of the flow characteristics of solid particles in fluids transported in pipelines is important in order to accurately predict solid particles deposition in pipelines. This research thesis presents the development of a three-dimensional (3D) Computational Fluids Dynamics (CFD) modelling technique for the prediction of liquid-solids multiphase flow in pipes, with special emphasis on the flow in V-inclined pipe bends. The Euler-Euler (two-fluid) multiphase modelling methodology has been adopted and the multiphase model equations and closure models describing the liquid-solids flow have been implemented and calculated using the finite volume method in a CFD code software. The liquid phase turbulence has been modelled using a two-equation k−ε turbulence model which contains additional terms to account for the effects of the solid-particles phase on the multiphase turbulence structure. The developed CFD numerical framework has been verified for the relevant forces and all the possible interaction mechanisms of the liquid-solids multiphase flow by investigating four different numerical frameworks, in order to determine the optimum numerical framework that captures the underlying physics and covers the interaction mechanisms that lead to sand deposition and the range of sand transport flow regimes in pipes. The flow of liquid-sand in pipe has been studied extensively and the numerical results of sand concentration distribution across pipe and other flow properties are in good agreement with published experimental data on validation. The numerical framework has been employed to investigate the multiphase flow in V-inclined pipe bends of ±4o−6o, seemingly small inclined bend angles. The predicted results which include the sand segregation, deposition velocity and flow turbulence modulation in the pipe bend show that the seemingly small pipe bends have significant effect on the flow differently from that of horizontal pipes. The pipe bend causes abrupt local change in the multiphase flow characteristic and formation of stationary sand deposit in the pipe at a relatively high flow velocity. The threshold velocity to keep sand entrained in liquid in pipe bends is significantly higher than that required for flow horizontal pipes. A critical implication of this is that the correlations for predicting sand deposition in pipelines must account for the effect of pipe bend on flow characteristics in order to provide accurate predictions of the critical sand transport velocity (MTV) in subsea petroleum flowlines, which V-inclined pipe bends are inevitable due to seabed topology.

Rising fuel prices have been inflating the operating costs of the fishing industry. Trawl doors are used to hold the fishing net open during trawling operations, and they have a great influence on the fuel consumption of vessels. Improvements in the design of trawl doors could therefore contribute significantly to increased fuel efficiency. An efficient optimization algorithm using two- and three-dimensional (2D and 3D) computational fluid dynamics (CFD) models is presented. Accurate CFD models, especially 3D, are computationally expensive. The direct use of traditional optimization algorithms, which often require a large number of evaluations, can therefore be prohibitive. The proposed method is iterative and uses low-order local response surface approximation models as surrogates for the expensive CFD model to reduce the number of iterations. The algorithm is applied to the design of two types of geometries: a typical modern trawl door, and a novel airfoil-shaped trawl door. The results from the 2D design optimization show that the hydrodynamic efficiency of the typical modern trawl door could be increased by 32%, and the novel airfoil-shaped trawl door by 13%. When the 2D optimum designs for the two geometries are compared, the novel airfoil-shaped trawl door results to be 320% more efficient than the optimized design of the typical modern trawl door. The 2D optimum designs were used as the initial designs for the 3D design optimization. The results from the 3D optimization show that the hydrodynamic efficiency could be increased by 6% for both the typical modern and novel airfoil-shaped trawl doors. Results from a 3D CFD analysis show that 3D flow effects are significant, where the values for drag are significantly underestimated in 2D CFD models.

<p> This thesis presents a computational analysis of the performance of the model of a hydropower barge, consisting of 3 pairs of paddle wheels, in generating power from the flow of a river. The barge dimensions used in this research are one-thirtieth scale of the actual barge dimensions.</p><p> These new dimensions were used in the simulations to predict the power that can be generated, having carried out similitude between the scaled model and its full scale size. The dimensionless parameters employed in achieving similitude are Strouhal number, power coefficient, and pressure coefficient. The efficiency of the paddle wheel conforms to that obtained from most power turbines. The power generated was improved by the addition of bottom fin under the paddle wheel.</p><p> Computational Fluid Dynamics using ANSYS Fluent software was employed to simulate and generate the results. The simulation results generated will be compared to an experimental model that will be performed in the future based on the new scaled dimensions. This will help validate the viability of the paddle wheel, as proposed here, in the generation of power prior to designing the full scale sized barge.</p>

This master’s thesis was provided by the Swedish Defence Research Agency, FOI. The task is to test several RANS (Reynolds-averaged Navier-Stokes) models on two different case geometries and compare the results with LES and experimental data. The first is two dimensional, constructed for flow separation at a sharp edge. The second is three dimensional and flow separation occurs at a smooth surface. The models tested are implemented in the open source CFD (Computational Fluid Dynamics) program, OpenFOAM. OpenFOAM uses the finite volume method and the SIMPLE algorithm as solution procedure. The main flow features evaluated is the shape, position and size of the flow separation. Most of the models tested have problems describing the complex dynamics of flow separation in these particular cases. In addition to the simulations, the RANS k-epsilon turbulence model is presented and the RANS equations and the equation for the turbulent kinetic energy are derived from the Navier-Stokes equations. The theory behind wall functions is described and these equations together with the equations in the k-epsilon model are compared with the equations implemented in OpenFOAM.

The study and modelling of two-phase flow, even the simplest ones such as the bubbly flow, remains a challenge that requires exploring the physical phenomena from different spatial and temporal resolution levels. CFD (Computational Fluid Dynamics) is a widespread and promising tool for modelling, but nowadays, there is no single approach or method to predict the dynamics of these systems at the different resolution levels providing enough precision of the results. The inherent difficulties of the events occurring in this flow, mainly those related with the interface between phases, makes that low or intermediate resolution level approaches as system codes (RELAP, TRACE, ...) or 3D TFM (Two-Fluid Model) have significant issues to reproduce acceptable results, unless well-known scenarios and global values are considered. Instead, methods based on high resolution level such as Interfacial Tracking Method (ITM) or Volume Of Fluid (VOF) require a high computational effort that makes unfeasible its use in complex systems. In this thesis, an open-source simulation framework has been designed and developed using the OpenFOAM library to analyze the cases from microescale to macroscale levels. The different approaches and the information that is required in each one of them have been studied for bubbly flow. In the first part, the dynamics of single bubbles at a high resolution level have been examined through VOF. This technique has allowed to obtain accurate results related to the bubble formation, terminal velocity, path, wake and instabilities produced by the wake. However, this approach has been impractical for real scenarios with more than dozens of bubbles. Alternatively, this thesis proposes a CFD Discrete Element Method (CFD-DEM) technique, where each bubble is represented discretely. A novel solver for bubbly flow has been developed in this thesis. This includes a large number of improvements necessary to reproduce the bubble-bubble and bubble-wall interactions, turbulence, velocity seen by the bubbles, momentum and mass exchange term over the cells or bubble expansion, among others. But also new implementations as an algorithm to seed the bubbles in the system have been incorporated. As a result, this new solver gives more accurate results as the provided up to date. Following the decrease on resolution level, and therefore the required computational resources, a 3D TFM have been developed with a population balance equation solved with an implementation of the Quadrature Method Of Moments (QMOM). The solver is implemented with the same closure models as the CFD-DEM to analyze the effects involved with the lost of information due to the averaging of the instantaneous Navier-Stokes equation. The analysis of the results with CFD-DEM reveals the discrepancies found by considering averaged values and homogeneous flow in the models of the classical TFM formulation. Finally, for the lowest resolution level approach, the system code RELAP5/MOD3 is used for modelling the bubbly flow regime. The code has been modified to reproduce properly the two-phase flow characteristics in vertical pipes, comparing the performance of the calculation of the drag term based on drift-velocity and drag coefficient approaches.<br>El estudio y modelado de flujos bifásicos, incluso los más simples como el bubbly flow, sigue siendo un reto que conlleva aproximarse a los fenómenos físicos que lo rigen desde diferentes niveles de resolución espacial y temporal. El uso de códigos CFD (Computational Fluid Dynamics) como herramienta de modelado está muy extendida y resulta prometedora, pero hoy por hoy, no existe una única aproximación o técnica de resolución que permita predecir la dinámica de estos sistemas en los diferentes niveles de resolución, y que ofrezca suficiente precisión en sus resultados. La dificultad intrínseca de los fenómenos que allí ocurren, sobre todo los ligados a la interfase entre ambas fases, hace que los códigos de bajo o medio nivel de resolución, como pueden ser los códigos de sistema (RELAP, TRACE, etc.) o los basados en aproximaciones 3D TFM (Two-Fluid Model) tengan serios problemas para ofrecer resultados aceptables, a no ser que se trate de escenarios muy conocidos y se busquen resultados globales. En cambio, códigos basados en alto nivel de resolución, como los que utilizan VOF (Volume Of Fluid), requirieren de un esfuerzo computacional tan elevado que no pueden ser aplicados a sistemas complejos. En esta tesis, mediante el uso de la librería OpenFOAM se ha creado un marco de simulación de código abierto para analizar los escenarios desde niveles de resolución de microescala a macroescala, analizando las diferentes aproximaciones, así como la información que es necesaria aportar en cada una de ellas, para el estudio del régimen de bubbly flow. En la primera parte se estudia la dinámica de burbujas individuales a un alto nivel de resolución mediante el uso del método VOF (Volume Of Fluid). Esta técnica ha permitido obtener resultados precisos como la formación de la burbuja, velocidad terminal, camino recorrido, estela producida por la burbuja e inestabilidades que produce en su camino. Pero esta aproximación resulta inviable para entornos reales con la participación de más de unas pocas decenas de burbujas. Como alternativa, se propone el uso de técnicas CFD-DEM (Discrete Element Methods) en la que se representa a las burbujas como partículas discretas. En esta tesis se ha desarrollado un nuevo solver para bubbly flow en el que se han añadido un gran número de nuevos modelos, como los necesarios para contemplar los choques entre burbujas o con las paredes, la turbulencia, la velocidad vista por las burbujas, la distribución del intercambio de momento y masas con el fluido en las diferentes celdas por cada una de las burbujas o la expansión de la fase gaseosa entre otros. Pero también se han tenido que incluir nuevos algoritmos como el necesario para inyectar de forma adecuada la fase gaseosa en el sistema. Este nuevo solver ofrece resultados con un nivel de resolución superior a los desarrollados hasta la fecha. Siguiendo con la reducción del nivel de resolución, y por tanto los recursos computacionales necesarios, se efectúa el desarrollo de un solver tridimensional de TFM en el que se ha implementado el método QMOM (Quadrature Method Of Moments) para resolver la ecuación de balance poblacional. El solver se desarrolla con los mismos modelos de cierre que el CFD-DEM para analizar los efectos relacionados con la pérdida de información debido al promediado de las ecuaciones instantáneas de Navier-Stokes. El análisis de resultados de CFD-DEM permite determinar las discrepancias encontradas por considerar los valores promediados y el flujo homogéneo de los modelos clásicos de TFM. Por último, como aproximación de nivel de resolución más bajo, se investiga el uso uso de códigos de sistema, utilizando el código RELAP5/MOD3 para analizar el modelado del flujo en condiciones de bubbly flow. El código es modificado para reproducir correctamente el flujo bifásico en tuberías verticales, comparando el comportamiento de aproximaciones para el cálculo del término d<br>L'estudi i modelatge de fluxos bifàsics, fins i tot els més simples com bubbly flow, segueix sent un repte que comporta aproximar-se als fenòmens físics que ho regeixen des de diferents nivells de resolució espacial i temporal. L'ús de codis CFD (Computational Fluid Dynamics) com a eina de modelatge està molt estesa i resulta prometedora, però ara per ara, no existeix una única aproximació o tècnica de resolució que permeta predir la dinàmica d'aquests sistemes en els diferents nivells de resolució, i que oferisca suficient precisió en els seus resultats. Les dificultat intrínseques dels fenòmens que allí ocorren, sobre tots els lligats a la interfase entre les dues fases, fa que els codis de baix o mig nivell de resolució, com poden ser els codis de sistema (RELAP,TRACE, etc.) o els basats en aproximacions 3D TFM (Two-Fluid Model) tinguen seriosos problemes per a oferir resultats acceptables , llevat que es tracte d'escenaris molt coneguts i se persegueixen resultats globals. En canvi, codis basats en alt nivell de resolució, com els que utilitzen VOF (Volume Of Fluid), requereixen d'un esforç computacional tan elevat que no poden ser aplicats a sistemes complexos. En aquesta tesi, mitjançant l'ús de la llibreria OpenFOAM s'ha creat un marc de simulació de codi obert per a analitzar els escenaris des de nivells de resolució de microescala a macroescala, analitzant les diferents aproximacions, així com la informació que és necessària aportar en cadascuna d'elles, per a l'estudi del règim de bubbly flow. En la primera part s'estudia la dinàmica de bambolles individuals a un alt nivell de resolució mitjançant l'ús del mètode VOF. Aquesta tècnica ha permès obtenir resultats precisos com la formació de la bambolla, velocitat terminal, camí recorregut, estela produida per la bambolla i inestabilitats que produeix en el seu camí. Però aquesta aproximació resulta inviable per a entorns reals amb la participació de més d'unes poques desenes de bambolles. Com a alternativa en aqueix cas es proposa l'ús de tècniques CFD-DEM (Discrete Element Methods) en la qual es representa a les bambolles com a partícules discretes. En aquesta tesi s'ha desenvolupat un nou solver per a bubbly flow en el qual s'han afegit un gran nombre de nous models, com els necessaris per a contemplar els xocs entre bambolles o amb les parets, la turbulència, la velocitat vista per les bambolles, la distribució de l'intercanvi de moment i masses amb el fluid en les diferents cel·les per cadascuna de les bambolles o els models d'expansió de la fase gasosa entre uns altres. Però també s'ha hagut d'incloure nous algoritmes com el necessari per a injectar de forma adequada la fase gasosa en el sistema. Aquest nou solver ofereix resultats amb un nivell de resolució superior als desenvolupat fins la data. Seguint amb la reducció del nivell de resolució, i per tant els recursos computacionals necessaris, s'efectua el desenvolupament d'un solver tridimensional de TFM en el qual s'ha implementat el mètode QMOM (Quadrature Method Of Moments) per a resoldre l'equació de balanç poblacional. El solver es desenvolupa amb els mateixos models de tancament que el CFD-DEM per a analitzar els efectes relacionats amb la pèrdua d'informació a causa del promitjat de les equacions instantànies de Navier-Stokes. L'anàlisi de resultats de CFD-DEM permet determinar les discrepàncies ocasionades per considerar els valors promitjats i el flux homogeni dels models clàssics de TFM. Finalment, com a aproximació de nivell de resolució més baix, s'analitza l'ús de codis de sistema, utilitzant el codi RELAP5/MOD3 per a analitzar el modelatge del fluxos en règim de bubbly flow. El codi és modificat per a reproduir correctament les característiques del flux bifàsic en canonades verticals, comparant el comportament d'aproximacions per al càlcul del terme de drag basades en velocitat de drift flux model i de les basades en coe<br>Peña Monferrer, C. (2017). Computational fluid dynamics multiscale modelling of bubbly flow. A critical study and new developments on volume of fluid, discrete element and two-fluid methods [Tesis doctoral no publicada]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/90493<br>TESIS

Efficient analysis and storage of data is an integral but often challenging task when working with computation fluid dynamics mainly due to the amount of data it can output. Methods centered around the proper orthogonal decomposition were used to analyze, compress, and model various simulation cases. Two different high-fidelity, time-accurate turbomachinery simulations were investigated to show various applications of the analysis techniques. The first turbomachinery example was used to illustrate the extraction of turbulent coherent structures such as traversing shocks, vortex shedding, and wake variation from deswirler and rotor blade passages. Using only the most dominant modes, flow fields were reconstructed and analyzed for error. The reconstructions reproduced the general dynamics within the flow well, but failed to fully resolve shock fronts and smaller vortices. By decomposing the domain into smaller, independent pieces, reconstruction error was reduced by up to 63 percent. A new method of data compression that combined an image compression algorithm and the proper orthogonal decomposition was used to store the reconstructions of the flow field, increasing data compression ratios by a factor of 40.The second turbomachinery simulation studied was a three-stage fan with inlet total pressure distortion. Both the snapshot and repeating geometry methods were used to characterize structures of static pressure fluctuation within the blade passages of the third rotor blade row. Modal coefficients filtered by frequencies relating to the inlet distortion pattern were used to produce reconstructions of the pressure field solely dependent on the inlet boundary condition. A hybrid proper orthogonal decomposition method was proposed to limit burdens on computational resources while providing high temporal resolution analysis.Parametric reduced order models were created from large databases of transient and steady conjugate heat transfer and airfoil simulations. Performance of the models were found to depend heavily on the range of the parameters varied as well as the number of simulations used to traverse that range. The heat transfer models gave excellent predictions for temperature profiles in heated solids for ambitious parameter ranges. Model development for the airfoil case showed that accuracy was highly dependent on modal truncation. The flow fields were predicted very well, especially outside the boundary layer region of the flow.

Fuel sprays play a major role in order to achieve the required combustion characteristics and pollutant emissions reduction on internal combustion engines, and thus, an accurate prediction of its behavior is required to perform reliable engine combustion and pollutant simulations. A great effort both on experimental and theoretical studies of spray atomization and dispersion has been performed in the latest years. As a result, Computational Fluid Dynamics (CFD) calculations have become a standard tool not only for spray physics understanding but also for design and optimization of engine spray systems. However, spray modeling in its different uses in the Internal Combustion Engine (ICE) context is still nowadays a challenging task due to the complex interrelated phenomena taking place, some of them still not fully understood. Primary atomization and secondary breakup, droplet collision, coalescence and vaporization, turbulent interactions between phases have to be solved under high Reynolds (so they are turbulent) and Weber numbers conditions due to the high speed (~500 m/s) and small nozzle diameter (~100 µm) imposed by current engine injection systems technologies. Moreover, Taylor numbers cover a wide range, according to the composition of the injected liquid. Those conditions make experimental observations quite challenging and probably insufficient, especially in the very near nozzle region, where primary atomization takes place. Most of the CFD spray models are currently based on the Discrete Droplet Method. The continuous liquid jet is discretized into 'blobs' or 'parcels', which consists in a number of droplets with the same characteristics. A Lagrangian method is applied to track the liquid phase parcels, which are subject to breakup according to atomization models mainly based on the linear instability theory proposed by Reitz and later extended by Huh and Gosman for liquid turbulence effects to be considered. This approach has been successfully applied b<br>Khuong ., AD. (2012). The Eulerian-Lagrangian Spray Atomization (ELSA) Model of the Jet Atomization in CFD Simulations: Evaluation and Validation [Tesis doctoral no publicada]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/17237<br>Palancia

The aim of the diploma thesis was creating a decarbonisation model of lime, implementation the model into CFD tool ANSYS Fluent and to test the decarbonisation model in a model of a real reactor. The required model was based on assumptions for a Shrinking Core Model (SCM). The main objective of this work was the non-catalytic conversion of substances and the search for the most used mathematical models for calcination. The CFD calculation, the sensitivity analysis and the Fluente parametric study were used. Data on the composition of gas flow, temperature, pressure and mass flow of limestone particles were selected for input variables. The particle model called Multiple Surface Reactions (MSR), which is a standard part of Fluent, was used at first. Subsequently, a UDF which was based on the SCM assumptions was written in the programming language C. The results of the CFD calculation were compared with the experimental values from the dissertation. It has been found that the MSR is sufficiently precise for calculation purposes but neglects the internal diffusion of CO2 through the CaO layer which forms behind the reaction front during calcination. It was found that it is possible to solve the flow with ongoing calcination without the need to know the parameters of the Arrhenian equation if the UDF is used. The created UDF incorporates the influence of intraparticular CO2 diffusion on the overall reaction rate.

The aim of this master thesis is to simulate, in a multiphase environment (air and water), a free surface wave by making use of the computational fluid dynamics software “OpenFoam” and evaluate the hydrodynamic forces and physical parameters acting on subsea structures. It will be discussed how the hydrodynamic load and other physical parameters vary, according to water depth, when the action of waves and current act together and impact structures of different cross-sections and geometry that lie on the seabed. A comparison between the obtained software results and the analytical solutions will be presented to assess the reliability of “OpenFoam” and the model created for such a problem, trying to give an explanation of the fluid dynamics around the elements when running the simulations. In particular, starting from the multiphase solver “InterFoam”, will be shown how the “wave” tutorial case, available in “OpenFoam V5.0”, has been modified and adapted in order to implement objects on the seabed and eventually determine the hydrodynamic load on submerged structures. For the implemented elements, it will be examined the horizontal and the lift force, trying to verify and validate the model, making a comparison between the theoretical hydrodynamic load, available from the literature and the model post-processed outcomes, extrapolated by OpenFoam object functions. The evaluation of the force coefficients (Drag, Lift, Added Mass) will be also made on the same structures to see which level of influence certain parameters have to respect with others. It will also be presented the Morison’s theory for the evaluation of waves induced load that represents not only the theoretical load with which the model results have been compared to testify the project but it also plays the key role concept of all the following aspects.

This document is a report summering the master thesis work dealing with the Computational Fluid Dynamic (CFD) study of the flow around a high-speed train. The model is a scaled 1:50 generic train with two cars, one inter-car gap and simplified bogies. A platform is set on the side of the train since one of the aim of the study is to look at the consequences of the phenomena in the wake on people or objects standing on the platform. The slipstream is one of this phenomena, it is due to the fact that the viscous air is dragged when the train is passing. If too strong, it can move or destabilize people or objects on the platform. In addition of the slipstream study, a velocity profile study, a drag and lift coefficients analyze as well as a Q-factor study and a frequency study have been realized. Some results of these different studies are compared with the ones obtained on the same model with a Delayed Detached Eddy Simulation (DDES). Since the flow is turbulent, for those different studies, the flow has been simulated with a Reynolds Averaged Navier-Stokes equation model (RANS) which is the k-ω SST model for the turbulence. The study of the slipstream allowed to calculate the Technical Specification for Interoperability (TSI) which must not be higher that the European Union requirement set at 15.5 m/s, the result obtained is 8.1 m/s which is then lower than the limit. The velocity profile shows similarities with the DDES results even though it is less detailed. The same conclusion is done for the Q-plot where is clearly visible the two counter-rotating vortices in the wake. Finally, a Fast Fourier Transform algorithm has been applied to instantaneous velocity results in the wake of the train in order to get the frequency of the aerodynamic phenomena in that wake. The main frequency is 25 Hz and corresponds to a Strouhal number of 0.1, quite closed to the results obtained with DDES which is 0.085. The results of the RANS and DDES are reasonably similar and by regarding at the large difference between the cell numbers (respectively 8 500 000 and 20 000 000) it can be conclude that in some ways the RANS model can be preferred at the DDES to save time for the computation but it does not contain the small scales resolved by the DDES.

A physical model and a computational fluid dynamic (cfd) model (HYDRO-3D) have been developed to simulate the effects of novel maturation pond configurations, inlet and outlet positions, channels and critical environmental factors (wind speed and direction) on the hydraulic efficiency of maturation ponds at a small sewage treatment works in southern England. The ponds form the final, polishing stage of an otherwise conventional percolating filter treatment works. The final effluent of the plant joins a stream that discharges directly to a bathing beach. Thus the overall objective of the project is to maximise the efficiency of pathogen indicator removal and meet European bathing water standards. The specific aim of the study reported here was to calibrate the physical model and to assess its reliability as a tool for assessing and predicting hydraulic performance of ponds. Work focused on the comparison of simulation outputs from the cfd model and full scale system performance in order to increase the reliability of both models as design tools. The full scale ponds under study comprised three parallel, tertiary stage maturation ponds. A Geopacks flow meter was used to calibrate the ratio of water and wind velocities in the full scale pond channels and the model. Additionally floats and novel drogues were used to measure low flow velocities within the physical model. Initial tracer experiments using salt and rhodamine-WT were conducted to define the hydraulic retention time (HRT) of the physical model in still air and with top and bottom inlets and outlets. In the physical model the effect of back-mixing and dead zones can be readily observed and photographed, whereas this is more difficult to display in the outputs of computational models. Potential sources of disagreement are discussed in the thesis, but their concurrence suggests they both have a valuable role to play in the future design of ponds aimed at maximising hydraulic efficiency. The study also confirmed field observations which indicated that wind tends to be damaging rather than beneficial to pond performance.

Tato diplomová práce se zaměřuje na modelování proudění tekutiny v paralelních distribučních systémech, kde hraje důležitou roli rovnoměrnost distribuce tekutin. Pro tento účel je vytvořen vlastní CFD kód. Kód je napsán v programovacím jazyce Java a používá ke zlepšení výkonu knihovny třetích stran, které se vyznačují přímým přístupem ke hardwarovým a systémovým prostředkům. Kód se také vyznačuje určitými zjednodušeními, u nichž se očekává, že sníží výpočetní časy. Vliv použitých zjednodušujících opatření je vyhodnocen porovnáním výsledků simulací proudění na několika geometriích s údaji získanými z podrobných modelů CFD. Geometrie použitých svazků trubek se odlišují různými uspořádáními toku a trubek a také různým počtem trubek.

<p>Nuclear reactors release small amounts of radioactivity during their normal operations. The most common method of calculating the dose to the public that results from such releases uses Gaussian Plume models. We are investigating these methods using CAP88-PC, a computer code developed for the Environmental Protection Agency (EPA) in the USA that calculates the concentration of radionuclides released from a stack using Pasquill stability classification. A buoyant or momentum driven part is also included. The uptake of the released radionuclide by plants, animals and humans, directly and indirectly, is then calculated to obtain the doses to the public. This method is well established but is known to suffer from many approximations and does not give answers that are accurate to be better than 50% in many cases. More accurate, though much more computer-intensive methods have been developed to calculate the movement of gases&nbsp<br>using fluid dynamic models. Such a model, using the code FLUENT can model complex terrains and will also be investigated in this work. This work is a preliminary study to compare the results of the traditional Gaussian plume model and a fluid dynamic model for a simplified case. The results indicate that Computational Fluid Dynamics calculations give qualitatively similar results with the possibility of including much more effects than the simple Gaussian plume model.</p>

The high turbine inlet temperature of modern gas turbines poses a challenge to the material used in the turbine blading of the primary stages. Mechanical failure mechanisms are more pronounced at these high temperatures, setting the lifetime of the blade. It is therefore crucial to obtain accurate local metal temperature predictions of the turbine blade. Accurately predicting the external heat transfer coefficient (HTC) distribution of the blade is therefore of uttermost importance. At present time, Siemens Energy uses the boundary layer code TEXSTAN for this purpose. The limitations coupled to such codes however make them less applicable for the complex flow physics involved in the hot gas path of turbine blading. The thesis therefore aims at introducing CFD for calculating the external HTC. This includes conducting an extensive literature study to find and validate a suitable methodology. The literature study was centered around RANS modeling, reviewing how the calculation of the HTC has evolved and the performance of some common turbulence and transition models. From the literature study, the SST k − ω model in conjunction with the γ − Reθ transition model, the v2 − f model and the Lag EB k − ε model were chosen for the investigation of a suitable methodology. The validation of the methodology was based on the extensively studied LS89 vane linear cascade of the von Karman Institute. In total 13 test cases of the cascade were chosen to represent a wide range of flow conditions. Both a periodic model and a model of the entire LS89 cascade were tested but there were great uncertainties whether or not the correct flow conditions were achieved with the model of the entire cascade. It was therefore abandoned and a periodic model was used instead. The decay of turbulence intensity is not known in the LS89 cascade. This made the case difficult to model since the turbulence boundary conditions then were incomplete. Two approaches were attempted to handle this deficiency, where one was ultimately found invalid. It was recognized that the Steelant-Dick postulation could be used in order to find a turbulent length scale which when specified at the inlet, lead to fairly good agreement with data of the HTC. The validation showed that the SST γ − Reθ model performs relatively well on the suction side and in transition onset predictions but worse on the pressure side for certain flow conditions. The v2 − f model performed better on the pressure side and on a small portion of the suction side. Literature emphasized the importance of obtaining proper turbulence characteristics around the vane for accurate HTC-predictions. It was found that the results of the validation step could be closely coupled to this statement and that further work is needed regarding this. Further research must also be done on the Steelant-Dick postulation to validate it as a reliable method in prescribing the inlet length scale.

The master´s thesis deals with the numerical modeling of flow in the emergency spillways of small dams. It is solved the flow rate capacity of a weir and of spillway called „duckbill-type“. Sharp crested weirs and long crested weirs are modeled using numerical simulations in ANSYS. Hydraulically complex spillway called „duckbill-type“ was modeled in the Flow-3D. The results of the discharge coefficients are compared with data in the hydraulic literature.

This diploma thesis discuss a flow in a vicinity of exhaust valve using computational of fluid dynamics. In a light of current state of the problem solution this approach is forward but very sophisticated. Intention of author is description of multidimensional boundary of characteristic variables which determinates the flow behaviour for given geometry and boundary conditions. Technical knowlegde of fluid flow in vicinity of exhaust valve allow to design geometry which provide more effective cylinder flush out berofe the exhaust stroke will become. This process reduce quantity of loss work of piston and improve its effective pointers.

In the interest of furthering the understanding of hemodynamics, this study has developed a method for modeling fluid mechanics behavior in individual human carotid arteries. A computational model was constructed from magnetic resonance (MR) data of a phantom carotid bifurcation model, and relevant flow conditions were simulated. Results were verified by comparison with previous in vitro experiments. The methodology was extended to create subject-specific carotid artery models from geometry data and fluid flow boundary conditions which were determined from MR and phase contrast MR (PCMR) scans of human subjects. The influence of subject-specific boundary conditions on the flow field was investigated by comparing a model based on measured velocity boundary conditions to a model based on the assumption of idealized velocity boundary conditions. It is shown that subject-specific velocity boundary conditions in combination with a subject-specific geometry and flow waveform influence fluid flow phenomena associated with plaque development. Comparing a model with idealized Womersley flow boundary conditions to a model with subject-specific velocity boundary conditions demonstrated the importance of employing inlet and flow division data obtained from individual subjects in order to predict accurate, clinically relevant, fluid flow phenomena such as low wall shear stress areas and negative axial velocity regions. This study also illustrates the influence of the bifurcation geometry, especially the flow divider position, with respect to the velocity distribution of the common carotid artery on the development of flow characteristics. Overall it is concluded that accurate geometry and velocity measurements are essential for modeling fluid mechanics in individual human carotid arteries for the purpose of understanding atherosclerosis in the carotid artery bifurcation.

A Computational Fluid Dynamics/Computational Fluid Dynamics (CFD/CSD) coupling interface was developed to obtain aeroelastic solutions of a morphing rotor. The methodology was implemented in Fully Unstructured Navier-Stokes (FUN3D) solver, which communicates aerodynamic forces on the blade surface to University of Michigan’s Nonlinear Active Beam Solver (UM/NLABS) and then imports structural deflections of the blade surface during each time step. Development of this methodology adds the capability to model elastic rotors with flexible airfoils. The method was validated through an aerodynamic work analysis, comparison of sectional blade loads and deflections with experimental data, and two-dimensional stability analyses for pitch/plunge flutter and camber flutter. Computational simulations were performed for a rotor in forward flight with the CFD/CSD solver and with a comprehensive CSD solver using finite-state (F-S) aerodynamics, and results were compared. Prescribed three-per-revolution camber deflections were then applied, and solutions of the CFD/CSD and comprehensive CSD computations indicated that three-per-revolution camber actuation has the potential to minimize hub forces and moments with deflections as small as 0.25%c. In anticipation of active rotor experiments inside enclosed facilities, the capability of CFD for accurately simulating flow inside enclosed volumes was examined. It was determined that URANS models are not suitable for rotor simulations in an enclosed facility, and components that are a distance of two to three rotor radii from the hub were also observed to have a large influence on recirculation and performance.

This study involves the modeling of small cold-gas (N2) thrusters nozzle and plume flows, their interactions with spacecraft surfaces and the induced pressure environment. These small cold-gas thrusters were used for pitch, yaw and roll control and were mounted on the bottom of the conical Environmental Monitor Payload (EMP) suborbital spacecraft. The pitch and yaw thrusters had 0.906 mm throat diameter and 4.826 mm exit diameter, while the roll thrusters had 1.6 mm throat diameter and 5.882 mm exit diameter. During thruster firing, at altitudes between 670 km and 1200 km, pressure measurements exhibited non-periodic pulses (Gatsonis et al., 1999). The pressure sensor was located inside the EMP and was connected to it's sidewall with a 0.1-m long, 0.022-m diameter tube and the pressure pulses appeared instantaneously with the firings for thrusters without a direct line-of-sight with the sensor entrance. Preliminary analysis showed that the plume of these small EMP thrusters undergoes transition from continuous to rarefied. Therefore, nozzle and plume simulations are performed using a combination of Navier-Stokes and Direct Simulation Monte Carlo codes. This study presents first a validation of the Navier-Stokes code Rampant used for the continuous EMP nozzle and plume simulations. The first Rampant validation example involves a two-dimensional axisymetric freejet expansion and is used to demonstrate the use of Bird's breakdown parameter. Results are compared favorably with those of Bird (1980) obtained through the method of characteristics. The second validation example involves three-dimensional plume simulations of a NASA thruster. This nitrogen nozzle has a throat diameter of 3.18 mm, an exit diameter of 31.8 mm, half-angle of 20 degrees, stagnation temperature of 699 K, stagnation pressure of 6,400 Pa. Simulation results are compared favorably with previous Navier-Stokes and Direct Simulation Monte Carlo numerical work. The third validation example involves three-dimensional simulations of Rothe's (1970) nozzle that has a throat diameter of 2.5 mm, an exit diameter of 20.3 mm, half-angle of 20 degrees, operating at stagnation temperature of 300 K and pressure of 1975 Pa. Numerical results also compared favorably to experimental data. The combined Navier-Stokes/DSMC approach and the EMP simulation results are presented and discussed. The continuous part of the EMP nozzle and plume flow is modeled using the three-dimensional Navier-Stokes Rampant code. The Navier-Stokes domain includes the geometry of the nozzle and the EMP base until transition of the continuous flow established by Bird's breakdown parameter. The rarefied part of the plume flow is modeled using the Direct Simulation Monte Carlo code DAC. Flowfield data obtained inside the breakdown surface from the Navier-Stokes simulation are used as inputs to the DSMC simulations. The DSMC domain includes the input surface and the EMP spacecraft geometry. The combined Navier-Stokes/DSMC simulations show the complex structure of the plume flow as it expands over the EMP surfaces. Plume reflection and backflow are demonstrated. The study also summarizes findings presented by Gatsonis et al. (2000), where the DSMC predictions at the entrance of the pressure sensor are used as inputs to a semi-analytical model to predict the pressure inside the sensor. It is shown that the pressure predictions for the pitch/yaw thrusters are close to the measurements. The plume of a pitch or yaw thruster reaches the pressure sensor after expanding on the EMP base. The pressure predicted for the roll thruster is larger that the measured. This is attributed to the uncertainty in the roll thruster location on the EMP base resulting, in the simulation, in a component of direct flow to the sensor.

Abstract In ladle metallurgy, gas stirring and the behaviour of the slag layer are very important for alloying and the homogenization of the steel. When gas is injected through a nozzle located at the bottom of the ladle into the metal bath, the gas jet exiting the nozzle breaks up into gas bubbles. The rising bubbles break the slag layer and create an open-eye. The size of the open-eye is very important as the efficiency of the metal-slag reactions depend on the interaction between the slag and steel created during the stirring process, and information about the position and size of the open-eye is important for effective alloying practice. Moreover, the open-eye has an effect on the energy balance since it increases heat losses. In this study, experimental measurements and numerical simulations were performed to study the effect of different operating parameters on the formation of the open-eye and mixing time in a water model and industrial ladle. Experimental measurements were performed to study the effect of the gas flow rate, slag layer thickness, slag layer densities and number of porous plugs in a 1/5 scale water model and in a 150-ton steelmaking ladle. For numerical modelling, a multi-phase volume of fluid (VOF) model was used to simulate the system including the behaviour of the slag layer. The numerical simulation of the open-eye size and mixing time was found to be in good agreement with the experimental data obtained from the water model and data obtained from the industrial measurements<br>Tiivistelmä Senkkametallurgiassa kaasuhuuhtelu ja kuonakerroksen käyttäytyminen ovat tärkeitä teräksen seostamisen ja homogenisoinnin näkökulmasta. Senkan pohjalla sijaitsevasta suuttimesta puhallettava kaasu hajoaa kupliksi, jotka rikkovat kuonakerroksen ja muodostavat avoimen silmäkkeen. Avoimen silmäkkeen koko on yhteydessä voimakkaampaan kuonan emulgoitumiseen, joka tehostaa metallisulan ja kuonan välisiä reaktioita. Tietoa avoimen silmäkkeen paikasta ja koosta tarvitaan myös tehokkaaseen seostuspraktiikkaan. Avoin silmäke vaikuttaa lisäksi prosessin energiataseeseen lisäten sen lämpöhäviöitä. Tässä tutkimuksessa tutkittiin kokeellisesti ja laskennallisesti erilaisten operointiparametrien vaikutusta avoimen silmäkkeen muodostumiseen vesimallissa ja terässenkassateollisessa senkassa. Kokeellisia mittauksia tehtiin kaasuhuuhtelun, kuonakerroksen paksuuden, ja suuttimien määrän vaikutuksen tutkimiseksi 1/5-mittakaavan vesimallissa ja 150 tonnin terässenkassa. Numeerisessa mallinnuksessa systeemin ja siihen lukeutuvan kuonakerroksen käyttäytymisen simuloimiseen käytettiin volume of fluid (VOF) –monifaasimenetelmää. Avoimen silmäkkeen kokoon ja sekoittumisaikaan liittyvien numeeristen simulointien havaittiin vastaavan hyvin vesimallista ja teollisista mittauksista saatua kokeellista aineistoa

The process of shock-induced ignition has been investigated both computa- tionally and experimentally, with particular emphasis on the concept of radical farming. The first component of the investigation contained Computational Fluid Dynamic (CFD) calculations of an ignition delay study, a 2D pre-mixed flow over flat plate at a constant angle to the freestream, and through a generic 2D scramjet model. The focal point of the investigation however examined the complex 3D flow through a generic scramjet model. Five experimental test conditions were ex- amined over flow enthalpies from 3.4 MJ/kg to 6.4 MJ/kg. All test conditions simulated flight at 21000 metres ([symbol=almost equal to] 70000 ft), while the equivalent flight Mach number varied from approximately 8.5 at the lowest enthalpy, to approximately Mach 12 at the highest enthalpy condition. The presence of H2 fuel injected in the intake caused a separated region to form on the lower surface of the model at the entrance to the combustor. A fraction of the total mass of fuel was entrained in this separated region, providing long residence times, hence increased time for the chemical reactions that lead to ignition to occur. In addition, extremely high temperatures were found to exist between each fuel jet. Both fuel and air are present in these regions, therefore the chance of ignition in these regions is high. Streamlines passing through the recirculation zone ignited within this zone, while streamlines passing between the fuel jets ignited soon after entry into the combustor. The first instance of a pressure rise from combustion was observed on the centreline of the model where the reflected bow shock around the fuel jets crossed the centreline of the combus- tor. Upstream of this location the static pressure of the flow was too low for the chemical reactions that release heat to occur. The comparison between the experimental and computational results was lim- ited due to inaccuracies in modelling the thermal state of the gas in the CFD calculations. The gas was modelled as being in a state of thermal equilibrium at all times, which incorrectly models the freestream flow from the nozzle of the shock tunnel, and also the flow downstream of oblique shock wave within the scramjet model. As a result combustion occurs sooner in the CFD calculations than in the experimental result.

Computational fluid dynamics (CFD) has become an attractive tool in understanding the characteristic of air flow in the human lungs. Inter-subject variations make subject-specific simulations essential for understanding structure-function relationship, assessing lung function and improving drug delivery. However, currently the subject-specific CFD analysis remains challenging due, in large part to, two issues: construction of realistic deforming airway geometry and imposition of physiological boundary conditions. To address these two issues, we develop subject-specific, dynamic lung models by utilizing two or multiple volume multi-detector row computed tomography (MDCT) data sets and image registrations in this thesis. A mass-preserving nonrigid image registration algorithm is first proposed to match a pair of three-dimensional (3D) MDCT data sets with large deformations. A novel similarity criterion, the sum of squared tissue volume difference (SSTVD), is introduced to account for changes in intensity with lung inflation. We then demonstrate the ability to develop dynamic lung models by using a pair of lung volumes to account for deformations of airway geometries and subject-specific boundary conditions. The deformation of the airway geometry is derived by the registration-derived deformation field and subject-specific boundary condition is estimated from regional ventilation in a 3D and one-dimensional (1D) coupled multi-scale framework. Improved dynamic lung models are then proposed from three lung volumes by utilizing nonlinear interpolations. The improved lung models account for nonlinear geometry motions and time-varying boundary conditions during breathing. The capability of the proposed dynamic lung model is expected to move the CFD-based interrogation of lung function to the next plateau.

A low-emissions power generator comprising a solid oxide fuel cell coupled to a gas turbine has been developed by Rolls-Royce Fuel Cell Systems. As part of the cycle, a fraction of the unreacted fuel (the off-gas) and oxidizer streams is reacted in a burner, which is the main source of pollutant formation. In this thesis a computational model of the burner has been developed which captures the formation of NOx and the oxidation of CO. This model gives accurate predictions at low computational cost, making it suitable for use as a design tool in future burner design optimization through parametric studies. A key factor in increasing computational efficiency was the development of a reduced H2/CO/N2 kinetic mechanism; from a starting mechanism of 30 species to 10 and 116 reactions to 6. The results of laminar opposed-flow diffusion flames have been used to validate the reduced mechanism. Several different turbulent combustion models have been evaluated by creating an interface between the reduced kinetic mechanism and the commercial CFD solver FLUENT. Comparison of model predictions with well-characterized turbulent syngas flames, which share a similar fuel composition to the experimental work conducted on the off-gas burner, shows acceptable agreement. These studies have demonstrated the sensitivity of modelling constants. Improved predictions were achieved by calibrating these constants and including radiative heat losses. Following suitable modification to reflect the predominantly laminar flow present in the current burner design, the relevant modelling approaches were applied to the off-gas burner. Comparison was made to previous detailed measurements, showing that the important trends of NOx and CO are captured in general. The model was extended to high pressure conditions, similar to those in the actual off-gas burner, with the emissions predictions within design limits. The outcome of this work is a fast, accurate design tool for CFD which has capabilities to simulate beyond the laminar burner studied here. It may be applied to more general types of off-gas/syngas burners where turbulence-chemistry interaction is expected to be more significant.

The premise of the work presented here is to use a common analytical tool, Computational Fluid dynamics (CFD), along with a difference turbulence models. Eddy viscosity models as well as state-of-the-art Large Eddy Simulation (LES) were used to study the flow past bluff bodies. A suitable CFD code (CFX5.6b) was selected and implemented. Simulation of turbulent transport for the gas through the gaps of the randomly distributed spherical fuel elements (pebbles) was performed. Although there are a number of numerical studies () on flows around spherical bodies, none of them use the necessary turbulence models that are required to simulate flow where strong separation exists. With the development of high performance computers built for applications that require high CPU time and memory; numerical simulation becomes one of the more effective approaches for such investigations and LES type of turbulence models can be used more effectively. Since there are objects that are touching each other in the present study, a special approach was applied at the stage of building computational domain. This is supposed to be a considerable improvement for CFD applications. Zero thickness was achieved between the pebbles in which fission reaction takes place. Since there is a strong pressure gradient as a result of high Reynolds Number on the computational domain, which strongly affects the boundary layer behavior, heat transfer in both laminar and turbulent flows varies noticeably. Therefore, noncircular curved flows as in the pebble-bed situatio n, in detailed local sense, is interesting to be investigated. Since a compromise is needed between accuracy of results and time/cost of effort in acquiring the results numerically, selection of turbulence model should be done carefully. Resolving all the scales of a turbulent flow is too costly, while employing highly empirical turbulence models to complex problems could give inaccurate simulation results. The Large Eddy Simulation (LES) method would achieve the requirements to obtain a reasonable result. In LES, the large scales in the flow are solved and the small scales are modeled. Eddy viscosity and Reynolds stress models were also be used to investigate the applicability of these models for this kind of flow past bluff bodies at high Re numbers.

Lithium ion (Li-ion) batteries have been identified as a promising solution to meet the increasing demands for alternative energy in electric vehicles (EVs) and hybrid electric vehicle (HEVs). This work describes experimental and numerical study of thermal management of battery module consisting of cylindrical Li-ion cells, with an emphasis on the use of active control to achieve optimal cooling performance with minimal parasitic power consumption. The major contribution from this work is the first experimental demonstration (based on our review of archival journal and conference literature) and the corresponding analysis of active thermal control of battery modules. The results suggest that the active control strategy, when combined with reciprocating cooling flow, can reduce the parasitic energy consumption and cooling flow amount substantially. Compared with results using passive control with unidirectional cooling flow, the parasitic energy consumption was reduced by about 80%. This contribution was achieved in three steps, which was detailed in this dissertation in chapters 2, 3, and 4, respectively. In the first step, an experimental facility and a corresponding CFD model were developed to capture the thermal behavior of multiple battery cells. Based on the experimental and CFD results, a reduced-order model (ROM) was then developed for active monitoring and control purposes. In the second step, the ROM was parameterized and an observer-based control strategy was developed to control the core temperature of battery cells. Finally, based on the experimental facility and the ROM model, the active control of a battery module was demonstrated. Each of these steps represents an important facet of the thermal management problem, and it is expected that the results and specifics documented in this dissertation lay the groundwork to facilitate further study.<br>Ph. D.

A busca por uma nova matriz energética intensificou significativamente o desenvolvimento de fontes renováveis de energia. Dentre essas, o álcool teve grande destaque, dado que sua utilização industrial já é intensa no Brasil. Em concordância com a tendência industrial, o comportamento do álcool em vários fenômenos físico-químicos tornou-se foco de intensa pesquisa e desenvolvimento, o que levou a parcerias entre faculdades e empresas. O presente trabalho objetiva ao desenvolvimento de simulações numéricas que consigam descrever com exatidão o comportamento de chamas de sprays turbulentos de etanol. Para a verificação da aderência da simulação a processos reais, serão utilizados dados experimentais. Esses dados serviram como base comparativa de validação do modelo numérico. A descrição do spray segue a abordagem Euler-Lagrange, na qual a fase líquida segue uma abordagem lagrangeana e a fase gasosa uma abordagem euleriana. Foi considerado o acoplamento completo entre as duas fases, considerando o impacto da fase líquida no regime turbulento da fase gasosa. A turbulência da fase gasosa é descrita pelo modelo RSM (Reynolds Stress Model), baseado na média temporal e decomposição da turbulência, enquanto o escoamento interno da fase líquida é negligenciado. A chama do spray caracterizada pelo modelo flamelet, o qual utilizará um mecanismo de cinética química reduzido da oxidação do etanol.<br>The research for a new energy source to substitute fossil fuels intensified the development of renewable energy sources. Among all, ethanol has been given great attention, once it has been widely used in many Brazilian industrial branches. In agreement with industrial trends, the behavior of ethanol became the center of innumerous researches and developments, which led to partnerships between academic and private companies. This thesis aims the development of numerical simulations that describe with accurately the behavior of ethanol turbulent spray flames. To verify correlation between model results and real process behavior, experimental data will be used for validation of the numerical model. The spray physical description is made through the use of the Euler- Lagrange approach, in which the liquid phase is described by the lagrangean approach and the gaseous phase is described by the eulerian approach. There is full coupling between both phases. The turbulent flow of the gaseous phase is modeled through the use of the Reynolds Stress Model (RSM), based on Reynolds decomposition and time averaging process. The liquid phase has no internal flow, the droplets behave as a solid sphere translating through the domain (gaseous phase). The spray flame is modeled by the flamelet model which will use a simplified chemical mechanism.

In recent years, the increasing impacts of landslide hazards on human lives and lifeline facilities worldwide has advanced the necessity to find out both economically acceptable and useful techniques to predict the occurrence and destructive power of landslides. Though many projects exist to attain this goal, the current investigation set out to establish an understanding of the initiation and propagation mechanisms of landslides via numerical simulations, so that mitigation strategies to reduce the long-term losses from landslide hazards can be made. In this research, the Discrete Element Method (DEM) and Computational Fluid Dynamics (CFD) have been used to investigate the mechanical and hydraulic behaviour of granular materials involved in landslides. The main challenge is to provide rational analyses of large scale landslides via small scale numerical simulations. To solve this problem, dimensional analyses have been performed on a simple granular column collapse model. The influence of governing dimensionless groups on the debris runout distance and deposit height has been studied for the terrestrial and submerged granular flows. 3D DEM investigations of granular flows in plane strain conditions have been performed in this research. The input parameters of the DEM model have been calibrated by the numerical triaxial tests, based on which, the relationships between the microscopic variables and the macroscopic soil strength properties are analysed. Using the simple granular column collapse model, the influences of column aspect ratio, characteristic strain, model size ratio and material internal friction angle on the runout distance and deposit height of granular materials have been examined. Additionally, the deformation and energy evolution of dry granular materials are also discussed. The DEM-CFD coupling model has been employed to study the mechanical and hydraulic behaviour of highly mobilized terrestrial / submarine landslides. This model has been validated via numerical simulations of fluid flow through a porous soil sample and grain batch sedimentations. The simulations of granular flows in the submerged environment have led to some meaningful insights into the flow mechanisms, such as the mobilization of sediments, the generation and dissipation of excess pore water pressures and the evolution of effective stresses. Overall, this study shows that the proposed numerical tools are capable of modelling the mechanical and hydraulic behaviour of terrestrial and submarine landslides.

Orientador: Milton Mori<br>Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Quimica<br>Made available in DSpace on 2018-08-14T18:54:23Z (GMT). No. of bitstreams: 1 Bastos_JaciCarloSchrammCamara_D.pdf: 4353350 bytes, checksum: 0de3633c8edcc4ae2ec7fc57a3b590d0 (MD5) Previous issue date: 2009<br>Resumo: Apresenta-se nesta pesquisa uma análise experimental e numérica do comportamento da fase dispersa em um jato circular bifásico confinado e uma comparação com jato circular bifásico livre. Nas análises experimentais, uma câmara pentagonal em acrílico foi utilizada como sistema de confinamento para a obtenção de perfis axiais e radiais de velocidade média, flutuação de velocidade (RMS) e intensidade de turbulência. Estes dados foram analisados a fim de desenvolver uma análise completa da região desenvolvida do jato. Três diferentes jatos foram utilizados para a alimentação da fase gás no topo da câmara, mas apenas o jato central foi carregado com partículas entre 60 e 90µm de diâmetro. Os outros dois foram utilizados para proporcionar uma maior interação entre as fases no interior da câmara. A técnica óptica Phase Doppler Anemometry (PDA), foi empregada na medição da velocidade instantânea da fase sólida e do diâmetro das partículas nas diferentes posições axiais a partir do bico do jato. Nove casos de estudo distintos são investigados individualmente e, em seguida, comparados entre si. Estes casos fornecem informações importantes sobre o comportamento e o efeito do confinamento dos jatos sobre o transporte macrocóspico e turbulento das partículas entre o centro e as regiões de contorno do jato. As análises numéricas tratam da modelagem matemática tridimensional, turbulenta e transiente do escoamento no jato bifásico confinado. O modelo trata as fases gás e sólida a partir de uma abordagem Euleriana. O fechamento das equações de transporte foi realizado utilizando o modelo de turbulência de duas equações k-e para a fase gás e modelos de turbulência de zero-equação para a fase sólida, e ainda em alguns casos esta última apenas sofreu efeitos turbulentos advindos da fase contínua. A acurácia das previsões do modelo em um jato de partículas confinadas com as características médias no tempo, assim como os coeficientes da correlação de turbulência foram avaliados. Perfis radiais de velocidade média e fração volumétrica das partículas foram capturados em quarenta e dois níveis, subdivididos em nove casos e comparados aos dados experimentais adquiridos. O diâmetro médio das partículas utilizado nas simulações foi de 75µm e as velocidades iniciais utilizadas variam entre 3 e 11m/s no jato central. O modelo matemático previu um escoamento desenvolvido semelhante ao que foi encontrado experimentalmente.<br>Abstract: It is presented in this research an experimental and numerical analisys of the dispersed phase behavior in a circular confined two-phase jet and a comparison with circular free two-phase jet. In the experimental analysis, a pentagonal plexiglass chamber was used as confined system for the axial and radial profiles investigation of mean velocity, fluctuation velocity known as RMS velocity and turbulence intensity. These data were analyzed in order to develop a complete analysis in the developed region of the jet. Three different nozzles were used to feed the gas phase at the top of the chamber, but just the central nozzle was loaded with particles between 60 and 90µm of diameter. The other two were used to increase the interaction between the phases in the chamber. An optical technique known as Phase Doppler Anemometry was used to measure the instantaneous velocity of the solid phase and particle diameter in different axial positions of the jet nozzle. Nine different cases of study are investigated individually and then compared among each other. These cases provide important information about the jets behavior and the confinement effect on the macrocospic and turbulent transport of particles between the jet center and the jet contour regions. The numerical analysis deals with three-dimensional, turbulent and transient mathematical modeling of a confined two-phase jet flow. The model treats the gas and the solid phases from an Eulerian approach. The closure of the transport equations have been accomplished by using the k-e turbulence model for the gas phase and the zero-equation turbulence model for the solid phase, and in some cases the latter suffered turbulent effects occuring only from the continuos phase. The accuracy of the model predictions in a particle-laden confined jet with the characteristics as well as turbulence correlation coefficients have been evaluated. Radial mean velocity profiles for the solid phase were computed on forty two axial levels, subdivided in nine cases and compared to the obtained experimental data. The mean particle diameter used in the simulations was 75µm and the initial velocities used vary between 3 and 11m/s. The mathematical model predicted a flow development similar to that found experimentally.<br>Doutorado<br>Desenvolvimento de Processos Químicos<br>Mestre em Engenharia Química