Die in dieser Arbeit vorgestellten Simulationen aufsteigender fluider Partikel wurden mit dem CFD-Programm FS3D durchgeführt, welches auf der Volume-of-Fluid (VoF) Methode basiert. Die Validierung des Codes erfolgt durch Vergleich der numerischen Lösungen für schleichende Strömungen mit analytischen Lösungen, wobei eine gute Übereinstimmung festgestellt wird. Im ersten Teil der Dissertation werden Simulationen für den freien Aufstieg von Öltropfen in Wasser mit experimentellen Beobachtungen hinsichtlich der Aufstiegsgeschwindigkeit, der Tropfenform und der Bewegungsbahn verglichen. Die Aufstiegsgeschwindigkeiten und Widerstandsbeiwerte sind vergleichbar, die simulierten Tropfen sind jedoch deutlich flacher. Dieser Unterschied kann durch Verunreinigungen der Grenzfläche im Experiment verursacht sein. Der Übergang von einem gradlinigen Aufstieg zu zickzack-förmigen Aufstiegsbahnen kann mit Hilfe der Simulationen auf Instabilitäten im Nachlauf der Blasen zurückgeführt werden, die zu einer periodischen Wirbelablösung führen. Im zweiten Teil der Dissertation wird der Aufstieg von Blasen in linearen Scherströmungen untersucht. Steigen die Blasen in einer vertikalen Scherströmung auf, so beobachtet man eine seitliche Migration. Diese seitliche Migration der Blasen wird durch die sogenannte Liftkraft verursacht, deren Vorzeichen und Betrag von der Blasengröße und den Stoffeigenschaften der Flüssigkeit abhängt. Die Simulationen zeigen, daß das Vorzeichen der Liftkraft für eher sphärische Blasen durch den Bernoulli-Effekt erklärt werden kann. An stark deformierten Blasen hingegen wirkt die Liftkraft in umgekehrter Richtung. Dieses Phänomen tritt auch in den Simulationen auf. Verschiedene Hypothesen für die Ursache dieses Phänomens werden überprüft. Die bekannteste experimentelle Korrelation für die Liftkraft von Tomiyama u.a. (2002) wird durch Simulation von realen Flüssigkeiten mit bekannten Stoffeigenschaften wie auch von Modellfluiden mit willkürlichen Stoffeigenschaften validiert und weitgehend bestätigt. Die Lift-Korrelation hat demnach hinsichtlich der Stoffeigenschaften der Flüssigkeit einen größeren Geltungsbereich, als bisher experimentell überprüft wurde. The simulations presented in this thesis were performed with the CFD code FS3D which is based on the Volume of Fluid method. The code is validated using analytical solutions for creeping flows and a good agreement is observed between simulation and analytical solution. In the first part of the thesis, the free rise of oil drops in water is simulated and compared with experimental observations. The results show that the rising velocities and the drag coefficients are similar in both cases, but the simulated drops are flatter (more oblate). This difference may be caused by impurities of the particle surface (surfactants) in the experiments. The simulations show that the transition from rectilinear to periodic trajectories is caused by instabilities in the wake, which lead to a periodic vortex shedding. In the second part of the thesis, the rise of bubbles in linear shear flows is investigated. If bubbles rise in a vertical shear flow, a lateral migration can be observed. This migration is caused by the so called lift force. Sign and magnitude of the lift force depend on the size of the bubble and the material properties of the liquid. The simulation results show that the sign of the lift force on spherical bubbles can be explained by the Bernoulli effect. However, the lift force on more distorted bubbles acts in the opposite direction. This phenomenon can also be observed in the simulation. In this work several hypotheses for the reason of this phenomenon are checked. Furthermore, most common correlation for the lift force (developed by Tomiyama et al. in 2002) is validated for fluids of known material and model fluids with arbitrary material data. The correlation is valid in a wider range of fluid material properties than proved experimentally up to now.

The problem of real-time simulation of ocean surface waves, ship movement and the coupling in between is tackled, and a number of different methods are covered and discussed. Among these methods, the finite volume method has been implemented in an attempt to solve the problem, along with the compressible Euler equations, an octree based staggered grid which allows for easy adaptive mesh refinement, the volume of fluid method and a variant of the Hyper-C advection scheme for compressible flows for advection of the phase fraction field. The process of implementing the methods that were chosen proved to be tricky in many ways, as they involve a large number of advanced topics, and the implementation that was implemented in this thesis work suffered from numerous issues. There were for example problems with keeping the interface intact, as well as a harsh restriction on the time step size due to the CFL condition. Improvements required to make the method sustainable for real-time applications are discussed, and a few suggestions on alternative approaches that are already in use for similar purposes are also given and discussed. Furthermore, a method for compensating for gain/loss of mass when solving the incompressible flow equations with an inaccurately solved pressure Poisson equation is presented and discussed. A momentum conservative method for transporting the velocity field on staggered grids without introducing unnecessary smearing is also presented and implemented. A simple, physically based illumination model for sea surfaces is derived, discussed and compared to the Blinn–Phong shading model, although it is never implemented. Finally, a two-dimensional partial differential equation in the spatial domain for simulating water surface waves for mildly varying bottom topography is derived and discussed, although it is deemed to be too slow for real-time purposes and is therefore never implemented. / <p>This publication differs from the printed version of the report in the sense that links are blue in this version and black in the printed version.</p>

The development of a fuel filler pipe is based solely on experience and physical experiment. The challenge lies in designing the pipe to fulfill the customer needs. In other words designing the pipe such as the fuel flow does not splash back on the fuel dispenser causing a premature shut off. To improve this “trial-and-error” based development a computational fluid dynamics (CFD) model of the refueling process is investigated. In this thesis a CFD model has been developed that can predict the fuel flow in the filler pipe. Worst case scenario of the refueling process is during the first second when the tank is partially filled. The most critical fluid is diesel due to the commercially high volume flow of 55 l/min. Due to limitations of computational resources the simulations are focused on the first second of the refueling process. The challenge in this project is creating a CFD model that is time efficient, thus require the least amount of computational resources necessary to provide useful information. A multiphase model is required to simulate the refueling process. In this project the implicit volume of fluid (VOF) has been used which has previously proven to be a suitable choice for refueling simulations. The project is divided into two parts. Part one starts with experiments and simulations of a simplified fuel system with water as acting liquid with a Reynolds number of 90 000. A short comparison between three different turbulence models has been investigated (LES, DES and URANS) where the most promising turbulence model is URANS, specifically the SST k-ω model. A sensitivity analysis was performed on the chosen turbulence model. Between the chosen mesh and the densest mesh the difference of streamwise velocity in the boundary layer was 2.6 %. The chosen mesh with 1.9 M cells and a time step of 1e-4 s was found to be the best correlating model with respect to the experiments. In part two a real fuel filling system was investigated both with experiments and simulations with the same computational model as the chosen one from part one. The change of fluid and geometry resulted in a lower Reynolds number of 12 000. Two different versions of the fuel system was investigated; with a bypass pipe and without a bypass pipe. Because of a larger volumetric region the resulting mesh had 3.7 M cells. The finished model takes about 230 h on a local workstation with 11 cores. On a cluster with 200 cores the same simulation takes 30 h. The resulting model suffered from interpolation errors at the inlet which resulted in a volume flow of 50 l/min as opposed to 55 l/min in the experiments. Despite the difference the model could capture the key flow characteristics. With the developed model a new filler pipe can be easily implemented and provide results in shorter time than a prototype filler pipe can be ordered. This will increase the chances of ordering one single prototype that fulfills all requirements. While the simulation model cannot completely replace verification by experiments it provides information that transforms the development of the filler pipe to knowledge based development.

abstract: The Volume-of-Fluid method is a popular method for interface tracking in Multiphase applications within Computational Fluid Dynamics. To date there exists several algorithms for reconstruction of a geometric interface surface. Of these are the Finite Difference algorithm, Least Squares Volume-of-Fluid Interface Reconstruction Algorithm, LVIRA, and the Efficient Least Squares Volume-of-Fluid Interface Reconstruction Algorithm, ELVIRA. Along with these geometric interface reconstruction algorithms, there exist several volume-of-fluid transportation algorithms. This paper will discuss two operator-splitting advection algorithms and an unsplit advection algorithm. Using these three interface reconstruction algorithms, and three advection algorithms, a comparison will be drawn to see how different combinations of these algorithms perform with respect to accuracy as well as computational expense. / Dissertation/Thesis / Masters Thesis Mechanical Engineering 2015

Le travail de cette thèse concerne l’étude des écoulements diphasiques dans les mini-canaux d’une pile à combustible. Ces canaux, de taille millimétrique, ont un double rôle d’alimentation de la pile en combustibles gazeux et d’évacuation de l’eau produite. Il a été montré expérimentalement que la configuration d’écoulement a un impact direct sur les performances de la pile. L’objectif est donc de progresser dans la compréhension des mécanismes physiques présents dans ces canaux. Pour cela, l’outil de simulation numérique JADIM est développé pour modéliser la physique de la ligne triple afin de permettre la description des écoulements diphasiques et les transitions entre configurations, où la mouillabilité joue un rôle déterminant. Deux phénomènes physiques sont essentiellement étudiés : la migration dans le canal de gouttelettes formées par l’introduction d’eau dans le canal, et l’évolution et la stabilité de la répartition spatiale des phases lors du remplissage progressif des canaux.

The aim of this master’s thesis is to determine the torque momentum of the jet deflector of the Pelton turbine using a two-phased CFD simulation. The reason for determining the moments is the inaccurate formula according to a standard, whose results are assumed to be too excessive. The first part is devoted to theory of Pelton turbine and its main components. The second part contains the design of deflector, explanation of the used k- turbulence method and the Volume of Fluid two-phased flow modeling method. The results of the simulation are presented at the end of the thesis.

Die in dieser Arbeit vorgestellten Simulationen aufsteigender fluider Partikel wurden mit dem CFD-Programm FS3D durchgeführt, welches auf der Volume-of-Fluid (VoF) Methode basiert. Die Validierung des Codes erfolgt durch Vergleich der numerischen Lösungen für schleichende Strömungen mit analytischen Lösungen, wobei eine gute Übereinstimmung festgestellt wird. Im ersten Teil der Dissertation werden Simulationen für den freien Aufstieg von Öltropfen in Wasser mit experimentellen Beobachtungen hinsichtlich der Aufstiegsgeschwindigkeit, der Tropfenform und der Bewegungsbahn verglichen. Die Aufstiegsgeschwindigkeiten und Widerstandsbeiwerte sind vergleichbar, die simulierten Tropfen sind jedoch deutlich flacher. Dieser Unterschied kann durch Verunreinigungen der Grenzfläche im Experiment verursacht sein. Der Übergang von einem gradlinigen Aufstieg zu zickzack-förmigen Aufstiegsbahnen kann mit Hilfe der Simulationen auf Instabilitäten im Nachlauf der Blasen zurückgeführt werden, die zu einer periodischen Wirbelablösung führen. Im zweiten Teil der Dissertation wird der Aufstieg von Blasen in linearen Scherströmungen untersucht. Steigen die Blasen in einer vertikalen Scherströmung auf, so beobachtet man eine seitliche Migration. Diese seitliche Migration der Blasen wird durch die sogenannte Liftkraft verursacht, deren Vorzeichen und Betrag von der Blasengröße und den Stoffeigenschaften der Flüssigkeit abhängt. Die Simulationen zeigen, daß das Vorzeichen der Liftkraft für eher sphärische Blasen durch den Bernoulli-Effekt erklärt werden kann. An stark deformierten Blasen hingegen wirkt die Liftkraft in umgekehrter Richtung. Dieses Phänomen tritt auch in den Simulationen auf. Verschiedene Hypothesen für die Ursache dieses Phänomens werden überprüft. Die bekannteste experimentelle Korrelation für die Liftkraft von Tomiyama u.a. (2002) wird durch Simulation von realen Flüssigkeiten mit bekannten Stoffeigenschaften wie auch von Modellfluiden mit willkürlichen Stoffeigenschaften validiert und weitgehend bestätigt. Die Lift-Korrelation hat demnach hinsichtlich der Stoffeigenschaften der Flüssigkeit einen größeren Geltungsbereich, als bisher experimentell überprüft wurde. The simulations presented in this thesis were performed with the CFD code FS3D which is based on the Volume of Fluid method. The code is validated using analytical solutions for creeping flows and a good agreement is observed between simulation and analytical solution. In the first part of the thesis, the free rise of oil drops in water is simulated and compared with experimental observations. The results show that the rising velocities and the drag coefficients are similar in both cases, but the simulated drops are flatter (more oblate). This difference may be caused by impurities of the particle surface (surfactants) in the experiments. The simulations show that the transition from rectilinear to periodic trajectories is caused by instabilities in the wake, which lead to a periodic vortex shedding. In the second part of the thesis, the rise of bubbles in linear shear flows is investigated. If bubbles rise in a vertical shear flow, a lateral migration can be observed. This migration is caused by the so called lift force. Sign and magnitude of the lift force depend on the size of the bubble and the material properties of the liquid. The simulation results show that the sign of the lift force on spherical bubbles can be explained by the Bernoulli effect. However, the lift force on more distorted bubbles acts in the opposite direction. This phenomenon can also be observed in the simulation. In this work several hypotheses for the reason of this phenomenon are checked. Furthermore, most common correlation for the lift force (developed by Tomiyama et al. in 2002) is validated for fluids of known material and model fluids with arbitrary material data. The correlation is valid in a wider range of fluid material properties than proved experimentally up to now.

Dam spillways and downstream areas are used to guide large flows of water during for example heavy rainfall. The large flows give way to turbulent pattern sand velocities that may damage the river banks or the dam structure. Investigation of these water patterns at certain flows are therefore done to examine at risk areas. In this study CFD simulations were performed for different flows with different boundary conditions for varying surface roughness level. Results were then compared to a previous model study from 2009. The ANSYS ecosystem was used in production of the 3D model, construction of mesh and running of simulations.The flow for the maximum discharge capacity of the sluices was simulated as well as the design flow which is the highest flow the dam is supposed to be able to withstand. In this report the flow has been modeled using RANS with the SST kω-model in a VOF transient setup. Results showed that for both the design flow and the maximum discharge capacity flow the energy conversion is functioning poorly and that a considerable backward circulation exists on the right riverside. This behavior could possibly injure the right dam structure by moving debris upwards against the stream.

No description available.

Wading is the phenomenon where a vehicle drives through water with a relatively deep water level. Sincea large portion of the vehicle is submerged in water it can affect the driveability and function of individualcomponents. Wading is therefore an important phenomenon to be aware of especially today where society moves towards alternative energy sources. This includes water sensitive components when contact with water can generate major consequences. Previous knowledge and experience of wading has been from performing physical tests, but using Computational Fluid Dynamics (CFD) to examine the phenomenon can accelerate the iterative design process. In this thesis, numerical method of wading simulations on complete heavy-duty vehicles using the software STAR-CCM+ are developed. Furthermore, the results from the numerical methods are validated against results from physical tests performed at Scania’s test facility in Södertälje. The numerical methods are divided into a simplified model of a Battery Electric Vehicle (BEV) and a detailed geometry of a gas-driven vehicle from Scania. Beside dividing the wading scenario into the geometries, two different methods are developed, Wave and Wading. The Wave-method includes the vehicle standing still while a water wave is fed in through the inlet of the domain, i.e. allowed to flush over the vehicle, with a velocity of 3.6 km/h and 8 km/h. This method is implemented for both a generic simplified BEV truck and a detailed real-life Scania truck. For the Wading-method, motion is applied to the vehicle where itis driving with a velocity of 3.6 km/h through a digital twin of the water trench available at the test facility. This method is further divided into two cases, Zero Gap and Floating, where the difference is the distance between the tires of the vehicle and ground of the domain. The Floating-case includes a 10 cm distance and the Zero Gap-case has no gap between the tires and ground. The Wading-method is only implemented for the simplified geometry due to the computational cost and complexity. All methods use the Volume of Fluid (VOF) method for multiphase modelling and the Zero Gap-case uses Overset Mesh for modelling motion. The validation of the simulations focuses on the water behaviour such as water surface topology and water flowing inside the vehicle while wading. The results for the Wave-method with both the simplified and detailed truck at 8 km/h shows similarities in the water surface topology between the numerical model and the physical test. The simulations of the Wading-method is not visualising any similarities since the visible wave pattern are few and unclear in the numerical model. An isosurface is used to visualise the surface of the water which generated a smooth topology since no other options, such as vector fields, are added. It is found that the water movement inside the vehicle will affect water sensitive areas, e.g. on the battery packs. It is concluded that the derived methods are a first draft and should be directed towards future development in optimising the methods to lower the computational cost, but also to improve the capturing of the interface between the two phases. Due to instability and computational cost the detailed geometry is not implemented in the Wading-method. The methods are adapted to use different vehicle types since the simplified and detailed geometry are a BEV and a gas-driven truck respectively.