Dispersion of passive particles by a quasi-two-dimensional turbulent flow

Using the experimental data of Paret and Tabeling [Phys. Rev. Lett. 79, 4162 (1997)] we consider in detail the dispersion of particle pairs by a two-dimensional turbulent flow and its relation to the kinematic properties of the velocity field. We show that the mean square separation of a pair of particles is governed by rather rare, extreme events and that the majority of initially close pairs are not dispersed by the flow. Another manifestation of the same effect is the fact that the dispersion of an initially dense cluster is not the result of homogeneously spreading the particles within the whole system. Instead it proceeds through a splitting into smaller but also dense clusters. The statistical nature of this effect is discussed.

Langevin approach to generate synthetic turbulent flows

We present an analytical scheme, easily implemented numerically, to generate synthetic Gaussian turbulent flows by using a linear Langevin equation, where the noise term acts as a stochastic stirring force. The characteristic parameters of the velocity field are well introduced, in particular the kinematic viscosity and the spectrum of energy. As an application, the diffusion of a passive scalar is studied for two different energy spectra. Numerical results are compared favorably with analytical calculations.

A multilevel multiscale finite volume method

The Multiscale Finite Volume (MsFV) method has been developed to efficiently solve reservoir-scale problems while conserving fine-scale details. The method employs two grid levels: a fine grid and a coarse grid. The latter is used to calculate a coarse solution to the original problem, which is interpolated to the fine mesh. The coarse system is constructed from the fine-scale problem using restriction and prolongation operators that are obtained by introducing appropriate localization assumptions. Through a successive reconstruction step, the MsFV method is able to provide an approximate, but fully conservative fine-scale velocity field. For very large problems (e.g. one billion cell model), a two-level algorithm can remain computational expensive. Depending on the upscaling factor, the computational expense comes either from the costs associated with the solution of the coarse problem or from the construction of the local interpolators (basis functions). To ensure numerical efficiency in the former case, the MsFV concept can be reapplied to the coarse problem, leading to a new, coarser level of discretization. One challenge in the use of a multilevel MsFV technique is to find an efficient reconstruction step to obtain a conservative fine-scale velocity field. In this work, we introduce a three-level Multiscale Finite Volume method (MlMsFV) and give a detailed description of the reconstruction step. Complexity analyses of the original MsFV method and the new MlMsFV method are discussed, and their performances in terms of accuracy and efficiency are compared.

A framework for hybrid simulations of two-phase flow in porous media

Rough a global coarse problem. Although these techniques are usually employed for problems in which the fine-scale processes are described by Darcy's law, they can also be applied to pore-scale simulations and used as a mathematical framework for hybrid methods that couples a Darcy and pore scales. In this work, we consider a pore-scale description of fine-scale processes. The Navier-Stokes equations are numerically solved in the pore geometry to compute the velocity field and obtain generalized permeabilities. In the case of two-phase flow, the dynamics of the phase interface is described by the volume of fluid method with the continuum surface force model. The MsFV method is employed to construct an algorithm that couples a Darcy macro-scale description with a pore-scale description at the fine scale. The hybrid simulations results presented are in good agreement with the fine-scale reference solutions. As the reconstruction of the fine-scale details can be done adaptively, the presented method offers a flexible framework for hybrid modeling.

New horizons for black holes and branes

We initiate a systematic scan of the landscape of black holes in any spacetime dimension using the recently proposed blackfold effective worldvolume theory. We focus primarily on asymptotically flat stationary vacuum solutions, where we uncover large classes of new black holes. These include helical black strings and black rings, black odd-spheres, for which the horizon is a product of a large and a small sphere, and non-uniform black cylinders. More exotic possibilities are also outlined. The blackfold description recovers correctly the ultraspinning Myers-Perry black holes as ellipsoidal even-ball configurations where the velocity field approaches the speed of light at the boundary of the ball. Helical black ring solutions provide the first instance of asymptotically flat black holes in more than four dimensions with a single spatial U(1) isometry. They also imply infinite rational non-uniqueness in ultraspinning regimes, where they maximize the entropy among all stationary single-horizon solutions. Moreover, static blackfolds are possible with the geometry of minimal surfaces. The absence of compact embedded minimal surfaces in Euclidean space is consistent with the uniqueness theorem of static black holes

Implementation of population balances in simulation of gas-liquid reactors using a CFD-PBM coupled approach

In bubbly flow simulations, bubble size distribution is an important factor in determination of hydrodynamics. Beside hydrodynamics, it is crucial in the prediction of interfacial area available for mass transfer and in the prediction of reaction rate in gas-liquid reactors such as bubble columns. Solution of population balance equations is a method which can help to model the size distribution by considering continuous bubble coalescence and breakage. Therefore, in Computational Fluid Dynamic simulations it is necessary to couple CFD and Population Balance Model (CFD-PBM) to get reliable distribution. In the current work a CFD-PBM coupled model is implemented as FORTRAN subroutines in ANSYS CFX 10 and it has been tested for bubbly flow. This model uses the idea of Multi Phase Multi Size Group approach which was previously presented by Sha et al. (2006) [18]. The current CFD-PBM coupled method considers inhomogeneous flow field for different bubble size groups in the Eulerian multi-dispersed phase systems. Considering different velocity field for bubbles can give the advantageof more accurate solution of hydrodynamics. It is also an improved method for prediction of bubble size distribution in multiphase flow compared to available commercial packages.

Numerical and ExperimentalModelling of Gas Flow and Heat Transfer in the Air Gap of An Electric Machine

This work deals with the cooling of high-speed electric machines, such as motors and generators, through an air gap. It consists of numerical and experimental modelling of gas flow and heat transfer in an annular channel. Velocity and temperature profiles are modelled in the air gap of a high-speed testmachine. Local and mean heat transfer coefficients and total friction coefficients are attained for a smooth rotor-stator combination at a large velocity range. The aim is to solve the heat transfer numerically and experimentally. The FINFLO software, developed at Helsinki University of Technology, has been used in the flow solution, and the commercial IGG and Field view programs for the grid generation and post processing. The annular channel is discretized as a sector mesh. Calculation is performed with constant mass flow rate on six rotational speeds. The effect of turbulence is calculated using three turbulence models. The friction coefficient and velocity factor are attained via total friction power. The first part of experimental section consists of finding the proper sensors and calibrating them in a straight pipe. After preliminary tests, a RdF-sensor is glued on the walls of stator and rotor surfaces. Telemetry is needed to be able to measure the heat transfer coefficients at the rotor. The mean heat transfer coefficients are measured in a test machine on four cooling air mass flow rates at a wide Couette Reynolds number range. The calculated values concerning the friction and heat transfer coefficients are compared with measured and semi-empirical data. Heat is transferred from the hotter stator and rotor surfaces to the coolerair flow in the air gap, not from the rotor to the stator via the air gap, althought the stator temperature is lower than the rotor temperature. The calculatedfriction coefficients fits well with the semi-empirical equations and precedingmeasurements. On constant mass flow rate the rotor heat transfer coefficient attains a saturation point at a higher rotational speed, while the heat transfer coefficient of the stator grows uniformly. The magnitudes of the heat transfer coefficients are almost constant with different turbulence models. The calibrationof sensors in a straight pipe is only an advisory step in the selection process. Telemetry is tested in the pipe conditions and compared to the same measurements with a plain sensor. The magnitudes of the measured data and the data from the semi-empirical equation are higher for the heat transfer coefficients than thenumerical data considered on the velocity range. Friction and heat transfer coefficients are presented in a large velocity range in the report. The goals are reached acceptably using numerical and experimental research. The next challenge is to achieve results for grooved stator-rotor combinations. The work contains also results for an air gap with a grooved stator with 36 slots. The velocity field by the numerical method does not match in every respect the estimated flow mode. The absence of secondary Taylor vortices is evident when using time averagednumerical simulation.

Kuumalanka-anemometri virtausmittauksissa

Tässä työssä selostetaan kuumalanka-anemometrin käyttö virtausmittauksissa. Kuumalanka-anemometrilla saadaan mitattua virtausnopeuden ja -suunnan lisäksi nopeusheilahteluja. Mittaustaajuus on tyypillisesti useita kymmeniä tuhansia mittauksia sekunnissa ja signaali on jatkuva. Nykytekniikalla pystytään helposti tallentamaan mittauslaitteistolta saatu viesti tietokoneelle ja muuntamaan se nopeudeksi. Hetkellisten nopeuksien avulla voidaan laskea turbulenttisen virtauksen ominaisuuksia, kuten turbulenssin intensiteetti ja spektri. Kuumalanka-anemometrissa lämmitetään sähköisesti ohutta lankaa, joka on mitattavassa virtauksessa. Langan sähköteho on suunnilleen yhtäsuuri kuin langasta konvektiolla siirtyvä lämpöteho. Tällöin on teoreettisesti mahdollista laskea virtausnopeus lämpötehosta lämmönsiirtokorrelaatioilla. Käytännössä laitteisto joudutaan kuitenkin erikseen kalibroimaan, mutta sähkötehon teoreettista riippuvuutta konvektiosta käytetään hyväksi. Kuumalangan lämmitettävä osuus on tyypillisesti halkaisijaltaan 5 µm ja pituudeltaan noin 1 mm. Sitä käytetään pääasiassa kaasuvirtausten mittaamiseen ja valtaosassa mittauksissa virtausaineena on ilma. Kuumalanka voi olla toteutettu kuumakalvotekniikalla, jossa halkaisijaltaan noin 50 - 70 µm paksuinen kuitu on päällystetty ohuella sähköä johtavalla kalvolla. Kuumakalvoanturin ei tarvitse olla muodoltaan sylinterimäinen, se voi olla mm. kartiomainen tai kiilamainen. Erikoispäällystetyllä kuumakalvoanturilla on mahdollista mitata myös nestevirtauksia. Mitattaessa kaasuvirtauksia kuumakalvon etuna on selvästi parempi kestävyys verrattuna kuumalankaan. Nimitystä kuumalanka-anemometri käytetään yleisesti molemmista anturityypeistä Tämän työn alussa käsitellään sylinterin yli tapahtuvaan virtaukseen liittyvää virtausmekaniikkaa ja lämmönsiirtoa. Anemometrin sähköinen osa, laitteisto ja sen kalibrointi käydään läpi. Langan suuntariippuvuuden laskentaan esitetään tarvittavat yhtälöt. Työssä esitellään kolme laitteistolla tehtyä perusmittausta: anturin kohtauskulman muuttaminen, pyörähdyssymmerisen suihkun nopeuskenttä ja tuulitunnelin rajakerros. Lisäksi esitellään yksi käytännöllinen ja vaativampi mittaus, jossa on mitattu nopeusprofiili radiaalikompressorin diffuusorin loppuosassa.

The Optimisation of the Draft Tube of a Hydro Power Plant at Off-Design Conditions

Several possible methods of increasing the efficiency and power of hydro power plants by improving the flow passages are investigated in this stydy. The theoretical background of diffuser design and its application to the optimisation of hydraulic turbine draft tubes is presented in the first part of this study. Several draft tube modernisation projects that have been carried out recently are discussed. Also, a method of increasing the efficiency of the draft tube by injecting a high velocity jet into the boundary layer is presented. Methods of increasing the head of a hydro power plant by using an ejector or a jet pump are discussed in the second part of this work. The theoretical principles of various ejector and jet pump types are presented and four different methods of calculating them are examined in more detail. A self-made computer code is used to calculate the gain in the head for two example power plants. Suitable ejector installations for the example plants are also discussed. The efficiency of the ejector power was found to be in the range 6 - 15 % for conventional head increasers, and 30 % for the jet pump at its optimum operating point. In practice, it is impossible to install an optimised jet pump with a 30 % efficiency into the draft tube as this would considerabely reduce the efficiency of the draft tube at normal operating conditions. This demonstrates, however, the potential for improvement which lies in conventional head increaser technology. This study is based on previous publications and on published test results. No actual laboratory measurements were made for this study. Certain aspects of modelling the flow in the draft tube using computational fluid dynamics are discussed in the final part of this work. The draft tube inlet velocity field is a vital boundary condition for such a calculation. Several previously measured velocity fields that have successfully been utilised in such flow calculations are presented herein.

Present-day deformation of the Pyrenees revealed by GPS surveying and earthquake focal mechanisms until 2011

The Pyrenean mountain range is a slowly deforming belt with continuous and moderate seismic activity. To quantify its deformation field, we present the velocity field estimated from a GPS survey of the Pyrenees spanning 18 yr. The PotSis and ResPyr networks, including a total of 85 GPS sites, were installed and first measured in 1992 and 1995 1997, respectively, and remeasured in 2008 and 2010. We obtain a deformation field with velocities less than 1 mm yr−1 across the range. The estimated velocities for individual stations do not differ significantly from zero with 95 per cent confidence. Even so, we estimate a maximum extensional horizontal strain rate of 2.0 ± 1.7 nanostrain per year in a N S direction in the western part of the range. We do not interpret the vertical displacements due to their large uncertainties. In order to compare the horizontal strain rates with the seismic activity, we analyse a set of 194 focal mechanisms using three methods: (i) the 'r' factor relating their P and T axes, (ii) the stress tensors obtained by fault slip inversion and (iii) the strain-rate tensors. Stress and strain-rate tensors are estimated for: (i) the whole data set, (ii) the eastern and western parts of the range separately, and (iii) eight zones, which are defined based on the seismicity and the tectonic patterns of the Pyrenees. Each of these analyses reveals a lateral variation of the deformation style from compression and extension in the east to extension and strike-slip in the west of the range. Although the horizontal components of the strain-rate tensors estimated from the seismic data are slightly smaller in magnitude than those computed from the GPS velocity field, they are consistent within the 2σ uncertainties. Furthermore, the orientations of their principal axes agree with the mapped active faults.

Kuumailmapolttolaitteiston numeerinen virtausmallinnus

Tämä diplomityö perustuu Lappeenrannan teknillisen yliopiston Uusiutuvien energiajärjestelmien laboratorion koelaitteistoon, jolla tutkitaan voimakkaan savukaasunkierrätyksen ja kuumailmapolton soveltuvuutta pienen kokoluokan energiantuotantoprosesseihin. Työn teoriaosassa esitellään tavanomaisesta palamisesta eroavaa kuumailmapolttoa ja tarkastellaan sen ominaisuuksia. Myös työssä käytetyn tutkimusmenetelmän, numeerisen virtauslaskennan, periaatteita ja ominaisuuksia tarkastellaan. Työssä tutkitaan numeerisella virtausmallinnuksella kuumailmapolttolaitteiston virtauskentän käyttäytymistä, kun takaisin tulipesään kierrätettävän savukaasun määrä sekä tulipesän lämpöhäviöiden suuruus vaihtelevat. Virtauskentän tarkastelu on tärkeää, sillä palamisilman ja kierrätetyn savukaasun täytyy sekoittua kuumailmapolton aikaansaamiseksi. Työn virtausmallinnus suoritettiin Finflo-virtausratkaisijalla kaksiulotteisena palamisreaktioita mallintamatta. Vaikka työssä käytetyt mallit olivat kaksiulotteisia ja niissä käytettiin yksinkertaistuksia, virtausten käyttäytymisestä tulipesässä saatiin olennaista tietoa, jota voidaan mahdollisesti hyödyntää jatkotutkimuksissa.

Utilization of Particle Image velocimetry in PPOOLEX condensation experiments

Particle Image Velocimetry, PIV, is an optical measuring technique to obtain velocity information of a flow in interest. With PIV it is possible to achieve two or three dimensional velocity vector fields from a measurement area instead of a single point in a flow. Measured flow can be either in liquid or in gas form. PIV is nowadays widely applied to flow field studies. The need for PIV is to obtain validation data for Computational Fluid Dynamics calculation programs that has been used to model blow down experiments in PPOOLEX test facility in the Lappeenranta University of Technology. In this thesis PIV and its theoretical background are presented. All the subsystems that can be considered to be part of a PIV system are presented as well with detail. Emphasis is also put to the mathematics behind the image evaluation. The work also included selection and successful testing of a PIV system, as well as the planning of the installation to the PPOOLEX facility. Already in the preliminary testing PIV was found to be good addition to the measuring equipment for Nuclear Safety Research Unit of LUT. The installation to PPOOLEX facility was successful even though there were many restrictions considering it. All parts of the PIV system worked and they were found out to be appropriate for the planned use. Results and observations presented in this thesis are a good background to further PIV use.

Instabilities in electrochemical systems with a rotating disc electrode

Polarization curves experimentally obtained in the electro-dissolution of iron in a 1 M H2SO4 solution using a rotating disc as the working electrode present a current instability region within the range of applied voltage in which the current is controlled by mass transport in the electrolyte. According to the literature (Barcia et. al., 1992) the electro-dissolution process leads to the existence of a viscosity gradient in the interface metal-solution, which leads to a velocity field quantitatively different form the one developed in uniform viscosity conditions and may affect the stability of the hydrodynamic field. The purpose of this work is to investigate whether a steady viscosity profile, depending on the distance to the electrode surface, affects the stability properties of the classic velocity field near a rotating disc. Two classes of perturbations are considered: perturbations monotonically varying along the radial direction, and perturbations periodically modulated along the radial direction. The results show that the hydrodynamic field is always stable with respect to the first class of perturbations and that the neutral stability curves are modified by the presence of a viscosity gradient in the second case, in the sense of reducing the critical Reynolds number beyond which perturbations are amplified. This result supports the hypothesis that the current oscillations observed in the polarization curve may originate from a hydrodynamic instability.

Simulation de flammes interactives en temps réel

Vidéos et images des résultats disponible à : http://www.iro.umontreal.ca/labs/infographie/theses/fatnasss/

Approche cartésienne pour le calcul du vent en terrain complexe avec application à la propagation des feux de forêt

La méthode de projection et l'approche variationnelle de Sasaki sont deux techniques permettant d'obtenir un champ vectoriel à divergence nulle à partir d'un champ initial quelconque. Pour une vitesse d'un vent en haute altitude, un champ de vitesse sur une grille décalée est généré au-dessus d'une topographie donnée par une fonction analytique. L'approche cartésienne nommée Embedded Boundary Method est utilisée pour résoudre une équation de Poisson découlant de la projection sur un domaine irrégulier avec des conditions aux limites mixtes. La solution obtenue permet de corriger le champ initial afin d'obtenir un champ respectant la loi de conservation de la masse et prenant également en compte les effets dûs à la géométrie du terrain. Le champ de vitesse ainsi généré permettra de propager un feu de forêt sur la topographie à l'aide de la méthode iso-niveaux. L'algorithme est décrit pour le cas en deux et trois dimensions et des tests de convergence sont effectués.