Prosessihyötysuhteen parantamiskohteiden kartoitus painevesireaktorityyppisessä ydinvoimalaitoksessa

Työn tavoitteena on kartoittaa painevesireaktorityyppisen ydinvoimalaitoksen prosessihyötysuhteen parantamiskohteita. Aluksi kirjallisuudesta etsitään hyötysuhteen parantamiskeinoja ideaalisessa höyryvoimalaitosprosessissa. Näistä valitaan sopivimmat tarkastelun kohteeksi todellisessa voimalaitoksessa: syöttöveden esilämmityksen tehostaminen väliottohöyryvirtausta kasvattamalla ja syöttöveden esilämmittimen lämmönsiirtopintaa lisäämällä. Tarkastelussa pyritään löytämään paras mahdollinen hyötysuhde väliottohöyrylinjojen putkikokoa sekä esilämmittimien putkien lukumäärää muuttamalla. Diskreetin optimoinnin iteraatioaskel määritetään hyötysuhteen osittaisderivaattojen avulla. Tehtäviä muutoksia simuloidaan APROS-simulointiohjelmalla, jossa käytetään Loviisan voimalaitoksesta tehtyä mallia VVER-440. Työssä havaittiin, että pelkkiä väliottohöyrylinjojen putkikokoja – ja massavirtaa – kasvattamalla Loviisan voimalaitoksen hyötysuhdetta voidaan parantaa parhaimmillaan 32,75%:sta 32,85%:iin. Syöttöveden esilämmittimien lämmönsiirtopintaa lisäämällä saadaan suurempi parannus hyötysuhteeseen: 32,75%:sta 32,99%:iin. Näissä tapauksissa muutettiin kaikkia väliottohöyrylinjoja tai syöttöveden esilämmittimien lämpöpintoja. Työssä tarkasteltiin myös joitakin pienempiä muutoskohteita, joista paras hyötysuhteen kasvu saatiin korkeapaine-esilämmittimien lämmönsiirtopintaa kasvattamalla sekä toisen väliottohöyrylinjan (RD12) ja sitä vastaavan syöttöveden esilämmittimen muutosten yhteisvaikutuksena.

Simulation of Boundary Layers and Heat Transfer in the Entrance Region of Pipes

The study of fluid flow in pipes is one of the main topic of interest for engineers in industries. In this thesis, an effort is made to study the boundary layers formed near the wall of the pipe and how it behaves as a resistance to heat transfer. Before few decades, the scientists used to derive the analytical and empirical results by hand as there were limited means available to solve the complex fluid flow phenomena. Due to the increase in technology, now it has been practically possible to understand and analyze the actual fluid flow in any type of geometry. Several methodologies have been used in the past to analyze the boundary layer equations and to derive the expression for heat transfer. An integral relation approach is used for the analytical solution of the boundary layer equations and is compared with the FLUENT simulations for the laminar case. Law of the wall approach is used to derive the empirical correlation between dimensionless numbers and is then compared with the results from FLUENT for the turbulent case. In this thesis, different approaches like analytical, empirical and numerical are compared for the same set of fluid flow equations.

Ash Particle Erosion on Steam Boiler Convective Section

In this study, equations for the calculation of erosion wear caused by ash particles on convective heat exchanger tubes of steam boilers are presented. Anew, three-dimensional test arrangement was used in the testing of the erosion wear of convective heat exchanger tubes of steam boilers. When using the sleeve-method, three different tube materials and three tube constructions could be tested. New results were obtained from the analyses. The main mechanisms of erosionwear phenomena and erosion wear as a function of collision conditions and material properties have been studied. Properties of fossil fuels have also been presented. When burning solid fuels, such as pulverized coal and peat in steam boilers, most of the ash is entrained by the flue gas in the furnace. In bubbling andcirculating fluidized bed boilers, particle concentration in the flue gas is high because of bed material entrained in the flue gas. Hard particles, such as sharp edged quartz crystals, cause erosion wear when colliding on convective heat exchanger tubes and on the rear wall of the steam boiler. The most important ways to reduce erosion wear in steam boilers is to keep the velocity of the flue gas moderate and prevent channelling of the ash flow in a certain part of the cross section of the flue gas channel, especially near the back wall. One can do this by constructing the boiler with the following components. Screen plates can beused to make the velocity and ash flow distributions more even at the cross-section of the channel. Shield plates and plate type constructions in superheaters can also be used. Erosion testing was conducted with three types of tube constructions: a one tube row, an inline tube bank with six tube rows, and a staggered tube bank with six tube rows. Three flow velocities and two particle concentrations were used in the tests, which were carried out at room temperature. Three particle materials were used: quartz, coal ash and peat ash particles. Mass loss, diameter loss and wall thickness loss measurements of the test sleeves were taken. Erosion wear as a function of flow conditions, tube material and tube construction was analyzed by single-variable linear regression analysis. In developing the erosion wear calculation equations, multi-variable linear regression analysis was used. In the staggered tube bank, erosion wear had a maximum value in a tube row 2 and a local maximum in row 5. In rows 3, 4 and 6, the erosion rate was low. On the other hand, in the in-line tube bank the minimum erosion rate occurred in tube row 2 and in further rows the erosion had an increasing value, so that in a six row tube bank, the maximum value occurred in row 6.

Lämmöntalteenottokattilan dynaaminen simulointimalli

Työssä mallinnettiin kombivoimalaitoksen lämmöntalteenottokattila Apros-simulointiohjelmalla. Simulointimalli valmistettiin vastaamaan Helsingin Energian Vuosaari B:n voimalaitoksen lämmöntalteenottokattilaa, joka toimii kahdella painetasolla. Kattila on Foster Wheelerin valmistama. Ennen mallinnuksen aloittamista tutustuttiin laitoksen termodynamiikkaan, jolloin saatiin riittävä teoreettinen tieto koko laitoksen toiminnasta. Kattilan reunaehtoina ovat kaasuturbiiniprosessi ja laitoksen höyrykierto. Kaasuturbiini korvattiin laskentayhtälöillä, jotka antavat alkuarvot mm. savukaasun massavirralle ja lämpötilalle ennen kattilaa kaasuturbiinin tehon funktiona. Kattila liitetään höyrykiertoon tuorehöyry- ja syöttövesilinjasta, jolloin reunaehtoina annetaan lämpötilat ja paineet massavirroille. Valmistettua mallia testattiin ylösajo- ja kuormanmuutostilanteessa. Ylösajotilanteessa saatuja laskentatuloksia verrattiin todellisen laitoksen mittaustuloksiin, jolloin varmistuttiin simulointimallin oikeasta fysikaalisesta toiminnasta. Kuormanmuutostilanteissa kaasuturbiinin tehoa muutettiin ja samalla seurattiin kattilan reagointia muutostilanteessa. Kuormanmuutosmittauksessa varmistettiin vielä, että kattila reagoi kuormanmuutokseen oikealla tavalla, eikä muutos aiheuta kattilan toiminnalle haitallista värähtelyä.

Fluid flow in T-junction of pipes

The aim of this work is to study flow properties at T-junction of pipe, pressure loss suffered by the flow after passing through T-junction and to study reliability of the classical engineering formulas used to find head loss for T-junction of pipes. In this we have compared our results with CFD software packages with classical formula and made an attempt to determine accuracy of the classical formulas. In this work we have studies head loss in T-junction of pipes with various inlet velocities, head loss in T-junction of pipes when the angle of the junction is slightly different from 90 degrees and T-junction with different area of cross-section of the main pipe and branch pipe. In this work we have simulated the flow at T-junction of pipe with FLUENT and Comsol Multiphysics and observed flow properties inside the T-junction and studied the head loss suffered by fluid flow after passing through the junction. We have also compared pressure (head) losses obtained by classical formulas by A. Vazsonyi and Andrew Gardel and formulas obtained by assuming T-junction as combination of other pipe components and observations obtained from software experiments. One of the purposes of this study is also to study change in pressure loss with change in angle of T-junction. Using software we can have better view of flow inside the junction and study turbulence, kinetic energy, pressure loss etc. Such simulations save a lot of time and can be performed without actually doing the experiment. There were no real life experiments made, the results obtained completely rely on accuracy of software and numerical methods used.

Modelling of the PWR PACTEL vertical steam generator with the TRACE code

The behavior of the nuclear power plants must be known in all operational situations. Thermal hydraulics computer applications are used to simulate the behavior of the plants. The computer applications must be validated before they can be used reliably. The simulation results are compared against the experimental results. In this thesis a model of the PWR PACTEL steam generator was prepared with the TRAC/RELAP Advanced Computational Engine computer application. The simulation results can be compared against the results of the Advanced Process Simulator analysis software in future. Development of the model of the PWR PACTEL vertical steam generator is introduced in this thesis. Loss of feedwater transient simulation examples were carried out with the model.

Modeling of water/steam circulation in Circulating Fluidized Bed boiler

A distinctive design feature of steam boiler with natural circulation is the presence of the steam drum which plays a role of the separator of vapor from the flow of water-and-steam mixture coming into steam drum from the furnace tubes. Steam drum with unheated downcomer tubes, deducing from it, and riser (screen/furnace tubes) inside the furnace is a closed circulation loop in which movement of water (downcomer tubes) and water-and-steam mixture (riser tubes) is organized. The movement of the working fluid is appears due to occurrence of the natural pressure, determined by the difference in hydrostatic pressure and the mass of water and water-and-steam mixtures in downcomer and riser tubes and called the driving pressure of the natural circulation:

Modeling of water/steam circulation in Circulating Fluidized Bed boiler

A distinctive design feature of steam boiler with natural circulation is the presence of the steam drum which plays a role of the separator of vapor from the flow of water-and-steam mixture coming into steam drum from the furnace tubes. Steam drum with unheated downcomer tubes, deducing from it, and riser (screen/furnace tubes) inside the furnace is a closed circulation loop in which movement of water (downcomer tubes) and water-and-steam mixture (riser tubes) is organized. The movement of the working fluid is appears due to occurrence of the natural pressure, determined by the difference in hydrostatic pressure and the mass of water and water-and-steam mixtures in downcomer and riser tubes and called the driving pressure of the natural circulation: S drive = H steam (ρ down + ρ mix) g where: ρ down - density of water in downcomer tubes; ρ mix - density of water in riser tubes; H steam - height of steam content section; g - acceleration of gravity. In steam boilers with natural circulation the circulation rate is usually between 10 and 30. Thus, consumption of water in the circulation circuit “circulation rate times” more than steam output of the boiler. There are two aspects of the design of natural water circulation loops. One is to ensure a sufficient mass flux of circulating water to avoid burnout of evaporator tubes. The other is to avoid tube wall temperature fluctuation and tube vibration due to oscillation of circulation velocity. The design criteria are therefore reduced, in principle, to those of critical heat flux, critical flow rate for burnout, and flow instability. In practical design, however, the circulation velocity and the void fraction at the evaporator tube outlet are used as the design criteria (Seikan I., et. al., 1999). This study has been made with assumption that the heat flux in the furnace of the boiler even all the time. The target of the study was to define the circulation rate of the boiler, thus average heat flux do not change it. I would like to acknowledge professionals from “Foster Wheeler” company for good and comfortable cooperation.

Vesihöyrypiirin kiertolaskenta kuplapetikattilassa

Työssä selvitetään kiertolaskennan periaatteita kuplapetikattilassa, esitetään lyhyesti kattilan toimintaperiaate ja paneudutaan alan laskentaohjelmistoihin. Luonnonkierto kattilan vesihöyrypiirissä on seurausta hydrostaattisesta paineesta, joka aiheutuu tiheyserosta nousu- ja laskuputkien välillä. Kiertolaskennassa on huomioitava kaksifaasivirtauksen ominaispiirteet. Höyry ja neste virtaavat putkistossa eri nopeuksilla, jolloin esimerkiksi painehäviön määrityksessä käytetään erilaisia korrelaatioita ja käyrästöjä. Kaksifaasivirtauksen laskennassa tarvitaan kolmea eri taseyhtälöä: energiatasetta, massatasetta ja liikemäärätasetta. Luonnonkiertokattiloissa höyrykierron suunnittelussa on kaksi pääasiallista ehtoa. Ensimmäiseksi tulee varmistaa riittävä kiertoveden massavirta, jotta vältetään höyrystinputkien puhki palaminen. Toiseksi tulee välttää höyrystinputken pinnan lämpötilavaihteluita ja värähtelyitä. Alustavassa luonnonkierron mitoituksessa turvaudutaan kokemusperäiseen tietoon. Myöhemmässä tarkemmassa mitoituksessa käytettäviä ohjelmistoja ovat NOWA sekä kaupalliset PPSD ja Apros. Laskenta lähtee liikkeelle siitä, että lasketaan ensin vesihöyrykierron massavirrat erilaisilla lämpökuormilla ja höyryntuotanto määritellään painehäviöiden perusteella. NOWA- ja PPSD- ohjelmistoilla tehtyjen esimerkkilaskelmien perusteella voidaan sanoa, että tulokset riippuvat käytetystä laskentamallista.

Integration of Torrefaction with Steam Power Plant

Torrefaction is one of the pretreatment technologies to enhance the fuel characteristics of biomass. The efficient and continuous operation of a torrefaction reactor, in the commercial scale, demands a secure biomass supply, in addition to adequate source of heat. Biorefinery plants or biomass-fuelled steam power plants have the potential to integrate with the torrefaction reactor to exchange heat and mass, using available infrastructure and energy sources. The technical feasibility of this integration is examined in this study. A new model for the torrefaction process is introduced and verified by the available experimental data. The torrefaction model is then integrated in different steam power plants to simulate possible mass and energy exchange between the reactor and the plants. The performance of the integrated plant is investigated for different configurations and the results are compared.

Comparison of energy balances of steam explosion and torrefaction technologies

The present study introduce two pretreatment technologies which are torrefaction and steam explosion, and compare energy balance for both technologies to investigate and compare the use of these technologies to improve pelletization. In this research, torrefaction and steam explosion pretreatments were accomplished on the mixed small diameter wood (70%) with moisture content of 40 %, and logging residues (30%) with moisture content of 45 % at temperature 230 ̊C, and treatment duration 10 min. Competing methods were evaluated, and the results showed higher volumetric energy for steam explosion pellet than torrefied pellet.

Droplet Deposition in the Last Stage of Steam Turbine

During the expansion of steam in turbine, the steam crosses the saturation line and hence subsequent turbine stages run under wet condition. The stages under wet condition run with low efficiency as compared to stages running with supersaturated steam and the life of the last stage cascade is reduced due to erosion. After the steam crosses the saturation line it does not condense immediately but instead it becomes supersaturated which is a meta-stable state and reversion of equilibrium results in the formation of large number of small droplets in the range of 0.05 - 1 μm. Although these droplets are small enough to follow the stream lines of vapor however some of the fog droplets are deposited on the blade surface. After deposition they coagulate into films and rivulets which are then drawn towards the trailing edge of the blade due to viscous drag of the steam. These large droplets in the range of radius 100 μm are accelerated by steam until they impact on the next blade row causing erosion. The two phenomenon responsible for deposition are inertial impaction and turbulent-diffusion. This work shall discuss the deposition mechanism in steam turbine in detail and numerically model and validate with practical data.

Computational modelling of non-equilibrium condensing steam flows in low-pressure steam turbines

The steam turbines play a significant role in global power generation. Especially, research on low pressure (LP) steam turbine stages is of special importance for steam turbine man- ufactures, vendors, power plant owners and the scientific community due to their lower efficiency than the high pressure steam turbine stages. Because of condensation, the last stages of LP turbine experience irreversible thermodynamic losses, aerodynamic losses and erosion in turbine blades. Additionally, an LP steam turbine requires maintenance due to moisture generation, and therefore, it is also affecting on the turbine reliability. Therefore, the design of energy efficient LP steam turbines requires a comprehensive analysis of condensation phenomena and corresponding losses occurring in the steam tur- bine either by experiments or with numerical simulations. The aim of the present work is to apply computational fluid dynamics (CFD) to enhance the existing knowledge and understanding of condensing steam flows and loss mechanisms that occur due to the irre- versible heat and mass transfer during the condensation process in an LP steam turbine. Throughout this work, two commercial CFD codes were used to model non-equilibrium condensing steam flows. The Eulerian-Eulerian approach was utilised in which the mix- ture of vapour and liquid phases was solved by Reynolds-averaged Navier-Stokes equa- tions. The nucleation process was modelled with the classical nucleation theory, and two different droplet growth models were used to predict the droplet growth rate. The flow turbulence was solved by employing the standard k-ε and the shear stress transport k-ω turbulence models. Further, both models were modified and implemented in the CFD codes. The thermodynamic properties of vapour and liquid phases were evaluated with real gas models. In this thesis, various topics, namely the influence of real gas properties, turbulence mod- elling, unsteadiness and the blade trailing edge shape on wet-steam flows, are studied with different convergent-divergent nozzles, turbine stator cascade and 3D turbine stator-rotor stage. The simulated results of this study were evaluated and discussed together with the available experimental data in the literature. The grid independence study revealed that an adequate grid size is required to capture correct trends of condensation phenomena in LP turbine flows. The study shows that accurate real gas properties are important for the precise modelling of non-equilibrium condensing steam flows. The turbulence modelling revealed that the flow expansion and subsequently the rate of formation of liquid droplet nuclei and its growth process were affected by the turbulence modelling. The losses were rather sensitive to turbulence modelling as well. Based on the presented results, it could be observed that the correct computational prediction of wet-steam flows in the LP turbine requires the turbulence to be modelled accurately. The trailing edge shape of the LP turbine blades influenced the liquid droplet formulation, distribution and sizes, and loss generation. The study shows that the semicircular trailing edge shape predicted the smallest droplet sizes. The square trailing edge shape estimated greater losses. The analysis of steady and unsteady calculations of wet-steam flow exhibited that in unsteady simulations, the interaction of wakes in the rotor blade row affected the flow field. The flow unsteadiness influenced the nucleation and droplet growth processes due to the fluctuation in the Wilson point.

A Numerical study of fluid flow in peristaltic pumps and polymeric elastic pipes

This thesis work deals with a mathematical description of flow in polymeric pipe and in a specific peristaltic pump. This study involves fluid-structure interaction analysis in presence of complex-turbulent flows treated in an arbitrary Lagrangian-Eulerian (ALE) framework. The flow simulations are performed in COMSOL 4.4, as 2D axial symmetric model, and ABAQUS 6.14.1, as 3D model with symmetric boundary conditions. In COMSOL, the fluid and structure problems are coupled by monolithic algorithm, while ABAQUS code links ABAQUS CFD and ABAQUS Standard solvers with single block-iterative partitioned algorithm. For the turbulent features of the flow, the fluid model in both codes is described by RNG k-ϵ. The structural model is described, on the basis of the pipe material, by Elastic models or Hyperelastic Neo-Hookean models with Rayleigh damping properties. In order to describe the pulsatile fluid flow after the pumping process, the available data are often defective for the fluid problem. Engineering measurements are normally able to provide average pressure or velocity at a cross-section. This problem has been analyzed by McDonald's and Womersley's work for average pressure at fixed cross section by Fourier analysis since '50, while nowadays sophisticated techniques including Finite Elements and Finite Volumes exist to study the flow. Finally, we set up peristaltic pipe simulations in ABAQUS code, by using the same model previously tested for the fl uid and the structure.

Vaneri- ja kertopuutehtaiden kuivatuskoneiden lämmöntalteenottojen järjestelmäselvitys

Vanerin tai kertopuun valmistusprosessissaviilun kuivaukseen käytetään suurin osa koko valmistusprosessin primäärienergiasta. Viilunkuivauskoneessa viilun sisältämä vesi siirretään tyypillisesti prosessihöyryllä lämmitettyyn viilunkuivaajan kiertoilmaan höyrystämällä ja poistetaanviilunkuivaajasta poistoilman mukana. Viilunkuivaajan poistoilma on lämmintä jaerittäin suuren kosteuspitoisuutensa takia sisältää runsaasti energiaa. Tyypillisellä viilunkuivaajalla poistoilmaan sitoutunut lämpöteho vaihtelee prosessiolosuhteista riippuen välillä 2,7-5,7 MW. Diplomityössä tutkittiin viilunkuivaajan poistoilman sisältämän lämmön talteenottoa laitteistolla, johon kuuluu lämmöntalteenottopesuri, jossa poistoilmalla lämmitetään tuotantolaitoksen tukkipuun hautomon kiertovettä sekä ilma-ilma-lämmönsiirrin, jolla lämmitetään pesurista poistuvan ilman jäännöslämmöllä ulkoilmaa tehdassalin tuloilmakäyttöön. Työn tavoitteena oli kehittää lämmöntalteenottojärjestelmän suunnittelua, mitoitusta ja ajotapoja. Työssä analysoitiin teoreettisesti pesuria ja ilmalämmönsiirrintä, kehitettiin lämmöntalteenottopesurin simulointimenetelmä ja mitattiin toiminnassa olevia talteenottolaitteistoja. Tutkimuksessa todettiin lämmöntalteenottohyötysuhteen vaihtelevan lämmityskaudella välillä 50-70 %. Lämmöntalteenottolaitteiston pesurin veteen saatava teho riippuu ensisijaisesti viilunkuivaajan poistoilman lämpösisällöstä, joka on enimmäkseen kosteusriippuvainen ja ilmanvaihtoilmaan saatava teho ulkolämpö-tilan määräämästä tehontarpeesta. Pesurin vesijärjestelmän vaikutusmekanismit pesurin suorituskykyyn tunnistettiin ja niiden pohjalta annetaan suositukset mitoitukseen ja ajotapaan. Lämmöntalteenottolaitteiston lämpötehon tasapainottamiseen pesurin ja ilma-ilma-lämmönsiirtimen välillä mitoituksen avulla esitellään työkalut.