京都大学 / 0048 / 新制・課程博士 / 博士(理学) / 甲第18084号 / 理博第3962号 / 新制||理||1571(附属図書館) / 30942 / 京都大学大学院理学研究科地球惑星科学専攻 / (主査)准教授 齊藤 昭則, 教授 余田 成男, 准教授 藤 浩明 / 学位規則第4条第1項該当 / Doctor of Science / Kyoto University / DGAM

In urban areas, discharging wastewater into rivers is a common way to dispose of contaminants, and it is usually the most economical. Accurate information about how effluents are distributed in the receiving water body is desirable when designing industrial plants. Flow structures will be influenced by an effluent’s dilution processes during the mixing. Meanwhile, the cross-stream motions resulting from the streamline curvature can redistribute both the velocity and the shear stress, which favors the mixing behavior compared to a straight channel. However, the interactions between jet mixing behavior and the bend flow requires further investigation. In the present study, jets with different densities were discharged horizontally into a laboratory flume with a 135-degree open channel bend, and both the main and secondary flow behaviors in the bend were observed after the introduction of effluents. The acquired three- dimensional velocity data were used to validate numerical models of the effluent-bend flow. Numerical turbulence models such as the standard k-ε eddy viscosity model, non-linear k-ε model (Shih quadratic k-ε), and the k-ω SST (shear stress transport) model were employed to evaluate their accuracy. OpenFOAM was selected in the analysis for proposing better numerical models since it gives high-quality results to individualized complex fluid flows, and as an open source CFD software it can be beneficial to further develop and maintain. The first part of this study presents the implementation of the physical modelling of the proposed problems. Detailed descriptions of the experimental process were elaborated. Specifically, the three velocity components at four cross-sectional planes in the bend section were measured with and without saltwater jets by using the stereo Particle Image Velocimetry (PIV) ii technique in the laboratory flume. The experimental results show that the more pronounced effects with the jets were found at the beginning and exit of the bend. Although the jets had little effect on the maximum streamwise velocity, it was found that the occurrence of the negatively buoyant jets would affect the patterns and properties of the secondary flow in the bend. The second part of this study investigated the mechanisms underlying the two cells system, particularly when interacting with a discharged effluent jet. Detailed experimental data were used in interpreting the large center-region cell as well as small structures in the 135-degree open channel bend. A term-by-term analysis of the downstream vorticity equation was executed to investigate the various mechanisms underlying these cross-stream flow motions considering the influence of the negatively buoyant jets. The results indicated the generation and the dissipation of the streamwise vorticity with the effective terms of the vorticity equation. The third part of this study evaluated the performance of three different turbulence models with the experimental measurements. It can be concluded that fully 3D numerical models are capable of simulating the primary flow pattern in a strongly curved channel with the presence of a negatively buoyant jet. The comparison also shows that, although the outer bank cell was not predicted, the k-omega SST model can satisfactorily predict some of the smaller flow features in bend flow, such as the inner bank circulation cell and the overall form of the vorticity distributions. The results enable more reliable predictions for the characteristics and development of jets in a bend.

No / Vegetation remains to be an important factor that can hinder the river flow. It needs innovative management scheme, in order to adapt these changes and ensure sustainability of their multiple usages. Vegetation plays an important role in floods and droughts adaptation within river system to alleviate any flood that may propagates from river to its surrounding. Vegetation within river can also retard its flow to cause building-up of deposition, and further adding to uncertainty of water use under extreme droughts. Due to these, it is important to study and understand vegetation drag behaviour toward flow in order to prevent flood risk and water security with hydrological drought in the basin and any other negative impact caused by it. In this study, an analytical approach for river flooding has been studied by improved representation of drag coefficient CD in flow velocity distribution modelling. The analysis of flow parameters, i.e. Reynolds number, on the drag coefficient CD has been conducted. The presented model has been used and analysed in open channel flows with flexible vegetation. In modelling, the flexible vegetated channel layers were divided into vegetation, top of vegetation and water layer zones in the model. The balance of forces for each layer has been established by validation using different reported measured data. The modelling results showed reasonably corresponding prediction of velocity profile in flows with flexible vegetation.

No description available.

Hydropower is a more environmental friendly way of producing electric power than many other alternatives today. Though, the effects of constructing mega dams are much tangible for the local eco systems in addition to changing many people’s lives forever. In order to prevent floods, riverbank erosions or landslides, proper investigations of the environmental impact from dam constructions must be performed. One of the key parameters in such investigations is the flow discharge velocity. This master thesis treats experimental measurements and numerical simulations of the velocity downstream a model of Xiluodu dam. The Xiluodu dam is a mega dam under construction in China and will have a total capacity of 12 600 MW when completed. The model is in scale 1:100 and the experiments have been performed at Department of Hydraulic Engineering, Tsinghua University, Beijing, China. The velocity profile shows that the velocity in the middle of the river is larger than the velocity at the surface and near the riverbank. The comparison between the measured and the simulated velocities shows a difference of less than 20 percent in almost all points which can be considered as a good result. In those points where the difference is more than 20 percent, this is believed to be due to the position of these points. Some of them were located near a vortex and others very close to the bottom. This is a problem when sparsely measured topography in combination with linear interpolation makes the boundaries of the simulations incorrect. In order to perform better simulations, more densely topography data and better flow boundary conditions should be used. More measuring points of the velocity could also improve the result.

Abstract This study aims to measure the flow field distribution in a microchannel with different heights adjusted. Two different materials, PDMS and Coverglass, were used to observe the flow velocity distribution change resulting from the difference in Zeta potential. The velocity distribution data were also obtained. In the experiment, 1¡Ñ TBE buffer solution with viscosity of 1 cp was used with the electric field intensity controlled under 5, 7.5 and 10 kV/m, respectively. Micrometer resolution Particle Image Velocimetry (£gPIV) was used to measure partial velocity distribution in order to explore the hydrodynamic stretch effect on DNA molecules when the microchannel, where the solution was placed, was adjusted to different heights. This study also statistically analyzed the stretch length distribution of DNA molecules in the microchannel and calculated the time of DNA molecule deformation and stress relaxation time in order to understand the stretch condition under different heights as well as the stretch and deformation of DNA molecules in microchannels.

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.

The purpose of this study was to investigate the velocity distribution through an annular bed packed randomly with equal sized spheres. Extensive research has been conducted on the velocity distribution inside packed beds packed with equal sized spheres, different sized spheres, deformed spheres, cylinders and Raschig-rings. A majority of these experimental and numerical studies focused on the cylindrical packed bed. These studies and numerical models are all confined to the velocity profile once the fluid flow is fully developed. The development of the velocity through the inlet region of the bed and the fluid flow redistribution in the outlet of the bed is thus neglected. The experimental investigation into the velocity distribution down stream of the annular packed bed of the HTTU indicated that the velocity profile was independent of the mass flow rate for a particle Reynolds number range of 439 £ Re £ 3453 . These velocity profiles did not represent the distribution of the axial velocity due to shortcomings associated with the single sensor hot wire anemometry system used to measure the velocity distribution. A numerical investigation, using the RANS CFD code STAR-CCM+®, into the velocity distribution downstream of an explicitly modelled bed of spheres indicated that the axial velocity distribution could be extracted from the experimental velocity profiles by using an adjustment factor of 0.801. This adjusted velocity profile was used in the verification of the implicit bed simulation model. The implicit bed simulation model was developed in STAR-CCM+®. The resistance of the spheres was modelled using the KTA (1981) pressure drop correlation and the structure of the bed was modelled using the porosity correlation proposed by Martin (1978), while the effective viscosity model of Giese et al. (1998), adjusted by a factor of 0.8, was used to model the velocity distribution in the near wall region. It was found that the structure in the inlet region of the bed, where two walls disturb the packing structure, can be modelled as the weighted average of the radial and axial porosity while the structure in the outlet regions can be modelled by letting the radial porosity increase linearly to unity. The basic shape of the velocity profile is established immediately when the fluid enters the bed. The amplitude of the velocity peaks however increase in magnitude until the velocity profile is fully developed at a distance approximately of five sphere diameters from the bed inlet. The profile remains constant throughout the bed until the outlet region of the bed is reached. In the outlet region a significant amount of fluid redistribution is observed. The amplitude of the velocity peaks is reduced and the position of the velocity peaks is shifted inwards towards the centre of the annular region. The fully developed velocity profile, predicted by the simulation model is in good agreement with profiles presented by amongst others Giese et al. (1998). The current model however also offers insight into the development of the profile through the inlet of the bed and the fluid redistribution, which occurs in the outlet region of the bed. / Thesis (M.Ing. (Mechanical Engineering))--North-West University, Potchefstroom Campus, 2010.

The purpose of this study was to investigate the velocity distribution through an annular bed packed randomly with equal sized spheres. Extensive research has been conducted on the velocity distribution inside packed beds packed with equal sized spheres, different sized spheres, deformed spheres, cylinders and Raschig-rings. A majority of these experimental and numerical studies focused on the cylindrical packed bed. These studies and numerical models are all confined to the velocity profile once the fluid flow is fully developed. The development of the velocity through the inlet region of the bed and the fluid flow redistribution in the outlet of the bed is thus neglected. The experimental investigation into the velocity distribution down stream of the annular packed bed of the HTTU indicated that the velocity profile was independent of the mass flow rate for a particle Reynolds number range of 439 £ Re £ 3453 . These velocity profiles did not represent the distribution of the axial velocity due to shortcomings associated with the single sensor hot wire anemometry system used to measure the velocity distribution. A numerical investigation, using the RANS CFD code STAR-CCM+®, into the velocity distribution downstream of an explicitly modelled bed of spheres indicated that the axial velocity distribution could be extracted from the experimental velocity profiles by using an adjustment factor of 0.801. This adjusted velocity profile was used in the verification of the implicit bed simulation model. The implicit bed simulation model was developed in STAR-CCM+®. The resistance of the spheres was modelled using the KTA (1981) pressure drop correlation and the structure of the bed was modelled using the porosity correlation proposed by Martin (1978), while the effective viscosity model of Giese et al. (1998), adjusted by a factor of 0.8, was used to model the velocity distribution in the near wall region. It was found that the structure in the inlet region of the bed, where two walls disturb the packing structure, can be modelled as the weighted average of the radial and axial porosity while the structure in the outlet regions can be modelled by letting the radial porosity increase linearly to unity. The basic shape of the velocity profile is established immediately when the fluid enters the bed. The amplitude of the velocity peaks however increase in magnitude until the velocity profile is fully developed at a distance approximately of five sphere diameters from the bed inlet. The profile remains constant throughout the bed until the outlet region of the bed is reached. In the outlet region a significant amount of fluid redistribution is observed. The amplitude of the velocity peaks is reduced and the position of the velocity peaks is shifted inwards towards the centre of the annular region. The fully developed velocity profile, predicted by the simulation model is in good agreement with profiles presented by amongst others Giese et al. (1998). The current model however also offers insight into the development of the profile through the inlet of the bed and the fluid redistribution, which occurs in the outlet region of the bed. / Thesis (M.Ing. (Mechanical Engineering))--North-West University, Potchefstroom Campus, 2010.