The works aimed at the further development and validation of models for CFD codes. For this reason, the new thermal-hydraulic test facility TOPFLOW was erected and equipped with wire-mesh sensors with high spatial and time resolution. Vertical test sections with nominal diameters of DN50 and DN200 operating with air-water as well as steam-water two-phase flows provided results on the evaluation of flow patterns, on the be¬haviour of the interfacial area as well as on interfacial momentum and heat transfer. The validation of the CFD-code for complex geometries was carried out using 3D void fraction and velocity distributions obtained in an experiment with an asymmetric obstacle in the large DN200 test section. With respect to free surface flows, stratified co- and counter-current flows as well as slug flows were studied in two horizontal test channels made from acrylic glass. Post-test calculations of these experiments succeeded in predicting the slug formation process. Corresponding to the main goal of the project, the experimental data was used for the model development. For vertical flows, the emphasis was put on lateral bubble forces (e.g. lift force). Different constitutive laws were tested using a Multi Bubble Size Class Test Solver that has been developed for this purpose. Basing on the results a generalized inhomogeneous Multiple Size Group (MUSIG) Model has been proposed and implemented into the CFD code CFX (ANSYS). Validation calculations with the new code resulted in the conclusion that particularly the models for bubble coalescence and fragmentation need further optimisation. Studies of single effects, like the assessment of turbulent dissipation in a bubbly flow and the analysis of trajectories of single bubbles near the wall, supplied other important results of the project.

Ziel der Arbeiten war die Weiterentwicklung und Validierung von Modellen in CFD-Codes. Hierzu wurde am FZD die thermohydraulische Versuchsanlage TOPFLOW aufgebaut und mit räumlich und zeitlich hochauflösenden Gittersensoren ausgestattet. Vertikale Teststrecken mit Nenndurchmessern von DN50 bzw. DN200 für Luft/Wasser- sowie Dampf/Wasser-Strömungen lieferten Ergebnisse zur Entwicklung von Strömungsformen, zum Verhalten der Zwischenphasengrenzfläche sowie zum Wärme- und Impulsaustausch zwischen den Phasen. Die Validierung des CFD-Codes in komplexen Geometrien erfolgte anhand von 3D Gasgehalts- und Geschwindigkeitsfeldern, die bei Umströmung eines asymmetrischen Hindernisses auftreten, das in der Teststrecke DN200 eingebaut war. Im Hinblick auf Strömungen mit freier Oberfläche untersuchte das FZD in zwei horizontalen Acrylglas-Kanälen geschichtete Zweiphasenströmungen im Gleich- bzw. Gegenstrom sowie Schwallströmungen. Bei den Nachrechnungen dieser Versuche gelang die Simulation der Schwallentstehung. Entsprechend des Projektziels wurden die experimentellen Ergebnisse zur Modellentwicklung genutzt. Bei vertikalen Strömungen stand die Wirkung der lateralen Blasenkräfte (z.B. Liftkraft) im Vordergrund. Zum Test unterschiedlicher Modellansätze wurde hierzu ein Mehrblasenklassen-Testsolver entwickelt und genutzt. Darauf aufbauend wurde ein neues Konzept für ein Mehrblasenklassenmodell, das Inhomogene MUSIG Modell erarbeitet und in den kommerziellen CFD Code CFX (ANSYS) implementiert. Bei Validierungsrechnungen zeigte sich, dass vor allem die Blasenkoaleszenz- und -zerfallsmodelle weiter optimiert werden müssen. Untersuchungen zu Einzeleffekten, wie z.B. die Abschätzung von Turbulenzkoeffizienten und die Analyse der Trajektoren von Einzelblasen in unmittelbarer Wandnähe, lieferten weitere wichtige Ergebnisse des Projekts.

Die in diesem Bericht vorgestellten Experimente betreffen die Kondensation von Dampf in horizontalen bzw. leicht geneigten Rohren bei hohen auftretenden Temperaturdifferenzen bis zu über 200 K. Weitere Besonderheiten sind die detaillierte Untersuchung des transienten Verhaltens eines nichtkondensierbaren Gases mit einer neuartigen Messtechnik sowie die Ermittlung des Einflusses des Gases auf den Kondensationsvorgang. Beim Experiment wurden schnelle Übergangsvorgänge ausgelöst, indem ein in einer Kühlwanne liegendes, leicht geneigtes Wärmeübertragerrohr plötzlich mit dem Dampfraum eines unter Druck stehenden Kessels verbunden wurde. Dabei wurden im Rohr unterschiedliche Anfangsbedingungen hinsichtlich der Vorlage von nichtkondensierbaren Gasen (in diesem Falle Luft) eingestellt. Es wurden Versuche mit Atmosphärendruck, mit erhöhtem Druck, aber auch mit vorheriger Evakuierung des Versuchsrohrs durchgeführt. Durch eine Instrumentierung mit neuartigen Nadelsonden, die eine Phasendetektion kombiniert mit einer schnellen lokalen Temperaturmessung ermöglichen, konnte die Umverteilung von Dampf, Kondensat und nichtkondensierbarem Gas als Funktion der Zeit beobachtet werden. Damit bieten die erhaltenen Daten die Möglichkeit, insbesondere die in den Thermohydraulikprogrammen vorhandenen Optionen zur Berechnung der Ausbreitung von nichtkondensierbaren Gasen unter transienten Bedingungen zu validieren.

In order to investigate the two-phase flow behaviour in a complex reactor-typical geometry and to supply suitable data for CFD code validation, a model of the hot leg of a pressurised water reactor was built at FZD. The hot leg model is operated in the pressure chamber of the TOPFLOW test facility, which is used to perform high-pressure experiments under pressure equilibrium with the inside atmosphere of the chamber. This technique makes it possible to visualise the two-phase flow through large windows, also at reactor-typical pressure levels. In order to optimise the optical observation possibilities, the test section was designed with a rectangular cross-section. Experiments were performed with air and water at 1.5 and 3.0 bar at room temperature as well as with steam and water at 15, 30 and 50 bar and the corresponding saturation temperature (i.e. up to 264°C). The total of 194 runs are divided into 4 types of experiments covering stationary co-current flow, counter-current flow, flow without water circulation and transient counter-current flow limitation (CCFL) experiments. This report provides a detailed documentation of the experiments including information on the experimental setup, experimental procedure, test matrix and on the calibration of the measuring devices. The available data is described and data sheets were arranged for each experiment in order to give an overview of the most important parameters. For the cocurrent flow experiments, water level histograms were arranged and used to characterise the flow in the hot leg. In fact, the form of the probability distribution was found to be sensitive to the boundary conditions and, therefore, is useful for the CFD comparison. Furthermore, the flooding characteristics of the hot leg model plotted in terms of the classical Wallis parameter or Kutateladze number were found to fail to properly correlate the data of the air/water and steam/water series. Therefore, a modified Wallis parameter is proposed, which takes the effect of viscosity into account.

<p>Consider a mixture of two immiscible, incompressible fluids e.g. oil and water. Since the fluids do not mix, an interface between the two fluids will form and move in time. The motion of the two fluids can be modelled by the incompressible Navier-Stokes equations for two phase flow with surface tension together with a representation of the moving interface. The parameters in the Navier-Stokes equations will depend on the position and other properties of the interface. The interface should move with the velocity of the flow at the interface. Since the fluids are incompressible, the density of each fluid is constant. Mass conservation then implies that the volume occupied by each of the two fluids should not change with time. The object of this thesis has been to develop a new numerical method to simulate incompressible two phase flow accurately that conserves mass and volume of each fluid correctly.</p><p>Numerical simulations of incompressible two phase flow with surface tension have been a challenge for many years. Several methods have been developed and used prior to the work presented in this thesis. The two most commonly used methods are volume of fluid methods and level set methods. There are advantages and disadvantages of both of the methods.</p><p>In volume of fluid methods the interface is represented by a discontinuity of a globally defined function. Because of the discontinuity it is hard both to move the interface as well as to calculate properties of the interface such as curvature. Specially designed methods have to be used, and all these methods are low order accurate. Volume of fluid methods do however conserve the volumes of the two fluids correctly.</p><p>In level set methods the interface is represented by the zero contour of the globally defined signed distance function. This function is smooth across the interface. Since the function is smooth, standard methods for partial differential equations can be used to advect the interface accurately. A reinitialization is however needed to make sure that the level set function remains a signed distance function. During this process the zero contour might move slightly. Because of this, the volume conservation of the method becomes poor.</p><p>In this thesis we present a new level set method. The method is designed such that the volume of each fluid is conserved, at least approximately. The interface is represented by the 0.5 contour of a regularized characteristic function. As for standard level set methods, the interface is moved first by an advective step, and then reinitialized. Unlike traditional level set methods, we can formulate the reinitialization as a conservation law. Conservative methods can then be used to move and to reinitialize the level set function numerically. Since the level set function is a regularized characteristic function, we can expect good conservation of the volume bounded by the interface.</p><p>The method is discretized using both finite differences and finite elements. Uniform and adaptive grids are used in both two and three space dimensions. Good convergence as well as volume conservation is observed. Theoretical studies are performed to investigate the conservation and the computational time needed for reinitialization.</p>

When a lander approaches a dusty surface, the plume from the descent engine impinges on the ground and entrains loose regolith into a high velocity spray. This problem exhibits a wide variety of complex phenomena such as highly under-expanded plume impingement, transition from continuum to free molecular flow, erosion, coupled gas-dust motions, and granular collisions for a polydisperse distribution of aerosolized particles. The focus of this work is to identify and model the important physical phenomena and to characterize the dust motion that would result during typical lunar landings. A hybrid continuum-kinetic solver is used, but most of the complex physics are simulated using the direct simulation Monte Carlo method. A descent engine of comparable size and thrust to the Lunar Module Descent Engine is simulated because it allows for direct comparison to Apollo observations. Steady axisymmetric impingement was first studied for different thrust engines and different hovering altitudes. The erosion profiles are obtained from empirically derived scaling relationships and calibrated to closely match the net erosion observed during the Apollo missions. Once entrained, the dust motion is strongly influenced by particle-particle collisions and the collision elasticity. The effects of two-way coupling between the dust and gas motions are also studied. Small particles less than 1 µm in diameter are accelerated to speeds that exceed 1000 m/s. The larger particles have more inertia and are accelerated to slower speeds, approximately 350 m/s for 11 µm grains, but all particle sizes tend obtain their maximum speed within approximately 40 m from the lander. The maximum particle speeds and erosion rates tend to increase as the lander approaches the lunar surface. The erosion rates scale linearly with engine thrust and the maximum particle speed increases for higher thrust engines. Dust particles are able to travel very far from the lander because there is no background atmosphere on the moon to inhibit their motion. The far field deposition is obtained by using a staged calculation, where the first stages are in the near field where the flow is quasi-steady and the outer stages are unsteady. A realistic landing trajectory is approximated by a set of discrete hovering altitudes which range from 20 m to 3 m. Larger particles are accelerated to slower speeds and are deposited closer to the lander than smaller particles. Many of the gas molecules exceed lunar escape speed, but some gas molecules become trapped within the dust cloud and remain on the moon. The high velocity particulate sprays can be damaging to nearby structures, such as a lunar outpost. One way of mitigating this damage is to use a berm or fence to shield nearby structures from the dust spray. This work attempts to predict the effectiveness of such a fence. The effects of fence height, placement, and angle as well as the model sensitivity to the fence restitution coefficient are discussed. The expected forces exerted on fences placed at various locations are computed. The pressure forces were found to be relatively small at fences placed at practical distances from the landing site. The trajectories of particles that narrowly avoid the fence were not significantly altered by the fence, suggesting that the dust motion is weakly coupled to the gas in the near vicinity of the fence. Future landers may use multi-engine configurations that can form 3-dimensional gas and dust flows. There are multiple plume-plume and plume-surface interactions that affect the erosion rates and directionality of the dust sprays. A 4-engine configuration is simulated in this work for different hovering altitudes. The focusing of dust along certain trajectories depends on the lander hovering altitude, where at lower altitudes the dust particles focus along symmetry planes while at higher altitudes the sprays are more uniform. The surface erosion and trenching behavior for a 4-engine lander are also discussed. / text

The diffusion absorption refrigeration (DAR) cycle is a promising technology for fully thermally driven cooling. It is well suited to applications in medicine refrigeration and air-conditioning in off-grid settings. However, design and engineering knowhow for the technology is limited; therefore, system development has historically been an iterative and expensive process. Additionally, conventional system designs require high-grade energy input for operation, and are unsuitable for low-temperature solar- or waste-heat activated applications. In the present effort, component- and system-level DAR engineering analyses are performed. Detailed bubble-pump generator (BPG) component models are developed, and are validated experimentally and with direct simulations. Investigations into the BPG focus on the Taylor flow pattern in the intermediate Bond number regime, which has not yet been thoroughly characterized in the literature, and has numerous industry applications, including nuclear fuel processing and well dewatering. A coupling-fluid heated BPG design is also investigated experimentally for low-source-temperature operation. Phase-change simulation methodologies are developed to rigorously study the continuously developing flow pattern in this BPG configuration. Detailed component-level models are also formulated for all of the other DAR heat and mass exchangers, and are integrated to yield a complete system-level model. Results from these modeling studies are applied to develop a novel fully passive low-source-temperature (110 - 130°C) DAR system that delivers refrigeration grade cooling. This design achieves operation at target conditions through the use of alternate working fluids (NH3-NaSCN-He), the coupling-fluid heated BPG, and a novel absorber configuration. The complete DAR system is demonstrated experimentally, and evaluated over a range of operating conditions. Experimental results are applied to assess and refine component- and system- level models.

The aim of the work is to investigate heat transfer of the staggered tube bank to downward statically stable foam flow.

As a result of the problem of hydrate in multiphase pipelines in offshore production is becoming more and more severe with the increasing of the water depth, the study on oil-gas-water-hydrate has became a hot point of multiphase flow. In this paper, the hydrate particle and liquid phase was treated as pseudo-fluid, the steady hydraulic, thermodynamical and phase equilibrium calculation method of gas-hydrate slurry was developed. Comparison was carried out between calculated data and experimental data from flow loop in our laboratory. With strict flash calculation the following items were determined: the amount of hydrate; phase number; the location that hydrate appeared; flowrate and molar component of gas phase and liquid phase. Then thermodynamic quantities were carried out with proper relational expression. When Compositional model is used to simulate two phase flow, it is required to couple mass, momentum, energy equation and equation of state. In the other word, the parameters in these four equations are interacted. However they are all the functions of p, T and z. In steady condition, it’s assumed that the composition of fluid is unchangeable along the pipeline and the flow can be described by pressure and temperature. In this paper, calculation method of gas-liquid two phase flow which respectively was improved. Liquid holdup and pressure drop were calculated by momentum equation. Enthalpy balance equation was substituted by explicit formulation of temperature calculation which meant that the loop of temperature was not required.

Bubble nucleation at subcooled flow boiling condition in a horizontal annular channel with a square cross section by the use of high-speed camera is investigated. The channel represents a scaled-down version of a single rod of CANDU reactor core. The experiments were performed by the use of water at pressures between 1-3 atm, constant heat flux of 0.124 MW/m2, liquid bulk subcooling of 32-1oC and mean flow velocities of 0.3-0.4 m/s. Bubble lift-off diameters were obtained from direct high speed videography. The developed model for the bubble lift-off diameter was obtained by analyzing the forces acting on a bubble. Furthermore, a model for the bubble growth rate constant was suggested. The proposed model was then compared to experimental data and it has shown a good agreement with the experimental data. Additionally, the effects of liquid bulk subcooling, liquid pressure and mean flow velocity on bubble lift-off diameter were investigated.