Ce travail est consacré à l'étude expérimentale des écoulements diphasiques et triphasiques d'inclusions (bulles, gouttes, sphères) en milieux tant newtoniens que non newtoniens à l'échelle microscopique et mésoscopique, en utilisant la visualisation par une caméra rapide, la vélocimétrie par images des particules (PIV) ainsi que la micro-vélocimétrie par images des particules. Des bulles et des gouttes ont été étudiées expérimentalement depuis leur formation, en passant par leur déformation jusqu'à leur coalescence. La formation de bulles dans des micro-mélangeurs a été étudiée et caractérisée par l'obtention de champs de vitesses. Différents paramètres, tels que le cisaillement, la géométrie de la zone de formation, les débits ou encore les propriétés physiques ont été testés afin de développer des lois d'échelles. La traversée d'une interface liquide-liquide par une inclusion a été abordée par des expériences originales, permettant de décrire la dynamique du phénomène, de définir des nombres adimensionnels et de mettre en évidence des instabilités interfaciales. L'effet Weissenberg a également été étudié aux différentes échelles afin de comprendre les phénomènes conduisant à son amplification. Enfin, en milieu viscoélastique et rhéofluidifiant, nous avons caractérisé l'écoulement autour d'une inclusion isolée solide par l'obtention de champs de vitesses. Ces résultats ont permis de confirmer l'origine viscoélastique du sillage négatif et de prédire ses caractéristiques / The present work was devoted to the experimental study of the multiphase flow around inclusions in both Newtonian and non-Newtonian media at respectively microscopic and mesoscopic scales, by means of the Particle Image Velocimetry (both PIV and µ-PIV) and fast camera visualization. Bubbles and drops were experimentally studied starting from their formation, the rising and up to their coalescence and fragmentation. Bubble formation in micro-mixers was also investigated and characterized by measuring the liquid velocity fields. Different parameters, such as the shear rate, the geometry of the micro-mixer, the flow rates or the physical properties were tested to develop correlations of power-law kind. The deformation of a liquid-liquid interface due to the passage of an inclusion was investigated by original experiments in order to observe and describe the dynamics of the phenomenon, to define the reliable dimensionless numbers and also to highlight several interfacial instabilities. The Weissenberg effect was also studied at different scales in various non-Newtonian fluids in comparison with Newtonian fluids to understand the amplification phenomena under the combining effects of instability and normal forces. Finally, the flow around an isolated solid inclusion was characterized by performing the measurements of velocity fields in viscoelastic and shear-thinning fluids. These results were used to confirm the viscoelastic origin of the negative wake behind the solid sphere and to model its characteristics

This thesis reports on investigations carried out to study of the effect of horizontal wall jets on rough, fixed and mobile beds in open channel flow. Experimental tests were carried out, using fixed and mobile sediment beds. Computer simulation models for the flow within the jet and resulting sediment transport were developed and their results analysed in this study. In the experimental phase, tests were carried out with both fixed and mobile sediment beds. The shape of the water surface, numerous point velocity measurements and measurements of the evolving scour hole shape were made. Detailed descriptions of the turbulent flow field over a fixed rough bed and for scour holes at equilibrium were obtained for a range of initial jet conditions. Fully turbulent, multiphase flow was modelled using the Fluent Computational Fluid Dynamics software. This was used to analyze the flow caused by a jet in a rectangle open-channel with a rough bed, and also the flow pattern in a channel with a local scour hole. The volume of fluid (VOF) multiphase method and K- model was used to model the fluid flow in both cases. The model predictions of velocity and shear stress were compared against experimental observations. The experimental data was used to develop new empirical relationships to describe the pattern of boundary shear stress caused by a wall jet over fixed beds and in equilibrium scour holes. These relationships were linked with existing bed-load transport rate models in order to predict the temporal evolution of scour holes. An analytical model describing the relationship between the wall jet flow and the development of a local scour hole shape was reported and its predictions compared with experimental data.

Ce manuscrit s'intéresse au devenir de contaminants organiques dans le sol, et plus précisément celui des composés organo-halogénés volatils (COHV).Il propose des outils d'aide à la décision en utilisant la modélisation numérique appliquée à des problématiques environnementales portant sur le traitement de sites et sols pollués par des COHV. Il présente, à ma connaissance, la première modélisation numérique en 3D, grâce au simulateur TMVOC, d'une technique de dépollution physique, à savoir le pompage réussi au sein d'un aquifère d'un liquide en phase non-aqueuse plus dense que l'eau ou Dense Non-Aqueous Phase Liquid (DNAPL). Les très bons résultats de cette simulation permettent d’envisager l’optimisation d’un système de pompage asservi pour dépolluer un site contaminé aux COHV. Cette thèse s'intéresse aussi à une méthode de d'évaluation, à la fois qualitative et quantitative, de l'efficacité du pompage : un test de traçage utilisant des traceurs bisolubles à coefficients de partage variables (partitioning interwell tracer test – PITT). Ce PITT permet de connaître à la fois la répartition spatiale, au sein d’un aquifère, d’une bulle de DNAPL et aussi d’en évaluer sa saturation et donc son volume. Ce manuscrit met à disposition des outils élaborés et adaptés au traitement d'un DNAPL dans un aquifère. La combinaison des deux techniques présentées, à savoir la modélisation d'un pompage de DNAPL et un PITT, sont parfaitement reproductibles dans des environnements similaires et à l'échelle industrielle. Enfin, ces méthodes permettent de réduire considérablement les coûts de caractérisation (PITT) et d'exploitation par l'optimisation de systèmes de pompage / This manuscript deals with the fate of organic contaminants in soil, more precisely of volatile organo-chlorinated compounds (VOHC) and offers some decision making techniques and tools using numerical modelling applied to environmental issues about the treatment of soils contaminated by VOHC. It presents, to the best knowledge of the auhor, the first 3D numerical modelling, with the simulator TMVOC, of a physical treatment technique, namely the successful pumping within an aquifer, of a Dense Non-Aqueous Phase Liquid (DNAPL). The very good results of this simulation give the opportunity to design a controlled and automated pumping system to remediate a polluted site. This thesis also deols with a tracer test assessment method, both qualitative and quantitative, for the efficiency of this pumping : a partitioning interwell tracer test (PITT). The PITT allows us to determine the spatial repartion of the DNAPL and also to assess its saturation, hence its volume. This manuscript offers elaborated tools adapted to the remediation of a DNAPL in an aquifer. The combination of these two techniques, namely the DNAPL pumping and the PITT, are perfectly reproducible in similar environments up to an industrial scale. Finally, exploitation and characterisation costs for DNAPL remediation can be extremely reduced by numerical modelling and optimisation

In the present work, the flow behavior of an oil-in-water emulsion around and inside the measurement chamber of an in-line optical sensor probe is studied. The emulsion consists of a metalworking fluid, with oil droplets diameter ranging from 100 nanometers to 100 micrometers. The design of the UV-Vis light spectrometer probe is in accordance with the concept proposed in the research project named EPM (Emulsion Process Monitor in Metalworking Processes), carried out within the scope of the BRAGECRIM program between the University of São Paulo and the University of Bremen. This study is based on the numerical simulation of the interactions between the emulsion and the measurement system using computational fluid dynamic techniques, and is aimed at evaluating the effects of the probe geometry, its position relative to the flow field, and fluid properties on the expected spectrometer readings. Such effects are correlated with changes in droplet concentration and or droplet size segregation inside the measurement chamber of the optical probe, which can cause changes in the scattered light intensity readings. Segregation effects due to flow disturbances around the probe can be neglected under normal measurement conditions, with the probe facing the slit area inlet against the flow stream. Based on the simulation results, even if the probe is misplaced, the effect on the measurements is still insignificant. In-situ measurements carried out in a laboratory set up installed in the injection tube of a drilling machine support the simulation results, since no segregation effect related to the measurement system was observed. In addition to the study, the possibility of bacteria attachment on the internal glass walls of the probe was evaluated and it was found that when the flow velocity is large enough to produce a wall shear stress of about 3-5 Pa the bacterial contamination may be avoided. Changes in the probe geometry are proposed in order to attain an isokinetic condition for the flow around and inside the probe, resulting in a higher wall shear stress for lower inlet flow velocities. An additional study was performed using a tracking particle model to understand the relevance of the individual particles behavior under different flow conditions. The results do not indicate any significant effect on the measurements inside the probe, although additional studies should be carried out in this topic by considering a population balance model for the oil droplets. / O presente trabalho tem como objetivo o estudo do comportamento do escoamento de uma emulsão do tipo óleo - em água que flui no interior de câmara/duto de medição e que tem como obstáculo em seu caminho uma sonda de um sensor óptico, sensor óptico este que deve avaliar em tempo real a estabilidade da emulsão onde está inserido. A emulsão é constituída por um fluido de corte para usinagem, com gotículas de óleo de diâmetro variando de 100 nanômetros para 100 micrometros. A sonda utilizada junto ao espectrômetro de luz UV- Vis está de acordo com o conceito proposto no projeto de pesquisa chamado EPM (Emulsion Process Monitor in Metalworking Fluid), realizado no âmbito do programa BRAGECRIM entre a Universidade de São Paulo e a Universidade de Bremen. Este estudo baseia-se na simulação numérica das interações entre a emulsão e o sistema de medição proposto, utilizando técnicas de Fluido Dinâmica Computacional (CFD), e tem por objetivo avaliar os efeitos da geometria da sonda, a sua posição em relação ao campo do escoamento, e propriedades do fluido, em especial as propriedades a serem medidas pelo espectrômetro. Tais efeitos estão correlacionadas com alterações na concentração de gotas e a segregação ou o tamanho das gotas dentro da câmara de medição da sonda óptica, o que pode causar mudanças nas leituras de intensidade de luz difusa. Efeitos de segregação devido a perturbações do escoamento em torno da sonda podem ser negligenciados, em condições normais de medição, com a sonda voltada para frente e sua área de entrada contra a corrente do escoamento. Com base nos resultados de simulação, mesmo que a sonda seja deslocada, o efeito sobre as medições ainda é insignificante. Medições foram efetuadas em laboratório e também foram realizadas medições in-situ utilizando um adaptador de medição acoplado diretamente no tubo de injeção de fluido de corte da máquina de perfuração, estes testes foram feitos para de validar os resultados obtidos por simulação, uma vez que não se observou qualquer efeito de segregação relacionada com o sistema de medição. Além disso, foi acrescentada ao estudo, a possibilidade de fixação de bactérias nas paredes de vidro internas da sonda e verificou-se que quando a velocidade de escoamento é suficientemente grande para produzir uma tensão de cisalhamento de cerca de 3-5 Pa a fixação de bacteriana pode ser evitado. Aproveitando os resultados do estudo para prevenção de contaminação por bactéria, mudanças na geometria da sonda foram propostas a fim de atingir uma condição isocinética para o escoamento ao redor e dentro da sonda, resultando em uma maior tensão de cisalhamento para baixas velocidades de fluxo de entrada. Por fim um estudo adicional foi realizado utilizando um modelo de rastreamento de partículas para compreender a relevância do comportamento individual de cada partícula no escoamento da emulsão. Os resultados não indicam qualquer efeito significativo sobre as medições no interior da sonda, embora estudos adicionais devem ser realizados neste tópico, considerando um modelo de balanço populacional para as gotículas de óleo.

Uma das dificuldades relacionadas com a exploração de petróleo é o transporte de óleo pesado, que devido a sua alta viscosidade, acarreta em uma elevada perda de carga no sistema. Para proporcionar economia de energia aplica-se o método do Core Annular Flow (CAF) onde é utilizado um escoamento anular bifásico em que a água escoa na periferia da tubulação para redução do gasto energético. O presente trabalho visou compreender e reproduzir este fenômeno, desenvolvendo simulações em CFD através do pacote comercial ANSYS FLUENT considerando o escoamento 3D, turbulento, isotérmico e incompressível para casos estacionários e transientes. A interface entre a água e o óleo foi adequadamente reproduzida em diferentes geometrias como tubo reto e com curva. O método LES para simulação de grandes escalas provou ser o melhor método de turbulência dentre os testados, como k-epsilon e modelo de tensores de Reynolds, de forma com que a interface fosse representada corretamente. O modelo para sistema multifásico adotado foi o Volume de Fluido (VOF), comparado com o comportamento experimental e com dados da literatura. Os fenômenos de swirl observados experimentalmente também foram reproduzidos de forma satisfatória. / One of the difficulties related to oil exploration is the transportation of heavy oil, which due to its high viscosity, causes a high pressure drop in the system. In order to provide energy savings, the Core Annular Flow (CAF) method applies where a two-phase annular flow occurs in which water flows at the periphery adjacent to the pipe to reduce energy expenditure. The present work aimed to understand and reproduce this phenomenon, developing CFD simulations through the commercial package ANSYS FLUENT considering flow as 3D, turbulent, isothermal and incompressible for stationary and transient cases. The interface between water and oil has been properly reproduced in different geometries such as straight pipe and pipe with a curve. The LES method for large scale simulation proved to be the best turbulence method among the tested, such as k-epsilon and Reynolds stress model, so that the interface was correctly represented. The model for the multiphase system adopted was the Volume of Fluid (VOF), compared to the experimental behavior and with data from the literature. The experimentally observed swirl phenomena were also reproduced satisfactorily.

Simulating the exact chaotic turbulent flow field about any geometry is a dilemma between accuracy and computational resources, which has been continuously studied for just over a hundred years. This thesis is a complete walk-through of the entire process utilized to approximate the flow ingested by a Sevik-type rotor based on solutions to the Reynolds Averaged Navier-Stokes equations (RANS). The Multiple Reference Frame fluid model is utilized by the code of ANSYS-FLUENT and results are validated by experimental wake data. Three open rotor configurations are studied including a uniform inflow and the rotor near a plate with and without a thick boundary layer. Furthermore, observations are made to determine the variation in velocity profiles of the ingested turbulent flow due to varying flow conditions. / by Felipe Ferreira Lachowski. / Thesis (M.S.C.S.)--Florida Atlantic University, 2013. / Includes bibliography. / Mode of access: World Wide Web. / System requirements: Adobe Reader.

This thesis presents an investigation on multiphase reactor hydrodynamics using magnetic resonance imaging (MRI). The study demonstrates experimental techniques by which computational and quasi-analytical fluid models may be validated. Three types of industrially-important multiphase reaction vessels are considered: a co-current upflow gas-liquid-solid bed, a co-current downward trickle bed (gas, liquid, solid), and a gas-solid fluidised bed. These reactors were selected as they commonly demonstrate local hydrodynamic anisotropy which affects the global performance of industrial units. MRI was used to obtain 2D velocity images of the gas and liquid phases in the packed beds, and of the gas and the solid phases in the fluidised bed. This study reports the first spatially resolved velocity measurements of both the gas and liquid phases in a co-current upflow bed, and the gas and solid phases of an isolated bubble in a fluidised bed. The experimental vessels were: 52 mm in diameter using 5 mm glass spheres in the upflow bed at 8 bara, 27 mm with 5 mm glass spheres in the trickle bed at 6.75 bara, and 52 mm using 1.2 mm poppy seeds as the fluidised particles at 8.5 bara. The experiments were conducted at a laboratory temperature of 25.0 ± 3.0 °C. In the upflow bed, time-averaged velocity images were acquired over a 2.5 h experimental time. This was done to capture the steady state behaviour of the vessel operating in the pulsing flow regime. The temporally-stable trickle flow state in the trickle bed was imaged over 15-100 minutes. In both packed beds, severe spatial anisotropy in the distribution of flow between pores was revealed. Furthermore, the data were used to determine classical design features such as catalyst wetting and liquid holdup which compared well with literature models. The trickle bed data were further analysed using a morphological algorithm which unambiguously identified the gas-liquid and liquid-solid interfaces. The interfacial flow fields were found to be similar to the bulk flow, with most voxels exhibiting static behaviour. The amount of interaction between the phases was found to be minimal, which is typical of the low interaction regime. A single bubble injection system was employed in the fluidised bed which allowed the injection of isolated bubbles into the incipiently fluidised bed. It also enabled the triggered acquisition of NMR data at precise time intervals. The bubble was found to be an indented ellipsoidal shape, which rose with atypical behaviour which caused it to collapse. Rise velocity was found to be consistent with theory, and the injected bubbles were sufficiently spatially reproducible to acquire 2D velocity images using single-point imaging. These velocity images showed flow behaviour characteristic of a 'fast' rising bubble, with a gas recirculation cloud 37 mm in diameter. The particle field was shown to have very high flow in the bubble wake, revealing the mechanism of bubble collapse. The flow data were compared to classical two-phase fluidisation theory, which revealed noteworthy differences in the division of flow between the particulate and bubbling regions.

Uma das dificuldades relacionadas com a exploração de petróleo é o transporte de óleo pesado, que devido a sua alta viscosidade, acarreta em uma elevada perda de carga no sistema. Para proporcionar economia de energia aplica-se o método do Core Annular Flow (CAF) onde é utilizado um escoamento anular bifásico em que a água escoa na periferia da tubulação para redução do gasto energético. O presente trabalho visou compreender e reproduzir este fenômeno, desenvolvendo simulações em CFD através do pacote comercial ANSYS FLUENT considerando o escoamento 3D, turbulento, isotérmico e incompressível para casos estacionários e transientes. A interface entre a água e o óleo foi adequadamente reproduzida em diferentes geometrias como tubo reto e com curva. O método LES para simulação de grandes escalas provou ser o melhor método de turbulência dentre os testados, como k-epsilon e modelo de tensores de Reynolds, de forma com que a interface fosse representada corretamente. O modelo para sistema multifásico adotado foi o Volume de Fluido (VOF), comparado com o comportamento experimental e com dados da literatura. Os fenômenos de swirl observados experimentalmente também foram reproduzidos de forma satisfatória. / One of the difficulties related to oil exploration is the transportation of heavy oil, which due to its high viscosity, causes a high pressure drop in the system. In order to provide energy savings, the Core Annular Flow (CAF) method applies where a two-phase annular flow occurs in which water flows at the periphery adjacent to the pipe to reduce energy expenditure. The present work aimed to understand and reproduce this phenomenon, developing CFD simulations through the commercial package ANSYS FLUENT considering flow as 3D, turbulent, isothermal and incompressible for stationary and transient cases. The interface between water and oil has been properly reproduced in different geometries such as straight pipe and pipe with a curve. The LES method for large scale simulation proved to be the best turbulence method among the tested, such as k-epsilon and Reynolds stress model, so that the interface was correctly represented. The model for the multiphase system adopted was the Volume of Fluid (VOF), compared to the experimental behavior and with data from the literature. The experimentally observed swirl phenomena were also reproduced satisfactorily.

Hydrodynamic vortex separators (HDVSs) are used in wastewater treatment to separate solids from wastewater. The aim of this research is to devise a CFD-based methodology that optimises their performance through modification of their design. A validation study is performed to assess whether OpenFOAM can be used to reliably model the flow of water in an HDVS. The results of the simulations are compared with experimental readings, showing a good fit when the appropriate boundary layer height and turbulence model are used. The continuous adjoint method is employed to derive the adjoint equations, associated with the drift flux equations used to model the flow of wastewater. They are specialised to the typical boundary conditions of ducted flows and are coded using OpenFOAM. An optimal design is found for boundary conditions, corresponding to typical values used in practice, and is shown to improve the performance of a simplified initial design by 40%. This optimal design is subsequently subjected to a different hydraulic loading rate and dispersed-phase volume fraction at the inlet, to assess the performance variation in these circumstances. Though the optimal design removes all the solids when the dispersed-phase fraction is reduced at the inlet, initial results suggest that the design is sensitive to hydraulic loading rate and further tests are recommended before drawing more explicit conclusions. This is the first time the adjoint drift flux equations have been derived. It is also the first time they have been coded and applied to an HDVS to optimise its performance. The methodology developed in this thesis could be applied to any device that separates solids from liquid or two immiscible liquids, in order to optimise its performance.

A new measurement technique called fiber optic reflectometer (FOR) was developed to investigate multiphase flows. The principle and setup of the FOR technique were introduced and applied to various experiments. Based on the coherently mixed signal between the Fresnel reflection off the fiber-liquid interface and the scattered signal off the object, such as a gas bubble, and a solid particle, this single probe technique is capable of simultaneously measuring the velocity of the object with a high accuracy and the phase of the fluid. In addition, bubble diameter, velocity, and void fraction were measured directly. By means of a simple modification of the FOR technique, solute concentration and refractive index change were measured with a greatly improved accuracy. This modified technique was used for measuring of a NaCl concentration in deionized water to validate a new normalization technique. In the second part of this thesis, a plunging breaking wave in deep water has been studied. Using the wave focusing method, a strong plunging breaker was generated with accuracy in the deep water condition in a two-dimensional wave tank. It was possible to describe the breaking process in detail using a high speed camera with a frame rate of 500 or 1000 fps. Four kinds of experimental techniques were employed or developed to investigate the plunging breaker. Bubble image velocimetry (BIV) and particle image velocimetry (PIV) were used to measure the velocity fields. The velocity fields of the highly aerated region were obtained from the BIV measurements. In addition, the modified PIV technique is capable of measuring the velocities in the entire flow field including the aerated region. Mean and turbulent properties were obtained by the ensemble average. The mean velocity, mean vorticity, and mean kinetic energy were examined over the entire flow field. In addition, the Reynolds stresses and turbulent kinetic energy were calculated with high temporal and spatial resolutions. Free surface elevation was obtained from wave gauge measurements. BIV and PIV images were also used to obtain the free surface elevation and the boundary of the aerated region for more accurate results. The FOR technique was used to obtain the void ratio at each splash-up region. Compressibility of the plunging breaker was considered. Mass flux, momentum flux, kinetic energy, and Reynolds stresses at each FOR station were recalculated using the void ratio obtained from the FOR measurements. All terms at the first splash-up region were highly overestimated more than 100 percent unless the void ratio was applied to the calculation of fluxes and energies. Compared with the fully developed first splash-up region, the overestimation at the second and third splash-up was less significant. However, most terms were overestimated by 20~30 percent when the void ratio was not considered.