A Development and Application of a Three-Dimensional Model for a CFB Gasifier

The developing energy markets and rising energy system costs have sparked the need to find new forms of energy production and increase the self-sufficiency of energy production. One alternative is gasification, whose principles have been known for decades, but it is only recently when the technology has become a true alternative. However, in order to meet the requirements of modern energy production methods, it is necessary to study the phenomenon thoroughly. In order to understand the gasification process better and optimize it from the viewpoint of ecology and energy efficiency, it is necessary to develop effective and reliable modeling tools for gasifiers. The main aims of this work have been to understand gasification as a process and furthermore to develop an existing three-dimensional circulating fluidized bed modeling tool for modeling of gasification. The model is applied to two gasification processes of 12 and 50 MWth. The results of modeling and measurements have been compared and subsequently reviewed. The work was done in co-operation with Lappeenranta University of Technology and Foster Wheeler Energia Oy.

Reactive dissolution of sedimentary rocks in flue gas desulfurization : modeling and experimental investigation

Svavel förekommer i kol och olja och oxideras vid förbränning till svaveldioxid (SO2). Årligen utsläpps stora mängder svaveldioxid som åstadkommer sura regn, minskning av stratosfäriskt ozon och sjukdomar. Av dessa orsaker är det nödvändigt att minska utsläppen av svaveldioxid. Den teknologi som rör planering av reaktorer och processer för rökgasavsvavling (FGD) har utvecklats kraftigt och idag använder man olika typer av lösningar. De vanligaste alternativen för FGD är våtskrubber- och semitorra skrubbersystem (Spray Dry Scrubbers, SDS) och injektionsprocesser för absorbenter. SDS-processer har studerats bland annat av Ahlbeck [4] och Klingspor [5] och injektionsprocesser av Alvfors [6]. Kalksten, som i huvusak består av kalciumkarbonat, används i rökgasavsvavling på grund av sin förmåga att binda svavel i form av sulfatsalter. Den vanligaste rökgasavsvavlingsmetoden är våtskrubbning, där det sedimentära stenmaterialets upplösningshastighet är en av de faktorer som påverkar resultatet mest. Utvärdering av kalkstensreaktivitet är därför speciellt viktig vid planering och drift av anläggningar för rökgasavsvavling. Målsättningen med detta arbete var att modellera upplösningen av olika typers kalksten för att få en kvantitativ utvärdering av kvaliteten på de analyserade proverna. Därtill testades även karbonatbiprodukter från stålindustrin för att utvärdera möjligheter att använda andra råmaterial. Det transienta förloppet har analyserats, varvid upplösningshastigheten modellerades bl.a. i avseende på tid och pH. Under arbetets gång har antalet empiriska korrelationer minskats till fördel för fysikaliska modeller av diffusiva och konvektiva masstransportfenomen. En målsättning var att skapa en effektiv och snabb metod för att testa olika absorbenter för rökgasavsvavling under transienta förlopp. I arbetet användes PSD-analys, gjordes pH-mätningar och andra utvärderingar av de fysikaliska parametrar som ingår i beräkningarna. On-line mätningar för de icke-stationära variablerna tid och pH ger möjlighet att eliminera osäkerheter. Vissa modeller kan vara komplicerade. En modell för upplösningshastigheten med mer detaljerad utvärdering av parametrar och färre approximationer är därför nödvändig då man vill utvärdera reaktionshastigheten för fasta partiklar i sur miljö. Arbetet utfördes under fyra år och fem peer review-artiklar ingår i avhandlingen.

Molecular modeling of asymmetric induction in heterogeneously catalyzed hydrogenation of the C=O bond

Modifiering av metallytor med starkt adsorberade kirala organiska molekyler är eventuellt den mest relevanta teknik man vet i dag för att skapa kirala ytor. Den kan utnyttjas i katalytisk produktion av enantiomeriskt rena kirala föreningar som behövs t.ex. som läkemedel och aromkemikalier. Trots många fördelar av asymmetrisk heterogen katalys jämfört med andra sätt för att få kirala föreningar, har den ändå inte blivit ett allmänt verktyg för storskaliga tillämpningar. Detta beror t.ex. på brist på djupare kunskaper i katalytiska reaktionsmekanismer och ursprunget för asymmetrisk induktion. I denna studie användes molekylmodelleringstekniker för att studera asymmetriska, heterogena katalytiska system, speciellt hydrering av prokirala karbonylföreningar till motsvarande kirala alkoholer på cinchona-alkaloidmodifierade Pt-katalysatorer. 1-Fenyl-1,2-propandion (PPD) och några andra föreningar, som innehåller en prokiral C=O-grupp, användes som reaktanter. Konformationer av reaktanter och cinchona-alkaloider (som kallas modifierare) samt vätebundna 1:1-komplex mellan dem studerades i gas- och lösningsfas med metoder som baserar sig på vågfunktionsteori och täthetsfunktionalteori (DFT). För beräkningen av protonaffiniteter användes också högst noggranna kombinationsmetoder såsom G2(MP2). Den relativa populationen av modifierarnas konformationer varierade som funktion av modifieraren, dess protonering och lösningsmedlet. Flera reaktant–modifierareinteraktionsgeometrier beaktades. Slutsatserna på riktning av stereoselektivitet baserade sig på den relativa termodynamiska stabiliteten av de diastereomeriska reaktant–modifierare-komplexen samt energierna hos π- och π*-orbitalerna i den reaktiva karbonylgruppen. Adsorption och reaktioner på Pt(111)-ytan betraktades med DFT. Regioselektivitet i hydreringen av PPD och 2,3-hexandion kunde förklaras med molekyl–yta-interaktioner. Storleken och formen av klustret använt för att beskriva Pt-ytan inverkade inte bara på adsorptionsenergierna utan också på de relativa stabiliteterna av olika adsorptionsstrukturer av en molekyl. Populationerna av modifierarnas konformationer i gas- och lösningsfas korrelerade inte med populationerna på Pt-ytan eller med enantioselektiviteten i hydreringen av PPD på Pt–cinchona-katalysatorer. Vissa modifierares konformationer och reaktant–modifierare-interaktionsgeometrier var stabila bara på metallytan. Teoretiskt beräknade potentialenergiprofiler för hydrering av kirala α-hydroxiketoner på Pt implicerade preferens för parvis additionsmekanism för väte och selektiviteter i harmoni med experimenten. De uppnådda resultaten ökar uppfattningen om kirala heterogena katalytiska system och kunde därför utnyttjas i utvecklingen av nya, mera aktiva och selektiva kirala katalysatorer.

CRM system implementation for managing customer relationships in business-to-business markets

Tutkielman tavoitteena oli selvittää, kuinka CRM-järjestelmä otetaan käyttöön asiakkuuksienhallinnan tueksi yritysten välisillä markkinoilla. Tutkielma keskittyy erityisesti käyttöönoton varhaisiin vaiheisiin selvittämällä, mikä on järjestelmien rooli asiakkuudenhallinnassa ja miten ne tukevat erilaisia käyttäjätasoja, jotka työssä jaettiin strategiseksi, operatiiviseksi ja analyyttiseksi tasoksi. Lisäksi työ esittelee käyttöönottoon yleisimmin liittyviä sudenkuoppia. Tutkielma on laadullinen tutkimus ja siinä on kuvattu CRM-järjestelmän käyttöönottoon liittyviä kysymyksiä yritysten välisillä markkinoilla toimivan caseyrityksen näkökulmasta. Tutkimuksessa havaittiin, että CRM-järjestelmän onnistunut käyttöönotto vaatii mahdollisimman yksinkertaista ja helppokäyttöistä järjestelmää. Käyttöönoton ei myöskään tulisi vaatia monimutkaista koulutusta, vaikkakin koulutusten tärkeys tunnistettiin. Tutkimus osoitti, että vaikka tämänhetkiset CRM-järjestelmät keskittyvät enemmän yksityiskohtaisten asiakasanalyysien tekoon, järjestelmä voi myös palvella asiakkaisiin liittyvän tiedon yhteisenä tallennus- ja jakamispaikkana. Lisäksi tutkimuksessa todettiin, että useimmiten potentiaaliset sudenkuopat olivat luonteeltaan hyvin käytännönläheisiä, kuten järjestelmän käyttämättömyys sekä huono motivointi ja sitouttaminen.

Applications of graph transformation in tools for domain-specific modeling languages

The use of domain-specific languages (DSLs) has been proposed as an approach to cost-e ectively develop families of software systems in a restricted application domain. Domain-specific languages in combination with the accumulated knowledge and experience of previous implementations, can in turn be used to generate new applications with unique sets of requirements. For this reason, DSLs are considered to be an important approach for software reuse. However, the toolset supporting a particular domain-specific language is also domain-specific and is per definition not reusable. Therefore, creating and maintaining a DSL requires additional resources that could be even larger than the savings associated with using them. As a solution, di erent tool frameworks have been proposed to simplify and reduce the cost of developments of DSLs. Developers of tool support for DSLs need to instantiate, customize or configure the framework for a particular DSL. There are di erent approaches for this. An approach is to use an application programming interface (API) and to extend the basic framework using an imperative programming language. An example of a tools which is based on this approach is Eclipse GEF. Another approach is to configure the framework using declarative languages that are independent of the underlying framework implementation. We believe this second approach can bring important benefits as this brings focus to specifying what should the tool be like instead of writing a program specifying how the tool achieves this functionality. In this thesis we explore this second approach. We use graph transformation as the basic approach to customize a domain-specific modeling (DSM) tool framework. The contributions of this thesis includes a comparison of di erent approaches for defining, representing and interchanging software modeling languages and models and a tool architecture for an open domain-specific modeling framework that e ciently integrates several model transformation components and visual editors. We also present several specific algorithms and tool components for DSM framework. These include an approach for graph query based on region operators and the star operator and an approach for reconciling models and diagrams after executing model transformation programs. We exemplify our approach with two case studies MICAS and EFCO. In these studies we show how our experimental modeling tool framework has been used to define tool environments for domain-specific languages.

Ballistisen simuloinnin liikeidean kehittämismahdollisuudet

Asevaikutusta laukaisusta kohteeseen simuloiva integroitu laskentaketju koostuu sisä-, ulko- ja maaliballistiikan malleista. Ulkoballistiikka kattaa laskentamallit radanlaskennan, sääkorjauksen ja ammusaerodynamiikan alueilla. Graafisella käyttöliittymällä toteutetulla, fysikaalisesti tarkkaan mallinnukseen perustuvalla ja kokonaisuuden kattavalla laskentajärjestelmällä on kasvavaa tarvetta teknisiä ja koulutuksellisia tarkoituksia varten. Erikoisesti, jos laskentaketjuun lisätään räjähdysvaikutuksen mallintaminen, voidaan simuloida asejärjestelmien vaikutusta kohteessa käyttäjien arvostamalla tavalla. Tietointensiiviset ballistiikan laskentamallit ovat välttämättömiä työkaluja teknisen suunnitteluosaamisen kattamiseksi ja kilpailuedun luomiseksi verkostoituneessa yritysympäristössä. Yliopistotutkimuksen tuottamien laskennallisten menetelmien hyötykäyttö yritysten suunnittelujärjestelmissä syventää teknistä osaamista, jolla on myös henkilöstöä motivoiva vaikutus teknisesti vaikeutuvilla markkinoilla. Työssä arvioidaan toimialaa analysoimalla eri käyttötarpeita samoille tietokantoihin tukeutuville laskentamalleille. Tarkastellaan teknisiä perusteita, käyttöympäristöjä ja markkinoita liiketoimintamahdollisuuksien tunnistamiseksi. Työn tuloksena syvennetään näkemystä ydinosaamisista ja visioidaan liikeidean erottumista kilpailijoista, markkinoita ja sen kehittämistä.

Knowledge modeling for innovative companies: case of business incubator

The paper presents a study which is aimed at building a knowledge model for a case company – business incubator “Ingria” (St. Petersburg, Russia). The business incubator is one of its kind organization in St. Petersburg, and one of the few in Russia, providing services for innovative entrepreneurial companies at an international level. Business incubation impact is deeply researched from the point of view of knowledge engineering. The paper also provides a broad analysis of various knowledge engineering tools used for visualization of knowledge, as well as knowledge modeling techniques.

Charge transport in disordered materials : simulations, theory, and numerical modeling of hopping transport and electron-hole recombination

Min avhandling behandlar hur oordnade material leder elektrisk ström. Bland materialen som studeras finns ledande polymerer, d.v.s. plaster som leder ström, och mer allmänt organiska halvledare. Av de här materialen har man kunnat bygga elektroniska komponenter, och man hoppas på att kunna trycka hela kretsar av organiska material. För de här tillämpningarna är det viktigt att förstå hur materialen själva leder elektrisk ström. Termen oordnade material syftar på material som saknar kristallstruktur. Oordningen gör att elektronernas tillstånd blir lokaliserade i rummet, så att en elektron i ett visst tillstånd är begränsad t.ex. till en molekyl eller ett segment av en polymer. Det här kan jämföras med kristallina material, där ett elektrontillstånd är utspritt över hela kristallen (men i stället har en väldefinierad rörelsemängd). Elektronerna (eller hålen) i det oordnade materialet kan röra sig genom att tunnelera mellan de lokaliserade tillstånden. Utgående från egenskaperna för den här tunneleringsprocessen, kan man bestämma transportegenskaperna för hela materialet. Det här är utgångspunkten för den så kallade hopptransportmodellen, som jag har använt mig av. Hopptransportmodellen innehåller flera drastiska förenklingar. Till exempel betraktas elektrontillstånden som punktformiga, så att tunneleringssannolikheten mellan två tillstånd endast beror på avståndet mellan dem, och inte på deras relativa orientation. En annan förenkling är att behandla det kvantmekaniska tunneleringsproblemet som en klassisk process, en slumpvandring. Trots de här grova approximationerna visar hopptransportmodellen ändå många av de fenomen som uppträder i de verkliga materialen som man vill modellera. Man kan kanske säga att hopptransportmodellen är den enklaste modell för oordnade material som fortfarande är intressant att studera. Man har inte hittat exakta analytiska lösningar för hopptransportmodellen, därför använder man approximationer och numeriska metoder, ofta i form av datorberäkningar. Vi har använt både analytiska metoder och numeriska beräkningar för att studera olika aspekter av hopptransportmodellen. En viktig del av artiklarna som min avhandling baserar sig på är att jämföra analytiska och numeriska resultat. Min andel av arbetet har främst varit att utveckla de numeriska metoderna och applicera dem på hopptransportmodellen. Därför fokuserar jag på den här delen av arbetet i avhandlingens introduktionsdel. Ett sätt att studera hopptransportmodellen numeriskt är att direkt utföra en slumpvandringsprocess med ett datorprogram. Genom att föra statisik över slumpvandringen kan man beräkna olika transportegenskaper i modellen. Det här är en så kallad Monte Carlo-metod, eftersom själva beräkningen är en slumpmässig process. I stället för att följa rörelsebanan för enskilda elektroner, kan man beräkna sannolikheten vid jämvikt för att hitta en elektron i olika tillstånd. Man ställer upp ett system av ekvationer, som relaterar sannolikheterna för att hitta elektronen i olika tillstånd i systemet med flödet, strömmen, mellan de olika tillstånden. Genom att lösa ekvationssystemet fås sannolikhetsfördelningen för elektronerna. Från sannolikhetsfördelningen kan sedan strömmen och materialets transportegenskaper beräknas. En aspekt av hopptransportmodellen som vi studerat är elektronernas diffusion, d.v.s. deras slumpmässiga rörelse. Om man betraktar en samling elektroner, så sprider den med tiden ut sig över ett större område. Det är känt att diffusionshastigheten beror av elfältet, så att elektronerna sprider sig fortare om de påverkas av ett elektriskt fält. Vi har undersökt den här processen, och visat att beteendet är väldigt olika i endimensionella system, jämfört med två- och tredimensionella. I två och tre dimensioner beror diffusionskoefficienten kvadratiskt av elfältet, medan beroendet i en dimension är linjärt. En annan aspekt vi studerat är negativ differentiell konduktivitet, d.v.s. att strömmen i ett material minskar då man ökar spänningen över det. Eftersom det här fenomenet har uppmätts i organiska minnesceller, ville vi undersöka om fenomenet också kan uppstå i hopptransportmodellen. Det visade sig att det i modellen finns två olika mekanismer som kan ge upphov till negativ differentiell konduktivitet. Dels kan elektronerna fastna i fällor, återvändsgränder i systemet, som är sådana att det är svårare att ta sig ur dem då elfältet är stort. Då kan elektronernas medelhastighet och därmed strömmen i materialet minska med ökande elfält. Elektrisk växelverkan mellan elektronerna kan också leda till samma beteende, genom en så kallad coulombblockad. En coulombblockad kan uppstå om antalet ledningselektroner i materialet ökar med ökande spänning. Elektronerna repellerar varandra och ett större antal elektroner kan leda till att transporten blir långsammare, d.v.s. att strömmen minskar.

Use data from permanent sampling points in growth and yield modeling

This study was carried to evaluate the efficiency of the Bitterlich method in growth and yield modeling of the even-aged Eucalyptus stands. 25 plots were setup in Eucalyptus grandis cropped under a high bole system in the Central Western Region of Minas Gerais, Brazil. The sampling points were setup in the center of each plot. The data of four annual mesurements were colleted and used to adjust the three model types using the age, the site index and the basal area as independent variables. The growths models were fitted for volume and mass of trees. The efficiency of the Bitterlich method was confirmed for generating the data for growth and yield modeling.

Hipsometric relationship modeling using data sampled in tree scaling and inventory plots

This work evaluated eight hypsometric models to represent tree height-diameter relationship, using data obtained from the scaling of 118 trees and 25 inventory plots. Residue graphic analysis and percent deviation mean criteria, qui-square test precision, residual standard error between real and estimated heights and the graybill f test were adopted. The identity of the hypsometric models was also verified by applying the F(Ho) test on the plot data grouped to the scaling data. It was concluded that better accuracy can be obtained by using the model prodan, with h and d1,3 data measured in 10 trees by plots grouped into these scaling data measurements of even-aged forest stands.

Modeling and Analysis of Noise and Interconnects for On-Chip Communication Link Design

This thesis considers modeling and analysis of noise and interconnects in onchip communication. Besides transistor count and speed, the capabilities of a modern design are often limited by on-chip communication links. These links typically consist of multiple interconnects that run parallel to each other for long distances between functional or memory blocks. Due to the scaling of technology, the interconnects have considerable electrical parasitics that affect their performance, power dissipation and signal integrity. Furthermore, because of electromagnetic coupling, the interconnects in the link need to be considered as an interacting group instead of as isolated signal paths. There is a need for accurate and computationally effective models in the early stages of the chip design process to assess or optimize issues affecting these interconnects. For this purpose, a set of analytical models is developed for on-chip data links in this thesis. First, a model is proposed for modeling crosstalk and intersymbol interference. The model takes into account the effects of inductance, initial states and bit sequences. Intersymbol interference is shown to affect crosstalk voltage and propagation delay depending on bus throughput and the amount of inductance. Next, a model is proposed for the switching current of a coupled bus. The model is combined with an existing model to evaluate power supply noise. The model is then applied to reduce both functional crosstalk and power supply noise caused by a bus as a trade-off with time. The proposed reduction method is shown to be effective in reducing long-range crosstalk noise. The effects of process variation on encoded signaling are then modeled. In encoded signaling, the input signals to a bus are encoded using additional signaling circuitry. The proposed model includes variation in both the signaling circuitry and in the wires to calculate the total delay variation of a bus. The model is applied to study level-encoded dual-rail and 1-of-4 signaling. In addition to regular voltage-mode and encoded voltage-mode signaling, current-mode signaling is a promising technique for global communication. A model for energy dissipation in RLC current-mode signaling is proposed in the thesis. The energy is derived separately for the driver, wire and receiver termination.

From Fluid Dynamics to Human Psychology.What Drives Financial Markets Towards Extreme Events

For decades researchers have been trying to build models that would help understand price performance in financial markets and, therefore, to be able to forecast future prices. However, any econometric approaches have notoriously failed in predicting extreme events in markets. At the end of 20th century, market specialists started to admit that the reasons for economy meltdowns may originate as much in rational actions of traders as in human psychology. The latter forces have been described as trading biases, also known as animal spirits. This study aims at expressing in mathematical form some of the basic trading biases as well as the idea of market momentum and, therefore, reconstructing the dynamics of prices in financial markets. It is proposed through a novel family of models originating in population and fluid dynamics, applied to an electricity spot price time series. The main goal of this work is to investigate via numerical solutions how well theequations succeed in reproducing the real market time series properties, especially those that seemingly contradict standard assumptions of neoclassical economic theory, in particular the Efficient Market Hypothesis. The results show that the proposed model is able to generate price realizations that closely reproduce the behaviour and statistics of the original electricity spot price. That is achieved in all price levels, from small and medium-range variations to price spikes. The latter were generated from price dynamics and market momentum, without superimposing jump processes in the model. In the light of the presented results, it seems that the latest assumptions about human psychology and market momentum ruling market dynamics may be true. Therefore, other commodity markets should be analyzed with this model as well.

Modeling of electrolyte sorption – from phase equilibria to dynamic separation systems

In this thesis, general approach is devised to model electrolyte sorption from aqueous solutions on solid materials. Electrolyte sorption is often considered as unwanted phenomenon in ion exchange and its potential as an independent separation method has not been fully explored. The solid sorbents studied here are porous and non-porous organic or inorganic materials with or without specific functional groups attached on the solid matrix. Accordingly, the sorption mechanisms include physical adsorption, chemisorption on the functional groups and partition restricted by electrostatic or steric factors. The model is tested in four Cases Studies dealing with chelating adsorption of transition metal mixtures, physical adsorption of metal and metalloid complexes from chloride solutions, size exclusion of electrolytes in nano-porous materials and electrolyte exclusion of electrolyte/non-electrolyte mixtures. The model parameters are estimated using experimental data from equilibrium and batch kinetic measurements, and they are used to simulate actual single-column fixed-bed separations. Phase equilibrium between the solution and solid phases is described using thermodynamic Gibbs-Donnan model and various adsorption models depending on the properties of the sorbent. The 3-dimensional thermodynamic approach is used for volume sorption in gel-type ion exchangers and in nano-porous adsorbents, and satisfactory correlation is obtained provided that both mixing and exclusion effects are adequately taken into account. 2-Dimensional surface adsorption models are successfully applied to physical adsorption of complex species and to chelating adsorption of transition metal salts. In the latter case, comparison is also made with complex formation models. Results of the mass transport studies show that uptake rates even in a competitive high-affinity system can be described by constant diffusion coefficients, when the adsorbent structure and the phase equilibrium conditions are adequately included in the model. Furthermore, a simplified solution based on the linear driving force approximation and the shrinking-core model is developed for very non-linear adsorption systems. In each Case Study, the actual separation is carried out batch-wise in fixed-beds and the experimental data are simulated/correlated using the parameters derived from equilibrium and kinetic data. Good agreement between the calculated and experimental break-through curves is usually obtained indicating that the proposed approach is useful in systems, which at first sight are very different. For example, the important improvement in copper separation from concentrated zinc sulfate solution at elevated temperatures can be correctly predicted by the model. In some cases, however, re-adjustment of model parameters is needed due to e.g. high solution viscosity.

Developing markets of energy biomass – local and global perspectives

The thesis explores global and national-level issues related to the development of markets for biomass for energy. The thesis consists of five separate papers and provides insights on selected issues. The aim of Paper I was to identify methodological and statistical challenges in assessing international solid and liquid biofuels trade and provide an overview of the Finnish situation with respect to the status of international solid and liquid biofuels trade. We found that, for the Finnish case, it is possible to qualify direct and indirect trade volumes of biofuels. The study showed that indirect trade of biofuels has a highly significant role in Finland and may be a significant sector also in global biofuels trade. The purpose of Paper II was to provide a quantified insight into Finnish prospects for meeting the national 2020 renewable energy targets and concurrently becoming a largescale producer of forest-biomass-based second-generation biofuels for feeding increasing demand in European markets. We found that Finland has good opportunities to realise a scenario to meet 2020 renewable energy targets and for large-scale production of wood-based biofuels. The potential net export of transport biofuels from Finland in 2020 would correspond to 2–3% of European demand. Paper III summarises the global status of international solid and liquid biofuels trade as illuminated by several separate sources. International trade of biofuels was estimated at nearly 1 EJ for 2006. Indirect trade of biofuels through trading of industrial roundwood and material by-products comprises the largest proportion of the trading, with a share of about two thirds. The purpose of Paper IV was to outline a comprehensive picture of the coverage of various certification schemes and sustainability principles relating to the entire value-added chain of biomass and bioenergy. Regardless of the intensive work that has been done in the field of sustainability schemes and principles concerning use of biomass for energy, weaknesses still exist. The objective of Paper V was to clarify the alternative scenarios for the international biomass market until 2020 and identify the underlying steps needed toward a wellfunctioning and sustainable market for biomass for energy purposes. An overall conclusion drawn from this analysis concerns the enormous opportunities related to the utilisation of biomass for energy in the coming decades.