Experimental analysis of the flow between stay and guide vanes of a pump-turbine in pumping mode

An experimental study was conducted in a pump-turbine model in pumping mode, in order to characterize the flow field structure in the region between stay and guide vanes, using mainly the laser-Doppler anemometry in a two-color and back-scattered light-based system. The structure of the steady and unsteady flow was analyzed. The measurements were performed at three operation points. The obtained data provide appropriate boundary conditions and a good base of validation for numerical codes, and for the understanding of main loss mechanisms of this complex flow.

Gas-solid two-phase flow in the riser of circulating fluidized beds: mathematical modelling and numerical simulation

A mathematical model is developed for gas-solids flows in circulating fluidized beds. An Eulerian formulation is followed based on the two-fluids model approach where both the fluid and the particulate phases are treated as a continuum. The physical modelling is discussed, including the formulation of boundary conditions and the description of the numerical methodology. Results of numerical simulation are presented and discussed. The model is validated through comparison to experiment, and simulation is performed to investigate the effects on the flow hydrodynamics of the solids viscosity.

Surface Tension Implementation for Gensmac2d

In the present work we describe a method which allows the incorporation of surface tension into the GENSMAC2D code. This is achieved on two scales. First on the scale of a cell, the surface tension effects are incorporated into the free surface boundary conditions through the computation of the capillary pressure. The required curvature is estimated by fitting a least square circle to the free surface using the tracking particles in the cell and in its close neighbors. On a sub-cell scale, short wavelength perturbations are filtered out using a local 4-point stencil which is mass conservative. An efficient implementation is obtained through a dual representation of the cell data, using both a matrix representation, for ease at identifying neighbouring cells, and also a tree data structure, which permits the representation of specific groups of cells with additional information pertaining to that group. The resulting code is shown to be robust, and to produce accurate results when compared with exact solutions of selected fluid dynamic problems involving surface tension.

Numerical study of wedge supported oblique shock wave-oblique detonation wave transitions

The results of a numerical study of premixed Hydrogen-air flows ignition by an oblique shock wave (OSW) stabilized by a wedge are presented, in situations when initial and boundary conditions are such that transition between the initial OSW and an oblique detonation wave (ODW) is observed. More precisely, the objectives of the paper are: (i) to identify the different possible structures of the transition region that exist between the initial OSW and the resulting ODW and (ii) to evidence the effect on the ODW of an abrupt decrease of the wedge angle in such a way that the final part of the wedge surface becomes parallel to the initial flow. For such a geometrical configuration and for the initial and boundary conditions considered, the overdriven detonation supported by the initial wedge angle is found to relax towards a Chapman-Jouguet detonation in the region where the wedge surface is parallel to the initial flow. Computations are performed using an adaptive, unstructured grid, finite volume computer code previously developed for the sake of the computations of high speed, compressible flows of reactive gas mixtures. Physico-chemical properties are functions of the local mixture composition, temperature and pressure, and they are computed using the CHEMKIN-II subroutines.

Modeling of fireside deposit formation in two industrial furnaces

Fireside deposits can be found in many types of utility and industrial furnaces. The deposits in furnaces are problematic because they can reduce heat transfer, block gas paths and cause corrosion. To tackle these problems, it is vital to estimate the influence of deposits on heat transfer, to minimize deposit formation and to optimize deposit removal. It is beneficial to have a good understanding of the mechanisms of fireside deposit formation. Numerical modeling is a powerful tool for investigating the heat transfer in furnaces, and it can provide valuable information for understanding the mechanisms of deposit formation. In addition, a sub-model of deposit formation is generally an essential part of a comprehensive furnace model. This work investigates two specific processes of fireside deposit formation in two industrial furnaces. The first process is the slagging wall found in furnaces with molten deposits running on the wall. A slagging wall model is developed to take into account the two-layer structure of the deposits. With the slagging wall model, the thickness and the surface temperature of the molten deposit layer can be calculated. The slagging wall model is used to predict the surface temperature and the heat transfer to a specific section of a super-heater tube panel with the boundary condition obtained from a Kraft recovery furnace model. The slagging wall model is also incorporated into the computational fluid dynamics (CFD)-based Kraft recovery furnace model and applied on the lower furnace walls. The implementation of the slagging wall model includes a grid simplification scheme. The wall surface temperature calculated with the slagging wall model is used as the heat transfer boundary condition. Simulation of a Kraft recovery furnace is performed, and it is compared with two other cases and measurements. In the two other cases, a uniform wall surface temperature and a wall surface temperature calculated with a char bed burning model are used as the heat transfer boundary conditions. In this particular furnace, the wall surface temperatures from the three cases are similar and are in the correct range of the measurements. Nevertheless, the wall surface temperature profiles with the slagging wall model and the char bed burning model are different because the deposits are represented differently in the two models. In addition, the slagging wall model is proven to be computationally efficient. The second process is deposit formation due to thermophoresis of fine particles to the heat transfer surface. This process is considered in the simulation of a heat recovery boiler of the flash smelting process. In order to determine if the small dust particles stay on the wall, a criterion based on the analysis of forces acting on the particle is applied. Time-dependent simulation of deposit formation in the heat recovery boiler is carried out and the influence of deposits on heat transfer is investigated. The locations prone to deposit formation are also identified in the heat recovery boiler. Modeling of the two processes in the two industrial furnaces enhances the overall understanding of the processes. The sub-models developed in this work can be applied in other similar deposit formation processes with carefully-defined boundary conditions.

OL1/OL2-suojarakennuksen jäähdytteenmenetysonnettomuusanalyysit APROS 5.09 -ohjelmistolla

Tässä työssä on tutkittu OL1/OL2-suojarakennuksen käyttäytymistä jäähdytteenmenetysonnettomuuden eli LOCA:n aikana. Onnettomuuden simulointiin on kehitetty suojarakennusmalli suomalaiseen APROS 5.09 - ohjelmistoon sisältyvällä LP-koodilla (Lumped Parameter Code). Työssä on keskitytty suojarakennuksen kannalta oleellisimpien suureiden: kaasutilavuuksien paineen sekä lämpötilan ja lauhdutusaltaan lämpötilan ja pinnankorkeuden ajalliseen käyttäytymiseen. Mallinnetut LOCA:t ovat päähöyrylinjan ja sammutetun reaktorin jäähdytysjärjestelmän putkikatkoksia. Simulointeja on tehty laitoksen täyden tehon ja kuumavalmiuden lähtötiloissa ja tarkasteltavien suureiden käyttäytymistä on tutkittu erikseen valituilla konservatiivisilla oletuksilla. Laskentatapaukset rajoittuvat lyhyeen aikaväliin 27.78 tuntia kuvitellusta putkirikosta eteenpäin. Tuloksia on verrattu OL1/OL2- laitostoimittajan, Westinghouse Electric Sweden AB:n, tekemiin lisensiointianalyyseihin. Työssä on myös kuvattu APROS-laskennan epävarmuustekijöitä. Suojarakennuksen on todettu käyttäytyvän fysikaalisesti yhtenevästi lisensiointianalyyseissä ja APROS:lla tehdyissä vertailulaskuissa.

Höyryturbiinin kanavistojen numeerinen virtauslaskenta

Diplomityössä tutkittiin Fortumin Loviisan ydinvoimalaitoksen ulosvirtauskanaviston ja suurnopeuskosteudenerottimen toimintaa, sekä selvitettiin taustalla olevaa teoriaa ja aiemmin tehtyjä tutkimuksia. Tavoitteena oli ymmärtää ja esittää laitteiden toimintaa, sekä tutkia voiko ulosvirtauskanaviston suorituskykyä parantaa geometrian muutoksilla. Työssä luotiin tutkittaville kohteille geometriat ja laskentahilat, joiden avulla simuloitiin niiden toimintaa eri käyttötilanteissa numeerisen virtauslaskennan avulla. Laskennan reunaehdot saatiin olemassa olevasta prosessimallista ja aiemmista turbiiniselvityksistä. Ulosvirtauskanaviston suorituskyky laskettiin kolmella eri lauhdutinpaineella neljällä eri geometrialla. Geometrian muutokset vaikuttivat selkeästi ulosvirtauskanaviston suorituskykyyn ja sitä saatiin parannettua. Kaksi kolmesta muutoksesta, lisäkanavat ja oikaistu vesilippa, pa-ransivat suorituskykyä. Lokinsiipien poistaminen heikensi ulosvirtauskanaviston toi-mintaa. Suurnopeuskosteudenerottimen mallintaminen jäi lähtötietojen ja ajan puutteen takia hieman tavoitteesta. Sekä ulosvirtauskanaviston että suurnopeuskosteudenerotti-men jatkotutkimusta ja mahdollisia toimenpiteitä varten saatiin arvokasta uutta tietoa.

Large-eddy simulation of wind flows over complex terrains for wind energy applications

Wind energy has obtained outstanding expectations due to risks of global warming and nuclear energy production plant accidents. Nowadays, wind farms are often constructed in areas of complex terrain. A potential wind farm location must have the site thoroughly surveyed and the wind climatology analyzed before installing any hardware. Therefore, modeling of Atmospheric Boundary Layer (ABL) flows over complex terrains containing, e.g. hills, forest, and lakes is of great interest in wind energy applications, as it can help in locating and optimizing the wind farms. Numerical modeling of wind flows using Computational Fluid Dynamics (CFD) has become a popular technique during the last few decades. Due to the inherent flow variability and large-scale unsteadiness typical in ABL flows in general and especially over complex terrains, the flow can be difficult to be predicted accurately enough by using the Reynolds-Averaged Navier-Stokes equations (RANS). Large- Eddy Simulation (LES) resolves the largest and thus most important turbulent eddies and models only the small-scale motions which are more universal than the large eddies and thus easier to model. Therefore, LES is expected to be more suitable for this kind of simulations although it is computationally more expensive than the RANS approach. With the fast development of computers and open-source CFD software during the recent years, the application of LES toward atmospheric flow is becoming increasingly common nowadays. The aim of the work is to simulate atmospheric flows over realistic and complex terrains by means of LES. Evaluation of potential in-land wind park locations will be the main application for these simulations. Development of the LES methodology to simulate the atmospheric flows over realistic terrains is reported in the thesis. The work also aims at validating the LES methodology at a real scale. In the thesis, LES are carried out for flow problems ranging from basic channel flows to real atmospheric flows over one of the most recent real-life complex terrain problems, the Bolund hill. All the simulations reported in the thesis are carried out using a new OpenFOAM® -based LES solver. The solver uses the 4th order time-accurate Runge-Kutta scheme and a fractional step method. Moreover, development of the LES methodology includes special attention to two boundary conditions: the upstream (inflow) and wall boundary conditions. The upstream boundary condition is generated by using the so-called recycling technique, in which the instantaneous flow properties are sampled on aplane downstream of the inlet and mapped back to the inlet at each time step. This technique develops the upstream boundary-layer flow together with the inflow turbulence without using any precursor simulation and thus within a single computational domain. The roughness of the terrain surface is modeled by implementing a new wall function into OpenFOAM® during the thesis work. Both, the recycling method and the newly implemented wall function, are validated for the channel flows at relatively high Reynolds number before applying them to the atmospheric flow applications. After validating the LES model over simple flows, the simulations are carried out for atmospheric boundary-layer flows over two types of hills: first, two-dimensional wind-tunnel hill profiles and second, the Bolund hill located in Roskilde Fjord, Denmark. For the twodimensional wind-tunnel hills, the study focuses on the overall flow behavior as a function of the hill slope. Moreover, the simulations are repeated using another wall function suitable for smooth surfaces, which already existed in OpenFOAM® , in order to study the sensitivity of the flow to the surface roughness in ABL flows. The simulated results obtained using the two wall functions are compared against the wind-tunnel measurements. It is shown that LES using the implemented wall function produces overall satisfactory results on the turbulent flow over the two-dimensional hills. The prediction of the flow separation and reattachment-length for the steeper hill is closer to the measurements than the other numerical studies reported in the past for the same hill geometry. The field measurement campaign performed over the Bolund hill provides the most recent field-experiment dataset for the mean flow and the turbulence properties. A number of research groups have simulated the wind flows over the Bolund hill. Due to the challenging features of the hill such as the almost vertical hill slope, it is considered as an ideal experimental test case for validating micro-scale CFD models for wind energy applications. In this work, the simulated results obtained for two wind directions are compared against the field measurements. It is shown that the present LES can reproduce the complex turbulent wind flow structures over a complicated terrain such as the Bolund hill. Especially, the present LES results show the best prediction of the turbulent kinetic energy with an average error of 24.1%, which is a 43% smaller than any other model results reported in the past for the Bolund case. Finally, the validated LES methodology is demonstrated to simulate the wind flow over the existing Muukko wind farm located in South-Eastern Finland. The simulation is carried out only for one wind direction and the results on the instantaneous and time-averaged wind speeds are briefly reported. The demonstration case is followed by discussions on the practical aspects of LES for the wind resource assessment over a realistic inland wind farm.

Prosessilaitteen teräsrakenteen optimointi

Diplomityön tarkoituksena oli optimoida puukentän nosturin teräsrakenteiden kustannuksia ja selvittää halvemman valmistusmaan vaikutus nosturin kokonaiskustannukseen. Tavoitteena oli saada 15 % kokonaiskustannussäästö. Työn alussa perehdytään optimoitavan nosturin rakenteeseen ja mahdollisiin rakenteen optimoitaviin osiin. Teräsrakenteen optimointi alkoi kustannusrakenteen tarkalla selvittämisellä. Näin saatiin kuva nosturin kokonaiskustannuksesta ja mahdollisista säästökohteista. Uusien ideoiden kehittämisessä käytettiin hyväksi osatoimintojen pisteytystä ja ideamatriisia. Ideamatriisin lopputuloksena uudet rakenteet tarkastettiin nosturin kuormituksien suhteen ja valittiin kustannusten kannalta paras ratkaisu. Teräsrakenteeseen tehtiin muutoksia myös merirahdin konttien asettamien rajaehtojen mukaan. Kiinassa valmistettavalle teräsrakenteelle tehtiin rahtisuunnitelma, jossa selvitettiin optimaalinen rahtimuoto ja tarvittavien konttien lukumäärä. Tarjouskilpailun avulla saatiin halvin kuljetushinta valmistavalta tehtaalta työmaalle. Lopuksi nosturin optimoidun teräsrakenteen kustannuksia verrattiin työn alussa selvitettyyn kustannusrakenteeseen. Työn lopputuloksena päästiin 12 % potentiaaliseen kustannussäästöön, joka jäi kolme prosenttiyksikköä kustannussäästötavoitteesta.

Simulation and material calibration of ultra high strength steel (UHSS) S960 welded joint under static tensile test utilizing finite element method (FEM)

The aim of this work was to calibrate the material properties including strength and strain values for different material zones of ultra-high strength steel (UHSS) welded joints under monotonic static loading. The UHSS is heat sensitive and softens by heat due to welding, the affected zone is heat affected zone (HAZ). In this regard, cylindrical specimens were cut out from welded joints of Strenx® 960 MC and Strenx® Tube 960 MH, were examined by tensile test. The hardness values of specimens’ cross section were measured. Using correlations between hardness and strength, initial material properties were obtained. The same size specimen with different zones of material same as real specimen were created and defined in finite element method (FEM) software with commercial brand Abaqus 6.14-1. The loading and boundary conditions were defined considering tensile test values. Using initial material properties made of hardness-strength correlations (true stress-strain values) as Abaqus main input, FEM is utilized to simulate the tensile test process. By comparing FEM Abaqus results with measured results of tensile test, initial material properties will be revised and reused as software input to be fully calibrated in such a way that FEM results and tensile test results deviate minimum. Two type of different S960 were used including 960 MC plates, and structural hollow section 960 MH X-joint. The joint is welded by BöhlerTM X96 filler material. In welded joints, typically the following zones appear: Weld (WEL), Heat affected zone (HAZ) coarse grained (HCG) and fine grained (HFG), annealed zone, and base material (BaM). Results showed that: The HAZ zone is softened due to heat input while welding. For all the specimens, the softened zone’s strength is decreased and makes it a weakest zone where fracture happens while loading. Stress concentration of a notched specimen can represent the properties of notched zone. The load-displacement diagram from FEM modeling matches with the experiments by the calibrated material properties by compromising two correlations of hardness and strength.

The solutions and tunneling properties for a linear potential and an inverted harmonic oscillator in an infinite well

In this work we look at two different 1-dimensional quantum systems. The potentials for these systems are a linear potential in an infinite well and an inverted harmonic oscillator in an infinite well. We will solve the Schrödinger equation for both of these systems and get the energy eigenvalues and eigenfunctions. The solutions are obtained by using the boundary conditions and numerical methods. The motivation for our study comes from experimental background. For the linear potential we have two different boundary conditions. The first one is the so called normal boundary condition in which the wave function goes to zero on the edge of the well. The second condition is called derivative boundary condition in which the derivative of the wave function goes to zero on the edge of the well. The actual solutions are Airy functions. In the case of the inverted oscillator the solutions are parabolic cylinder functions and they are solved only using the normal boundary condition. Both of the potentials are compared with the particle in a box solutions. We will also present figures and tables from which we can see how the solutions look like. The similarities and differences with the particle in a box solution are also shown visually. The figures and calculations are done using mathematical software. We will also compare the linear potential to a case where the infinite wall is only on the left side. For this case we will also show graphical information of the different properties. With the inverted harmonic oscillator we will take a closer look at the quantum mechanical tunneling. We present some of the history of the quantum tunneling theory, its developers and finally we show the Feynman path integral theory. This theory enables us to get the instanton solutions. The instanton solutions are a way to look at the tunneling properties of the quantum system. The results are compared with the solutions of the double-well potential which is very similar to our case as a quantum system. The solutions are obtained using the same methods which makes the comparison relatively easy. All in all we consider and go through some of the stages of the quantum theory. We also look at the different ways to interpret the theory. We also present the special functions that are needed in our solutions, and look at the properties and different relations to other special functions. It is essential to notice that it is possible to use different mathematical formalisms to get the desired result. The quantum theory has been built for over one hundred years and it has different approaches. Different aspects make it possible to look at different things.

The solutions and tunneling properties for a linear potential and an inverted harmonic oscillator in an infinite well

In this work we look at two different 1-dimensional quantum systems. The potentials for these systems are a linear potential in an infinite well and an inverted harmonic oscillator in an infinite well. We will solve the Schrödinger equation for both of these systems and get the energy eigenvalues and eigenfunctions. The solutions are obtained by using the boundary conditions and numerical methods. The motivation for our study comes from experimental background. For the linear potential we have two different boundary conditions. The first one is the so called normal boundary condition in which the wave function goes to zero on the edge of the well. The second condition is called derivative boundary condition in which the derivative of the wave function goes to zero on the edge of the well. The actual solutions are Airy functions. In the case of the inverted oscillator the solutions are parabolic cylinder functions and they are solved only using the normal boundary condition. Both of the potentials are compared with the particle in a box solutions. We will also present figures and tables from which we can see how the solutions look like. The similarities and differences with the particle in a box solution are also shown visually. The figures and calculations are done using mathematical software. We will also compare the linear potential to a case where the infinite wall is only on the left side. For this case we will also show graphical information of the different properties. With the inverted harmonic oscillator we will take a closer look at the quantum mechanical tunneling. We present some of the history of the quantum tunneling theory, its developers and finally we show the Feynman path integral theory. This theory enables us to get the instanton solutions. The instanton solutions are a way to look at the tunneling properties of the quantum system. The results are compared with the solutions of the double-well potential which is very similar to our case as a quantum system. The solutions are obtained using the same methods which makes the comparison relatively easy. All in all we consider and go through some of the stages of the quantum theory. We also look at the different ways to interpret the theory. We also present the special functions that are needed in our solutions, and look at the properties and different relations to other special functions. It is essential to notice that it is possible to use different mathematical formalisms to get the desired result. The quantum theory has been built for over one hundred years and it has different approaches. Different aspects make it possible to look at different things.

Opistorakennusten energiatehokkuuden parantaminen

Tämän diplomityön tavoitteena on selvittää Kymenlaakson Opiston energiatehokkuuden parantamista ja tutkia onko olemassa selvästi taloudellisempi sekä ekologisempi tapa kattaa Opiston lämmitystarve verrattuna nykyisin käytössä olevaan kaukolämpöön. Työn teoriaosuudessa tehdään katsaus rakennusten energiatehokkuuteen vaikuttaviin seikkoihin, lähienergian tuotantoon ja energiatehokkuuden parantamiseen liittyvään lainsäädäntöön ja säädöksiin. Useasta rakennuksesta koostuva kansanopisto tarjoaa mielenkiintoisen pohjan selvitystyölle ja suuri lämmitystehontarve yhdistettynä monille saneerauskohteille tyypillisiin ahtaisiin teknisiin tiloihin asettaa rajoituksia lämmitysjärjestelmän suunnittelulle. Soveltavassa osuudessa määritellään reunaehdot mahdolliselle kaukolämmön korvaavalle lämmitysratkaisulle. Tutkitaan vesistölämmön hyödyntämisen mahdollisuutta ja lasketaan aurinkosähkön ja -lämmön tuotantopotentiaalia. Maalämpöjärjestelmän mitoituksessa ja taloudellisessa vertailussa käytettiin apuna maalämpöjärjestelmiä toimittavia yrityksiä. Työssä saatujen tulosten perusteella maalämpöjärjestelmä on taloudellisesti kannattava isossa kohteessa, tosin järjestelmän asennukseen liittyy ahtaiden tilojen johdosta ongelmia. Maalämpö on myös selvästi ekologisempi, kuin nykyisin käytössä oleva kaukolämpö. Aurinkosähkön tuotannolle on Kymenlaakson Opistolla hyvä potentiaali ja sähkön tuotanto kohtaa hyvin sähkön käytön.

Työkoneen levypyörän kantavuus

Työkoneen vanteen kantavuuslaskentaan käytettävän taulukon epäiltiin tuottavan todellista huonompia kantavuuksia etenkin suurten offset-arvon vanteiden tapauksessa. Tutkimuksessa yhdestä tyypillistä kokoa edustavasta vannetyypistä valmistettiin offset-säädettävä erikoisvanne. Vanteen käytön aikaisia rasituksia seurattiin venymäliuska-antureilla koeajoradalla tavanomaisissa ajotilanteissa eri offset-arvoilla traktorilla vedettävän kuormitetun testivaunun avulla. Saatuja tuloksia vertailtiin vanteen elementtimallin tuloksiin ja käytettyjä kuormitusolettamia muutettiin. Tutkimuksen tuloksena kantavuuden laskentaan käytettävän taulukon todettiin toimivan pääosin mittauksen mukaisesti. Vanteen offset-arvolla on laskentataulukon mukaisesti merkittävä vaikutus vanteen keskiön väsymiskestävyyteen. Suurilla negatiivisilla offset-arvoilla vanteen kestoiän määrä ajo taisella alustalla, kun taas suurten positiivisten vanteiden väsymisvaurio aiheutuu pääosin kaltevalla pinnalla ajosta. Mittauksella pystyttiin osoittamaan laskentataulukon konservatiivinen olettama etenkin suurten negatiivisten offset-arvojen vanteille, jolla kyseenomaisten vanteiden väsymiskestoikää voidaan parantaa merkittävästi. Lisäksi laskentaohjelman tekemisessä käytettyjen kuormitus- ja reunaehto-olettamien vahvistettiin toimivan keskiön kannalta riittävällä tarkkuudella. Havainnolla on tärkeä rooli tulevissa elementtimenetelmään perustuvissa kantavuuslaskelmissa.

When job dissatisfaction leads to customer-oriented citizenship behaviors

This thesis places boundary conditions on the withdrawal model in the frontline setting of service organizations by considering continuance commitment and supervisory support as moderators of the relationship between job dissatisfaction and customer-oriented citizenship behaviors (COCBs). Departing from traditional research in the areas of the service-profit chain and employee withdrawal, the author advances our understanding of conditions that may lead frontline service employees who are dissatisfied to deposit COCBs into the organizational system. Specifically, based on principles derived from social exchange theory, high continuance commitment and high supervisory support are expected to lead to COCBs, because under this condition the benefits of performing such behaviors are increased (i.e., promotion-based, reciprocity-based), while the costs are decreased (i.e., opportunity costs). Utilizing a sample of 127 frontline employees from both the financial services and travel agency industries, the hypothesized relationships are empirically supported using moderated hierarchical regression analysis. To conclude discussion, implications of the results for both academics and p