Modeling of non-isothermalvapor membrane separation with thermodynamic models and generalized mass transfer equations

A rigorous unit operation model is developed for vapor membrane separation. The new model is able to describe temperature, pressure, and concentration dependent permeation as wellreal fluid effects in vapor and gas separation with hydrocarbon selective rubbery polymeric membranes. The permeation through the membrane is described by a separate treatment of sorption and diffusion within the membrane. The chemical engineering thermodynamics is used to describe the equilibrium sorption of vapors and gases in rubbery membranes with equation of state models for polymeric systems. Also a new modification of the UNIFAC model is proposed for this purpose. Various thermodynamic models are extensively compared in order to verify the models' ability to predict and correlate experimental vapor-liquid equilibrium data. The penetrant transport through the selective layer of the membrane is described with the generalized Maxwell-Stefan equations, which are able to account for thebulk flux contribution as well as the diffusive coupling effect. A method is described to compute and correlate binary penetrant¿membrane diffusion coefficients from the experimental permeability coefficients at different temperatures and pressures. A fluid flow model for spiral-wound modules is derived from the conservation equation of mass, momentum, and energy. The conservation equations are presented in a discretized form by using the control volume approach. A combination of the permeation model and the fluid flow model yields the desired rigorous model for vapor membrane separation. The model is implemented into an inhouse process simulator and so vapor membrane separation may be evaluated as an integralpart of a process flowsheet.

Rikkiyhdisteiden laskenta simuloinneissa

Diplomityön tarkoituksena on tutkia eri laskentamenetelmien soveltuvuutta kevyiden rikkiyhdisteiden laskentaan ja kuinka mitatusta kaasu-neste tasapainotiedoista sovitetut binääriset vuorovaikutusparametrit parantavat kaasu-neste tasapainojen laskentaa simuloinneissa. Kirjallisuusosassa paneudutaan kevyisiin rikkiyhdisteisiin ja niiden aineominaisuuksiin. Lisäksi käsitellään öljynjalostuksessa nykyisin käytettäviä ja uusia kehitteillä olevia rikinpoistomenetelmiä.Kokeellisessa osassa tarkastellaan eri laskentamenetelmien soveltuvuutta rikkiyhdisteiden ja kevyiden hiilivetyjen kaasu-neste tasapainon laskentaan. Mitatusta rikkiyhdisteiden ja hiilivetyjen kaasu-neste tasapainoista sovitetaan binäärisiä vuorovaikutusparametrejä tarkentamaan käytettäviä laskentamenetelmiä. Osassa verrataan binääristen seosten mittaustuloksia eri laskentamenetelmillä saatuihin simulointituloksiin. Tarkasteluiden perusteella tehdään johtopäätöksiä laskentamenetelmien soveltuvuudesta kevyiden hiilivetyjen ja rikkiyhdisteiden laskentaan. Tarkastellaan kahden prosessin (rikkivetystripperi ja butaaninpoistokolonni) rikkiyhdisteiden laskentaa. Prosesseille tehdään taseajot, joista saatuja analyysituloksia verrataan simulointien antamiin tuloksiin. Työssä tarkastellaan myös veden liukoisuuden laskentaa ja mahdollisten laskentamenetelmien käytön vaikutusta rikkiyhdisteiden laskentaan.

Micronization of Pharmaceuticals and Food Ingredients using Supercritical Fluid Techniques

Micronization techniques based on supercritical fluids (SCFs) are promising for the production of particles with controlled size and distribution. The interest of the pharmaceutical field in the development of SCF techniques is increasing due to the need for clean processes, reduced consumption of energy, and to their several possible applications. The food field is still far from the application of SCF micronization techniques, but there is increasing interest mainly for the processing of products with high added value. The aim of this study is to use SCF micronization techniques for the production of particles of pharmaceuticals and food ingredients with controlled particle size and morphology, and to look at their production on semi-industrial scale. The results obtained are also used to understand the processes from the perspective of broader application within the pharmaceutical and food industries. Certain pharmaceuticals, a biopolymer and a food ingredient have been tested using supercritical antisolvent micronization (SAS) or supercritical assisted atomization (SAA) techniques. The reproducibility of the SAS technique has been studied using physically different apparatuses and on both laboratory and semi-industrial scale. Moreover, a comparison between semi-continuous and batch mode has been performed. The behaviour of the system during the SAS process has been observed using a windowed precipitation vessel. The micronized powders have been characterized by particle size and distribution, morphology and crystallinity. Several analyses have been performed to verify if the SCF process modified the structure of the compound or caused degradation or contamination of the product. The different powder morphologies obtained have been linked to the position of the process operating point with respect to the vapour-liquid equilibrium (VLE) of the systems studied, that is, mainly to the position of the mixture critical point (MCP) of the mixture. Spherical micro, submicro- and nanoparticles, expanded microparticles (balloons) and crystals were obtained by SAS. The obtained particles were amorphous or with different degrees of crystallinity and, in some cases, had different pseudo-polymorphic or polymorphic forms. A compound that could not be processed using SAS was micronized by SAA, and amorphous particles were obtained, stable in vials at room temperature. The SCF micronization techniques studied proved to be effective and versatile for the production of particles for several uses. Furthermore, the findings of this study and the acquired knowledge of the proposed processes can allow a more conscious application of SCF techniques to obtain products with the desired characteristics and enable the use of their principles for broader applications.

Synteesikaasun liukoisuus ja aineensiirto fermentointiliuokseen

Työssä tutkittiin synteesikaasun komponenttien: hiilimonoksidin, vedyn ja hiilidioksidin liukoisuutta ja aineensiirtonopeutta fermentointiliuokseen. Kirjallisuusosassa käsitellään kaasujen liukoisuuksiin ja kaasu-nesteaineensiirtoon vaikuttavia tekijöitä ja esitellään korrelaatioita, jotka on kehitetty volumetrisen aineensiirtokertoimen ennustamiseen sekoitetussa fermentorissa. Kirjallisuus-osassa esitetään myös synteesikaasun komponenttien liukoisuudet veteen, etanoliin ja etikkahappoon 37 ºC lämpötilassa ja esitellään Flowbat-simulointiohjelman MHV2-mallin käyttöä kaasu-neste tasapainojen mallin-nuksessa. Työn kokeellisessa osassa tutkittiin synteesikaasun komponenttien liukoisuuksia veteen ja kasvatusalustaan sekä kehitettiin mittausmenetelmä kaasu-neste tasapainojen mittaukseen. Tasapainomittauksissa tutkittiin etanolin ja etikkahapon konsentraatioiden vaikutusta synteesikaasun liukoisuuteen. Lisäksi mallinnettiin kaasu-neste tasapainoja monikomponenttisysteemeissä MHV2-mallin avulla. Kokeellisen osan aineensiirtomittauksissa tutkittiin sekoitusnopeuden ja kaasun volumetrisen syöttönopeuden vaikutusta hiilimonoksidin ja vedyn volumetriseen aineensiirtokertoimeen kLa kahden litran tilavuuksisessa laboratoriofermentorissa. Mittaustulosten perusteella kasvatusalustan komponentit vaikuttavat merkittävästi hiilidioksidin liukoisuuteen. Lisäys etanolin ja etikkahapon konsentraatioissa parantaa hiilimonoksidin liukoisuutta kasvatusalustaan. Kaasun volumetrinen syöttönopeus ja sekoitusnopeus vaikuttavat voimakkaasti volumetrisen aineensiirtokertoimen arvoon. Tutkitussa systeemissä korkein teoreettinen solutiheys, joka voitiin saavuttaa suurimmalla hiilimonoksidin aineensiirto-nopeudella, oli 3 g/L. Tämä on kaksinkertainen verrattuna aiemmissa VTT:n kokeissa saavutettuihin solutiheyksiin.

Nesteseosten liukoisuuksien mittaaminen

Diplomityön tarkoituksena on tarkastella isobuteenin dimeroitumisprosessin neste-nestetasapainoja. Tavoitteena on määrittää prosessissa esiintyvien komponenttien väliset neste-nestetasapainot. Työn kirjallisuusosassa on tarkasteltu neste-nestetasapainojen teoriaa. Erityisesti on tarkasteltu mittausmenetelmiä sekä kirjallisuudesta löytyneitä laitteistoja binääristen ja ternääristen systeemien neste-nestetasapainojen määritystä varten. Menetelmät ja laitteistot on esitetty erikseen matalassa ja korkeassa paineessa suoritettaville mittauksille. Lisäksi on tarkasteltu näytteenottoa sekä näytteiden analysointimenetelmiä. Kirjallisuusosassa on myös sivuttu kaasu-neste-nestetasapainojen määritystä, mutta työn varsinainen kohde on neste-nestetasapainojen määritys. Työn kokeellisessa osassa määritettiin iso-oktaaniprosessissa esiintyvien komponenttien välisiä binäärisiä ja ternäärisiä neste-nestetasapainoja. Mitattavien komponenttiparien määrää karsittiin ja jäljellejääneiden parien välillä suoritettavat mittaukset jaoteltiin matalassa ja korkeassa paineessa suoritettaviin määrityksiin. Ternääriset mittaukset tulivat kyseeseen sellaisten komponenttiparien kohdalla, joissa toisiinsa täysin liukenevien nesteiden systeemiin kolmatta komponenttia lisättäessä saatiin aikaiseksi kaksi nestefaasia. Tällaisesta mittaustiedosta voidaan määrittää neste-nestetasapainomallien parametrejä. Mittausten lisäksi kokeellisessa osassa tarkasteltiin näytteenottoa sekä näytteiden analysointia.

Influence of multi-phase phenomena on semibatch crystallization processes of aqueous solutions

Crystal properties, product quality and particle size are determined by the operating conditions in the crystallization process. Thus, in order to obtain desired end-products, the crystallization process should be effectively controlled based on reliable kinetic information, which can be provided by powerful analytical tools such as Raman spectrometry and thermal analysis. The present research work studied various crystallization processes such as reactive crystallization, precipitation with anti-solvent and evaporation crystallization. The goal of the work was to understand more comprehensively the fundamentals, phenomena and utilizations of crystallization, and establish proper methods to control particle size distribution, especially for three phase gas-liquid-solid crystallization systems. As a part of the solid-liquid equilibrium studies in this work, prediction of KCl solubility in a MgCl2-KCl-H2O system was studied theoretically. Additionally, a solubility prediction model by Pitzer thermodynamic model was investigated based on solubility measurements of potassium dihydrogen phosphate with the presence of non-electronic organic substances in aqueous solutions. The prediction model helps to extend literature data and offers an easy and economical way to choose solvent for anti-solvent precipitation. Using experimental and modern analytical methods, precipitation kinetics and mass transfer in reactive crystallization of magnesium carbonate hydrates with magnesium hydroxide slurry and CO2 gas were systematically investigated. The obtained results gave deeper insight into gas-liquid-solid interactions and the mechanisms of this heterogeneous crystallization process. The research approach developed can provide theoretical guidance and act as a useful reference to promote development of gas-liquid reactive crystallization. Gas-liquid mass transfer of absorption in the presence of solid particles in a stirred tank was investigated in order to gain understanding of how different-sized particles interact with gas bubbles. Based on obtained volumetric mass transfer coefficient values, it was found that the influence of the presence of small particles on gas-liquid mass transfer cannot be ignored since there are interactions between bubbles and particles. Raman spectrometry was successfully applied for liquid and solids analysis in semi-batch anti-solvent precipitation and evaporation crystallization. Real-time information such as supersaturation, formation of precipitates and identification of crystal polymorphs could be obtained by Raman spectrometry. The solubility prediction models, monitoring methods for precipitation and empirical model for absorption developed in this study together with the methodologies used gives valuable information for aspects of industrial crystallization. Furthermore, Raman analysis was seen to be a potential controlling method for various crystallization processes.

Kinetic modeling of mechanisms of industrially important organic reactions in gas and liquid phase

This dissertation is based on 5 articles which deal with reaction mechanisms of the following selected industrially important organic reactions: 1. dehydrocyclization of n-butylbenzene to produce naphthalene 2. dehydrocyclization of 1-(p-tolyl)-2-methylbutane (MB) to produce 2,6-dimethylnaphthalene 3. esterification of neopentyl glycol (NPG) with different carboxylic acids to produce monoesters 4. skeletal isomerization of 1-pentene to produce 2-methyl-1-butene and 2-methyl-2-butene The results of initial- and integral-rate experiments of n-butylbenzene dehydrocyclization over selfmade chromia/alumina catalyst were applied when investigating reaction 2. Reaction 2 was performed using commercial chromia/alumina of different acidity, platina on silica and vanadium/calcium/alumina as catalysts. On all catalysts used for the dehydrocyclization, major reactions were fragmentation of MB and 1-(p-tolyl)-2-methylbutenes (MBes), dehydrogenation of MB, double bond transfer, hydrogenation and 1,6-cyclization of MBes. Minor reactions were 1,5-cyclization of MBes and methyl group fragmentation of 1,6- cyclization products. Esterification reactions of NPG were performed using three different carboxylic acids: propionic, isobutyric and 2-ethylhexanoic acid. Commercial heterogeneous gellular (Dowex 50WX2), macroreticular (Amberlyst 15) type resins and homogeneous para-toluene sulfonic acid were used as catalysts. At first NPG reacted with carboxylic acids to form corresponding monoester and water. Then monoester esterified with carboxylic acid to form corresponding diester. In disproportionation reaction two monoester molecules formed NPG and corresponding diester. All these three reactions can attain equilibrium. Concerning esterification, water was removed from the reactor in order to prevent backward reaction. Skeletal isomerization experiments of 1-pentene were performed over HZSM-22 catalyst. Isomerization reactions of three different kind were detected: double bond, cis-trans and skeletal isomerization. Minor side reaction were dimerization and fragmentation. Monomolecular and bimolecular reaction mechanisms for skeletal isomerization explained experimental results almost equally well. Pseudohomogeneous kinetic parameters of reactions 1 and 2 were estimated by usual least squares fitting. Concerning reactions 3 and 4 kinetic parameters were estimated by the leastsquares method, but also the possible cross-correlation and identifiability of parameters were determined using Markov chain Monte Carlo (MCMC) method. Finally using MCMC method, the estimation of model parameters and predictions were performed according to the Bayesian paradigm. According to the fitting results suggested reaction mechanisms explained experimental results rather well. When the possible cross-correlation and identifiability of parameters (Reactions 3 and 4) were determined using MCMC method, the parameters identified well, and no pathological cross-correlation could be seen between any parameter pair.

Liquid transportation and distribution during (re)wetting : impact on coated papers in printing

In many industrial applications, such as the printing and coatings industry, wetting of porous materials by liquids includes not only imbibition and permeation into the bulk but also surface spreading and evaporation. By understanding these phenomena, valuable information can be obtained for improved process control, runnability and printability, in which liquid penetration and subsequent drying play important quality and economic roles. Knowledge of the position of the wetting front and the distribution/degree of pore filling within the structure is crucial in describing the transport phenomena involved. Although exemplifying paper as a porous medium in this work, the generalisation to dynamic liquid transfer onto a surface, including permeation and imbibition into porous media, is of importance to many industrial and naturally occurring environmental processes. This thesis explains the phenomena in the field of heatset web offset printing but the content and the analyses are applicable in many other printing methods and also other technologies where water/moisture monitoring is crucial in order to have a stable process and achieve high quality end products. The use of near-infrared technology to study the water and moisture response of porous pigmented structures is presented. The use of sensitive surface chemical and structural analysis, as well as the internal structure investigation of a porous structure, to inspect liquid wetting and distribution, complements the information obtained by spectroscopic techniques. Strong emphasis has been put on the scale of measurement, to filter irrelevant information and to understand the relationship between interactions involved. The near-infrared spectroscopic technique, presented here, samples directly the changes in signal absorbance and its variation in the process at multiple locations in a print production line. The in-line non-contact measurements are facilitated by using several diffuse reflectance probes, giving the absolute water/moisture content from a defined position in the dynamic process in real-time. The nearinfrared measurement data illustrate the changes in moisture content as the paper is passing through the printing nips and dryer, respectively, and the analysis of the mechanisms involved highlight the roles of the contacting surfaces and the relative liquid carrier properties of both non-image and printed image areas. The thesis includes laboratory studies on wetting of porous media in the form of coated paper and compressed pigment tablets by mono-, dual-, and multi-component liquids, and paper water/moisture content analysis in both offline and online conditions, thus also enabling direct sampling of temporal water/moisture profiles from multiple locations. One main focus in this thesis was to establish a measurement system which is able to monitor rapid changes in moisture content of paper. The study suggests that near-infrared diffuse reflectance spectroscopy can be used as a moisture sensitive system and to provide accurate online qualitative indicators, but, also, when accurately calibrated, can provide quantification of water/moisture levels, its distribution and dynamic liquid transfer. Due to the high sensitivity, samples can be measured with excellent reproducibility and good signal to noise ratio. Another focus of this thesis was on the evolution of the moisture content, i.e. changes in moisture content referred to (re)wetting, and liquid distribution during printing of coated paper. The study confirmed different wetting phases together with the factors affecting each phase both for a single droplet and a liquid film applied on a porous substrate. For a single droplet, initial capillary driven imbibition is followed by equilibrium pore filling and liquid retreat by evaporation. In the case of a liquid film applied on paper, the controlling factors defining the transportation were concluded to be the applied liquid volume in relation to surface roughness, capillarity and permeability of the coating giving the liquid uptake capacity. The printing trials confirmed moisture gradients in the printed sheet depending on process parameters such as speed, fountain solution dosage and drying conditions as well as the printed layout itself. Uneven moisture distribution in the printed sheet was identified to be one of the sources for waving appearance and the magnitude of waving was influenced by the drying conditions.

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

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

Equilibrium moisture content of flax/linseed and fibre hemp straw fractions

Selostus: Pellavan ja kuituhampun korren jakeiden tasapainokosteus

Equilibrium curves and growth models to deal with forests in transition to uneven-aged structure : application in two sample stands

Abstract

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.

Water vapor tightness of packages: methods and optimizing

Työn tarkoituksena oli testata jo tutkimuskeskuksella käytössä ollutta ja tutkimuskeskukselle tässä työssä kehitettyä pakkauksen vesihöyrytiiveyteen liittyvää mittausmenetelmää. Saatuja tuloksia verrattiin keskenään sekä materiaalista mitattuihin arvoihin. Elintarvikepakkauksia tutkittiin myös kosteussensoreiden, säilyvyyskokeen sekä kuljetussimuloinnin avulla. Optimoinnilla tutkittiin pakkauksen muodon vaikutusta vesihöyrytiiveyteen. Pakkauksen vesihöyrynläpäisyn mittaamiseen kehitetty menetelmä toimi hyvin ja sen toistettavuus oli hyvä. Verrattaessa sitä jo olemassa olleeseen menetelmään tulokseksi saatiin, että uusi menetelmä oli nopeampi ja vaati vähemmän työaikaa, mutta molemmat menetelmät antoivat hyviä arvoja rinnakkaisille näytteille. Kosteussensoreilla voitiin tutkia tyhjän pakkauksen sisällä olevan kosteuden muutoksia säilytyksen aikana. Säilyvyystesti tehtiin muroilla ja parhaan vesihöyrysuojan antoivat pakkaukset joissa oli alumiinilaminaatti- tai metalloitu OPP kerros. Kuljetustestauksen ensimmäisessä testissä pakkauksiin pakattiin muroja ja toisessa testissä nuudeleita. Kuljetussimuloinnilla ei ollutvaikutusta pakkausten sisäpintojen eheyteen eikä siten pakkausten vesihöyrytiiveyteen. Optimoinnilla vertailtiin eri muotoisten pakkausten tilavuus/pinta-ala suhdetta ja vesihöyrytiiveyden riippuvuutta pinta-alasta. Optimaalisimmaksi pakkaukseksi saatiin pallo, jonka pinta-ala oli pienin ja materiaalin sallima vesihöyrynläpäisy suurin ja vesihöyrybarrierin määrä pienin.

Solid and liquid biofuels markets in Finland - a study on international biofuels trade

This study considered the current situation of solid and liquid biofuels markets and international biofuels trade in Finland and identified the challenges ofthe emerging international biofuels markets for Finland. The fact that industryconsumes more than half of the total primary energy, widely applied combined heat and power production (CHP) and a high share of biofuels in the total energy consumption are specific to the Finnish energy system. One third of the electricity is generated in CHP plants. As much as 27% of the total energy consumption ismet by using wood and peat, which makes Finland the leading country in the use of biofuels. Finland has made a commitment to maintain greenhouse gas emissions at the 1990 level at the highest during the period 2008-2012. The Finnish energypolicy aims to achieve the target, and a variety of measures are taken to promote the use of renewable energy sources and especially wood fuels. In this study, the wooden raw material streams of the forest industry were included the international biofuels trade in addition to biomass streams that are traded for energy production. In 2004, as much as 45% of the raw wood importedinto Finland ended up in energy production. The total international trading of biofuels was evaluated at 72 PJ, of which the majority, 58 PJ, was raw wood. About 22% of wood based energy in Finland originated from imported raw wood. Tall oil and wood pellets composed the largest export streams of biofuels. The annual turnover of international biofuels trade was estimated at about ¤ 90 million fordirect trade and at about ¤ 190 million for indirect trade. The forest industryas the biggest user of wood, and the producer and user of wood fuels has a central position in biomass and biofuels markets in Finland. Lately, the international aspects of Finnish biofuels markets have been emphasised as the import of rawwood and the export of wood pellets have increased. Expanding the use of biofuels in the road transportation sector would increase the international streams ofbiofuels in Finland. In coming years, the international trading of biomass for energy purposes can be expected to continue growing.