Integrated techno-economic and environmental assessments of sixty scenarios for co-firing biomass with coal and natural gas

Displacement of fossil fuel-based power through biomass co-firing could reduce the greenhouse gas (GHG) emissions from fossil fuels. In this study, data-intensive techno-economic models were developed to evaluate different co-firing technologies as well as the configurations of these technologies. The models were developed to study 60 different scenarios involving various biomass feedstocks (wood chips, wheat straw, and forest residues) co-fired either with coal in a 500 MW subcritical pulverized coal (PC) plant or with natural gas in a 500 MW natural gas combined cycle (NGCC) plant to determine their technical potential and costs, as well as to determine environmental benefits. The results obtained reveal that the fully paid-off coal-fired power plant co-fired with forest residues is the most attractive option, having levelized costs of electricity (LCOE) of $53.12–$54.50/MW h and CO2 abatement costs of $27.41–$31.15/tCO2. When whole forest chips are co-fired with coal in a fully paid-off plant, the LCOE and CO2 abatement costs range from $54.68 to $56.41/MW h and $35.60 to $41.78/tCO2, respectively. The LCOE and CO2 abatement costs for straw range from $54.62 to $57.35/MW h and $35.07 to $38.48/tCO2, respectively.

A Schamel equation for ion acoustic waves in superthermal plasmas

An investigation of the propagation of ion acoustic waves in nonthermal plasmas in the presence of trapped electrons has been undertaken. This has been motivated by space and laboratory plasma observations of plasmas containing energetic particles, resulting in long-tailed distributions, in combination with trapped particles, whereby some of the plasma particles are confined to a finite region of phase space. An unmagnetized collisionless electron-ion plasma is considered, featuring a non-Maxwellian-trapped electron distribution, which is modelled by a kappa distribution function combined with a Schamel distribution. The effect of particle trapping has been considered, resulting in an expression for the electron density. Reductive perturbation theory has been used to construct a KdV-like Schamel equation, and examine its behaviour. The relevant configurational parameters in our study include the superthermality index κ and the characteristic trapping parameter β. A pulse-shaped family of solutions is proposed, also depending on the weak soliton speed increment u₀. The main modification due to an increase in particle trapping is an increase in the amplitude of solitary waves, yet leaving their spatial width practically unaffected. With enhanced superthermality, there is a decrease in both amplitude and width of solitary waves, for any given values of the trapping parameter and of the incremental soliton speed. Only positive polarity excitations were observed in our parametric investigation. 

An optimized process design for oxymethylene ether production from woody-biomass-derived syngas

The conversion of biomass for the production of liquid fuels can help reduce the greenhouse gas (GHG) emissions that are predominantly generated by the combustion of fossil fuels. Oxymethylene ethers (OMEs) are a series of liquid fuel additives that can be obtained from syngas, which is produced from the gasification of biomass. The blending of OMEs in conventional diesel fuel can reduce soot formation during combustion in a diesel engine. In this research, a process for the production of OMEs from woody biomass has been simulated. The process consists of several unit operations including biomass gasifi- cation, syngas cleanup, methanol production, and conversion of methanol to OMEs. The methodology involved the development of process models, the identification of the key process parameters affecting OME production based on the process model, and the development of an optimal process design for high OME yields. It was found that up to 9.02 tonnes day1 of OME3, OME4, and OME5 (which are suitable as diesel additives) can be produced from 277.3 tonnes day1 of wet woody biomass. Furthermore, an optimal combination of the parameters, which was generated from the developed model, can greatly enhance OME production and thermodynamic efficiency. This model can further be used in a techno- economic assessment of the whole biomass conversion chain to produce OMEs. The results of this study can be helpful for petroleum-based fuel producers and policy makers in determining the most attractive pathways of converting bio-resources into liquid fuels.

Novel furan polymer materials based on the reversible Diels-Alder reaction

A futura e inevitável escassez dos recursos fósseis, juntamente com o aumento imprevisível dos seus preços, levou, nas últimas décadas, a um aumento impressionante de iniciativas dedicadas não só à procura de fontes alternativas de fornecedores de energia, mas também de produtos químicos e polímeros a partir de fontes renováveis, em particular da biomassa vegetal. Entre estes, os polímeros derivados de monómeros furânicos constituem uma classe única de materiais cujas estruturas podem, em princípio, simular virtualmente os seus homólogos actualmente derivados de recursos fósseis. O anel furânico é uma estrutura heterocíclica com um carácter diénico pronunciado, o que torna-o um dieno particularmente apropriado para a reacção de Diels-Alder (DA) com dienófilos como a maleimida. Um dos aspectos mais relevantes da reacção de DA é a sua reversibilidade em função da temperatura, a qual permite que os aductos sejam facilmente revertidos nos seus precursores por aumento da temperatura (reacção de retro-DA). No caso específico da combinação furano-maleimida, a formação do aducto predomina até cerca de 60ºC, enquanto a reacção inversa é dominante acima de 100ºC. A combinação desta característica da reacção de DA com a química de compostos furânicos pode abrir um novo caminho para a preparação de materiais macromoleculares funcionais com base em fontes renováveis e com aplicações promissoras como auto-reparação e reciclabilidade. O principal objectivo desta Tese, é a síntese e caracterização de novos materiais poliméricos termo-reversíveis, aplicando a reacção de DA a monómeros complementares com estruturas dos tipos furânico (o dieno, designado por A) e de maleimida (o dienófilo, designado por B). A primeira etapa neste trabalho envolveu a síntese, purificação e caracterização de novos monómeros furânicos e de maleimida do tipo AA, A3, BB, B3, AB, AB2, cada um com diferentes grupos separadores das funções reactivas. Posteriormente, estes monómeros foram polimerizados e despolimerizados por ciclos de DA/retro-DA utilizando diferentes combinações. A formação e dissociação de todos os aductos de DA foram seguidas por ambas espectroscopias de UV e RMN de 1H. O primeiro sistema de DA estudado foi uma combinação modelo entre reagentes mono-funcionais (-A+-B), nomeadamente o acetato furfurílico (FA) e a N-metilmaleimida (MM), ambos comercialmente disponíveis. O objectivo desta abordagem foi estudar a cinética e o equilíbrio da formação/dissociação dos aductos de DA e obter indicações sobre as condições mais adequadas a serem usadas na preparação dos correspondentes novos materiais macromoleculares. Além disso, pretendia-se verificar a presença ou ausência de reacções secundárias que poderiam intervir em ambas as vias directa e inversa das reacções, mesmo após vários ciclos. A espectroscopia de UV forneceu informação quantitativa sobre a cinética de formação do aducto através da diminuição progressiva da absorvência máxima a 293 nm correspondente ao grupo maleimida, a diferentes temperaturas (35, 50, 65 ºC) Reciprocamente, a correspondente reacção de retro-DA foi seguida a 90 ºC através do aumento do mesmo pico. A reversibilidade destes sistemas foi verificada com sucesso após uma sequência de ciclos de DA/retro-DA. Adicionalmente, verificou-se que os espectros originaram um ponto isosbéstico, provando que estes sistemas não envolvem quaisquer reacções secundárias. Uma vez que foi usado um excesso de FA, as reacções de DA modelo apresentaram um comportamento cinético de pseudo-primeira ordem, com a constante de velocidade k mais alta (2.1x10-5 dm3mol-1s-1) para T=65 ºC. A correspondente energia de activação foi de 39.0 kJ.mol-1. A reacção de retro-DA seguiu um comportamento de primeira ordem, com constante de velocidade de 1.6x10-6 s-1. A evolução deste sistema por RMN de 1H a 65ºC deu-nos informações mais detalhadas sobre a sua evolução estrutural, ou seja, à medida que a intensidade dos picos atribuídos à formação do aducto aumentaram progressivamente ao longo do tempo, os pertencentes aos reagentes iniciais diminuiram proporcionalmente. O “rendimento final”, calculado após 20 dias à temperatura ambiente, foi de aproximadamente 70%. A reacção de retro-DA foi depois seguida a 90ºC, observando-se tal como na espectroscopia de UV, o deslocamento da reacção no sentido da regeneração dos reagentes de partida. A viabilidade de múltiplos ciclos de DA/retro-DA estabelecidos pela espectroscopia de UV foi igualmente confirmada por RMN de 1H. O passo seguinte envolveu o estudo de um sistema de policondensação linear baseado no crescimento gradual por reacção de DA entre um monómero bisfurânico A-A e um do tipo bismaleimida B-B, seguindo a mesma abordagem que no sistema modelo. O poliaducto linear foi obtido a partir de soluções equimolares dos monómeros, por reacção de DA a 65ºC. O progresso desta polimerização foi seguido por espectroscopia de UV e RMN de 1H e, mais qualitativamente, pelo aumento da viscosidade do meio. A reacção seguiu um comportamento de segunda ordem, com uma constante de velocidade de 9.4x10-6 dm3mol-1s-1, e observou-se novamente um ponto isosbéstico nos dados de UV. Os espectros de RMN apresentaram o padrão esperado, nomeadamente o aumento progressivo dos sinais associados ao aducto e a correspondente diminuição dos grupos furano e maleimida livres. A despolimerização do poliaducto através da reacção de retro-DA foi seguida a 110ºC usando as mesmas técnicas. Os dados de UV mostraram o retorno progressivo da absorção dos grupos de maleimida, seguindo um comportamento cinético de primeira ordem, com constante de velocidade de 2.5x10-6 s-1, até à completa regeneração de ambos os monómeros. Os espectros de RMN providenciaram mais uma vez informação estrutural sobre o progresso da despolimerização, a qual foi acompanhada por uma diminuição progressiva da viscosidade. Adicionalmente, para seguir a retro- DA, adicionou-se um excesso de composto furânico monofuncional, nomeadamente o 2,5-dimetilfurano (DMFu), ao sistema de modo a bloquear as funções maleimida complementares, evitando assim a repolimerização após arrefecimento. Os productos isolados foram então o monómero bisfurânico AA, DMFu que não reagiu e o bisaducto não-polimerizável de BB com DMFu. Este resultado indicou claramente que o polímero foi de facto revertido nos seus monómeros durante a reacção de retro-DA. O terceiro sistema estudado foi outra polimerização linear, seguindo as mesmas condições experimentais que os anteriores, mas com uma estratégia diferente de modo a contornar o problema clássico de assegurar a estequiometria exacta dos monómeros. As estruturas dos monómeros utilizados incorporam ambos os grupos reactivos, i.e, moléculas do tipo A-B. A polimerização prematura destes monómeros intrinsecamente reactivos foi evitada com a protecção do grupo maleimida na forma de um aducto de DA com furano, até a incorporação do substituinte furânico na outra extremidade. Portanto, a policondensação destes monómeros foi iniciada após a desprotecção in situ deste composto mediante aquecimento, seguido de arrefecimento até à temperatura adequada para polimerizar. Os resultados obtidos por UV e RMN sugerem que de facto o uso de monómeros do tipo A-B oferece um melhor sistema linear. Em seguida, foram estudados sistemas de policondensação não-linear por reacção de DA, entre monómeros (um ou ambos) com funcionalidade superior a dois, nomeadamente sistemas do tipo A3+B-B ou A-A+B3, seguindo mais uma vez as mesmas condições experimentais. Uma vez que utilizam monómeros complementares contendo, em média, mais de duas funcionalidades, estes sistemas conduzem a materiais reticulados. Nestes estudos, foram usadas três razões molares de [maleimida]/[furano], nomeadamente 1.0, 0.75 e 0.5, de modo a estudar ambas as situações de não-gelificação e reticulação. Ambos sistemas apresentaram um comportamento regular e boa reciclabilidade quer para gerar situações que possam conduzir à formação de redes a diferentes graus de conversão, ou que possam parar antes da sua obtenção, conforme previsto pela equação de Flory-Stockmayer. Como esperado, a utilização de grupos complementares em quantidades estequiométricas produziu o espessamento mais rápido e a reticulação quase completa; à medida que a quantidade relativa de monómero trifuncional decresceu, as reacções pararam antes da reticulação, ou seja, originaram meios altamente viscosos contendo polímeros solúveis altamente ramificados. As reacções de retro-DA a 110 ºC conduziram à gradual dissolução das partículas de gel (quando presentes), tendo sido comprovado pelos espectros de UV e de RMN de 1H, evidenciado a regeneração dos monómeros. Tal como no sistema do tipo A-A+B-B, a reacção de retro-DA foi seguida adicionando um excesso de DMFu ao sistema reaccional. Como esperado, os produtos finais foram os monómeros furânicos, o DMFu em excesso e o trisaducto ou o bisaducto maleimida-DMFu, o que confirma a eficiência da despolimerização com regeneração dos monómeros iniciais. O último sistema de policondensação por reacção de DA envolveu um monómero assimetricamente substituído do tipo AB2, capaz de originar estruturas macromoleculares hiper-ramificadas que não reticulam. Este estudo preliminar deste sistema foi seguido nas mesmas condições experimentais que os anteriores, apresentando um comportamento com as características esperadas.

Biomass fly ash incorporation in cement based materials

In recent years, pressures on global environment and energy security have led to an increasing demand on renewable energy sources, and diversification of Europe’s energy supply. Among these resources the biomass could exert an important role, since it is considered a renewable and CO2 neutral energy resource once the consumption rate is lower than the growth rate, and can potentially provide energy for heat, power and transports from the same installation. Currently, most of the biomass ash produced in industrial plants is either disposed of in landfill or recycled on agricultural fields or forest, and most times this goes on without any form of control. However, considering that the disposal cost of biomass ashes are raising, and that biomass ash volumes are increasing worldwide, a sustainable ash management has to be established. The main objective of the present study is the effect of biomass fly ashes in cement mortars and concretes in order to be used as a supplementary cementitious material. The wastes analyzed in the study were collected from the fluidized bed boilers and grate boilers available in the thermal power plants and paper pulp plants situated in Portugal. The physical as well as chemical characterisations of the biomass fly ashes were investigated. The cement was replaced by the biomass fly ashes in 10, 20 and 30% (weight %) in order to investigate the fresh properties as well as the hardened properties of biomass fly ash incorporated cement mortar and concrete formulations. Expansion reactions such as alkali silica reaction (ASR), sulphate attack (external and internal) were conducted in order to check the durability of the biomass fly ash incorporated cement mortars and concretes. Alternative applications such as incorporation in lime mortars and alkali activation of the biomass fly ashes were also attempted. The biomass fly ash particles were irregular in shape and fine in nature. The chemical characterization revealed that the biomass fly ashes were similar to a class C fly ash. The mortar results showed a good scope for biomass fly ashes as supplementary cementitious materials in lower dosages (<20%). The poor workability, concerns about the organic content, alkalis, chlorides and sulphates stand as the reasons for preventing the use of biomass fly ash in high content in the cement mortars. The results obtained from the durability tests have shown a clear reduction in expansion for the biomass fly ash mortars/concretes and the binder blend made with biomass fly ash (20%) and metakaolin (10%) inhibited the ASR reaction effectively. The biomass fly ash incorporation in lime mortars did not improve the mortar properties significantly though the carbonation was enhanced in the 15-20% incorporation. The biomass fly ash metakaolin blend worked well in the alkali activated complex binder application also. Portland cement free binders (with 30-40 MPa compressive strength) were obtained on the alkali activation of biomass fly ashes (60-80%) blended with metakaolin (20-40%).

Contribution of biomass combustion to air pollutant emissions

In Portugal, it was estimated that around 1.95 Mton/year of wood is used in residential wood burning for heating and cooking. Additionally, in the last decades, burnt forest area has also been increasing. These combustions result in high levels of toxic air pollutants and a large perturbation of atmospheric chemistry, interfere with climate and have adverse effects on health. Accurate quantification of the amounts of trace gases and particulate matter emitted from residential wood burning, agriculture and garden waste burning and forest fires on a regional and global basis is essential for various purposes, including: the investigation of several atmospheric processes, the reporting of greenhouse gas emissions, and quantification of the air pollution sources that affect human health at regional scales. In Southern Europe, data on detailed emission factors from biomass burning are rather inexistent. Emission inventories and source apportionment, photochemical and climate change models use default values obtained for US and Northern Europe biofuels. Thus, it is desirable to use more specific locally available data. The objective of this study is to characterise and quantify the contribution of biomass combustion sources to atmospheric trace gases and aerosol concentrations more representative of the national reality. Laboratory (residential wood combustion) and field (agriculture/garden waste burning and experimental wildland fires) sampling experiments were carried out. In the laboratory, after the selection of the most representative wood species and combustion equipment in Portugal, a sampling program to determine gaseous and particulate matter emission rates was set up, including organic and inorganic aerosol composition. In the field, the smoke plumes from agriculture/garden waste and experimental wildland fires were sampled. The results of this study show that the combustion equipment and biofuel type used have an important role in the emission levels and composition. Significant differences between the use of traditional combustion equipment versus modern equipments were also observed. These differences are due to higher combustion efficiency of modern equipment, reflecting the smallest amount of particulate matter, organic carbon and carbon monoxide released. With regard to experimental wildland fires in shrub dominated areas, it was observed that the largest organic fraction in the samples studied was mainly composed by vegetation pyrolysis products. The major organic components in the smoke samples were pyrolysates of vegetation cuticles, mainly comprising steradienes and sterol derivatives, carbohydrates from the breakdown of cellulose, aliphatic lipids from vegetation waxes and methoxyphenols from the lignin thermal degradation. Despite being a banned practice in our country, agriculture/garden waste burning is actually quite common. To assess the particulate matter composition, the smoke from three different agriculture/garden residues have been sampled into 3 different size fractions (PM2.5, PM2.5-10 and PM>10). Despite distribution patterns of organic compounds in particulate matter varied among residues, the amounts of phenolics (polyphenol and guaiacyl derivatives) and organic acids were always predominant over other organic compounds in the organosoluble fraction of smoke. Among biomarkers, levoglucosan, β-sitosterol and phytol were detected in appreciable amounts in the smoke of all agriculture/garden residues. In addition, inositol may be considered as an eventual tracer for the smoke from potato haulm burning. It was shown that the prevailing ambient conditions (such as high humidity in the atmosphere) likely contributed to atmospheric processes (e.g. coagulation and hygroscopic growth), which influenced the particle size characteristics of the smoke tracers, shifting their distribution to larger diameters. An assessment of household biomass consumption was also made through a national scale survey. The information obtained with the survey combined with the databases on emission factors from the laboratory and field tests allowed us to estimate the pollutant amounts emitted in each Portuguese district. In addition to a likely contribution to the improvement of emission inventories, emission factors obtained for tracer compounds in this study can be applied in receptor models to assess the contribution of biomass burning to the levels of atmospheric aerosols and their constituents obtained in monitoring campaigns in Mediterranean Europe.

Evaluation of thermochemical biomass conversion in fluidized bed

Dado o aumento acelerado dos preços dos combustíveis fósseis e as incertezas quanto à sua disponibilidade futura, tem surgido um novo interesse nas tecnologias da biomassa aplicadas à produção de calor, eletricidade ou combustíveis sintéticos. Não obstante, para a conversão termoquímica de uma partícula de biomassa sólida concorrem fenómenos bastante complexos que levam, em primeiro lugar, à secagem do combustível, depois à pirólise e finalmente à combustão ou gasificação propriamente ditas. Uma descrição relativamente incompleta de alguns desses estágios de conversão constitui ainda um obstáculo ao desenvolvimento das tecnologias que importa ultrapassar. Em particular, a presença de elevados conteúdos de matéria volátil na biomassa põe em evidência o interesse prático do estudo da pirólise. A importância da pirólise durante a combustão de biomassa foi evidenciada neste trabalho através de ensaios realizados num reator piloto de leito fluidizado borbulhante. Verificou-se que o processo ocorre em grande parte à superfície do leito com chamas de difusão devido à libertação de voláteis, o que dificulta o controlo da temperatura do reator acima do leito. No caso da gasificação de biomassa a pirólise pode inclusivamente determinar a eficiência química do processo. Isso foi mostrado neste trabalho durante ensaios de gasificação num reator de leito fluidizado de 2MWth, onde um novo método de medição permitiu fechar o balanço de massa ao gasificador e monitorizar o grau de conversão da biomassa. A partir destes resultados tornou-se clara a necessidade de descrever adequadamente a pirólise de biomassa com vista ao projeto e controlo dos processos. Em aplicações de engenharia há particular interesse na estequiometria e propriedades dos principais produtos pirolíticos. Neste trabalho procurou-se responder a esta necessidade, inicialmente através da estruturação de dados bibliográficos sobre rendimentos de carbonizado, líquidos pirolíticos e gases, assim como composições elementares e poderes caloríficos. O resultado traduziu-se num conjunto de parâmetros empíricos de interesse prático que permitiram elucidar o comportamento geral da pirólise de biomassa numa gama ampla de condições operatórias. Para além disso, propôs-se um modelo empírico para a composição dos voláteis que pode ser integrado em modelos compreensivos de reatores desde que os parâmetros usados sejam adequados ao combustível ensaiado. Esta abordagem despoletou um conjunto de ensaios de pirólise com várias biomassas, lenhina e celulose, e temperaturas entre os 600 e 975ºC. Elevadas taxas de aquecimento do combustível foram alcançadas em reatores laboratoriais de leito fluidizado borbulhante e leito fixo, ao passo que um sistema termo-gravimétrico permitiu estudar o efeito de taxas de aquecimento mais baixas. Os resultados mostram que, em condições típicas de processos de combustão e gasificação, a quantidade de voláteis libertada da biomassa é pouco influenciada pela temperatura do reator mas varia bastante entre combustíveis. Uma análise mais aprofundada deste assunto permitiu mostrar que o rendimento de carbonizado está intimamente relacionado com o rácio O/C do combustível original, sendo proposto um modelo simples para descrever esta relação. Embora a quantidade total de voláteis libertada seja estabelecida pela composição da biomassa, a respetiva composição química depende bastante da temperatura do reator. Rendimentos de espécies condensáveis (água e espécies orgânicas), CO2 e hidrocarbonetos leves descrevem um máximo relativamente à temperatura para dar lugar a CO e H2 às temperaturas mais altas. Não obstante, em certas gamas de temperatura, os rendimentos de algumas das principais espécies gasosas (e.g. CO, H2, CH4) estão bem correlacionados entre si, o que permitiu desenvolver modelos empíricos que minimizam o efeito das condições operatórias e, ao mesmo tempo, realçam o efeito do combustível na composição do gás. Em suma, os ensaios de pirólise realizados neste trabalho permitiram constatar que a estequiometria da pirólise de biomassa se relaciona de várias formas com a composição elementar do combustível original o que levanta várias possibilidades para a avaliação e projeto de processos de combustão e gasificação de biomassa.

Extraction of added-value products from biomass using ionic liquids

The main purpose of this thesis is to investigate the potential of ionic liquids (ILs) as a new class of extractive solvents for added-value products from biomass. These include phenolic compounds (vanillin, gallic, syringic and vanillic acids), alkaloids (caffeine) and aminoacids (L-tryptophan). The interest on these natural compounds relies on the wide variety of relevant properties shown by those families and further application in the food, cosmetic and pharmaceutical industries. Aiming at developping more benign and effective extraction/purification techniques than those used, a comprehensive study was conducted using aqueous biphasic systems (ABS) composed of ILs and inorganic/organic salts. In addition, ILs were characterized by a polarity scale, using solvatochromic probes, aiming at providing prior indications on the ILs affinity for particular added-value products. Solid-liquid (S-L) extractions from biomass and using aqueous solution of ILs were also investigated. In particular, and applying and experimental factorial design to optimize the operational conditions, caffeine was extracted from guaraná seeds and spent coffee. With both types of extractions it was found that it is possible to recover the high-value compounds and to recycle the IL and salt solutions. Finally, aiming at exploring the recovery of added-value compounds from biomass using a simpler and more suistainable technique, the solubility of gallic acid, vanillin and caffeine was studied in aqueous solutions of several ILs and common salts. With the gathered results it was possible to demonstrate that ILs act as hydrotropes and that water can be used as an adequate antisolvent. This thesis describes the use of ILs towards the development of more effective and sustainable processes.

New synthetic methodologies for the transformation of biomass derived intermediates to valuable molecules

Tese de doutoramento, Farmácia (Química Farmacêutica e Terapêutica), Universidade de Lisboa, Faculdade de Farmácia, 2014

Bacterial Biomass along the Kuroshio Extension

Senior thesis written for Oceanography 445

Logistics management and costs of biomass fuel supply

Discusses part of a project conducted by the authors into the logistics planning and management and costs of supplying biomass fuels to biomass-fired power stations in the UK. Defines biomass fuels and the reasons for the growth in interest in their use for electricity generation. The activities and parties involved in the biomass fuel supply chain are discussed together with the management of the chain in order to achieve smooth and consistent flow of biomass fuel to power stations. Explains the approach used to modelling the delivered costs of biomass fuels for four types of biomass fuel included in the project: forest fuel, short rotation coppice, straw and miscanthus. Comments are given on the environmental impacts of the fuel supply chains. The results indicate that straw supply systems are capable of producing the lowest delivered costs of the four fuels studied. Short rotation coppice and miscanthus, two new energy crops, are likely to have the highest delivered costs at present. This is due to the cost of growing these fuels and the financial incentives required by farmers to persuade them to grow these crops. Logistics costs (i.e. transport, storage and handling) are shown to represent a significant proportion of total delivered cost in biomass supply. Careful supply chain planning and logistics management will be of central importance to the success of the biomass industry.