Resistência ao escoamento em canais com vegetação submersa e rígida

Dissertação para a obtenção do grau de Mestre em Engenharia Civil - Perfil de Construção

A noninvasive method for measuring the velocity of diffuse hydrothermal flow by tracking moving refractive index anomalies

Diffuse flow velocimetry (DFV) is introduced as a new, noninvasive, optical technique for measuring the velocity of diffuse hydrothermal flow. The technique uses images of a motionless, random medium (e.g.,rocks) obtained through the lens of a moving refraction index anomaly (e.g., a hot upwelling). The method works in two stages. First, the changes in apparent background deformation are calculated using particle image velocimetry (PIV). The deformation vectors are determined by a cross correlation of pixel intensities across consecutive images. Second, the 2-D velocity field is calculated by cross correlating the deformation vectors between consecutive PIV calculations. The accuracy of the method is tested with laboratory and numerical experiments of a laminar, axisymmetric plume in fluids with both constant and temperaturedependent viscosity. Results show that average RMS errors are ∼5%–7% and are most accurate in regions of pervasive apparent background deformation which is commonly encountered in regions of diffuse hydrothermal flow. The method is applied to a 25 s video sequence of diffuse flow from a small fracture captured during the Bathyluck’09 cruise to the Lucky Strike hydrothermal field (September 2009). The velocities of the ∼10°C–15°C effluent reach ∼5.5 cm/s, in strong agreement with previous measurements of diffuse flow. DFV is found to be most accurate for approximately 2‐D flows where background objects have a small spatial scale, such as sand or gravel

Measurements of the bulk and interfacial velocity profiles in oscillating Newtonian and Maxwellian fluids

We present the dynamic velocity profiles of a Newtonian fluid (glycerol) and a viscoelastic Maxwell fluid (CPyCl-NaSal in water) driven by an oscillating pressure gradient in a vertical cylindrical pipe. The frequency range explored has been chosen to include the first three resonance peaks of the dynamic permeability of the viscoelastic-fluid¿pipe system. Three different optical measurement techniques have been employed. Laser Doppler anemometry has been used to measure the magnitude of the velocity at the center of the liquid column. Particle image velocimetry and optical deflectometry are used to determine the velocity profiles at the bulk of the liquid column and at the liquid-air interface respectively. The velocity measurements in the bulk are in good agreement with the theoretical predictions of a linear theory. The results, however, show dramatic differences in the dynamic behavior of Newtonian and viscoelastic fluids, and demonstrate the importance of resonance phenomena in viscoelastic fluid flows, biofluids in particular, in confined geometries.

Experimental investigation of the initial regime in fingering electrodeposition: Dispersion relation and velocity measurements.

Recently a fingering morphology, resembling the hydrodynamic Saffman-Taylor instability, was identified in the quasi-two-dimensional electrodeposition of copper. We present here measurements of the dispersion relation of the growing front. The instability is accompanied by gravity-driven convection rolls at the electrodes, which are examined using particle image velocimetry. While at the anode the theory presented by Chazalviel et al. [J. Electroanal. Chem. 407, 61 (1996)] describes the convection roll, the flow field at the cathode is more complicated because of the growing deposit. In particular, the analysis of the orientation of the velocity vectors reveals some lag of the development of the convection roll compared to the finger envelope.

Optimization of Snow Drifting Mitigation and Control Methods for Iowa Conditions Part II, TR-626, 2015

Blowing and drifting of snow is a major concern for transportation efficiency and road safety in regions where their development is common. One common way to mitigate snow drift on roadways is to install plastic snow fences. Correct design of snow fences is critical for road safety and maintaining the roads open during winter in the US Midwest and other states affected by large snow events during the winter season and to maintain costs related to accumulation of snow on the roads and repair of roads to minimum levels. Of critical importance for road safety is the protection against snow drifting in regions with narrow rights of way, where standard fences cannot be deployed at the recommended distance from the road. Designing snow fences requires sound engineering judgment and a thorough evaluation of the potential for snow blowing and drifting at the construction site. The evaluation includes site-specific design parameters typically obtained with semi-empirical relations characterizing the local transport conditions. Among the critical parameters involved in fence design and assessment of their post-construction efficiency is the quantification of the snow accumulation at fence sites. The present study proposes a joint experimental and numerical approach to monitor snow deposits around snow fences, quantitatively estimate snow deposits in the field, asses the efficiency and improve the design of snow fences. Snow deposit profiles were mapped using GPS based real-time kinematic surveys (RTK) conducted at the monitored field site during and after snow storms. The monitored site allowed testing different snow fence designs under close to identical conditions over four winter seasons. The study also discusses the detailed monitoring system and analysis of weather forecast and meteorological conditions at the monitored sites. A main goal of the present study was to assess the performance of lightweight plastic snow fences with a lower porosity than the typical 50% porosity used in standard designs of such fences. The field data collected during the first winter was used to identify the best design for snow fences with a porosity of 50%. Flow fields obtained from numerical simulations showed that the fence design that worked the best during the first winter induced the formation of an elongated area of small velocity magnitude close to the ground. This information was used to identify other candidates for optimum design of fences with a lower porosity. Two of the designs with a fence porosity of 30% that were found to perform well based on results of numerical simulations were tested in the field during the second winter along with the best performing design for fences with a porosity of 50%. Field data showed that the length of the snow deposit away from the fence was reduced by about 30% for the two proposed lower-porosity (30%) fence designs compared to the best design identified for fences with a porosity of 50%. Moreover, one of the lower-porosity designs tested in the field showed no significant snow deposition within the bottom gap region beneath the fence. Thus, a major outcome of this study is to recommend using plastic snow fences with a porosity of 30%. It is expected that this lower-porosity design will continue to work well for even more severe snow events or for successive snow events occurring during the same winter. The approach advocated in the present study allowed making general recommendations for optimizing the design of lower-porosity plastic snow fences. This approach can be extended to improve the design of other types of snow fences. Some preliminary work for living snow fences is also discussed. Another major contribution of this study is to propose, develop protocols and test a novel technique based on close range photogrammetry (CRP) to quantify the snow deposits trapped snow fences. As image data can be acquired continuously, the time evolution of the volume of snow retained by a snow fence during a storm or during a whole winter season can, in principle, be obtained. Moreover, CRP is a non-intrusive method that eliminates the need to perform man-made measurements during the storms, which are difficult and sometimes dangerous to perform. Presently, there is lots of empiricism in the design of snow fences due to lack of data on fence storage capacity on how snow deposits change with the fence design and snow storm characteristics and in the estimation of the main parameters used by the state DOTs to design snow fences at a given site. The availability of such information from CRP measurements should provide critical data for the evaluation of the performance of a certain snow fence design that is tested by the IDOT. As part of the present study, the novel CRP method is tested at several sites. The present study also discusses some attempts and preliminary work to determine the snow relocation coefficient which is one of the main variables that has to be estimated by IDOT engineers when using the standard snow fence design software (Snow Drift Profiler, Tabler, 2006). Our analysis showed that standard empirical formulas did not produce reasonable values when applied at the Iowa test sites monitored as part of the present study and that simple methods to estimate this variable are not reliable. The present study makes recommendations for the development of a new methodology based on Large Scale Particle Image Velocimetry that can directly measure the snow drift fluxes and the amount of snow relocated by the fence.

Development of Stage-Discharge Relations for Ungaged Bridge Waterways in Western Iowa, TR-567, 2010

Stage-discharge relations constitute a viable, alternative technique for estimating accurately flow for ungaged sites. In this research, we have utilized pressure transducers and Large Scale Particle Image Velocimetry techniques to develop stage-discharge relations at eleven sites in the Hungry Canyon Area (HCA) of southwestern Iowa under different hydrologic conditions. We have employed these data to calibrate and verify an established hydrologic model and then we have used this model to provide a stage discharge relation for different hydrologic conditions (i.e. rating curves). The benefits of the project are numerous including that the discharge data will be used for a number of purposes, including operational decision making in the HCA about the design of water-control and conveyance structures, input for hydraulic and hydrologic models, and calculation of sediment and other water quality constituents transport and “loads”, and for decision making. This project has also pointed out the difficulties in measuring flows in ungaged streams with ice jams, steep banks, erodible beds, and floating debris.

Measurement of a Structured Backflow in an Open Small Channel Induced by Surface-Tension Gradients

We present experiments in which the laterally confined flow of a surfactant film driven by controlled surface tension gradients causes the subtended liquid layer to self-organize into an inner upstream microduct surrounded by the downstream flow. The anomalous interfacial flow profiles and the concomitant backflow are a result of the feedback between two-dimensional and three-dimensional microfluidics realized during flow in open microchannels. Bulk and surface particle image velocimetry data combined with an interfacial hydrodynamics model explain the dependence of the observed phenomena on channel geometry.

Separation of Gas Bubbles from Circulation Waters and Pulp Suspensions in Open Channel Flow

The objective of this thesis was to study the removal of gases from paper mill circulation waters experimentally and to provide data for CFD modeling. Flow and bubble size measurements were carried out in a laboratory scale open gas separation channel. Particle Image Velocimetry (PIV) technique was used to measure the gas and liquid flow fields, while bubble size measurements were conducted using digital imaging technique with back light illumination. Samples of paper machine waters as well as a model solution were used for the experiments. The PIV results show that the gas bubbles near the feed position have the tendency to escape from the circulation channel at a faster rate than those bubbles which are further away from the feed position. This was due to an increased rate of bubble coalescence as a result of the relatively larger bubbles near the feed position. Moreover, a close similarity between the measured slip velocities of the paper mill waters and that of literature values was obtained. It was found that due to dilution of paper mill waters, the observed average bubble size was considerably large as compared to the average bubble sizes in real industrial pulp suspension and circulation waters. Among the studied solutions, the model solution has the highest average drag coefficient value due to its relatively high viscosity. The results were compared to a 2D steady sate CFD simulation model. A standard Euler-Euler k-ε turbulence model was used in the simulations. The channel free surface was modeled as a degassing boundary. From the drag models used in the simulations, the Grace drag model gave velocity fields closest to the experimental values. In general, the results obtained from experiments and CFD simulations are in good qualitative agreement.

Modelling of ultrasound assisted mixing in Newtonian and non-Newtonian liquids

The purpose of this work is to obtain a better understanding of behaviour of possible ultrasound appliance on fluid media mixing. The research is done in the regard to Newtonian and non-Newtonian fluids. The process of ultrasound appliance on liquids is modelled in COMSOL Multiphysics software. The influence of ultrasound using is introduced as waveform equation. Turbulence modelling is fulfilled by the k-ε model in Newtonian fluid. The modeling of ultrasound assisted mixing in non-Newtonian fluids is based on the power law. To verify modelling results two practical methods are used: Particle Image Velocimetry and measurements of mixing time. Particle Image Velocimetry allows capturing of velocity flow field continuously and presents detailed depiction of liquid dynamics. The second way of verification is the comparison of mixing time of homogeneity. Experimentally achievement of mixing time is done by conductivity measurements. In modelling part mixing time is achieved by special module of COMSOL Multiphysics – the transport of diluted species. Both practical and modelling parts show similar radial mechanism of fluid flow under ultrasound appliance – from the horn tip fluid moves to the bottom and along the walls goes back. Velocity profiles are similar in modelling and experimental part in the case of Newtonian fluid. In the case of non-Newtonian fluid velocity profiles do not agree. The development track of ultrasound-assisted mixing modelling is presented in the thesis.

Höyryturbiinin ulosvirtauskanaviston kokeellinen tutkimus pienoismallilla

Diplomityössä tutkittiin höyryturbiinin ulosvirtauskanavistojen kokeellisia tutkimusmenetelmiä ja suoritettiin käytännön mittauksia Fortum Oyj:n Loviisan ydinvoimalaitoksen höyryturbiinien huuvan pienoismallilla. Kirjallisuusselvityksen perusteella todettiin, että pienoismallitutkimuksella on ollut keskeinen asema ulosvirtauskanavistojen suunnittelussa. Kokeellisten menetelmien perusongelmana on höyryturbiinin ulosvirtausolosuhteiden jäljitteleminen. Käytetyt mittausmenetelmät perustuvat pääosin tavanomaisiin paine- ja nopeusmittauksiin. Lisäainepartikkeleihin ja laser-valaisuun perustuva PIV (particle image velocimetry) todettiin lupaavaksi menetelmäksi ulosvirtauskanavistojen tutkimuksen saralla. Työn käytännön osuudessa tehtiin mittauksia mittasuhteessa 1:8 rakennetulle höyryturbiinin huuvan pienoismallille. Mittauksilla tutkittiin virtausta mallin sisääntulo- ja ulostulotasoissa. Lisäksi mitattiin staattisen paineen jakauma huuvan sisällä. Kokonaispainetta mittaava kiel-putki todettiin käytännölliseksi työkaluksi huuvan virtauskentän tutkimuksessa. Tuloksista käy hyvin ilmi huuvan ulostuloon syntyvien pyörteiden muodostuminen ja ulostulon epätasainen nopeusjakauma. Staattinen paine huuvan sisällä havaittiin epätasaisesti jakautuneeksi. Ulostulotason ja staattisen paineen mittauksilla saadut tulokset sopivat hyvin yhteen kirjallisuudesta löytyvien tutkimustulosten kanssa ja tukevat Loviisan ulosvirtauskanavistosta aiemmin tehtyjä CFD-simulointeja.

Monitorointimenetelmät membraanin likaantumisen seurantaan suodatuksen aikana

Membraani on ohut kalvo, jossa on pieniä nanomittakaavan reikiä, jotka erottavat partikkelit ja liuenneet yhdisteet liuoksesta. Membraanisuodatuksen käyttö on lisääntynyt merkittävästi vedenpuhdistuksessa, johtuen lisääntyneestä puhtaan veden tarpeesta ja tiukentuneista ympäristövaatimuksista. Tässä työssä esitellään reaaliaikaisia mittausmenetelmiä membraanin likaantumisen seurantaan. Esiteltyjä menetelmiä ovat suora havainnointi pinnan läpi, lasertriangulometria, varjoanalyysi, taittokykymittaus, kuvakatkaisu-menetelmä, partikkelin nopeusmääritys, radioisotooppinen merkintä ja ydinmagneettinen resonanssispektrometria. Mittausmenetelmien avulla likakerroksen paksuutta ja sen leviämistä on mahdollista seurata reaaliaikaisesti. Mittausmenetelmien soveltuvuus olemassa oleviin prosesseihin on vielä epävarmaa. Suurin osa menetelmistä on rajoittunut tiettyyn membraanin materiaaliin, tietynlaiseen membraanisuodatusprosessin rakenteeseen tai tiettyihin olosuhteisiin. Vallitsevien prosessiolosuhteiden lisäksi mittausanturin tulisi kestää myös puhdistusolosuhteet. Lisätutkimuksia tarvitaan, jotta voidaan löytää toimiva laitekokonaisuus tarvittavan tiedon tuottamiseen.

Experimentelle Untersuchung der rotierenden Instabilität im Axialverdichtergitter

Die Strömung in Turboverdichtern ist geprägt von instationären Effekten. Ziel dieser experimentellen Arbeit ist die zeitlich aufgelöste Darstellung eines dieser Effekte: der rotierenden Instabilität. Hierbei handelt es sich um ein Phänomen, das dem Auftreten von rotierender Ablösung vorausgeht und bereits im stabilen Arbeitsbereich auftritt. Das bessere Verständnis der rotierenden Instabilität bietet somit auch einen Ansatz zur Klärung der Ursachen, die zu rotierender Ablösung führen. Zur Durchführung der zeitlich aufgelösten Untersuchung wird im Rahmen dieser Arbeit eine Messmethodik vorgestellt, welche erstmals bei der Strömungsuntersuchung in Turbomaschinen Anwendung findet. Die Methodik basiert auf simultanen Messungen mittels konventionellem PIV und zeitlich hochaufgelöster Hitzdrahtanemometrie. Über die Frequenzanalyse des Hitzdrahtsignals ist es im Postprocessing möglich, die einzelnen PIV-Bilder gemäß ihres Phasenwinkel zu ordnen. Dieser statistische Ansatz ermöglicht die Darstellung und Analyse der zeitlichen Entwicklung von rotierender Instabilität über einer Schwingungsperiode.

Dynamics of convectively driven banded jets in the laboratory

The banded organization of clouds and zonal winds in the atmospheres of the outer planets has long fascinated observers. Several recent studies in the theory and idealized modeling of geostrophic turbulence have suggested possible explanations for the emergence of such organized patterns, typically involving highly anisotropic exchanges of kinetic energy and vorticity within the dissipationless inertial ranges of turbulent flows dominated (at least at large scales) by ensembles of propagating Rossby waves. The results from an attempt to reproduce such conditions in the laboratory are presented here. Achievement of a distinct inertial range turns out to require an experiment on the largest feasible scale. Deep, rotating convection on small horizontal scales was induced by gently and continuously spraying dense, salty water onto the free surface of the 13-m-diameter cylindrical tank on the Coriolis platform in Grenoble, France. A “planetary vorticity gradient” or “β effect” was obtained by use of a conically sloping bottom and the whole tank rotated at angular speeds up to 0.15 rad s−1. Over a period of several hours, a highly barotropic, zonally banded large-scale flow pattern was seen to emerge with up to 5–6 narrow, alternating, zonally aligned jets across the tank, indicating the development of an anisotropic field of geostrophic turbulence. Using particle image velocimetry (PIV) techniques, zonal jets are shown to have arisen from nonlinear interactions between barotropic eddies on a scale comparable to either a Rhines or “frictional” wavelength, which scales roughly as (β/Urms)−1/2. This resulted in an anisotropic kinetic energy spectrum with a significantly steeper slope with wavenumber k for the zonal flow than for the nonzonal eddies, which largely follows the classical Kolmogorov k−5/3 inertial range. Potential vorticity fields show evidence of Rossby wave breaking and the presence of a “hyperstaircase” with radius, indicating instantaneous flows that are supercritical with respect to the Rayleigh–Kuo instability criterion and in a state of “barotropic adjustment.” The implications of these results are discussed in light of zonal jets observed in planetary atmospheres and, most recently, in the terrestrial oceans.

Ferramenta computacional para a medição de campos de velocidade utilizando processamento digital de imagens

A proposta deste trabalho, consiste na elaboração de uma ferramenta computacional para a medição de campos de velocidades em escoamentos com baixas velocidades (< 0,5 m/s) utilizando o processamento digital de imagens. Ao longo dos anos, inúmeras técnicas foram desenvolvidas com este objetivo. Para cada tipo de aplicação, uma técnica se aplica com maior ou menor eficiência do que outras. Para o caso de estudos em fluídos transparentes, onde o escoamento pode ser visualizado, técnicas que utilizam processamento digital de imagens vêm ganhando um grande impulso tecnológico nos últimos anos. Este impulso, é devido a fatores como: câmaras vídeo filmadoras de última geração, dispositivos de aquisição de imagens e componentes de processamento e armazenamento de dados cada vez mais poderosos. Neste contexto, está a velocimetria por processamento de imagens de partículas cuja sigla é PIV (particle image velocimetry). Existem várias formas de se implementar um sistema do tipo PIV. As variantes dependem, basicamente, do equipamento utilizado. Para sua implementação é necessário, inicialmente, um sistema de iluminação que incide em partículas traçadoras adicionadas ao fluido em estudo. Após, as partículas em movimento são filmadas ou fotografadas e suas imagens adquiridas por um computador através de dispositivos de captura de imagens. As imagens das partículas são então processadas, para a obtenção dos vetores velocidade. Existem diferentes formas de processamento para a obtenção das velocidades. Para o trabalho em questão, devido às características dos equipamentos disponíveis, optou-se por uma metodologia de determinação da trajetória de partículas individuais, que, apesar de limitada em termos de módulo de velocidade, pode ser aplicada a muitos escoamentos reais sob condições controladas Para validar a ferramenta computacional desenvolvida, imagens ideais de partículas foram simuladas como se estivessem em escoamento, através do deslocamento conhecido de vários pixels. Seguindo o objetivo de validação, foi utilizada ainda uma imagem real de partículas, obtida com o auxílio de um plano de iluminação de luz coerente (LASER) e câmaras de vídeo tipo CCD. O programa desenvolvido foi aplicado em situações de escoamento real e os resultados obtidos foram satisfatórios dentro da escala de velocidades inicialmente presumida.

Caracterização de motor aeronaútico utilizando misturas de gasolina de aviação e etanol: aspectos fluidodinâmicos, termodinâmicos e ecológicos

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)