Roal Oy:n käyttöhyödykejärjestelmän arviointi, optimointi ja raportointi

Prosessiteollisuudessa tarvitaan usein erilaisia apujärjestelmiä pääprosessin tueksi. Tyypillisiä tällaisia järjestelmiä ovat jäähdytys-, höyry- ja ilmajärjestelmät. Hyödykejärjestelmien kehitys jää helposti pääprosessin varjoon, joka usein johtaa tarpeettoman suuriin hyödykekustannuksiin ja järjestelmien teknisen tilan laskuun. Työn kirjallisuusosassa käsitellään Roal Oy:n fermentointiprosessin kannalta olennaisimpia hyödykevaatimuksia, niiden laskennallista ilmaisua ja vuorovaikutuksia sekä PINCH-menetelmää lämpöenergian hyötykäyttöön. Jäähdytysjärjestelmän osalta käydään läpi merkittävimmät laitetekniset ratkaisut, jäähdytystorni ja lämpöpumppu, toimintaperiaatteineen sekä luonnonvesien käyttö jäähdytykseen. Työn soveltavassa osassa seurattiin Aspergillus, Trichoderma ja Bacillus fermentointeja, joiden pohjalta luotiin kasvatuskohtainen empiirinen malli jäähdytystarpeen arviointiin perustuen sekoitustehoon, kasvatuksen hiilidioksidituottoon ja haihtumisen vaikutukseen. Kasvatuksien aikana seurattiin myös tilavuusperusteista lämmönsiirtokerrointa. Mitattujen lämmönsiirtokertoimien perustella laskettiin ominaislämmöntuottoon perustuva maksimilämpötila käytettävälle jäähdytysvedelle ja fermentorien maksimitilavuudet tunnetuilla kasvatusparametreilla eri lämpöisille jäähdytysvesille. Soveltavassa osassa käydään myös läpi Roalin höyry- ja kuumavesikulutukset ja tärkeimmät käyttökohteet. Mittaustulosten ja mallien perusteella on tehtiin kehitysehdotukset hyödykejärjestelmän optimoimiseksi.

Modelling of combustion and sorbent reactions in three-dimensional flow environment of a circulating fluidized bed furnace

This thesis presents a three-dimensional, semi-empirical, steady state model for simulating the combustion, gasification, and formation of emissions in circulating fluidized bed (CFB) processes. In a large-scale CFB furnace, the local feeding of fuel, air, and other input materials, as well as the limited mixing rate of different reactants produce inhomogeneous process conditions. To simulate the real conditions, the furnace should be modelled three-dimensionally or the three-dimensional effects should be taken into account. The only available methods for simulating the large CFB furnaces three-dimensionally are semi-empirical models, which apply a relatively coarse calculation mesh and a combination of fundamental conservation equations, theoretical models and empirical correlations. The number of such models is extremely small. The main objective of this work was to achieve a model which can be applied to calculating industrial scale CFB boilers and which can simulate all the essential sub-phenomena: fluid dynamics, reactions, the attrition of particles, and heat transfer. The core of the work was to develop the model frame and the required sub-models for determining the combustion and sorbent reactions. The objective was reached, and the developed model was successfully used for studying various industrial scale CFB boilers combusting different types of fuel. The model for sorbent reactions, which includes the main reactions for calcitic limestones, was applied for studying the new possible phenomena occurring in the oxygen-fired combustion. The presented combustion and sorbent models and principles can be utilized in other model approaches as well, including other empirical and semi-empirical model approaches, and CFD based simulations. The main achievement is the overall model frame which can be utilized for the further development and testing of new sub-models and theories, and for concentrating the knowledge gathered from the experimental work carried out at bench scale, pilot scale and industrial scale apparatus, and from the computational work performed by other modelling methods.

Design Aspects of Megawatt-Range Direct-Driven Permanent MagnetWind Generators

Direct-driven permanent magnet synchronous generator is one of the most promising topologies for megawatt-range wind power applications. The rotational speed of the direct-driven generator is very low compared with the traditional electrical machines. The low rotational speed requires high torque to produce megawatt-range power. The special features of the direct-driven generators caused by the low speed and high torque are discussed in this doctoral thesis. Low speed and high torque set high demands on the torque quality. The cogging torque and the load torque ripple must be as low as possible to prevent mechanical failures. In this doctoral thesis, various methods to improve the torque quality are compared with each other. The rotor surface shaping, magnet skew, magnet shaping, and the asymmetrical placement of magnets and stator slots are studied not only by means of torque quality, but also the effects on the electromagnetic performance and manufacturability of the machine are discussed. The heat transfer of the direct-driven generator must be designed to handle the copper losses of the stator winding carrying high current density and to keep the temperature of the magnets low enough. The cooling system of the direct-driven generator applying the doubly radial air cooling with numerous radial cooling ducts was modeled with a lumped-parameter-based thermal network. The performance of the cooling system was discussed during the steady and transient states. The effect of the number and width of radial cooling ducts was explored. The large number of radial cooling ducts drastically increases the impact of the stack end area effects, because the stator stack consists of numerous substacks. The effects of the radial cooling ducts on the effective axial length of the machine were studied by analyzing the crosssection of the machine in the axial direction. The method to compensate the magnet end area leakage was considered. The effect of the cooling ducts and the stack end area effects on the no-load voltages and inductances of the machine were explored by using numerical analysis tools based on the three-dimensional finite element method. The electrical efficiency of the permanent magnet machine with different control methods was estimated analytically over the whole speed and torque range. The electrical efficiencies achieved with the most common control methods were compared with each other. The stator voltage increase caused by the armature reaction was analyzed. The effect of inductance saturation as a function of load current was implemented to the analytical efficiency calculation.

Modeling and simulation of forced-air cooling of strawberries using variable convective coefficient

In the forced-air cooling process of fruits occurs, besides the convective heat transfer, the mass transfer by evaporation. The energy need in the evaporation is taken from fruit that has its temperature lowered. In this study it has been proposed the use of empirical correlations for calculating the convective heat transfer coefficient as a function of surface temperature of the strawberry during the cooling process. The aim of this variation of the convective coefficient is to compensate the effect of evaporation in the heat transfer process. Linear and exponential correlations are tested, both with two adjustable parameters. The simulations are performed using experimental conditions reported in the literature for the cooling of strawberries. The results confirm the suitability of the proposed methodology.

Technology for sugarcane agroindustry waste reuse as granulated organomineral fertilizer

Aiming to evaluate the use of sugarcane industry waste such as byproducts from vinasse concentration process, it was assessed the organomineral fertilizer BIOFOM (concentrated vinasse, filter cake, boiler ash, soot from chimneys and supplemented with mineral fertilizers). The study included characterization and agronomic potential analysis of a test plant (corn), by noting the differences between mineral fertilizers and BIOFOM fertilization until 45 days after sowing. The technology traditionally used to produce BIOFOM was based on vinasse evaporation with high heat transfer coefficients. It was observed that the technology, which can be formulated according to the needs of any crop, could be used in many cases as mineral fertilizer. Therefore, the use of this organomineral fertilizer reduces waste generation of sugarcane industry.

Ydinpolttoaineen lämmönsiirto

Kandidaatintyössä perehdytään ydinpolttoaineessa tapahtuvaan lämmönsiirtoon ja lämmönsiirron ilmiöihin. Lämmönsiirron tarkastelussa keskitytään erityisesti polttoainepelletissä tapahtuvaan lämmönsiirtoon, mutta työn edetessä esitellään myös lyhyesti lämmön siirtyminen polttoainepelletistä kaasunraon ja polttoainesauvan suojakuoren läpi jäähdytteeseen. Kandidaatintyössä tarkastellaan myös kiinteiden ydinpolttoaineiden lämmönsiirto-ominaisuuksia. Lämmönsiirto-ominaisuudet riippuvat materiaalien termodynaamisista ja kemiallisista ominaisuuksista. Lämmönsiirto-ominaisuuksien tunteminen on edellytys uusien, lämmönsiirrollisesti entistä parempien, polttoaineiden kehittämiselle.

Laskentasovelluksen kehittäminen modulaariselle vesiputkikattilalle

Modulaarisella vesiputkikattilalla tarkoitetaan täysin konepajalla valmistettavissa olevaa kattilaa, joka voidaan kuljettaa yhtenä tai muutamana suurena moduulina työmaalle. Tässä diplomityössä käsiteltiin modulaarisen vesiputkikattilan laskentasovelluksen ke-hittämistä KPA Unicon Oy:lle. Työn tavoitteena oli tarkastella modulaarisen vesiputkikattilan lämpöteknistä mitoitusta ja suunnittelua sekä kehittää laskentasovellus, jonka avulla voidaan arvioida kattilan mittoja ja painoa. Laskentasovellus laadittiin Microsoftin Excel-alustalle, josta se on myöhemmin mahdollista siirtää muille ohjelma-alustoille. Laskentasovelluksessa hyödynnetään lämmönsiirron ja virtaustekniikan laskentayhtälöitä sekä -menetelmiä. Sovellukseen valitut laskentayhtälöt sekä -menetelmät ovat yleisesti hyväksyttyjä ja käytännössä testattuja. Diplomityön tuloksena valmistui laskentasovellus, joka kykenee modulaarisen vesiputkikattilan lämpötekniseen mitoittamiseen. Sovelluksen avulla voidaan mitoittaa kattilan tulipesä, tulistimet, höyrystinpinnat sekä ekonomaiseri. Laskentasovellusta on tarkoitus hyödyntää yrityksen tarjousvaiheen projekteissa sekä mahdollisesti kattiloiden esisuunnittelussa. Laskentasovelluksen laatimista varten ei ollut mahdollista hyödyntää toiminnassa olevien kattiloiden prosessitietoja, koska ensimmäiset tämän tyyppiset kattilat ovat asenteilla. Sen sijaan sovelluksen antamia laskentatuloksia verrattiin toisen mitoitusohjelman antamiin tuloksiin, joiden perusteella laskentasovelluksen voidaan olettaa antavan oikeita tuloksia.

Lämmön- ja aineensiirto hiilipartikkelin palamisessa

Tässä tutkimusraportissa käsitellään yksittäisen hiilipartikkelin palamisen teoriaa. Tutkimusraportissa tullaan esittelemään palamisen kannalta olennaiset ilmiöt kuten lämmön- ja aineensiirto sekä palamiseen liittyvät reaktiot ja palamisen vaiheet. Lisäksi tarkastellaan eri aineominaisuuksien riippuvuutta eri parametreista. Käsitellyn teorian pohjalta luodaan palamista kuvaava malli, jonka avulla kuvataan aineensiirtoa ja lämpötilajakaumaa yksittäisessä partikkelissa. Tavoitteena on mallin avulla tarkastella partikkelikoon vaikutusta yllä esitettyihin ilmiöihin.

Production of Mg(OH)2 from Mg-silicate rock for CO2 mineral sequestration

Sequestration of carbon dioxide in mineral rocks, also known as CO2 Capture and Mineralization (CCM), is considered to have a huge potential in stabilizing anthropogenic CO2 emissions. One of the CCM routes is the ex situ indirect gas/sold carbonation of reactive materials, such as Mg(OH)2, produced from abundantly available Mg-silicate rocks. The gas/solid carbonation method is intensively researched at Åbo Akademi University (ÅAU ), Finland because it is energetically attractive and utilizes the exothermic chemistry of Mg(OH)2 carbonation. In this thesis, a method for producing Mg(OH)2 from Mg-silicate rocks for CCM was investigated, and the process efficiency, energy and environmental impact assessed. The Mg(OH)2 process studied here was first proposed in 2008 in a Master’s Thesis by the author. At that time the process was applied to only one Mg-silicate rock (Finnish serpentinite from the Hitura nickel mine site of Finn Nickel) and the optimum process conversions, energy and environmental performance were not known. Producing Mg(OH)2 from Mg-silicate rocks involves a two-staged process of Mg extraction and Mg(OH)2 precipitation. The first stage extracts Mg and other cations by reacting pulverized serpentinite or olivine rocks with ammonium sulfate (AS) salt at 400 - 550 oC (preferably < 450 oC). In the second stage, ammonia solution reacts with the cations (extracted from the first stage after they are leached in water) to form mainly FeOOH, high purity Mg(OH)2 and aqueous (dissolved) AS. The Mg(OH)2 process described here is closed loop in nature; gaseous ammonia and water vapour are produced from the extraction stage, recovered and used as reagent for the precipitation stage. The AS reagent is thereafter recovered after the precipitation stage. The Mg extraction stage, being the conversion-determining and the most energy-intensive step of the entire CCM process chain, received a prominent attention in this study. The extraction behavior and reactivity of different rocks types (serpentinite and olivine rocks) from different locations worldwide (Australia, Finland, Lithuania, Norway and Portugal) was tested. Also, parametric evaluation was carried out to determine the optimal reaction temperature, time and chemical reagent (AS). Effects of reactor types and configuration, mixing and scale-up possibilities were also studied. The Mg(OH)2 produced can be used to convert CO2 to thermodynamically stable and environmentally benign magnesium carbonate. Therefore, the process energy and life cycle environmental performance of the ÅAU CCM technique that first produces Mg(OH)2 and the carbonates in a pressurized fluidized bed (FB) were assessed. The life cycle energy and environmental assessment approach applied in this thesis is motivated by the fact that the CCM technology should in itself offer a solution to what is both an energy and environmental problem. Results obtained in this study show that different Mg-silicate rocks react differently; olivine rocks being far less reactive than serpentinite rocks. In summary, the reactivity of Mg-silicate rocks is a function of both the chemical and physical properties of rocks. Reaction temperature and time remain important parameters to consider in process design and operation. Heat transfer properties of the reactor determine the temperature at which maximum Mg extraction is obtained. Also, an increase in reaction temperature leads to an increase in the extent of extraction, reaching a maximum yield at different temperatures depending on the reaction time. Process energy requirement for producing Mg(OH)2 from a hypothetical case of an iron-free serpentine rock is 3.62 GJ/t-CO2. This value can increase by 16 - 68% depending on the type of iron compound (FeO, Fe2O3 or Fe3O4) in the mineral. This suggests that the benefit from the potential use of FeOOH as an iron ore feedstock in iron and steelmaking should be determined by considering the energy, cost and emissions associated with the FeOOH by-product. AS recovery through crystallization is the second most energy intensive unit operation after the extraction reaction. However, the choice of mechanical vapor recompression (MVR) over the “simple evaporation” crystallization method has a potential energy savings of 15.2 GJ/t-CO2 (84 % savings). Integrating the Mg(OH)2 production method and the gas/solid carbonation process could provide up to an 25% energy offset to the CCM process energy requirements. Life cycle inventory assessment (LCIA) results show that for every ton of CO2 mineralized, the ÅAU CCM process avoids 430 - 480 kg CO2. The Mg(OH)2 process studied in this thesis has many promising features. Even at the current high energy and environmental burden, producing Mg(OH)2 from Mg-silicates can play a significant role in advancing CCM processes. However, dedicated future research and development (R&D) have potential to significantly improve the Mg(OH)2 process performance.

Suurnopeus-turbokoneroottoreiden termodynaaminen ja mekaaninen mallinnus sekä rakenneanalyysi

Tässä väitöskirjassa tarkastellaan suurnopeustekniikan eri sovelluksissa ilmeneviä roottoreihin liittyviä rakenteellisia vaatimuksia ja haasteita. Tässä yhteydessä suurnopeustekniikalla tarkoitetaan järjestelyä, jossa sähkökone (moottori, generaattori) ja toimilaite (turbiini, kompressori, puhallin) on kytketty ilman vaihdetta suoraan mekaanisesti yhteen ja jossa yhteisen roottorin pyörimisnopeus on selvästi suurempi kuin 50/60 hertsin verkosta syötetyn kaksinapaisen vaihtovirtasähkökoneen tahtinopeus. Tyypillistä suurnopeuskoneen roottorille on suuri tehotiheys ja suuri mekaaninen kuormitus. Siksi esimerkiksi sähkökoneen jäähdytys on entistä haasteellisempaa kasvavien rautahäviöiden ja pienempien lämmönsiirtopinta-alojen vuoksi. Tämän työn tavoitteet voidaan jakaa kolmeen osaan:  Yhdistetyn sähkö- ja turbokoneen roottorin mekaanisen rakenteen tarkastelu, jonka tavoitteena on pienentää lämmönkehitystä ja tehostaa kriittisten kohtien jäähdytystä. Tähän liittyy sähkömagneettisten häviöiden keskittäminen jäähdytyksen kannalta edullisiin kohtiin  Yhdistetyn sähkö- ja turbokoneen roottorin mekaanisen rakenteen tarkastelu kriittisten ominaistaajuuksien kannalta  Yhdistetyn sähkö- ja turbokoneen roottorin mekaanisen rakenteen analysointi lujuustekniseltä kannalta. Tähän liittyvät mm. erilaiset ahdistussovitteet ja niiden säilyminen korkeilla pyörimisnopeuksilla sekä niiden roottoria jäykistävä vaikutus ja lämmön johtuminen kyseisissä liitospinnoissa. Tämän työn tieteellinen uutuusarvo on nimenomaan yhdistetyn sähkö- ja turbokoneen roottorin rakenteen analysointi ottamalla samanaikaisesti huomioon kaikki edellä mainitut näkökohdat: jäähtyminen erityisen kuumissa kohdissa, sähköisten häviöiden alentaminen ja niiden jakautuman huomioon ottaminen, roottorin jäykkyyden maksimointi, lujuusrasitusten hallinta ja rakenteen mekaaninen stabiliteetti sekä lämpöteknisten ylimenovastusten tarkastelu.

Theoretical Analysis and Numerical Simulation of Spectral Radiative Properties of Combustion Gases in Oxy/Air-Fired Combustion Systems

Energy efficiency is one of the major objectives which should be achieved in order to implement the limited energy resources of the world in a sustainable way. Since radiative heat transfer is the dominant heat transfer mechanism in most of fossil fuel combustion systems, more accurate insight and models may cause improvement in the energy efficiency of the new designed combustion systems. The radiative properties of combustion gases are highly wavelength dependent. Better models for calculating the radiative properties of combustion gases are highly required in the modeling of large scale industrial combustion systems. With detailed knowledge of spectral radiative properties of gases, the modeling of combustion processes in the different applications can be more accurate. In order to propose a new method for effective non gray modeling of radiative heat transfer in combustion systems, different models for the spectral properties of gases including SNBM, EWBM, and WSGGM have been studied in this research. Using this detailed analysis of different approaches, the thesis presents new methods for gray and non gray radiative heat transfer modeling in homogeneous and inhomogeneous H2O–CO2 mixtures at atmospheric pressure. The proposed method is able to support the modeling of a wide range of combustion systems including the oxy-fired combustion scenario. The new methods are based on implementing some pre-obtained correlations for the total emissivity and band absorption coefficient of H2O–CO2 mixtures in different temperatures, gas compositions, and optical path lengths. They can be easily used within any commercial CFD software for radiative heat transfer modeling resulting in more accurate, simple, and fast calculations. The new methods were successfully used in CFD modeling by applying them to industrial scale backpass channel under oxy-fired conditions. The developed approaches are more accurate compared with other methods; moreover, they can provide complete explanation and detailed analysis of the radiation heat transfer in different systems under different combustion conditions. The methods were verified by applying them to some benchmarks, and they showed a good level of accuracy and computational speed compared to other methods. Furthermore, the implementation of the suggested banded approach in CFD software is very easy and straightforward.

Validation of a 3D model for oxygen combustion in a circulating fluidized bed

In this thesis, a model called CFB3D is validated for oxygen combustion in circulating fluidized bed boiler. The first part of the work consists of literature review in which circulating fluidized bed and oxygen combustion technologies are studied. In addition, the modeling of circulating fluidized bed furnaces is discussed and currently available industrial scale three-dimensional furnace models are presented. The main features of CFB3D model are presented along with the theories and equations related to the model parameters used in this work. The second part of this work consists of the actual research and modeling work including measurements, model setup, and modeling results. The objectives of this thesis is to study how well CFB3D model works with oxygen combustion compared to air combustion in circulating fluidized bed boiler and what model parameters need to be adjusted when changing from air to oxygen combustion. The study is performed by modeling two air combustion cases and two oxygen combustion cases with comparable boiler loads. The cases are measured at Ciuden 30 MWth Flexi-Burn demonstration plant in April 2012. The modeled furnace temperatures match with the measurements as well in oxygen combustion cases as in air combustion cases but the modeled gas concentrations differ from the measurements clearly more in oxygen combustion cases. However, the same model parameters are optimal for both air and oxygen combustion cases. When the boiler load is changed, some combustion and heat transfer related model parameters need to be adjusted. To improve the accuracy of modeling results, better flow dynamics model should be developed in the CFB3D model. Additionally, more measurements are needed from the lower furnace to find the best model parameters for each case. The validation work needs to be continued in order to improve the modeling results and model predictability.

Thermal radiation in combustion systems

A numerical procedure for solving the nongray radiative transfer equation (RTE) in two-dimensional cylindrical participating media is presented. Nongray effects are treated by using a narrow-band approach. Radiative emission from CO, CO2, H2O, CH4 and soot is considered. The solution procedure is applied to study radiative heat transfer in a premixed CH4-O2, laminar, flame. Temperature, soot and IR-active species molar fraction distributions are allowed to vary in the axial direction of the flame. From the obtained results it is possible to quantify the radiative loss in the flame, as well as the importance of soot radiation as compared to gaseous radiation. Since the solution procedure is developed for a two-dimensional cylindrical geometry, it can be applied to other combustion systems such as furnaces, internal combustion engines, liquid and solid propellant combustion.

A law of the wall formulation for recirculating flows

This work presents a new law of the wall formulation for recirculating turbulent flows. An alternative expression for the internal length which can be applied in the separated region is also presented. The formulation is implemented in a numerical code which solves the k-epsilon model through a finite volume method. The theoretical results are compared with the experimental data of Vogel and Eaton (J. of Heat Transfer, Transactions of ASME, vol.107, pp. 922-929, 1985). The paper shows that the present formulation furnishes better results than the standard k-epsilon formulation.

Direct Simulation Monte Carlo Method Applied to Aerothermodynamics

An axisymmetric supersonic flow of rarefied gas past a finite cylinder was calculated applying the direct simulation Monte Carlo method. The drag force, the coefficients of pressure, of skin friction, and of heat transfer, the fields of density, of temperature, and of velocity were calculated as function of the Reynolds number for a fixed Mach number. The variation of the Reynolds number is related to the variation of the Knudsen number, which characterizes the gas rarefaction. The present results show that all quantities in the transition regime (Knudsen number is about the unity) are significantly different from those in the hydrodynamic regime, when the Knudsen number is small.