La présente thèse est consacrée à l'étude expérimentale et la simulation de la dynamique des gaz granulaires vibro-fluidisés. Les gaz granulaires sont caractérisés par une dissipation due aux collisions inélastiques. Pour maintenir cet état à l'équilibre mécanique (stationnarité), l'énergie est injectée en continu depuis les bords vibrant pour équilibrer la dissipation des vibrations. Ce système fournit une base d'étude de la physique des systèmes non-linéaires, hors équilibre thermodynamique et dissipatifs. Cette thèse insiste sur la nécessité d'intégrer, de comprendre et de rendre compte de la situation inhomogène de la distribution locale dans les gaz granulaires et permet la construction d'un nouveau modèle de gaz granulaires fluidisés par des vibrations. Cette approche inclut (i) des résultats expérimentaux 2d en micro-gravité dans l'Airbus A300 0-g de Novespace, des expériences 2d avec des cellules (et des vibrations) horizontales, des expériences 2d sur plan incliné (avec vibrations et cellules inclinées et avec une gravité effective variable), ainsi que des simulations de dynamique moléculaire par la méthode " event-driven" appliquée à chaque choc. Ces résultats confortent les simulations 3D de Liu et al. Les expériences en micro-gravité dans Airbus A380 (vol parabolique) permettent d'éviter les frottements avec les parois planes et éliminent l'effet de gravité. Les distributions locales de la vitesse dans la direction de vibrations sont asymétriques partout (à l'exception de la zone centrale de la cellule par raison de symétrie). L' équipartition de l'énergie n'est pas vérifiée dans la cellule, l'énergie est distribuée de manière inhomogène, anisotrope et directionnelle. La " température granulaire " n'est plus une mesure efficace pour décrire un tel système. On rend compte de ces résultats à l'aide d'une superposition de deux modèle gaussien pour décrire les profils locaux de vitesse asymétriques le long de la direction de vibration. Les résultats des simulations de dynamique moléculaire 2d en gravité nulle montrent les mêmes tendances et confortent les résultats l'expérimentaux (dissymétrie des distributions de vitesse locales). Cette dissymétrie est un effet à longue portée et est liée à la dissipation du système: Elle augmente si le coefficient de restitution billes-billes diminue ou lorsque le nombre de particules augmente. La dissymétrie disparaît lorsque les chocs billes-billes sont élastiques. Cet effet ne peut être ignoré et doit être traité comme la frontière d'une "nouvelle hydrodynamique". Dans les expériences de vibrations sur cellule 2D et plan inclinés parallèles aux vibrations, l'angle d'inclinaison a été modifié de façon systématique de l'horizontale à la verticale, pour simuler différentes gravités effectives. Les résultats confirment une dissymétrie locale des distributions de vitesse locales, à laquelle se rajoute une dissymétrie supplémentaire liée à la gravité, provoquant une densité différente en haut et en bas de la cellule. Ces études sont les prémices, nous le pensons, d’une nouvelle vision de la mécanique des gaz granulaires réels dissipatifs. / The present thesis is dedicated to the experimental and simulation study of vibro-fluidized granular gases dynamics. Granular gases are characterized by dissipation due to inelastic collisions. To keep a steady state, continuous energy is injected to balance dissipation by vibration. This system provides a platform to study the physics of non-linear, non-equilibrium and dissipative systems. This dissertation insisted on the necessity of understanding the local state in the granular gases and building a new model for vibration-fluidized granular gases. Research approach included experiments in micro-gravity, event-driven molecular dynamic simulation and experiments in tilted plane with various gravity. Micro-gravity experiments were performed on Airbus A380 (Parabolic flight) to avoid friction with the bottom and gravity field. A long range boundary effect is found to exist in 2D vibration granular gases. Local distributions of the velocity component in the vibration direction are asymmetric in the whole cell except for the center bin. In the system, energy equi-partition breaks down. ``Granular temperature" is not efficient to describe such a system. We proposed a superposition of two Gaussian model to describe the local asymmetric velocity profiles along the vibration direction. We demonstrated the performance of this model by the Airbus experimental data and others’ simulation works. Event-driven molecular dynamics simulation was utilized. Results showed support for experiment results. Furthermore, we found this long range boundary effect is related to the system dissipation. This effect becomes pronounced if the coefficient of restitution (e<1) decreased or the number of particles increased. For the elastic situation, there is no such effect. This effect cannot be ignored and treated only as a local boundary effect as in hydrodynamics. We studied a 2D vibration fluidized granular system in a tilted plane systematically. The inclined angle is changed from horizontal to vertical, changing the "effective gravity". These results also showed asymmetric local velocity distributions. Other than the number density profiles deviate from an exponential form, the spatial profiles of the number density of particles moving up and down are not equal, and non symetric from cell center.

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 separation of solids from biological wastewater treatment is an important step in the treatment process, as it has a significant impact on effluent water quality. The moving bed biofilm reactor (MBBR) technology is a proven upgrade or replacement wastewater treatment system for carbon and nitrogen removal. However, a challenge of this technology is the characteristics of the effluent solids that results in their poor settlement; with settling being the common method of solids removal. The main objective of this research is to understand and expand the current knowledge on the settling characteristics of MBBR produced solids and the parameters that influence them. In particular, in this dissertation, the impacts are studied of carrier types, biofilm thickness restraint design of carriers, and varying carbonaceous loading rates on MBBR performance, biofilm morphology, biofilm thickness, biofilm mass, biofilm density, biofilm detachment rate, solids production, particle size distribution (PSD) and particle settling velocity distribution (PSVD). With this aim, three MBBR reactors housing three different carrier types were operated with varying loading rates. In order to investigate the effect of carrier geometrical properties on the MBBR system, the conventional, cylindrically-shaped, flat AnoxK™ K5 carrier with protected voids was compared to two newly-designed, saddle-shaped Z-carriers with the fully exposed surface area. Moreover, the AnoxK™ Z-200 carrier was compared to the AnoxK™ Z-400 carrier to evaluate the biofilm thickness restraint design of these carriers, where the Z-200 carrier is designed for greater biofilm thickness-restraint. The Z-200 carrier is designed to limit the biofilm thickness to the level of 200 µm as opposed to 400 µm for the Z-400 carrier. Finally, to investigate the effects of varying carbonaceous loading rates on system removal performance, biofilm characteristics and solids characteristics, further analyses were performed at three different loading rates of 1.5 to 2.5 and 6.0 g-sBOD/m2·d in steady-state conditions. The PSD and the PSVD analyses were combined to relate these two properties. A settling velocity distribution analytical method, the ViCAs, was applied in combination with microscopy imaging and micro-flow imaging to investigate the relation of PSD and settling behaviour of MBBR produced particles. The obtained results have indicated that the carrier type significantly impacted the MBBR performance, biofilm, and particle characteristics. As such, the K5 carrier MBBR system demonstrated a statistically significantly higher carbonaceous removal rate and efficiency (3.8 ± 0.3 g-sBOD/m2·d and 59.9 ± 3.0% sBOD removal), higher biofilm thickness (281.1 ± 8.7 μm), higher biofilm mass per carrier (43.9 ± 1.0 mg), lower biofilm density (65.0 ± 1.5 kg/m3), lower biofilm detachment rate (1.7 ± 0.7 g-TSS/ m2·d) and hence lower solids production (0.7 ± 0.3 g-TSS/d) compared to the two Z-carriers. The Z-carriers' different shape exposes the biofilm to additional shear stress, which could explain why the Z-carriers have thinner and denser biofilm, resulting in higher solids production and lower system performance in comparison with K5. Moreover, the carrier type was also observed to impact the particle characteristics significantly. PSD analysis demonstrated a higher percentage of small particles in the Z-carrier system effluent and hence a significantly lower solids settling efficiency. Therefore, the solids produced in the K5 reactor have shown enhanced settling behaviour, consisting of larger particles with faster settling velocities compared to Z-carriers. This dissertation also investigated the effects of restraint biofilm thickness on MBBR performance by comparing the Z-200 biofilm thickness-restraint carrier to the Z-400 carrier. No significant difference was observed in removal efficiency, biofilm morphology, biofilm density, biofilm detachment rate, and solids production between the Z-200 to the Z-400 carriers. The PSD and the PSVD analyses did not illustrate any significant difference in the particles’ settling behaviour for these two biofilm thickness restraint carriers, indicating that the biofilm thickness-restraint carrier design was not a controlling factor in the settling potential of MBBR produced solids. Finally, this research studied the effect of varying loading rates and demonstrated a positive, strong linear correlation between the measured sBOD loading rate and the removal rate, indicating first-order BOD removal kinetics. The biofilm thickness, biofilm density and biofilm mass decreased when the surface area loading rate (SALR) was increased from 2.5 to 6.0 g-sBOD/m2·d. The solids retention time (SRT) was also shown to decrease by increasing the SALR, where the lowest SRT (1.7 ± 0.1 days) was observed at the highest SALR, with the highest cell viability (81.8 ± 1.7%). Significantly higher biofilm detachment rate and yield were observed at SALR 2.5, with the thickest biofilm and a higher percentage of dead cells. Consequently, a higher fraction of larger and rapidly settling particles was observed at SALR of 2.5 g-sBOD/m2·d, which leads to a significantly better settling behaviour of the MBBR effluent solids. This study expands the current knowledge of MBBR-produced particle characteristics and settling behaviour. A comprehensive understanding of the MBBR system performance and the potential influencing factors on the MBBR produced solids, particle characteristics, and their settleability will lead to optimized MBBR design for future pilot- and full-scale applications of the MBBR.

京都大学 / 0048 / 新制・課程博士 / 博士(工学) / 甲第18272号 / 工博第3864号 / 新制||工||1593(附属図書館) / 31130 / 京都大学大学院工学研究科原子核工学専攻 / (主査)教授 功刀 資彰, 教授 中部 主敬, 准教授 横峯 健彦 / 学位規則第4条第1項該当 / Doctor of Philosophy (Engineering) / Kyoto University / DGAM

In the advent of increasing the number of operable unmanned aerial systems (UAS) over the next years, a challenge exists in regard to the noise signature that these machines may generate. In this work, we perform advanced computational simulations to study the flow around an airfoil and the associated noise radiating to the near- and farield. The airfoil size and the freestream velocity are representative of a typical UAS. The study is aimed at investigating the characteristics of the aerodynamic noise radiating from an airfoil at various angles of attack, Reynolds number and Mach number. The numerical tool is a high-order compressible Navier-Stokes solver, using Runge-Kutta explicit time integration and dispersion-relation-preserving spatial discretization. Various results in terms of velocity and pressure distribution around the airfoil, and sound pressure level spectra calculated from different probe points located in the near- and farield are compared to each other and discussed.

No / This paper presents an integration procedure for the Reynolds-averaged Navier-Stokes equations for the determination of the distribution of the streamwise velocity using the vertical component. This procedure is dedicated to the outer region and central part of channels. The proposed model is applicable to both rough and smooth flow regimes, provided the velocity at the inner-outer boundary has been properly defined. To generate a simplified expansion, a number of hypotheses are proposed, focusing in particular on the analytical modeling of the vertical component by adopting a negligible viscosity. The proposed hypotheses are validated by the experimental data existing in the literature. The proposed simplified expansion is studied through a sensitivity analysis and proved consistent in regards to model experimental data. The proposed model seems capable of demonstrating different kinds of flows, including dip phenomenon flow patterns.

Yes / In this paper, an improved Shiono-Knight model (SKM) has been proposed to calculate the rectangular compound open channel flows by considering a Multi-Zonal (MZ) approach in modelling turbulence and secondary flows across lateral flow direction. This is an effort to represent natural flows with compound shape more closely. The proposed model improves the estimation of secondary flow by original SKM model to increase the accuracy of depthaveraged velocity profile solution formed within the transitional region between different sections (i.e. between main-channel and floodplain) of compound channel. This proposed MZ model works by sectioning intermediate zones between floodplain and main-channel for running computation in order to improve the modelling accuracy. The modelling results have been validated using the experimental data by national UK Flood Channel Facility (FCF). It has been proven to work reasonably well to model secondary flows within the investigated compound channel flow cases and hence produce better representation to their flow lateral velocity profile.

Possibilities of generating plasma turbulence to provide control of space weather processes have been of particular interest in recent years. Such turbulence can be created by chemical released into a magnetized background plasma. The released plasma clouds are heavy ions which have ring velocity distribution and large free energy to drive the turbulence. An electromagnetic hybrid (fluid electrons and particle ions) model incorporating electron inertia is developed to study the generation and nonlinear evolution of this turbulence. Fourier pseudo-spectral methods are combined with finite difference methods to solve the electron momentum equations. Time integration is accomplished by a 4th-order Runge-Kutta scheme or predicator-corrector method. The numerical results show good agreement with theoretical prediction as well as provide further insights on the nonlinear turbulence evolution. Initially the turbulence lies near harmonics of the ring plasma ion cyclotron frequency and propagates nearly perpendicular to the background magnetic field as predicted by the linear theory. If the amplitude of the turbulence is sufficiently large, the quasi-electrostatic short wavelength ion cyclotron waves evolve nonlinearly into electromagnetic obliquely propagating shear Alfven waves with much longer wavelength. The results indicate that ring densities above a few percent of the background plasma density may produce wave amplitudes large enough for such an evolution to occur. The extraction of energy from the ring plasma may be in the range of 10-15% with a generally slight decrease in the magnitude as the ring density is increased from a few percent to several 10's of percent of the background plasma density. Possibilities to model the effects of nonlinear processes on energy extraction by introducing electron anomalous resistivity are also addressed. Suitability of the nonlinearly generated shear Alfven waves for applications to scattering radiation belt particles is discussed. / Master of Science

Novel settling characteristic metrics were developed based on the fundamental mechanisms of coagulation, flocculation, and settling. The settling metrics determined parameters that are essential in monitoring and optimizing the activated sludge process without the need for expensive or specialized equipment. Current settling characteristic measurements that don't require specialized instruments such as sludge volume index (SVI) or initial settling velocity (ISV), have no fundamental basis in solid-liquid separation and only indicate whether settling is good or bad without providing information as to limitations present in a sludge matrix. Furthermore, the emergence of aerobic granulation as a potential pathway to mitigate solid-separation issues further stresses the need for new settling characteristic metrics to enable integration of the new technology with the current infrastructure. The granule or intrinsic aggregate fraction in different types was of sludge was quantified by simulating different surface overflow rates (SOR). The technique named Intrinsic Settling Classes (ISC) was able to separate granules and floc by simulating high SOR values due to the lack of a flocculation time needed for granules. The method had to be performed in a discrete settling environment to characterize a range of flocculation behavior and was able to classify the granular portion of five different types of sludge. ISC was proven to accurately (±2%) determine the granule fraction and discrete particle distribution. The major significance of the test is its ability to show if a system is producing particles that will eventually grow to become granules. This methodology proved to be very valuable in obtaining information as to the granular fraction of sludge and the granular production of a system. Flocculent settling (stokesian) was found to be predominant within ideally operating clarifiers, and the shift to 'slower' hindered settling (non-stokesian) causes both failure and poor effluent quality. Therefore, a new metric for settling characteristics was developed and classified as Limit of Stokesian Settling (LOSS). The technique consisted of determining the total suspended solids (TSS) concentration at which mixed liquor settling characteristics transition from stokesian to non–stokesian settling. An image analytical technique was developed with the aid of MATLAB to identify this transition. The MATLAB tool analyzed RGB images from video, and identified the presence of an interface by a dramatic shift in the Red indices. LOSS data for Secondary activated-sludge systems were analyzed for a period of 60 days at the Blue Plains Advanced Wastewater Treatment Plant. LOSS numbers collected experimentally were validated with the Takacs et al. (1991) settling model. When compared to flux curves with small changes in the sludge concentration matrix, LOSS was found to be faster at characterizing the hindered settling velocity and was less erratic. Simple batch experiments based on the critical settling velocity (CSV) selection were used as the basis for the development of two novel parameters: threshold of flocculation/flocculation limitation (TOF/a), and floc strength. TOF quantified the minimum solids concentration needed to form large flocs and was directly linked to collision efficiency. In hybrid systems, an exponential fitting on a CSV matrix was proposed to quantify the collision efficiency of flocs (a). Shear studies were conducted to quantify floc strength. The methods were applied to a wide spectrum of sludge types to show the broad applicability and sensitivity of the novel methods. Three different activated sludge systems from the Blue Plains AWWTP were monitored for a 1 year period to explore the relationship between effluent suspended solids (ESS) and activated sludge settling and flocculation behavior. Novel metrics based on the transitional solids concentration (TOF, and LOSS) were also collected weekly. A pilot clarifier and settling column were run and filmed to determine floc morphological properties. SVI was found to lose sensitivity (r < 0.20) when characterizing ISV above a hindered settling rate of 3 m h-1. ISV and LOSS had a strong correlation (r = 0.71), but ISV was subject to change, depending on the solids concentration. Two sludge matrix limitations influencing ESS were characterized by transition concentrations; pinpoint floc formation, and loose floc formation. Pinpoint flocs had TOF values above 400 mg TSS L-1; loose floc formation sludge had TOF and LOSS values below 400 mg TSS L-1 and 900 mg TSS L-1, respectively. TOF was found to correlate with the particle size distribution while LOSS correlated to the settling velocity distribution. The use of both TOF and LOSS is a quick and effective way to characterize limitations affecting ESS. / Ph. D.