Calculations of Boiling Two-Phase Flow Using a Porous Media Model

Boiling two-phase flow and the equations governing the motion of fluid in two-phase flows are discussed in this thesis. Disposition of the governing equations in three-dimensional complex geometries is considered from the perspective of the porous medium concept. The equations governing motion in two-phase flows were formulated, discretized and implemented in a subroutine for pressure-velocity solution utilizing the SIMPLE algorithm modified for two-phase flow. The subroutine was included in PORFLO, which is a three-dimensional 5-equation porous media model developed at VTT by Jaakko Miettinen. The development of two-phase flow and the resulting void fraction distribution was predicted in a geometry resembling a section of BWR fuel bundle in a couple of test cases using PORFLO.

CFD Simulation of Two-phase and Three-phase Flows in Internal-loop Airlift

Airlift reactors are pneumatically agitated reactors that have been widely used in chemical, petrochemical, and bioprocess industries, such as fermentation and wastewater treatment. Computational Fluid Dynamics (CFD) has become more popular approach for design, scale-up and performance evaluation of such reactors. In the present work numerical simulations for internal-loop airlift reactors were performed using the transient Eulerian model with CFD package, ANSYS Fluent 12.1. The turbulence in the liquid phase is described using κ- ε the model. Global hydrodynamic parameters like gas holdup, gas velocity and liquid velocity have been investigated for a range of superficial gas velocities, both with 2D and 3D simulations. Moreover, the study of geometry and scale influence on the reactor have been considered. The results suggest that both, geometry and scale have significant effects on the hydrodynamic parameters, which may have substantial effects on the reactor performance. Grid refinement and time-step size effect have been discussed. Numerical calculations with gas-liquid-solid three-phase flow system have been carried out to investigate the effect of solid loading, solid particle size and solid density on the hydrodynamic characteristics of internal loop airlift reactor with different superficial gas velocities. It was observed that averaged gas holdup is significantly decreased with increasing slurry concentration. Simulations show that the riser gas holdup decreases with increase in solid particle diameter. In addition, it was found that the averaged solid holdup increases in the riser section with the increase of solid density. These produced results reveal that CFD have excellent potential to simulate two-phase and three-phase flow system.

Two-phase flow measurements with spacer grid in a rod bundle geometry

Heat transfer effectiveness in nuclear rod bundles is of great importance to nuclear reactor safety and economics. An important design parameter is the Critical Heat Flux (CHF), which limits the transferred heat from the fuel to the coolant. The CHF is determined by flow behaviour, especially the turbulence created inside the fuel rod bundle. Adiabatic experiments can be used to characterize the flow behaviour separately from the heat transfer phenomena in diabatic flow. To enhance the turbulence, mixing vanes are attached to spacer grids, which hold the rods in place. The vanes either make the flow swirl around a single sub-channel or induce cross-mixing between adjacent sub-channels. In adiabatic two-phase conditions an important phenomenon that can be investigated is the effect of the spacer on canceling the lift force, which collects the small bubbles to the rod surfaces leading to decreased CHF in diabatic conditions and thus limits the reactor power. Computational Fluid Dynamics (CFD) can be used to simulate the flow numerically and to test how different spacer configurations affect the flow. Experimental data is needed to validate and verify the used CFD models. Especially the modeling of turbulence is challenging even for single-phase flow inside the complex sub-channel geometry. In two-phase flow other factors such as bubble dynamics further complicate the modeling. To investigate the spacer grid effect on two-phase flow, and to provide further experimental data for CFD validation, a series of experiments was run on an adiabatic sub-channel flow loop using a duct-type spacer grid with different configurations. Utilizing the wire-mesh sensor technology, the facility gives high resolution experimental data in both time and space. The experimental results indicate that the duct-type spacer grid is less effective in canceling the lift force effect than the egg-crate type spacer tested earlier.

Kaksifaasijäähdytysjärjestelmään liittyvien pumppujen vertailu ja kehittäminen sarjatuotannon tarpeisiin

Tehoelektroniikalta vaaditaan nykyään parempaa suorituskykyä entistä pienemmässä tilassa. Tämä luo haasteen riittävälle jäähdytykselle. Eräs ratkaisu on käyttää kaksifaasijäähdytystä, jolla aikaansaadaan tehokas lämmönsiirto komponenttien pinnalta. Lämmönsiirtonesteinä voidaan käyttää kylmäaineita tai muita alhaisessa lämpötilassa kiehuvia nesteitä. Tällaisille nesteille on tyypillistä alhainen höyrynpaine sekä matala viskositeetti. Nämä ominaisuudet tuovat haasteita nesteen pumppaukseen. Tässä työssä tarkastellaan R-134A:ta sekä Novec 7000:ta, perehdytään niiden fysikaalisiin ominaisuuksiin sekä materiaaliyhteensopivuuksiin ja näiden tietojen pohjalta etsitään sopivaa pumpputyyppiä kaksifaasijäähdytysjärjestelmään. Tehoelektroniikan jäähdytysjärjestelmän pumpun on oltava edullinen muuhun järjestelmään nähden. Tyypillinen kiertopumppu nestejäähdytysjärjestelmässä on pieni keskipakopumppu. Alhaisen kiehumispisteen vuoksi kavitointiriski kasvaa ja tämä voi vahingoittaa pumppua. Myös matala viskositeetti tuo haasteita vuotoherkkyyden kasvamisen myötä, joten mekaanisilla aksiaalitiivisteillä varustetut pumput eivät ole pitkäikäisiä. Kylmäainejärjestelmiin tarkoitetut pumput ovat arvokkaita, eikä näin ollen sovellu edullisiin jäähdytysjärjestelmiin. Tässä työssä käydään läpi erilaisia pumpputyyppejä, jotka voisivat soveltua pitkäikäiseen pumppaukseen ilman huoltotöitä. Näiden tietojen perusteella kehitetään edullista ja pitkäikäistä pumppua pieniin kaksifaasijäähdytysjärjestelmiin nesteiden fysikaaliset ominaisuudet huomioon ottaen. Kehitetyn pumpun ominaisuuksia ja kustannuksia vertaillaan kaupallisiin ratkaisuihin ottaen huomioon sarjavalmistus. Itse valmistettuna pienelle sisäryntöiselle hammaspyöräpumpulle jää hintaa alle kymmenesosa markkinoilta löytyviin kylmäaineille soveltuviin pumppuun.

Pulssitettu laserablaatio mikrofluidisten kanavistojen sarjatuotantomenetelmänä

Puolijohteiden yleistyttyä vuodesta 1948 alkaen, ovat elektroniset laitteet pienentyneet jatkuvasti tehojen kuitenkin kasvaessa. Kasvaneet tehotiheydet kuitenkin vaikeuttavat laitesuunnittelua, sillä puoljohdekomponenttien suorituskyvylle ja eliniälle on oleellista lämpötilojen ja lämpötilavaihteluiden minimointi. Perinteisen ilmajäähdytyksen lähestyessä rajojaan niin kokonaistehon kuin järkevän energiatehokkuudenkin suhteen, on parhaaksi seuraavaksi teknologiaksi ennustettu kaksifaasijäähdytystä, jonka suorituskyky ja energiatehokkuus ovat vaaditulla tasolla. Kaksifaasijäähdytyksen optimaaliselle toiminnalle tärkeää on hyvin suunniteltu ja tarkasti valmistettu lämmönsiirtopinta, jota kutsutaan mikrokanavistoksi. Pulssitettu laserkaiverrus on edistynyt valmistustekniikka, jonka tarkkuus ja luotettavuus sopisivat mikrokanavistojen valmistamiseen. Laserkaiverruksella saavutettavat lopputulokset vaihtelevat kuitenkin materiaalista riippuen ja kupari – jota käytetään yleisesti lämmönjohteena – on eräs huonoimmin lasertyöstöön reagoivista materiaaleista ja siksi on oleellista selvittää laser-kaiverruksen toimivuutta kuparisten mikrokanavistojen valmistuksessa. Pulssitetun laser-kaiverruksen eri variaatioista nanosekunti-luokan pulssinpituuksilla toimivat laitteet ovat jatkuvan tuotannon kannalta paras vaihtoehto niiden hyvän tuottavuuden, saatavuuden sekä kohtuullisen alkuinvestoinnin vuoksi. Käytännön kaiverruskokeiden perusteella selvisi, että menetelmä on laatunsa ja tarkkuutensa puolesta sopiva varsinaiseen tuotantoon. Kaiverruksen tehokkuus kuparia työstettäessä on kuitenkin ennakoituakin heikompi ja niin valmistus- kuin suunnitelu-prosessikin vaativat vielä jatkotutkimusta ja -kehitystä.

Analysis of Phase Separator Design Criteria Using Computational Fluid Dynamics

Gravitational phase separation is a common unit operation found in most large-scale chemical processes. The need for phase separation can arise e.g. from product purification or protection of downstream equipment. In gravitational phase separation, the phases separate without the application of an external force. This is achieved in vessels where the flow velocity is lowered substantially compared to pipe flow. If the velocity is low enough, the denser phase settles towards the bottom of the vessel while the lighter phase rises. To find optimal configurations for gravitational phase separator vessels, several different geometrical and internal design features were evaluated based on simulations using OpenFOAM computational fluid dynamics (CFD) software. The studied features included inlet distributors, vessel dimensions, demister configurations and gas phase outlet configurations. Simulations were conducted as single phase steady state calculations. For comparison, additional simulations were performed as dynamic single and two-phase calculations. The steady state single phase calculations provided indications on preferred configurations for most above mentioned features. The results of the dynamic simulations supported the utilization of the computationally faster steady state model as a practical engineering tool. However, the two-phase model provides more truthful results especially with flows where a single phase does not determine the flow characteristics.

Coolability of porous core debris beds : Effects of bed geometry and multi-dimensional flooding

This thesis addresses the coolability of porous debris beds in the context of severe accident management of nuclear power reactors. In a hypothetical severe accident at a Nordic-type boiling water reactor, the lower drywell of the containment is flooded, for the purpose of cooling the core melt discharged from the reactor pressure vessel in a water pool. The melt is fragmented and solidified in the pool, ultimately forming a porous debris bed that generates decay heat. The properties of the bed determine the limiting value for the heat flux that can be removed from the debris to the surrounding water without the risk of re-melting. The coolability of porous debris beds has been investigated experimentally by measuring the dryout power in electrically heated test beds that have different geometries. The geometries represent the debris bed shapes that may form in an accident scenario. The focus is especially on heap-like, realistic geometries which facilitate the multi-dimensional infiltration (flooding) of coolant into the bed. Spherical and irregular particles have been used to simulate the debris. The experiments have been modeled using 2D and 3D simulation codes applicable to fluid flow and heat transfer in porous media. Based on the experimental and simulation results, an interpretation of the dryout behavior in complex debris bed geometries is presented, and the validity of the codes and models for dryout predictions is evaluated. According to the experimental and simulation results, the coolability of the debris bed depends on both the flooding mode and the height of the bed. In the experiments, it was found that multi-dimensional flooding increases the dryout heat flux and coolability in a heap-shaped debris bed by 47–58% compared to the dryout heat flux of a classical, top-flooded bed of the same height. However, heap-like beds are higher than flat, top-flooded beds, which results in the formation of larger steam flux at the top of the bed. This counteracts the effect of the multi-dimensional flooding. Based on the measured dryout heat fluxes, the maximum height of a heap-like bed can only be about 1.5 times the height of a top-flooded, cylindrical bed in order to preserve the direct benefit from the multi-dimensional flooding. In addition, studies were conducted to evaluate the hydrodynamically representative effective particle diameter, which is applied in simulation models to describe debris beds that consist of irregular particles with considerable size variation. The results suggest that the effective diameter is small, closest to the mean diameter based on the number or length of particles.

Patient-perceived health-related quality of life during recovery after total hip arthroplasty – a 6-month follow-up study

Potilaiden käsitys terveyteen liittyvästä elämänlaadusta lonkan tekonivelleikkauksen jälkeisenä toipumisaikana – kuuden kuukauden seurantatutkimus Tässä kaksivaiheisessa seurantatutkimuksessa tarkasteltiin potilaiden käsitystä terveyteen liittyvästä elämänlaadusta lonkan tekonivelleikkauksen jälkeisenä toipumisaikana. Tutkimuksen ensimmäisessä vaiheessa tarkoituksena oli sekä kuvailla potilaiden kokemuksia potilaana olosta, saamastaan hoidosta ja terveyspalveluorganisaatiosta että analysoida aikaisempien tutkimusten perusteella leikkauksen tuloksia potilaan kannalta. Toisessa vaiheessa tarkoituksena oli arvioida potilaiden kokemaa elämänlaatua leikkauksen jälkeen, ja sitä vaikuttivatko primaaritulokset (fyysinen toimintakyky, kipu, ahdistus) tai taloudelliset seuraukset (potilaiden itsensämaksamat kustannukset, palvelujen käyttö) terveyteen liittyvään elämänlaatuun. Tutkimuksen tavoitteena oli löytää mahdolliset kriittiset ajankohdat tai tekijät, jotka saattavat hidastaa toipumista ja siten huonontaa potilaiden elämänlaatua. Tätä tietoa voidaan käyttää hoitotyössä kun suunnitellaan sopivaa hoitoa ja tukea toipumisajalle. Tutkimuksen ensimmäisessä vaiheessa primaarileikkaukseen tulevat potilaat (n = 17) kuvailivat teemahaastatteluissa kokemuksiaan kahdesti leikkauksen jälkeen. Haastatteluaineisto analysoitiin induktiivisella sisällönanalyysilla. Lisäksi 17 tutkimusartikkelista analysoitiin deduktiivisella sisällönanalyysilla leikkauksen tuloksia potilaalle, tuloksiin vaikuttavia tekijöitä ja käytetyt tutkimusmetodit. Toisessa vaiheessa primaari- tai revisioleikkaukseen tulevat potilaat (n = 100) arvioivat leikkauksen tuloksia kuuden kuukauden ajan leikkauksen jälkeen: terveyteen liittyvää elämänlaatua, primaarituloksia ja taloudellisia seurauksia. Aineisto kerättiin erilaisilla mittareilla: Sickness Impact Profile, Finnish Version, Stait-Trait Anxiety Inventory, ja Numeric Rating Scale. Lisäksi käytettiin tätä tutkimusta varten tehtyjä kyselylomakkeita: Fyysinen toimintakyky-mittari, Palvelujen käyttö-mittari ja Kustannusmittari. Tutkimuksen toiseen vaiheen tulokset analysoitiin tilastollisilla menetelmillä. Potilaiden terveyteen liittyvä elämänlaatu parani ja kipu lievittyi leikkauksen jälkeen ja fyysinen toimintakyky lisääntyi toipumisaikana. Positiivisista muutoksista huolimatta potilaat kokivat ahdistusta samassa määrin kuin ennen leikkaustakin. Palvelujen käyttö vaihteli toipumisajan kuluessa ja potilaiden maksamissa kustannuksissa oli suuria vaihteluita. Fyysisen toimintakyvyn lisääntyminen ja kivun lieveneminen paransivat terveyteen liittyvää elämänlaatua. Sen sijaan huonompi elämänlaatu toipumisaikana oli yhteydessä suurempaan palvelujen käyttöön, kun taas kustannuksilla ei ollut yhteyttä elämänlaatuun. Potilaiden ominaispiirteet tulisi ottaa enemmän huomioon suunniteltaessa sopivaa leikkauksenjälkeistä hoitoa ja tukea. Potilaat tarvitsevat yksilöllisiä ohjeita, sillä monet taustatekijät (esim. ikä, sukupuoli, preoperatiivinen kipu, siviilisääty, ja leikkaustyyppi) vaikuttavat toipumiseen.

Öljynjalostamon tyhjöyksikön esilämmitysuunin optimointi

Työn tarkoituksena oli vertailla kahta kaupallista prosessiuunisimulointiohjelmaa sekä selvittää onko toinen ohjelma soveltuva Neste Jacobs Oy:n käyttöön. Tavoitteena oli selvittää ohjelmien ominaisuuksien eroja ja rajoitteita sekä tuloksien yhdenmukaisuutta. Kirjallisuusosassa perehdyttiin prosessiuunien toimintaperiaatteisiin, yleisiin mitoitusperiaatteisiin ja säteilylämmönsiirron perusteisiin. Kirjallisuusosassa käsiteltiin myös prosessiuunien polttimien toimintaa ja NOx-päästöjä vähentäviä polttotekniikoita. Molemmilla laskentaohjelmilla laadittiin kuumaöljy- ja tyhjöuuni, joiden tuloksia vertailtiin tärkeimpien prosessiarvojen perusteella ja toisen ohjelman osalta myös laitevalmistajan ilmoittamiin arvoihin. Tällä ohjelmalla tutkittiin tyhjöuunin höyrynruiskutuksen ja koksaantumisen vaikutuksia prosessiin sekä säteilyosan leveyden muutoksen vaikutusta säteilylämmönsiirtoon. Ohjelmien ominaisuuksien keskeisimmät erot liittyivät kaksifaasivirtauksen laskentaan, prosessiaineen massavirran rajoituksiin, ilman esilämmittimen käyttöön ja laskentamallien paineen kontrollointiin. Tulokset erosivat savukaasun ja tuubien lämpötilojen ja painehäviön laskennassa. Tuloksien ja ominaisuuksien perusteella toisen ohjelman ei todettu olevan täysin soveltuva Neste Jacobs Oy:n käyttöön.

Virtauslaskennankäyttö keskipakopumpun eroosiokulumisen mallintamisessa

Diplomityössä kehitetään malli tiheän kaksifaasivirtauksen aiheuttaman eroosiokulumi-sen mallintamiseksi, ratkaistaan virtauskenttä kahdessa erilaisessa keskipakopumpussa, sovelletaan kehitettyä mallia, sekävertaillaan mallin antamia tuloksia käytännön kokeissa saavutettuihin tuloksiin. Työssä on erityisenä mielenkiinnon kohteena savukaasupesurin pumppu. Työn alkuosa sisältää tarkemman kuvauksen savukaasupesurin toi-minnasta. Numeerinen ratkaisu ja laskentahilan generointi suoritetaan ANSYS CFX- ja Turbo-Grid-ohjelmistoilla. Laskennassa virtauksen Navier-Stokesin yhtälöt on aikakeskiarvo-tettu ja ratkaistu käyttäen kontrollitilavuusmenetelmää. Tiheiden kaksifaasivirtausten eroosiokulumista on mallinnettu tekijän kehittämällä mallilla, jonka käytännön toteutus-ta ei kuitenkaan saateta julkiseksi, koodin kehittämisessä käytetty teoria on kuitenkin esitetty työssä. Työn piirissä tehtiin myös kulumiskokeita Sulzer Pumps Finlandin Karhulan tehtailla, sekä vertailtiin simuloinnin tuloksia aikaisemmissa kokeissa saavutettuihin. Koejärjestelyt kuvataan työssä.

Typpeä sisältävien jätevesien käsittely 2-n-PRO menetelmällä

Työssä tutkittiin kaksivaiheisen typenpoistoprosessin (2-N-PRO) soveltuvuutta Joutsenon Kukkuroinmäen aluejätekeskuksen kompostointilaitoksen jätevesille pilot-kokein 12.1.- 5.4.2006. Kompostilaitoksella on jätevesien esikäsittelytarve korkeista ammoniumtyppipitoisuuksista johtuen. Pilot-laitteisto koostuu sekoitussäiliöstä, strippaustornista ja katalyyttipolttimesta. Käsiteltävän jäteveden pH nostetaan korkealle tasolle, jolloin ammoniumtyppi muuttuu ammoniakiksi. Vesi johdetaan strippaustorniin, jossa se sadetetaan tornin pohjalle. Ammoniakki erottuu sadetuksessa ilmaan, joka imetään katalyyttipolttimelle. Katalyyttinen poltin käsittelee ammoniakkia typpikaasuksi. Pilot-kokeet suoritettiin jatkuvatoimisesti. Laitteisto pystyy erottamaan jätevedestä ammoniumtyppeä ammoniakiksi ja käsittelemään ammoniakin pääosin typpikaasuksi. Lisäksi suoritettiin panoskoe, jonka tulokset tukevat jatkuvatoimisesta käytöstä saatuja tuloksia.

Gas-liquid mass transfer in bubbly flow: Estimation of mass transfer, bubble size and reactor performance in various applications

Gas-liquid mass transfer is an important issue in the design and operation of many chemical unit operations. Despite its importance, the evaluation of gas-liquid mass transfer is not straightforward due to the complex nature of the phenomena involved. In this thesis gas-liquid mass transfer was evaluated in three different gas-liquid reactors in a traditional way by measuring the volumetric mass transfer coefficient (kLa). The studied reactors were a bubble column with a T-junction two-phase nozzle for gas dispersion, an industrial scale bubble column reactor for the oxidation of tetrahydroanthrahydroquinone and a concurrent downflow structured bed.The main drawback of this approach is that the obtained correlations give only the average volumetric mass transfer coefficient, which is dependent on average conditions. Moreover, the obtained correlations are valid only for the studied geometry and for the chemical system used in the measurements. In principle, a more fundamental approach is to estimate the interfacial area available for mass transfer from bubble size distributions obtained by solution of population balance equations. This approach has been used in this thesis by developing a population balance model for a bubble column together with phenomenological models for bubble breakage and coalescence. The parameters of the bubble breakage rate and coalescence rate models were estimated by comparing the measured and calculated bubble sizes. The coalescence models always have at least one experimental parameter. This is because the bubble coalescence depends on liquid composition in a way which is difficult to evaluate using known physical properties. The coalescence properties of some model solutions were evaluated by measuring the time that a bubble rests at the free liquid-gas interface before coalescing (the so-calledpersistence time or rest time). The measured persistence times range from 10 msup to 15 s depending on the solution. The coalescence was never found to be instantaneous. The bubble oscillates up and down at the interface at least a coupleof times before coalescence takes place. The measured persistence times were compared to coalescence times obtained by parameter fitting using measured bubble size distributions in a bubble column and a bubble column population balance model. For short persistence times, the persistence and coalescence times are in good agreement. For longer persistence times, however, the persistence times are at least an order of magnitude longer than the corresponding coalescence times from parameter fitting. This discrepancy may be attributed to the uncertainties concerning the estimation of energy dissipation rates, collision rates and mechanisms and contact times of the bubbles.

Development of a PC Tool for Estimating Pressure Drops in Forced-convection Boiler Tubes

Kaksifaasivirtauksen kuvaamiseen käytettävät mallit, ja menetelmät kaksifaasivirtauksen painehäviön määrittämiseksi kehittyvät yhä monimutkaisimmiksi. Höyrystinputkissa tapahtuvien painehäviöiden arvioinnin vaatiman laskennan suorittamiseksi tietokoneohjelman kehittäminen on välttämätöntä. Tässä työssä on kehitetty itsenäinen PC-ohjelma painehäviöiden arvioimiseksi pakotetulle konvektiovirtaukselle pystysuorissa höyrykattilan höyrystinputkissa. Veden ja vesihöyryn aineominaisuuksien laskentaan käytetään IAPWS-IF97 –yhtälökokoelmaa sekä muita tarvittavia IAPWS:n suosittelemia yhtälöitä. Höyrystinputkessa kulloinkin vallitsevan virtausmuodon määrittämiseen käytetään sovelluskelpoisia virtausmuotojen välisiä rajoja kuvaavia yhtälöitä. Ohjelmassa käytetään painehäviön määritykseen kirjallisuudessa julkaistuja yhtälöitä, virtausmuodosta riippuen, alijäähtyneelle virtaukselle, kupla-, tulppa- ja rengasvirtaukselle sekä tulistetun höyryn virtaukselle. Ohjelman laskemia painehäviöarvioita verrattiin kirjallisuudesta valittuihin mittaustuloksiin. Laskettujen painehäviöiden virhe vaihteli välillä –19.5 ja +23.9 %. Virheiden itseisarvojen keskiarvo oli 12.8 %.

ADVANCED NUMERICAL METHODS FOR SIMULATING MICROBUBBLE IN AN AERATED MIXING TANK

Fluid mixing in mechanically agitated tanks is one of the major unit operations in many industries. Bubbly flows have been of interest among researchers in physics, medicine, chemistry and technology over the centuries. The aim of this thesis is to use advanced numerical methods for simulating microbubble in an aerated mixing tank. Main components of the mixing tank are a cylindrical vessel, a rotating Rushton turbine and the air nozzle. The objective of Computational Fluid Dynamics (CFD) is to predict fluid flow, heat transfer, mass transfer and chemical reactions. The CFD simulations of a turbulent bubbly flow are carried out in a cylindrical mixing tank using large eddy simulation (LES) and volume of fluid (VOF) method. The Rushton turbine induced flow is modeled by using a sliding mesh method. Numerical results are used to describe the bubbly flows in highly complex liquid flow. Some of the experimental works related to turbulent bubbly flow in a mixing tank are briefly reported. Numerical simulations are needed to complete and interpret the results of the experimental work. Information given by numerical simulations has a major role in designing and scaling-up mixing tanks. The results of this work have been reported in the following scientific articles: ·Honkanen M., Koohestany A., Hatunen T., Saarenrinne P., Zamankhan P., Large eddy simulations and PIV experiments of a two-phase air-water mixer, in Proceedings of ASME Fluids Engineering Summer Conference (2005). ·Honkanen M., Koohestany A., Hatunen T., Saarenrinne P., Zamankhan P., Dynamical States of Bubbling in an Aerated Stirring Tank, submitted to J. Computational Physics.