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.

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.

Structural and operational optimisation : applications in energy systems

In this work mathematical programming models for structural and operational optimisation of energy systems are developed and applied to a selection of energy technology problems. The studied cases are taken from industrial processes and from large regional energy distribution systems. The models are based on Mixed Integer Linear Programming (MILP), Mixed Integer Non-Linear Programming (MINLP) and on a hybrid approach of a combination of Non-Linear Programming (NLP) and Genetic Algorithms (GA). The optimisation of the structure and operation of energy systems in urban regions is treated in the work. Firstly, distributed energy systems (DES) with different energy conversion units and annual variations of consumer heating and electricity demands are considered. Secondly, district cooling systems (DCS) with cooling demands for a large number of consumers are studied, with respect to a long term planning perspective regarding to given predictions of the consumer cooling demand development in a region. The work comprises also the development of applications for heat recovery systems (HRS), where paper machine dryer section HRS is taken as an illustrative example. The heat sources in these systems are moist air streams. Models are developed for different types of equipment price functions. The approach is based on partitioning of the overall temperature range of the system into a number of temperature intervals in order to take into account the strong nonlinearities due to condensation in the heat recovery exchangers. The influence of parameter variations on the solutions of heat recovery systems is analysed firstly by varying cost factors and secondly by varying process parameters. Point-optimal solutions by a fixed parameter approach are compared to robust solutions with given parameter variation ranges. In the work enhanced utilisation of excess heat in heat recovery systems with impingement drying, electricity generation with low grade excess heat and the use of absorption heat transformers to elevate a stream temperature above the excess heat temperature are also studied.

Heat transfer inside internal combustion engine: modelling and comparison with experimental data

The paper is devoted to study specific aspects of heat transfer in the combustion chamber of compression ignited reciprocating internal combustion engines and possibility to directly measure the heat flux by means of Gradient Heat Flux Sensors (GHFS). A one – dimensional single zone model proposed by Kyung Tae Yun et al. and implemented with the aid of Matlab, was used to obtain approximate picture of heat flux behavior in the combustion chamber with relation to the crank angle. The model’s numerical output was compared to the experimental results. The experiment was accomplished by A. Mityakov at four stroke diesel engine Indenor XL4D. Local heat fluxes on the surface of cylinder head were measured with fast – response, high – sensitive GHFS. The comparison of numerical data with experimental results has revealed a small deviation in obtained heat flux values throughout the cycle and different behavior of heat flux curve after Top Dead Center.

Heat flux measurement inside internal combustion engine with gradient heat flux sensor.

This master’s thesis is devoted to study different heat flux measurement techniques such as differential temperature sensors, semi-infinite surface temperature methods, calorimetric sensors and gradient heat flux sensors. The possibility to use Gradient Heat Flux Sensors (GHFS) to measure heat flux in the combustion chamber of compression ignited reciprocating internal combustion engines was considered in more detail. A. Mityakov conducted an experiment, where Gradient Heat Flux Sensor was placed in four stroke diesel engine Indenor XL4D to measure heat flux in the combustion chamber. The results which were obtained from the experiment were compared with model’s numerical output. This model (a one – dimensional single zone model) was implemented with help of MathCAD and the result of this implementation is graph of heat flux in combustion chamber in relation to the crank angle. The values of heat flux throughout the cycle obtained with aid of heat flux sensor and theoretically were sufficiently similar, but not identical. Such deviation is rather common for this type of experiment.

Modeling of residual stresses and distortion due to welding in fillet welds

Welding has a growing role in modern world manufacturing. Welding joints are extensively used from pipes to aerospace industries. Prediction of welding residual stresses and distortions is necessary for accurate evaluation of fillet welds in relation to design and safety conditions. Residual stresses may be beneficial or detrimental, depending whether they are tensile or compressive and the loading. They directly affect the fatigue life of the weld by impacting crack growth rate. Beside theoretical background of residual stresses this study calculates residual stresses and deformations due to localized heating by welding process and subsequent rapid cooling in fillet welds. Validated methods are required for this purpose due to complexity of process, localized heating, temperature dependence of material properties and heat source. In this research both empirical and simulation methods were used for the analysis of welded joints. Finite element simulation has become a popular tool of prediction of welding residual stresses and distortion. Three different cases with and without preload have been modeled during this study. Thermal heat load set is used by calculating heat flux from the given heat input energy. First the linear and then nonlinear material behavior model is modeled for calculation of residual stresses. Experimental work is done to calculate the stresses empirically. The results from both the methods are compared to check their reliability. Residual stresses can have a significant effect on fatigue performance of the welded joints made of high strength steel. Both initial residual stress state and subsequent residual stress relaxation need to be considered for accurate description of fatigue behavior. Tensile residual stresses are detrimental and will reduce the fatigue life and compressive residual stresses will increase it. The residual stresses follow the yield strength of base or filler material and the components made of high strength steel are typically thin, where the role of distortion is emphasizing.

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.

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.

Biogas production in regional integrated biodegradable waste treatment – Possibilities for improving energy performance and reducing GHG emissions

The greatest threat that the biodegradable waste causes on the environment is the methane produced in landfills by the decomposition of this waste. The Landfill Directive (1999/31/EC) aims to reduce the landfilling of biodegradable waste. In Finland, 31% of biodegradable municipal waste ended up into landfills in 2012. The pressure of reducing disposing into landfills is greatly increased by the forthcoming landfill ban on biodegradable waste in Finland. There is a need to discuss the need for increasing the utilization of biodegradable waste in regional renewable energy production to utilize the waste in a way that allows the best possibilities to reduce GHG emissions. The objectives of the thesis are: (1) to find important factors affecting renewable energy recovery possibilities from biodegradable waste, (2) to determine the main factors affecting the GHG balance of biogas production system and how to improve it and (3) to find ways to define energy performance of biogas production systems and what affects it. According to the thesis, the most important factors affecting the regional renewable energy possibilities from biodegradable waste are: the amount of available feedstock, properties of feedstock, selected utilization technologies, demand of energy and material products and the economic situation of utilizing the feedstocks. The biogas production by anaerobic digestion was seen as the main technology for utilizing biodegradable waste in agriculturally dense areas. The main reason for this is that manure was seen as the main feedstock, and it can be best utilized with anaerobic digestion, which can produce renewable energy while maintaining the spreading of nutrients on arable land. Biogas plants should be located close to the heat demand that would be enough to receive the produced heat also in the summer months and located close to the agricultural area where the digestate could be utilized. Another option for biogas use is to upgrade it to biomethane, which would require a location close to the natural gas grid. The most attractive masses for biogas production are municipal and industrial biodegradable waste because of gate fees the plant receives from them can provide over 80% of the income. On the other hand, directing gate fee masses for small-scale biogas plants could make dispersed biogas production more economical. In addition, the combustion of dry agricultural waste such as straw would provide a greater energy amount than utilizing them by anaerobic digestion. The complete energy performance assessment of biogas production system requires the use of more than one system boundary. These can then be used in calculating output–input ratios of biogas production, biogas plant, biogas utilization and biogas production system, which can be used to analyze different parts of the biogas production chain. At the moment, it is difficult to compare different biogas plants since there is a wide variation of definitions for energy performance of biogas production. A more consistent way of analyzing energy performance would allow comparing biogas plants with each other and other recovery systems and finding possible locations for further improvement. Both from the GHG emission balance and energy performance point of view, the energy consumption at the biogas plant was the most significant factor. Renewable energy use to fulfil the parasitic energy demand at the plant would be the most efficient way to reduce the GHG emissions at the plant. The GHG emission reductions could be increased by upgrading biogas to biomethane and displacing natural gas or petrol use in cars when compared to biogas CHP production. The emission reductions from displacing mineral fertilizers with digestate were seen less significant, and the greater N2O emissions from spreading digestate might surpass the emission reductions from displacing mineral fertilizers.

Mixing in rotating drums

Tämän diplomityön tarkoitus on parantaa meesauunin toiminnallista tehokkuutta tehostamalla lämmönsiirtoa. Lämmönsiirron parantamiseksi kehitetään erilaisia nostinratkaisuja. Kokeita suoritetaan käyttäen eri sekoitinratkaisuja ja erilaisia prosessiparametreja. Työn kirjallisuusosassa esitetään meesauuni sekä rumpumaisten uunien toiminta. Työssä selvitetään myös sekoituksen analysointiin käytettäviä tapoja ja laskukaavoja. Kirjallisuusosassa keskitytään myös rummussa tapahtuviin fysikaalisiin ilmiöihin sekä erilaisten fluidien reologiaan. Työn kokeellisessa osassa käytettiin LUT Kemiantekniikalla suunniteltua pilot -kokoluokan rumpu-uunia, jolla kokeitaan suoritettiin, käyttäen erilaisia sekoitinratkaisuja ja sekoitusprosessiparametreja. Kokeissa käytettiin myös eri viskositeetin omaavia materiaaleja. Valitut materiaalit olivat vesi, CMC (karboksimetyyliselluloosa) ja kiinteä meesa. Kokeiden tuloksena löydettiin nostinratkaisuja, joilla sekoittumista ja lämmönsiirtoa pystytään parantamaan sekä pidentämään viipymäaikaa.

Working fluid selection and design of small-scale waste heat recovery systems based on organic Rankine cycles

Demand for the use of energy systems, entailing high efficiency as well as availability to harness renewable energy sources, is a key issue in order to tackling the threat of global warming and saving natural resources. Organic Rankine cycle (ORC) technology has been identified as one of the most promising technologies in recovering low-grade heat sources and in harnessing renewable energy sources that cannot be efficiently utilized by means of more conventional power systems. The ORC is based on the working principle of Rankine process, but an organic working fluid is adopted in the cycle instead of steam. This thesis presents numerical and experimental results of the study on the design of small-scale ORCs. Two main applications were selected for the thesis: waste heat re- covery from small-scale diesel engines concentrating on the utilization of the exhaust gas heat and waste heat recovery in large industrial-scale engine power plants considering the utilization of both the high and low temperature heat sources. The main objective of this work was to identify suitable working fluid candidates and to study the process and turbine design methods that can be applied when power plants based on the use of non-conventional working fluids are considered. The computational work included the use of thermodynamic analysis methods and turbine design methods that were based on the use of highly accurate fluid properties. In addition, the design and loss mechanisms in supersonic ORC turbines were studied by means of computational fluid dynamics. The results indicated that the design of ORC is highly influenced by the selection of the working fluid and cycle operational conditions. The results for the turbine designs in- dicated that the working fluid selection should not be based only on the thermodynamic analysis, but requires also considerations on the turbine design. The turbines tend to be fast rotating, entailing small blade heights at the turbine rotor inlet and highly supersonic flow in the turbine flow passages, especially when power systems with low power outputs are designed. The results indicated that the ORC is a potential solution in utilizing waste heat streams both at high and low temperatures and both in micro and larger scale appli- cations.

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

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

Kasvihuonekaasujen päästökaupan ohjausvaikutus sellutehtaan prosessien ja laitteiden kehityksessä

EU:n päästökaupan ensimmäinen jakso alkoi 1.1.2005. Päästökauppa on aiheuttanut sen piiriin kuuluvalle teollisuudelle monia haasteita ja riskejä kasvavien kustannusten muodossa. Massa- ja paperiteollisuus on päästökaupanpiiriin kuuluva teollisuudenala, johon päästökaupan kustannukset vaikuttavat haitallisesti globaalin hinnoittelun vuoksi. Massa- ja paperiteollisuudelle päästökaupasta voi koitua kustannuksia päästöoikeuksien ostamisesta, sähkön, polttoaineiden ja kemikaalien hinnan noususta sekä raaka-ainehuollon vaikeutumisesta. Toisaalta tehtaat voivat hyötyä päästökaupasta alittaessaan päästöoikeutensa tai myydessään sähköä ulkopuoliseen verkkoon. Massa- ja paperiteollisuudessa sähköä kuluu enimmäkseen pumppauksiin eli massan siirtoon ja mekaanisen massan valmistukseen. Suurimpia sähköenergian kuluttajia sellun valmistuksessa ovat soodakattila, puunkäsittely, valkaisu ja lajittelu. Lämpöä tarvitaan haihdutus-, kuivaus- ja keittoprosesseissa. Kemikaaleista klooridioksidin valmistuksessa käytettävä natriumkloraatti on kustannusten kannalta merkittävin kemikaali. Tässä työssä tutkittiin päästökaupan aiheuttamien kustannusten vähentämismahdollisuuksia kohdetehtaalla. Suurin potentiaali liittyy meesauunissa poltettavan maakaasun korvaamiseen mäntyöljyllä tai biomassan kaasutuskaasulla. Kemikaalikulutuksen osalta happidelignifiointi on merkittävin mahdollisuuskustannusten alentamiseksi. Lisäksi päästökaupan kustannuksia voidaan alentaa muun muassa oikealla mitoituksella ja sekundäärilämpöjen optimaalisella käytöllä.

Lämmöntalteenotto kosteasta poistoilmasta prosessiteollisuudessa

Tämän diplomityön päämääränä on tehdä prosessiteollisuuden tarpeisiin Excel-taulukkolaskentaohjelmassa toimiva putkilämmönsiirtimen mitoitusohjelma. Prosessiteollisuudessa lämmönvaihtimien toimintaympäristöt ja olosuhteet vaihtelevat merkittävästi ja niinpä jokaisen vaihtimen suunnittelu ja mitoitus on toteutettava tapauskohtaisesti. Työssä käsitellään rekuperatiivisen ristivirtaputkilämmönvaihtimen yleinen lämpötekninen mitoitus sisältäen putken sisäpinnalle tapahtuvan mahdollisen lauhtumisen. Mitoitettava vaihdinkoostuu pystysuorista putkista, joissa lämmin ja kostea ilma virtaa putkien sisäpuolella ja kylmä kuiva ilma vaippapuolella vaakasuoraan. Vaihdinmateriaalina käytetään ruostumatonta AISI 304 -tai haponkestävää AISI 316 terästä. Kuuman ilman tila vaihtelee tarkasteltavan kohteen mukaan. Paperiteollisuuden kuivausyksiköiltä poistuva ilma on usein lämmintä ja kosteaa, ja infrakuivaimilta poistuva ilma on kuumaa. Mitoitettavalle lämmönvaihtimelle tulevan kuuman ilman lämpötila tapauksesta riippuen voi vaihdella 30°C, maksimissaan +300°C:een saakka, vesisisällön ollessa välillä 0,010...0,200 kg/kg ki tai jopa tämän ylikin. Vaihtimen mitoitus perustuu energiataseyhtälöiden käyttöön. Laskennassa määritetään vaihtimen pintalämpötila sekä mahdollinen kostean ilman lauhtuminen putken sisäpinnalle. Lisäksi teoria käsittää molempien virtausten tilanmuutosten laskennan. Työssä on esitetty esimerkkilaskelma, jossa on laskettu ilma- kostea ilma lämmönsiirrinkonstruktio. Esimerkissä on tarkasteltu vaihtimen hyötysuhdetta, virtausten lämpö- ja kosteuskäyttäytymistä ulkoilman lämpötilan funktiona. Ohjelmasta saadaan tulostettua mitoitettavanvaihtimen dimensiot; putkien lukumäärät syvyys- ja pituussuunnassa sekä kokonaisputkilukumäärä, putkien väliset etäisyydet toisiinsa nähden sekä syvyys, että pituussuunnassa, putkipituus ja putken sisä- ja ulkohalkaisijat. Nämä tiedot suunnittelija itse syöttää ohjelmalle alkuarvoina. Laskettuna tietona ohjelma antaa molempien virtausten poistolämpötilat, kuuman ilman poistuvan absoluuttisen kosteuden, kondenssivesimäärän, vaihtimen tehon ja painehäviöt vaippa- ja putkipuolelle. Lisäksi ohjelma laskee kuuman ilman ominaisentalpiat vaihtimen sisään- ja ulostulossa. Tämä mahdollistaa ilman tilapisteiden piirtämisen Mollier-piirrokseen.

Sähkön omatuotannon kannattavuustarkastelu

Tässä diplomityössä tutkitaan sähkön omatuotannon kannattavuutta M-realin Simpeleen tehtaiden voimalaitoksella. Erityisesti työssä arvioidaan lauhdesähköntuotannon kannattavuutta polttoaine- ja päästökustannuksista muodostuvien marginaalikustannusten osalta. Koska voimalaitoksen rakennusaste on varsin alhainen,sähköntuotannon kannattavuutta on tarkasteltu arvioimalla sähkön ja lämmön yhteistuotannon kustannuksia ja jakamalla syntyneet kustannukset suhdemenetelmän avulla eri tuotteille. Diplomityössä etsitään kustannustehokkain seospolttosuhde annettujen reunaehtojen puitteissa muodostamalla polttoaineista aiheutuvista kustannuksista laskentamalli, jota optimoidaan Microsoft Excelin Solver-toiminnolla. Lauhdesähköntuotannon marginaalikustannuksia verrataan Nord Poolin SPOT-tuntihintaan. Lauhdesähköntuotanto voimalaitoksella on kannattavaa, mikäli SPOT-tuntihinnan vuorokautinen keskiarvo ylittää lauhdesähköntuotannon marginaalikustannukset.