University, Diss., 2000--Heidelberg.

In a pressurized water reactor (PWR), pressurized water flows over fuel rods containing radioactive uranium. Potential failure of these nuclear fuel rods is a primary concern, as fuel rod failure typically results in power generation losses and reactor downtime. Thermal parameters such as critical heat flux have traditionally been utilized as performance metrics to ensure that the reactor core remains stable even during failure events. Recently, fuel leaking events have occurred which have resulted in excess debris buildup on fuel rods and fuel grid array mixing devices. Understanding the flow field surrounding these nuclear fuel rods is critical in predicting where crud could deposit. Although CFD simulations have been conducted to characterize the fluid flow around fuel rod bundles, limited experimental data characterizing the mechanics of this fluid flow exists in the current literature. This study will present experimental data collected detailing the fluid flow around a rod bundle geometry using a novel matched refractive index particle tracking velocimetry (PTV) technique over a 3D volume cross section of a prototypical nuclear fuel rod bundle. Velocimetry tracking will be performed in order to characterize the mechanics of the fluid flow. Using optical distortion mitigation techniques and various image processing methods, data from multiple cameras was used to assemble 3-dimensional velocity information of a turbulent fluid region. Results are compared to the solution of a k-epsilon unsteady RANS numerical simulation.

A particle tracking velocimetry (PTV) system is used to investigate the turbulent properties at the free surface of shallow shear flows and a shallow vortex street (VS) wake flow. The resolution of the PTV system enables information to be gathered regarding the large-scale turbulent structure of these flows, and also enables analysis to proceed in both the temporal and spatial domains. Statistical tools such as the probability density function (PDF), autocorrelation and power spectral density (PSD) are utilised to characterise the turbulent properties at the flow surface. Two supercritical flows and one subcritical shallow shear flow are analysed. Taylor's frozen turbulence hypothesis is shown to be valid for these flows, and the integral length scales indicate that 2D isotropic structures with scales larger than the flow depth are present at the free surface. Such large-scale structures at the free surface are consistent with observations from dye visualisation experiments and with "spiral eddies" identified by Kumar, et al (1998). The longitudinal extent of near and intermediate wake fields for the shallow VS wake flow is well defined by the integral wake length scale specified by v.Carmer (2005). The near wake region is characterised by high rates of exchange between the mean flow and large-scale 2D coherent structures (2DCS). In the intermediate field, the rate of decay of the turbulent stress components greatly diminishes as the 2DCS are stabilised and dissipated under the action of bed friction. Multiple peaks are observed in the power spectral density of the turbulent fluctuations. The periodic shedding of 2DCS behind the circular cylinder is characterised by an energy peak at a Strouhal number of 0.21, and further energy peaks are observed in the near-wake region. The PSD estimates are consistent with the findings of v.Carmer (2005) in which a -5/3 decay law to high frequencies is observed, and no evidence of an inverse energy cascade is present.

Heidelberg, Univ., Diplomarb., 2000.

Dans le cadre du design des assemblages combustibles, les écoulements turbulents induits par les grilles de maintien provoquent de fortes fluctuations de pression dues aux niveaux de vitesses du fluide induisant la vibration des crayons combustibles. Or ces vibrations sont une source de vieillissement accéléré et d’usure prématurée des assemblages combustibles dont l’origine convient d’être mieux analysée.Dans cette l'optique, des simulations numériques ont été utilisées pour prédire les niveaux de vitesses transverses ainsi que le niveau de turbulence en aval de la grille. Cependant, les codes CFD doivent être validés avec des mesures expérimentales conduisant ainsi à une meilleure compréhension des structures mis en jeu dans le processus de mélange. Or la caractéristique de ce type d'écoulement repose sur son caractère tridimensionnel entre les faisceaux de crayons où l'accès optique est limité et à l'heure actuelle, aucune mesure de vitesse 3D n'a été réalisée. C'est pourquoi nous nous proposons de présenter dans cette étude l'application d'une technique de mesure 3D résolue en temps dans une boucle hydraulique transparente. L’approche consiste à utiliser une stratégie de suivi, dans lequel des particules traceuses individuelles sont d'abord détectées dans l’écoulement et ensuite suivies dans le temps – PTV 3D.Cette thèse présente les résultats obtenus dans une veine d’essais hydraulique en utilisant une approche de suivi à l’aide de deux caméras. Dans un premier temps, une nouvelle extension de la technique de défocalisation permettant de récupérer les positions des particules dans le temps avec une caméra est présentée et la méthodologie pour récupérer les déplacements individuels est décrite. Des cas d’application et de validation de la technique sont présentés afin de mieux quantifier les erreurs de mesures. Ensuite, la boucle hydraulique est introduite et son écoulement caractérisé par des techniques établies de mesure de vitesse. Finalement, les mesures 3D à deux caméras sont réalisées et discutées. / In order to validate the numerical simulation of the mixing phenomena downstream a spacer grid in the reactor core, reaching the 3rd velocity component through experimental studies is of major importance since mechanical structures responses are linked to transverse velocities. Nevertheless the main difficulty relies on applying a non-intrusive velocimetry technique around obstacles composed of a vertical set of rods. So far only 1D and 2D temporal studies have been performed.Hence, numerical methods have been widely used to predict correctly those transverse velocities and the turbulence level downstream the grid. However, CFD codes are to be validated with experimental measurements leading to a better understanding of the detailed flow structure in the mixing process. Therefore we propose to present in this study the application of a 3D time resolved velocity measurement technique to a hydraulic test facility. The approach is to use a tracking strategy, in which individual particles are first detected then followed in time - 3D time resolved Particle Tracking Velocimetry.This PhD Thesis presents results obtained in a hydraulic test section using a tracking based approach with two cameras. At first, a new extension of the defocussing technique to recover particles location in time with one camera is presented and the methodology o get individual velocity vector is then described. Applications and validations of the technique in dedicated flows allow to determinate and quantify measurement uncertainties. Then the hydraulic test section is introduced together with preliminary flow characterization using Laser Doppler Velocimetry or 2D time resolved PIV analysis. Finally, two cameras measurements are reported and post-processing techniques are discussed.

The analysis of turbulent flows provides a unique scientific challenge whose solution remains central to unraveling the fundamental nature of all fluid dynamics. Measuring and predicting turbulent flows becomes even more difficult when considering a two-phase flow, which is a commonly encountered engineering problem across many disciplines. One such example, the ingestion of foreign debris into a gas turbine engine, provided the impetus for this study. Despite more than 40 years of research, operation with a particle-laden inlet flow remains a significant problem for modern turbomachines. The purpose, therefore, is to develop experimental methods for investigating multi-phase flows relevant to the cooling of gas turbine components. Initially, several generic components representing turbine cooling designs were evaluated with a particle-laden flow using a special high temperature test facility. The results of this investigation revealed that blockage was highly sensitive to the carrier flowfield as defined by the cooling geometry. A second group of experiments were conducted in one commonly used cooling design using a Time Resolved Digital Particle Image Velocimetry (TRDPIV) system that directly investigated both the carrier flowfield and particle trajectories. Traditional PIV processing algorithms, however, were unable to resolve the particle motions of the two-phase flow with sufficient fidelity. To address this issue, a new Particle Tracking Velocimetry (PTV) algorithm was developed and validated for both single-phase and two-phase flows. The newly developed PTV algorithm was shown to outperform other published algorithms as well as possessing a unique ability to handle particle laden two-phase flows. Overall, this work demonstrates several experimental methods that are well suited for the investigation of wall-bounded turbulent two-phase flows, with a special emphasis on a turbine cooling method. The studies contained herein provide valuable information regarding the previously unknown fluid and particle dynamics within the turbine cooling system. / Ph. D.

Ces travaux de thèse s’intègrent au sein du projet Biocarb lancé par le Commissariat à l’Énergie Atomique dont l’objectif est de développer des procédés de production de carburants liquides ou gazeux à partir de gaz de synthèse riche en H2 et CO obtenu par gazéification de la biomasse lignocellulosique. L’objectif de cette étude est d’étudier le comportement de particules de biomasse millimétriques lors de la pyrolyse dans des conditions types de gazéifieurs industriels tels que les réacteurs à lit fluidisé ou à flux entraîné, qui fonctionnent pour des flux de chaleur élevés (105 – 106 W⋅m-2) et pour de hautes températures (>800°C). Tout d’abord, des expériences de pyrolyse sont menées à 800 et 950°C dans un four à chute de laboratoire sur des particules de bois entre 350 et 800 μm. Les résultats montrent que dans les conditions de l’étude, l’augmentation de la taille de la particule augmente seulement la durée de la pyrolyse mais ne modifie pas les rendements ou la composition du solide et du gaz au cours de la pyrolyse. Par ailleurs, des mesures basées sur la technique de PTV (Particle Tracking Velocimetry) sont réalisées à température ambiante pour caractériser la taille et la densité des particules de bois brut et de résidu, et valider une corrélation donnant le coefficient de traînée qui sert à calculer le temps de séjour des particules dans le réacteur. On constate à la fin de la pyrolyse une diminution de la densité comprise entre 70 et 80% ainsi qu’une diminution de la taille des particules entre 25 et 40%. Les résultats montrent également que la vitesse de glissement de la particule et l’évolution de ses propriétés doivent être prises en compte lors du calcul de sa vitesse. Enfin, à partir des résultats expérimentaux, un modèle unidimensionnel à coeur rétrécissant est développé pour décrire le comportement d’une particule de bois lors de sa pyrolyse. Le modèle est capable de prévoir l’évolution du rendement en solide, en gaz total et en goudrons au cours de la pyrolyse ainsi que la vitesse de glissement de la particule et son temps de séjour dans le réacteur.L’analyse de sensibilité du modèle montre que même pour des particules millimétriques, une connaissance précise de la chaleur de réaction associée à la pyrolyse, de la densité du bois et de la conductivité thermique du résidu solide est essentielle / The present work is part of a project of the French energy research centre Commissariat à l’Energie Atomique. The goal of the project is to develop processes of production of gaseous or liquid fuel from synthesis gas obtained by gasification of lignocellulosic biomass. The objective of the present work is to study the pyrolysis behaviour of millimetric biomass particles under the operating conditions encountered in fluidized bed or entrained flow gasifiers, namely high external heat flux (105 – 106 W⋅m-2) and high temperature (> 800°C). First, pyrolysis experiments are conducted at 800 and 950°C in a lab-scale drop tube reactor on wood particles between 350 and 800 μm. The results show that under the explored conditions, the increase of the particle size only increases the time required for pyrolysis but does not affect the product distribution during pyrolysis. Since in the pyrolysis experiments, the particle residence time cannot be directly measured, PTV (Particle Tracking Velocimetry) measurements are performed at room temperature to characterize the evolution of the particle size and density along pyrolysis and to validate a drag coefficient correlation for the particle residence time calculation. The optical measurements show that at the end of pyrolysis there is a decrease of particle density of 70 – 80% and of particle size of 25 – 40%. It is also proven that the particle slip velocity cannot be neglected and that the change of these particle properties must be taken into account for the calculation of the particle slip velocity and residence time. Finally, based on these experimental results, a 1D shrinking-core model is developed that is able to predict the solid/gas/tar yields and the residence time of a single particle along pyrolysis in the drop tube reactor. It is validated on both the pyrolysis and optical experiments. The model sensitivity analysis shows that even for millimetric particles, the accurate knowledge of the heat of pyrolysis, of the wood density and of the char thermal conductivity is essential

La turbulence est connue pour sa capacité à disperser efficacement de la matière, que ce soit des polluantes dans les océans ou du carburant dans les moteurs à combustion. Deux considérations essentielles s’imposent lorsqu’on considère de telles situations. Primo, l’écoulement sous-jacente pourrait avoir une influence non-négligeable sur le comportement des particules. Secundo, la concentration locale de la matière pourrait empêcher le transport ou l’augmenter. Pour répondre à ces deux problématiques distinctes, deux dispositifs expérimentaux ont été étudiés au cours de cette thèse. Un premier dispositif a été mis en place pour étudier l’écoulement de von Kàrmàn, qui consiste en une enceinte fermé avec de l’eau forcé par deux disques en contra-rotation. Cette écoulement est connu pour être très turbulent, inhomogène, et anisotrope. Deux caméras rapides ont facilité le suivi Lagrangien des particules isodenses avec l’eau et petites par rapport aux échelles de la turbulence. Ceci a permis une étude du bilan d’énergie cinétique turbulente qui est directement relié aux propriétés de transport. Des particules plus lourdes que l’eau ont aussi été étudiées et montrent le rôle de l’anisotropie de l’écoulement dans la dispersion des particules inertielles. Un deuxième dispositif, un écoulement de soufflerie ensemencé avec des gouttelettes d’eau micrométriques a permis une étude de l’effet de la concentration locale de l’eau sur la vitesse de chute des gouttelettes grâce à une montage préexistant. Un modèle basé sur des méthodes théorique d'écoulements multiphasiques a été élaboré enfin de prendre en compte les effets collectifs de ces particules sedimentant dans un écoulement turbulent. Les résultats théoriques et expérimentaux mettent en évidence le rôle de la polydispersité et du couplage entre les deux phases dans l’augmentation de la sédimentation des gouttelettes. / Turbulence is well known for its ability to efficiently disperse matter, whether it be atmospheric pollutants or gasoline in combustion motors. Two considerations are fundamental when considering such situations. First, the underlying flow may have a strong influence of the behavior of the dispersed particles. Second, the local concentration of particles may enhance or impede the transport properties of turbulence. This dissertation addresses these points separately through the experimental study of two different turbulent flows. The first experimental device used is the so-called von K\'arm\'an flow which consists of an enclosed vessel filled with water that is forced by two counter rotating disks creating a strongly inhomogeneous and anisotropic turbulence. Two high-speed cameras permitted the creation a trajectory data base particles that were both isodense and heavier than water but were smaller than the smallest turbulent scales. The trajectories of this data base permitted a study of the turbulent kinetic energy budget which was shown to directly related to the transport properties of the turbulent flow. The heavy particles illustrate the role of flow anisotropy in the dispersive dynamics of particles dominated by effects related to their inertia. The second flow studied was a wind tunnel seeded with micrometer sized water droplets which was used to study the effects of local concentration of the settling velocities of these particles. A model based on theoretical multi-phase methods was developed in order to take into account the role of collective effects on sedimentation in a turbulent flow. The theoretical results emphasize the role of coupling between the underlying flow and the dispersed phase.

Inclined negatively buoyant jets are commonly used to dispose brine effluent produced by desalination plants. Desalination and associated research has expanded in recent years due to the continued depletion and degradation of natural potable water sources. Desalination plants are the preferred option for meeting water demand deficits in many countries around the world. Inclined negatively buoyant jets are produced when the brine is discharged at an upward inclined angle via an offshore pipeline and diffuser system. Previous experimental studies have focused on the rapid mixing and dilution achieved by these discharges, as well as geometric parameters. Dilution measurements between these experimental studies vary significantly, which is possibly due to variations in the location of a lower boundary on observed flow behaviour. In the present study, velocity field information is experimentally measured for inclined negatively buoyant jets and compared to integral model predictions. Experiments are conducted with and without a lower boundary influencing observed flow behaviour, thus allowing the effects of a lower boundary to be determined. The particle tracking velocimetry experimental technique is employed to measure near field velocities of these discharges. Firstly, discharges with source angles between 15\degree and 75\degree are investigated without boundary influence in stationary ambient conditions. The source was a minimum of 655 mm above the bottom of the experimental tank to ensure there was no lower boundary influence on observed behaviour. Time-averaged and fluctuating data are extracted along the trajectory of discharges. All non-dimensionalised geometric and centreline velocity parameters are found to collapse. Empirical coefficients are compared to previous experimental studies and integral model predictions. A new detrainment model is developed to predict the behaviour of inclined negatively buoyant jets without boundary influence. The model further develops recent attempts to allow for buoyancy flux reduction along the flow path. The reduction in buoyancy flux is dependent on the local parameters of the flow and simulates experimentally observed detrainment. Dilution, geometric, and velocity predictions are found to be improved over previous models when compared to experimental data. Finally, a raised platform was placed inside the experimental tank to determine the influence of a lower boundary on inclined negatively buoyant jets. Source angles of 30\degree, 45\degree, and 60\degree are investigated at three different non-dimensional source heights. The lower boundary is horizontal and ambient conditions are again stationary. Discharges impinge the lower boundary before forming a radially spreading layer along the boundary. Geometric and velocity data are compared to the first set of experiments in this study to determine the influence of the lower boundary on observed flow behaviour. Empirical coefficients at maximum height are similar with and without the influence of the boundary, whereas coefficients are substantially influenced at the return point when the boundary is present.

Tsunami are a very dangerous natural hazard, as highlighted in recent years by the Indian Ocean Tsunami of 2004 and the Japan Tsunami of 2011. In the last decade, tsunami have claimed hundreds of thousands of lives, and caused billions of dollars in damage around the world. The hazard posed to coastal communities by tsunami is expected to increase in the future, due to population growth, intensification of coastal development and sea level rise due to climate change. Tsunami may be generated by a number of different source mechanisms. One such source mechanism is a submarine landslide, which can occur in a number of marine environments containing significant sediment accumulation on a sloping seafloor. The high amplitudes and rapid celerities of landslide-generated tsunami make them very dangerous to communities in the vicinity of the landslide, although these waves do not possess the potential for transoceanic devastation. The objectives of this research project are to carry out a series of two-dimensional physical experiments investigating the waves generated by a rigid block landslide moving along a horizontal boundary. The use of a horizontal boundary has the advantage that waves propagating in the offshore and onshore directions may be measured (unlike previous studies using sloping boundaries). The landslide motion is provided by a mechanical system, allowing testing of a broad range of motion, and isolation of the wavemaking properties of different phases of landslide motion. Experiments are carried out in a 14.66 m long flume, with width 0.25 m and working depth 0.50 m. A false floor installed in the flume provides the sliding surface for the landslide motion, and houses the mechanical system. A series of preliminary particle tracking velocimetry experiments confirm the ability of the mechanical system to achieve its velocity targets to within 5% or better, depending on the parameters of the landslide motion. Full spatial and temporal resolution of the wave field is achieved using a laser-induced fluorescence technique to identify the air-water interface to sub-pixel accuracy. The measurements obtained using laser-induced fluorescence are validated against measurements from a resistance wave gauge, with sub-millimetre agreement. In an additional experiment, the particle tracking velocimetry technique provides measurements of the subsurface velocity field. The landslide motion during all experiments consists of an initial period of constant acceleration, followed by a period of constant velocity, followed by a deceleration to rest (at the same rate as the initial acceleration). The landslide acceleration generates two dispersive packets of waves, travelling in the offshore and onshore directions. The offshore-propagating wave packet contains a leading crest and the onshore-propagating wave packet contains a leading trough, with both waves approaching the shallow water limit. A free surface depression forms above the landslide during its constant-velocity motion, and its amplitude may be predicted to within approximately 20% using standard hydraulic theory (considering a frame of reference moving with the landslide). The offshore-propagating waves passing over the landslide cause the amplitude of this depression to fluctuate over time. The deceleration of the landslide generates two additional packets of waves with the opposite polarity to the waves generated by the landslide acceleration. The full spatial and temporal resolution of the generated wave field allows the calculation of the potential energy within the wave field. Additionally, the energy (and mass) within the onshore- and offshore-propagating wave packets may be estimated by calculating these quantities within the onshore and offshore regions of the experimental domain. The wave packets generated by the initial landslide acceleration transport positive mass in the offshore direction, and negative mass in the onshore direction. This mass transport is balanced by the waves generated during the deceleration of the landslide. The nondimensional landslide acceleration, landslide Froude number and submergence depth are varied during the physical experiments. The landslide Froude number has the greatest effect on the behaviour of the generated wave field. At low Froude numbers, the wave field is dominated by the waves generated by the acceleration and deceleration of the landslide. As the Froude number increases, the onshore-propagating waves become negligible in amplitude compared to the offshore-propagating waves. Additionally, the free surface depression increases in amplitude and a group of short-wavelength waves become trapped behind the landslide. These waves exhibit highly nonlinear behaviour at landslide Froude numbers greater than 0.5. The simple experimental geometry allows comparison between the measured wave fields with the predictions of three mathematical models. Two inviscid-irrotational models, differing in their treatment of the bottom boundary condition, provide comparisons over the entire parameter space. These models under-predict the amplitudes of the generated waves, and fail to correctly predict the ongoing interaction between the landslide and the offshore-propagating waves. The inclusion of bottom boundary nonlinearity improves the predictions of the amplitude of the leading waves, and the potential energy within the wave field. However, both of the inviscid models do not predict the extent of wave trapping behaviour behind the landslide observed in the experiments. A viscous model, formulated in the DNS solver Gerris, improves the predictions of wave trapping (and amplitude in general) in one experiment. Although the model still under-predicts the amplitudes of the generated waves, it correctly predicts the amplification of the waves behind the landslide during its constant-velocity motion. The failure of the inviscid models to predict the amplitudes of these waves can be mostly attributed to the linearised free surface condition used by both models. The presence of the turbulent wake may also have a secondary effect on these predictions. An extension of the linear inviscid-irrotational model to three dimensions allows the effect of the landslide width on the amplitudes of the generated waves to be determined. As the width increases, the behaviour of the waves approaches the two-dimensional limiting case.