Red clover grown in a mixture with grasses : yield, persistence and dynamics of quality characteristics

Selostus: Puna-apilan pysyvyys apila-heinänurmessa sekä seosnurmen satoisuus ja laadun muutokset erilaisissa kasvuoloissa

Implementation of population balances in simulation of gas-liquid reactors using a CFD-PBM coupled approach

In bubbly flow simulations, bubble size distribution is an important factor in determination of hydrodynamics. Beside hydrodynamics, it is crucial in the prediction of interfacial area available for mass transfer and in the prediction of reaction rate in gas-liquid reactors such as bubble columns. Solution of population balance equations is a method which can help to model the size distribution by considering continuous bubble coalescence and breakage. Therefore, in Computational Fluid Dynamic simulations it is necessary to couple CFD and Population Balance Model (CFD-PBM) to get reliable distribution. In the current work a CFD-PBM coupled model is implemented as FORTRAN subroutines in ANSYS CFX 10 and it has been tested for bubbly flow. This model uses the idea of Multi Phase Multi Size Group approach which was previously presented by Sha et al. (2006) [18]. The current CFD-PBM coupled method considers inhomogeneous flow field for different bubble size groups in the Eulerian multi-dispersed phase systems. Considering different velocity field for bubbles can give the advantageof more accurate solution of hydrodynamics. It is also an improved method for prediction of bubble size distribution in multiphase flow compared to available commercial packages.

Effect of Shielding Gases in Laser Welding of Aluminium Alloys

High reflectivity to laser light, alloying element evaporation during high power laser welding makes aluminium alloys highly susceptibility to weld defects such as porosity, cracking and undercutting. The dynamic behaviour of the keyhole, due to fluctuating plasma above the keyhole and the vaporization ofthe alloying elements with in the keyhole, is the key problem to be solved for the improvement of the weld quality and stabilization of the keyhole dynamics isperhaps the single most important development that can broaden the application of laser welding of aluminium alloys. In laser welding, the shielding gas is commonly used to stabilize the welding process, to improve the welded joint features and to protect the welded seam from oxidation. The chemicalcomposition of the shielding gas is a key factor in achieving the final qualityof the welded joints. Wide range of shielding gases varying from the pure gasesto complex mixtures based on helium, argon, nitrogen and carbon dioxide are commercially available. These gas mixtures should be considered in terms of their suitability during laser welding of aluminium alloys to produce quality welds. The main objective of the present work is to study the effect of the shielding gascomposition during laser welding of aluminium alloys. Aluminium alloy A15754 was welded using 3kW Nd-YAG laser (continuous wave mode). The alloy samples were butt welded with different shielding gases (pure and mixture of gases) so that high quality welds with high joint efficiencies could be produced. It was observed that the chemical composition of the gases influenced the final weld quality and properties. In general, the mixture gases, in correct proportions, enabled better utilisation of the properties of the mixing gases, stabilized the welding process and produced better weld quality compared to the pure shielding gases.

Plasmapolttoprosessin kehittäminen ja optimointi

Ydinvoimaloidenprimaarivesikierron puhdistukseen käytetään ioninvaihtohartsia. Käytönjälkeen ioninvaihtohartsi luokitellaan matalaja keskiaktiivisiin jätteisiin. Plasmakäsittelyllä käytetyn ioninvaihtohartsin tilavuutta voidaan pienentää sekä sen orgaaninen luonne poistaa. Plasmakäsittelyn tarkoituksena on hapettaa orgaaninen aines oksideiksi, jotka poistuvat prosessista savukaasuina. Epäorgaaninen aines, joka sisältää radioaktiivisen aineksen, on tarkoitus hapettaa oksideiksi ja sulfideiksi, jotka voidaan kerätä talteen tuhkana. Tässä diplomityössä käsitellään käytetyn ioninvaihtohartsin käsittelyyn suunnitellun plasmapolttoprosessin kehittämistä ja optimointia. Ioninvaihtohartsin plasmakäsittelyssä syntyvien reaktiotuotteiden selvittäminen suoritettiin tarkastelemalla ainetaseita sekä aihetta käsitteleviä tutkimuksia. Näiden perusteella parannettiin jäähdytystä, suunniteltiin jatkuvatoiminen syöttömenetelmä sekä laadittiin toimintaalueen reunaehdot laitteistolle. Koelaitteistossa 6,5 kW:n rfteho syötetään sovitinpiirin ja kuparisen induktiokelan kautta plasmaan. Plasmakaasuna on käytetty hapenja argonin seoskaasua. Plasmapolttoa on seurattu massaspektrometrilla, optisella emissiospektrometrilla, lämpösekä painemittareilla. Laskennan ja kokeiden pohjalta selvitettiin optimaalinen seossuhde plasmakaasulle, paineen ja tehon noston vaikutus hartsin polttonopeuteen sekä jatkuvatoimisen syöttömenetelmän edut panostoimiseen syöttöön. Rfgeneraattorin teho rajoitti jatkuvatoimisen polttonopeuden 130 g/h ja hetkellisen polttonopeuden 175 g/h. Radioaktiivisten aineiden pidätys oli 93,5 % cesiumin osalta. Tulosten perusteella 4 kg/h ioninvaihtohartsia polttavan laitteiston tehon lähteeksi tarvitaan 65 kW rfgeneraattori. Palamattoman hartsin ja tuhkan kulkeutuminen partikkelisuodattimille sekä reaktiotuotteena syntyvien rikinoksidien käsittely vaatii vielä jatkotutkimusta.

Typpeä sisältävien jätevesien käsittely 2-n-PRO menetelmällä

Työssä tutkittiin kaksivaiheisen typenpoistoprosessin (2-N-PRO) soveltuvuutta Joutsenon Kukkuroinmäen aluejätekeskuksen kompostointilaitoksen jätevesille pilot-kokein 12.1.- 5.4.2006. Kompostilaitoksella on jätevesien esikäsittelytarve korkeista ammoniumtyppipitoisuuksista johtuen. Pilot-laitteisto koostuu sekoitussäiliöstä, strippaustornista ja katalyyttipolttimesta. Käsiteltävän jäteveden pH nostetaan korkealle tasolle, jolloin ammoniumtyppi muuttuu ammoniakiksi. Vesi johdetaan strippaustorniin, jossa se sadetetaan tornin pohjalle. Ammoniakki erottuu sadetuksessa ilmaan, joka imetään katalyyttipolttimelle. Katalyyttinen poltin käsittelee ammoniakkia typpikaasuksi. Pilot-kokeet suoritettiin jatkuvatoimisesti. Laitteisto pystyy erottamaan jätevedestä ammoniumtyppeä ammoniakiksi ja käsittelemään ammoniakin pääosin typpikaasuksi. Lisäksi suoritettiin panoskoe, jonka tulokset tukevat jatkuvatoimisesta käytöstä saatuja tuloksia.

Reduction of Greenhouse Gas Emissions in North-West Russia - Finnish Business Opportunities

Tutkimus suomalaisten yritysten liiketoimintamahdollisuuksista hiilidoksidipäästöjen vähentämisen parissa Luoteis-Venäjällä.

Computational Analysis and Optimization of Real Gas Flow In Small Centrifugal Compressors

Small centrifugal compressors are more and more widely used in many industrialsystems because of their higher efficiency and better off-design performance comparing to piston and scroll compressors as while as higher work coefficient perstage than in axial compressors. Higher efficiency is always the aim of the designer of compressors. In the present work, the influence of four partsof a small centrifugal compressor that compresses heavy molecular weight real gas has been investigated in order to achieve higher efficiency. Two parts concern the impeller: tip clearance and the circumferential position of the splitter blade. The other two parts concern the diffuser: the pinch shape and vane shape. Computational fluid dynamics is applied in this study. The Reynolds averaged Navier-Stokes flow solver Finflo is used. The quasi-steady approach is utilized. Chien's k-e turbulence model is used to model the turbulence. A new practical real gas model is presented in this study. The real gas model is easily generated, accuracy controllable and fairly fast. The numerical results and measurements show good agreement. The influence of tip clearance on the performance of a small compressor is obvious. The pressure ratio and efficiency are decreased as the size of tip clearance is increased, while the total enthalpy rise keeps almost constant. The decrement of the pressure ratio and efficiency is larger at higher mass flow rates and smaller at lower mass flow rates. The flow angles at the inlet and outlet of the impeller are increased as the size of tip clearance is increased. The results of the detailed flow field show that leakingflow is the main reason for the performance drop. The secondary flow region becomes larger as the size of tip clearance is increased and the area of the main flow is compressed. The flow uniformity is then decreased. A detailed study shows that the leaking flow rate is higher near the exit of the impeller than that near the inlet of the impeller. Based on this phenomenon, a new partiallyshrouded impeller is used. The impeller is shrouded near the exit of the impeller. The results show that the flow field near the exit of the impeller is greatly changed by the partially shrouded impeller, and better performance is achievedthan with the unshrouded impeller. The loading distribution on the impeller blade and the flow fields in the impeller is changed by moving the splitter of the impeller in circumferential direction. Moving the splitter slightly to the suction side of the long blade can improve the performance of the compressor. The total enthalpy rise is reduced if only the leading edge of the splitter ismoved to the suction side of the long blade. The performance of the compressor is decreased if the blade is bended from the radius direction at the leading edge of the splitter. The total pressure rise and the enthalpy rise of thecompressor are increased if pinch is used at the diffuser inlet. Among the fivedifferent pinch shape configurations, at design and lower mass flow rates the efficiency of a straight line pinch is the highest, while at higher mass flow rate, the efficiency of a concave pinch is the highest. The sharp corner of the pinch is the main reason for the decrease of efficiency and should be avoided. The variation of the flow angles entering the diffuser in spanwise direction is decreased if pinch is applied. A three-dimensional low solidity twisted vaned diffuser is designed to match the flow angles entering the diffuser. The numerical results show that the pressure recovery in the twisted diffuser is higher than in a conventional low solidity vaned diffuser, which also leads to higher efficiency of the twisted diffuser. Investigation of the detailed flow fields shows that the separation at lower mass flow rate in the twisted diffuser is later than in the conventional low solidity vaned diffuser, which leads to a possible wider flow range of the twisted diffuser.

Numerical and ExperimentalModelling of Gas Flow and Heat Transfer in the Air Gap of An Electric Machine

This work deals with the cooling of high-speed electric machines, such as motors and generators, through an air gap. It consists of numerical and experimental modelling of gas flow and heat transfer in an annular channel. Velocity and temperature profiles are modelled in the air gap of a high-speed testmachine. Local and mean heat transfer coefficients and total friction coefficients are attained for a smooth rotor-stator combination at a large velocity range. The aim is to solve the heat transfer numerically and experimentally. The FINFLO software, developed at Helsinki University of Technology, has been used in the flow solution, and the commercial IGG and Field view programs for the grid generation and post processing. The annular channel is discretized as a sector mesh. Calculation is performed with constant mass flow rate on six rotational speeds. The effect of turbulence is calculated using three turbulence models. The friction coefficient and velocity factor are attained via total friction power. The first part of experimental section consists of finding the proper sensors and calibrating them in a straight pipe. After preliminary tests, a RdF-sensor is glued on the walls of stator and rotor surfaces. Telemetry is needed to be able to measure the heat transfer coefficients at the rotor. The mean heat transfer coefficients are measured in a test machine on four cooling air mass flow rates at a wide Couette Reynolds number range. The calculated values concerning the friction and heat transfer coefficients are compared with measured and semi-empirical data. Heat is transferred from the hotter stator and rotor surfaces to the coolerair flow in the air gap, not from the rotor to the stator via the air gap, althought the stator temperature is lower than the rotor temperature. The calculatedfriction coefficients fits well with the semi-empirical equations and precedingmeasurements. On constant mass flow rate the rotor heat transfer coefficient attains a saturation point at a higher rotational speed, while the heat transfer coefficient of the stator grows uniformly. The magnitudes of the heat transfer coefficients are almost constant with different turbulence models. The calibrationof sensors in a straight pipe is only an advisory step in the selection process. Telemetry is tested in the pipe conditions and compared to the same measurements with a plain sensor. The magnitudes of the measured data and the data from the semi-empirical equation are higher for the heat transfer coefficients than thenumerical data considered on the velocity range. Friction and heat transfer coefficients are presented in a large velocity range in the report. The goals are reached acceptably using numerical and experimental research. The next challenge is to achieve results for grooved stator-rotor combinations. The work contains also results for an air gap with a grooved stator with 36 slots. The velocity field by the numerical method does not match in every respect the estimated flow mode. The absence of secondary Taylor vortices is evident when using time averagednumerical simulation.

Theoretical Analysis and Simulations of Vertically Vibrated Granular Materials

The dynamical properties ofshaken granular materials are important in many industrial applications where the shaking is used to mix, segregate and transport them. In this work asystematic, large scale simulation study has been performed to investigate the rheology of dense granular media, in the presence of gas, in a three dimensional vertical cylinder filled with glass balls. The base wall of the cylinder is subjected to sinusoidal oscillation in the vertical direction. The viscoelastic behavior of glass balls during a collision, have been studied experimentally using a modified Newton's Cradle device. By analyzing the results of the measurements, using numerical model based on finite element method, the viscous damping coefficient was determinedfor the glass balls. To obtain detailed information about the interparticle interactions in a shaker, a simplified model for collision between particles of a granular material was proposed. In order to simulate the flow of surrounding gas, a formulation of the equations for fluid flow in a porous medium including particle forces was proposed. These equations are solved with Large Eddy Simulation (LES) technique using a subgrid-model originally proposed for compressible turbulent flows. For a pentagonal prism-shaped container under vertical vibrations, the results show that oscillon type structures were formed. Oscillons are highly localized particle-like excitations of the granular layer. This self-sustaining state was named by analogy with its closest large-scale analogy, the soliton, which was first documented by J.S. Russell in 1834. The results which has been reportedbyBordbar and Zamankhan(2005b)also show that slightly revised fluctuation-dissipation theorem might apply to shaken sand, which appears to be asystem far from equilibrium and could exhibit strong spatial and temporal variations in quantities such as density and local particle velocity. In this light, hydrodynamic type continuum equations were presented for describing the deformation and flow of dense gas-particle mixtures. The constitutive equation used for the stress tensor provides an effective viscosity with a liquid-like character at low shear rates and a gaseous-like behavior at high shear rates. The numerical solutions were obtained for the aforementioned hydrodynamic equations for predicting the flow dynamics ofdense mixture of gas and particles in vertical cylindrical containers. For a heptagonal prism shaped container under vertical vibrations, the model results were found to predict bubbling behavior analogous to those observed experimentally. This bubbling behavior may be explained by the unusual gas pressure distribution found in the bed. In addition, oscillon type structures were found to be formed using a vertically vibrated, pentagonal prism shaped container in agreement with computer simulation results. These observations suggest that the pressure distribution plays a key rolein deformation and flow of dense mixtures of gas and particles under vertical vibrations. The present models provide greater insight toward the explanation of poorly understood hydrodynamic phenomena in the field of granular flows and dense gas-particle mixtures. The models can be generalized to investigate the granular material-container wall interactions which would be an issue of high interests in the industrial applications. By following this approach ideal processing conditions and powder transport can be created in industrial systems.

Modeling dynamic behavior in tilting pad gas journal bearings

Thisresearch deals with the dynamic modeling of gas lubricated tilting pad journal bearings provided with spring supported pads, including experimental verification of the computation. On the basis of a mathematical model of a film bearing, a computer program has been developed, which can be used for the simulation of a special type of tilting pad gas journal bearing supported by a rotary spring under different loading conditions time dependently (transient running conditions due to geometry variations in time externally imposed). On the basis of literature, different transformations have been used in the model to achieve simpler calculation. The numerical simulation is used to solve a non-stationary case of a gasfilm. The simulation results were compared with literature results in a stationary case (steady running conditions) and they were found to be equal. In addition to this, comparisons were made with a number of stationary and non-stationary bearing tests, which were performed at Lappeenranta University of Technology using bearings designed with the simulation program. A study was also made using numerical simulation and literature to establish the influence of the different bearing parameters on the stability of the bearing. Comparison work was done with literature on tilting pad gas bearings. This bearing type is rarely used. One literature reference has studied the same bearing type as that used in LUT. A new design of tilting pad gas bearing is introduced. It is based on a stainless steel body and electron beam welding of the bearing parts. It has good operation characteristics and is easier to tune and faster to manufacture than traditional constructions. It is also suitable for large serial production.

Gas-liquid mass transfer in bubbly flow: Estimation of mass transfer, bubble size and reactor performance in various applications

Gas-liquid mass transfer is an important issue in the design and operation of many chemical unit operations. Despite its importance, the evaluation of gas-liquid mass transfer is not straightforward due to the complex nature of the phenomena involved. In this thesis gas-liquid mass transfer was evaluated in three different gas-liquid reactors in a traditional way by measuring the volumetric mass transfer coefficient (kLa). The studied reactors were a bubble column with a T-junction two-phase nozzle for gas dispersion, an industrial scale bubble column reactor for the oxidation of tetrahydroanthrahydroquinone and a concurrent downflow structured bed.The main drawback of this approach is that the obtained correlations give only the average volumetric mass transfer coefficient, which is dependent on average conditions. Moreover, the obtained correlations are valid only for the studied geometry and for the chemical system used in the measurements. In principle, a more fundamental approach is to estimate the interfacial area available for mass transfer from bubble size distributions obtained by solution of population balance equations. This approach has been used in this thesis by developing a population balance model for a bubble column together with phenomenological models for bubble breakage and coalescence. The parameters of the bubble breakage rate and coalescence rate models were estimated by comparing the measured and calculated bubble sizes. The coalescence models always have at least one experimental parameter. This is because the bubble coalescence depends on liquid composition in a way which is difficult to evaluate using known physical properties. The coalescence properties of some model solutions were evaluated by measuring the time that a bubble rests at the free liquid-gas interface before coalescing (the so-calledpersistence time or rest time). The measured persistence times range from 10 msup to 15 s depending on the solution. The coalescence was never found to be instantaneous. The bubble oscillates up and down at the interface at least a coupleof times before coalescence takes place. The measured persistence times were compared to coalescence times obtained by parameter fitting using measured bubble size distributions in a bubble column and a bubble column population balance model. For short persistence times, the persistence and coalescence times are in good agreement. For longer persistence times, however, the persistence times are at least an order of magnitude longer than the corresponding coalescence times from parameter fitting. This discrepancy may be attributed to the uncertainties concerning the estimation of energy dissipation rates, collision rates and mechanisms and contact times of the bubbles.

Data collection and criticality analysis for cost effective maintenance during the life cycle of a flue gas line

Työn tarkoituksena oli kerätä käyttövarmuustietoa savukaasulinjasta kahdelta suomalaiselta sellutehtaalta niiden käyttöönotosta aina tähän päivään asti. Käyttövarmuustieto koostuu luotettavuustiedoista sekä kunnossapitotiedoista. Kerätyn tiedon avulla on mahdollista kuvata tarkasti laitoksen käyttövarmuutta seuraavilla tunnusluvuilla: suunnittelemattomien häiriöiden lukumäärä ja korjausajat, laitteiden seisokkiaika, vikojen todennäköisyys ja korjaavan kunnossapidon kustannukset suhteessa savukaasulinjan korjaavan kunnossapidon kokonaiskustannuksiin. Käyttövarmuustiedon keräysmetodi on esitelty. Savukaasulinjan kriittisten laitteiden määrittelyyn käytetty metodi on yhdistelmä kyselytutkimuksesta ja muunnellusta vian vaikutus- ja kriittisyysanalyysistä. Laitteiden valitsemiskriteerit lopulliseen kriittisyysanalyysiin päätettiin käyttövarmuustietojen sekä kyselytutkimuksen perusteella. Kriittisten laitteiden määrittämisen tarkoitus on löytää savukaasulinjasta ne laitteet, joiden odottamaton vikaantuminen aiheuttaa vakavimmat seuraukset savukaasulinjan luotettavuuteen, tuotantoon, turvallisuuteen, päästöihin ja kustannuksiin. Tiedon avulla rajoitetut kunnossapidon resurssit voidaan suunnata oikein. Kriittisten laitteiden määrittämisen tuloksena todetaan, että kolme kriittisintä laitetta savukaasulinjassa ovat molemmille sellutehtaille yhteisesti: savukaasupuhaltimet, laahakuljettimet sekä ketjukuljettimet. Käyttövarmuustieto osoittaa, että laitteiden luotettavuus on tehdaskohtaista, mutta periaatteessa samat päälinjat voidaan nähdä suunnittelemattomien vikojen todennäköisyyttä esittävissä kuvissa. Kustannukset, jotka esitetään laitteen suunnittelemattomien kunnossapitokustannusten suhteena savukaasulinjan kokonaiskustannuksiin, noudattelevat hyvin pitkälle luotettavuuskäyrää, joka on laskettu laitteen seisokkiajan suhteena käyttötunteihin. Käyttövarmuustiedon keräys yhdistettynä kriittisten laitteiden määrittämiseen mahdollistavat ennakoivan kunnossapidon oikean kohdistamisen ja ajoittamisen laitteiston elinaikana siten, että luotettavuus- ja kustannustehokkuusvaatimukset saavutetaan.

Improvements in Air Intake on Frame 6 Gas Turbines

Työssä tarkastellaan kahta kaasuturbiinin imuilman sisäänottojärjestelmän kehitysmenetelmää, imuilman jäähdytystä ja sähköstaattista suodatusta. Imuilman jäähdytysmenetelmien tarkastelussa käytettiin kahta kaasuturbiinin tehonlisäystekniikoiden laskentatyökalua. Arviointi kohdistettiin Glanford Brigg Generating Station -voimalaitoksen kaasuturbiinityyppiin ja paikallisiin englantilaisiin ilmasto-olosuhteisiin. Tarkastelussa olivat kostutusjäähdytys ja overspray. Tuloksia vertailtiin keskenään ja näiden perusteella arvioitiin menetelmien vaikutuksia tehoon, hyötysuhteeseen sekä veden kulutukseen. Sähköstaattisen suodattimen prototyyppi oli rakennettu Briggin voimalaitokselle. Järjestelmää kehitetään kaupalliseksi tuotteeksi ja tätä varten kerättiin tekninen dokumentaatio kokonaisuudeksi, jota voitiin hyödyntää tuotteistusprosessissa. Imuilman jäähdyttämisellä voidaan saavuttaa merkittävä tehonlisäys ilmasto-olosuhteista riippuen. Menetelmällä voidaan myös tasata lämpötilan vuorokausi-vaihtelusta aiheutuvia tehoeroja. Sähköstaattisen suodattimen prototyyppi saavutti kehitysvaiheelle asetetut tavoitteet. Sähköstaattinen suodatus tarjoaa useita etuja perinteiseen mekaaniseen suodatukseen verrattuna.

Kestävän kehityksen näkökulmien huomioiminen kunnan energiaratkaisussa

Diplomityön tavoitteena on tuottaa informaatiota kunnalliseen päätöksentekoon, jonka avulla kestävän kehityksen näkökulmia voidaan huomioida kunnan energiaratkaisusta päätettäessä. Yhtenä työn lähtökohtana on ollut myös uusi EU-direktiivi, jonka mukaan ympäristönäkökohtia voidaan huomioida julkisten hankintojen tarjouspyyntömenettelyssä valintaperusteena. Tarkastelun kohteena oli kokoluokaltaan 0,5–3 MW:n aluelämpölaitokset sekä polttoaineiden tuotantoketjut. Työssä vertailtavat polttoaineet olivat metsähake, raskas polttoöljy, kevyt polttoöljy ja turve. Diplomityössä on perehdytty kestävän kehityksen käsitteeseen ja muodostettu sen mukaan ekologiselle, sosiaaliselle ja taloudelliselle näkökulmalle kunnallisen energiaratkaisun indikaattoreita. Empiirisessä osassa käsitellään kestävän kehityksen näkökulmien muodostumista Enon energiaosuuskunnan toimintaan perustuen. Käytettävät kestävän kehityksen näkökulmien mukaiset indikaattorit ovat polttoaineen tuotannosta ja käytöstä aiheutuvat kasvihuonekaasupäästöt, polttoaineen tuotannon työllisyysvaikutukset sekä energian hinnan muodostuminen osuuskunnan asiakkaille. Tässä diplomityössä tarkastelluilla kestävän kehityksen indikaattoreilla mitattuna, metsähakkeen käytöllä energiantuotannossa on positiivinen vaikutus niin kunnan kasvihuonekaasutaseessa, työllisyystilanteessa sekä myös enemmän kuluttajaystävällinen asema, lämmön hinnan vakauden ansiosta, kuin muilla työssä käsiteltävillä polttoaineilla. Polttoaineen tuotantoketjun osalta metsähakkeelle saatiin tuotannon ja käytön aiheuttamaksi kasvihuonekaasupäästöksi 2,9–4,2 g CO2-ekv/MJ. Tulos perustuu Enon energiaosuuskunnan polttoaineen hankinnassa käytössä oleviin keskimääräisiin etäisyyksiin metsäkuljetuksessa (250 m) ja kaukokuljetuksessa (15 km). Tuotannon ja käytön aiheuttamat kasvihuonekaasupäästöt olivat raskaalla polttoöljyllä 88,2 g CO2-ekv/MJ, kevyellä polttoöljyllä 85,0 g CO2-ekv/MJ ja turpeella 104,0–108,1 g CO2-ekv/MJ. Metsähakkeen osalta polttoaineen tuotannon osuus koko tarkastellun energiaketjun kasvihuonekaasupäästöistä oli noin 43–57 %. Enon energiaosuuskunnan tapauksessa vuoden 2005 odotetulla toiminta-asteella metsähakkeella tuotetun lämmön tuotantoketjun kasvihuonekaasupäästöt ovat noin 160 t CO2-ekv. Kevyellä polttoöljyllä tuotetun lämmön tuotantoketjun kasvihuonekaasupäästöt olisivat noin 3700 t CO2-ekv sekä turpeen (50 %) ja metsähakkeen (50 %) seoskäyttöön perustuvalla ketjulla noin 2300–2400 t CO2-ekv. Samaisella toiminta-asteella työllisyysvaikutukset ovat käytettäessä metsähaketta 2,2–8,6 htv, raakaöljyä 0,12 htv ja turvetta 1,4–1,6 htv. Metsähakkeen käyttö aluelämpölaitosten pääpolttoaineena takaa myös vakaan hintakehityksen osuuskunnan asiakkaille.