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.

Substitution-Newton-Raphson method applied to the modeling of a vapour compression refrigeration system using different representations of the thermodynamic properties of R-134a

This paper gives a detailed presentation of the Substitution-Newton-Raphson method, suitable for large sparse non-linear systems. It combines the Successive Substitution method and the Newton-Raphson method in such way as to take the best advantages of both, keeping the convergence features of the Newton-Raphson with the low requirements of memory and time of the Successive Substitution schemes. The large system is solved employing few effective variables, using the greatest possible part of the model equations in substitution fashion to fix the remaining variables, but maintaining the convergence characteristics of the Newton-Raphson. The methodology is exemplified through a simple algebraic system, and applied to a simple thermodynamic, mechanical and heat transfer modeling of a single-stage vapor compression refrigeration system. Three distinct approaches for reproducing the thermodynamic properties of the refrigerant R-134a are compared: the linear interpolation from tabulated data, the use of polynomial fitted curves and the use of functions derived from the Helmholtz free energy.

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.

Model based analysis of the post-combustion calcium looping process for carbon dioxide capture

This thesis presents a one-dimensional, semi-empirical dynamic model for the simulation and analysis of a calcium looping process for post-combustion CO2 capture. Reduction of greenhouse emissions from fossil fuel power production requires rapid actions including the development of efficient carbon capture and sequestration technologies. The development of new carbon capture technologies can be expedited by using modelling tools. Techno-economical evaluation of new capture processes can be done quickly and cost-effectively with computational models before building expensive pilot plants. Post-combustion calcium looping is a developing carbon capture process which utilizes fluidized bed technology with lime as a sorbent. The main objective of this work was to analyse the technological feasibility of the calcium looping process at different scales with a computational model. A one-dimensional dynamic model was applied to the calcium looping process, simulating the behaviour of the interconnected circulating fluidized bed reactors. The model incorporates fundamental mass and energy balance solvers to semi-empirical models describing solid behaviour in a circulating fluidized bed and chemical reactions occurring in the calcium loop. In addition, fluidized bed combustion, heat transfer and core-wall layer effects were modelled. The calcium looping model framework was successfully applied to a 30 kWth laboratory scale and a pilot scale unit 1.7 MWth and used to design a conceptual 250 MWth industrial scale unit. Valuable information was gathered from the behaviour of a small scale laboratory device. In addition, the interconnected behaviour of pilot plant reactors and the effect of solid fluidization on the thermal and carbon dioxide balances of the system were analysed. The scale-up study provided practical information on the thermal design of an industrial sized unit, selection of particle size and operability in different load scenarios.

Ilmatislaamon vesijäähdytysjärjestelmän lämmönvaihtajien tehokkuuden vaikutus tislausprosessin hyötysuhteeseen

Sähköenergiankulutuksella on suuri merkitys kryogeenisessa ilman tislausprosessissa. Sähköenergiankulutus määrittää pääosan tuotetun kaasun tuotantokustannuksista, koska raaka-aine eli ympäröivä ilmakehä ei aiheuta taloudellisia kuluja. Vesijäähdytysjärjestelmä vaikuttaa tislausprosessin tehokkuuteen merkittävästi. Jäähdytysveden lämpötilalla on oleellinen merkitys tuotettujen kaasujen massavirtoihin ja kulutettavaan energiaan. Lisäksi jäähdytysveden lämpötila vaikuttaa laitoksen käytettävyyteen. Jäähdytysjärjestelmän tehokkuuteen vaikuttavat useat eri tekijät. Lähtökohtana lämmönsiirtoon vaikuttavista tekijöistä on saatavan jäähdytysveden lämpötila. Lämmönsiirtimien mitoituksella, likaantumisella ja virtausnopeuksilla on myös oleellinen merkitys ilman tislausprosessin lämmönsiirrossa. Tässä raportissa on kuvattu kryogeeninen ilman tislausprosessi ja syvennytty prosessin vesijäähdytysjärjestelmään. Raportissa lämmönsiirto pohjautuu ainoastaan laitoksen vesilämmönsiirtimiin ulkoisessa ja sisäisessä kierrossa. Raportti sisältää prosessikuvauksen, tutkimustietoa jäähdytysjärjestelmästä ja lämmönsiirrosta, koejärjestelyjen tulokset ja päätelmiä laitoksen jäähdytysjärjestelmän ja tuotannon yhteyksistä. Raportti perustuu kattavaan lähdemateriaaliin.

Computational studies for the design of process equipment with complex geometries

This study combines several projects related to the flows in vessels with complex shapes representing different chemical apparata. Three major cases were studied. The first one is a two-phase plate reactor with a complex structure of intersecting micro channels engraved on one plate which is covered by another plain plate. The second case is a tubular microreactor, consisting of two subcases. The first subcase is a multi-channel two-component commercial micromixer (slit interdigital) used to mix two liquid reagents before they enter the reactor. The second subcase is a micro-tube, where the distribution of the heat generated by the reaction was studied. The third case is a conventionally packed column. However, flow, reactions or mass transfer were not modeled. Instead, the research focused on how to describe mathematically the realistic geometry of the column packing, which is rather random and can not be created using conventional computeraided design or engineering (CAD/CAE) methods. Several modeling approaches were used to describe the performance of the processes in the considered vessels. Computational fluid dynamics (CFD) was used to describe the details of the flow in the plate microreactor and micromixer. A space-averaged mass transfer model based on Fick’s law was used to describe the exchange of the species through the gas-liquid interface in the microreactor. This model utilized data, namely the values of the interfacial area, obtained by the corresponding CFD model. A common heat transfer model was used to find the heat distribution in the micro-tube. To generate the column packing, an additional multibody dynamic model was implemented. Auxiliary simulation was carried out to determine the position and orientation of every packing element in the column. This data was then exported into a CAD system to generate desirable geometry, which could further be used for CFD simulations. The results demonstrated that the CFD model of the microreactor could predict the flow pattern well enough and agreed with experiments. The mass transfer model allowed to estimate the mass transfer coefficient. Modeling for the second case showed that the flow in the micromixer and the heat transfer in the tube could be excluded from the larger model which describes the chemical kinetics in the reactor. Results of the third case demonstrated that the auxiliary simulation could successfully generate complex random packing not only for the column but also for other similar cases.

Modeling of temperature regimes in a furnace of bubbling fluidized bed boiler

The work is mainly focused on the technology of bubbling fluidized bed combustion. Heat transfer and hydrodynamics of the process were examined in the work in detail. Special emphasis was placed on the process of heat exchange in a freeboard zone of bubbling fluidized bed boiler. Operating mode of bubbling fluidized bed boiler depends on many parameters. To assess the influence of some parameters on a temperature regime inside the furnace a simplified method of zonal modeling was used in the work. Thus, effects of bed material fineness, excess air ratio and changes in boiler load were studied. Besides the technology of combustion in bubbling fluidized bed, other common technologies of solid fuels combustion were reviewed. In addition, brief survey of most widely used types of solid fuel was performed in the work.

Jätteen kaasutus ja energiakäyttö

Kiinteiden polttoaineiden muuntaminen kaasumaiseen muotoon eli kaasuttaminen herätti jo vuosia sitten tutkijoiden kiinnostusta. Suomessa yleinen mielenkiinto on viimeaikoina keskittynyt edullisempien polttoaineiden kaasuttamiseen ja saatavan tuotekaasun polttamiseen. Lahti Energian uudessa Kymijärvi 2 voimalaitoksessa kaasutetaan puuta ja kierrätyspolttoainetta (REF). Saatava tuotekaasu jäähdytetään ja suodatetaan, jotta epäpuhtaudet saadaan pois tuotekaasusta. Puhdistettu tuotekaasu poltetaan kaasukattilassa korkeilla höyrynarvoilla. Tuotekaasun jäähdyttämistä ei ole paljoa tehty eikä siitä ole juurikaan tieteellistä tutkimusta eikä kokemusta. Tuotekaasun jäähtyessä partikkelit tarttuvat lämpöpinnoille aiheuttaen kerrostumia. Kerrostumat heikentävät lämmönsiirtoa olennaisesti. Tämän työn tarkoitus on tutkia kaasutusprosessia, tuotekaasun jäähdyttimen likaantumista sekä antaa lisätietoja likaantumiseen vaikuttavista tekijöistä.