Numerical Simulation of Combustion and Gasification of Biomass Particles

Transcription

1 Numerical Simulation of Combustion and Gasification of Biomass Particles DOCTORAL DISSERTATION Hesameddin Fatehi Division of Fluid Mechanics, Department of Energy Sciences Faculty of Engineering LTH, Lund University Lund, Sweden, 2014

2 Hesameddin Fatehi Division of Fluid Mechanics Department of Energy Sciences Faculty of Engineering Lund University Box 118 SE Lund University Sweden ISBN (printed) ISBN (pdf) ISSN ISRN LUTMDN/TMHP-14/1102-SE Printed in Sweden by Media-Tryck, Lund University Lund 2014 ii

3 Abstract In this thesis, a numerical approach is adopted to study biomass thermochemical conversion. Detailed physical and chemical processes involved in the thermochemical conversion of biomass are considered. The aims are to improve the understanding of the physical and chemical processes involved and to develop and validate mathematical models for numerical simulation of the biomass conversion process. The main focus of the thesis work has been on large biomass particles under fixed bed conditions. The thermochemical conversion of single particle is first considered. A comprehensive detailed model is developed and evaluated; the results provide valuable insight into the underlying physics of thermochemical conversion of biomass. Based on the comprehensive model a simplified model is proposed which takes into account some of the detailed reaction and transport process inside the particle at high computational efficiency. A two-dimensional model for the fuel bed of biomass furnace is developed and validated. Mathematical description of various sub-processes involved in the conversion process is presented taking into account the main features of biomass conversion. Various approaches and assumptions for modeling the different sub-processes are assessed by comparing the numerical results with experimental measurements. Specifically, the different moisture evaporation models and assumptions regarding moisture diffusion and vapor re-condensation are investigated. The kinetic scheme and rate constants of devolatilization process are studied. A systematical approach is presented for the evaluation of the heat of pyrolysis which ensures elemental mass and energy balance. The effects of shrinkage on the particle and the change of porosity and biomass density are considered in the mathematical modeling. The formation of ash around the particle, the ash melting at high temperature and the consequences on the particle conversion are investigated. By means of a joint numerical study and advanced experimental measurements, a mechanism is proposed for the release of alkali metals from low chlorine biomass and the corresponding kinetic rate constants during devolatilization and char reaction stages are obtained. It is shown that the proposed model is able to predict the release behavior of the alkali metal during biomass gasification at various stages of the biomass conversion. The heterogeneous char reactions at the regime II reaction is affected by both the intra-particle mass transfer and the chemical kinetic rates. The evolution of iii

4 char porous structure can affect the conversion rate of the char. A model based on the classical capillary pores is developed taking into account different conversion rates for pores having different radii. This model is used to examine the contribution of each group of pores (micro, meso and macro-pores) in the conversion of biomass char. iv

5 Popular Science Description Most of us are familiar with using plant-derived material such as wood for producing heat for warming the cottages or cooking food when we are camping. Apart from that, plant-derived materials such as trees, forest residues, grasses, agricultural crops, agricultural residues and straw have great potential in producing energy in a larger scale. This group of material is referred to as biomass. Sweden has a great source of biomass as more than half of its area is covered with forest. Among the European countries, Sweden has the top share of using renewable energy sources which provide around 40% of the total energy of the country. The main reasons of recent attraction towards using renewable energy sources are the impact of fossil fuels in climate changes and concerns about permanency of the fossil fuels sources. In contrast to fossil fuels, biomass is renewable and can be CO 2 -natural if produced in a sustainable manner. The CO 2 produced from combustion of biomass will be absorbed by the plants that will replace biomass in a shorter term. In this sense, there is much less increase in CO 2 compared with the fossil fuels. As one of the green-house gases, CO 2 is believed to have a great impact on climate changes and global warming and therefore it raises concerns about its environmental impacts. Although the replacement of fossil fuel such as coal with renewable sources like biomass can help in producing green energy, there are challenges and problems that should be addressed. Emission of harmful pollutants, deposition of remaining ash in the biomass combustion systems, and corrosion of these systems need to be understood in order to improve the efficiency and decrease the cost of energy production. The aim of this thesis is to improve the understanding of the underlying physics of conversion of biomass to thermal energy in order to provide suggestions and solutions to improve the energy production system. The questions we are trying to answer are (a) how a biomass particle evolves during the conversion to thermal energy, (b) how various forms of pollutant are formed and released during biomass conversion, and (c) how to mathematically model the process with a balance between the computational time and accuracy of the results. Mathematical equations representing the physical phenomena are solved using numerical methods with the help of advanced computational programs. A biomass conversion model is developed which is comprehensive and can take into account the processes such as evaporation of moisture, decomposition of biomass, combustion of gases extracts from biomass decomposition, surface reactions of char and changes in the thermo-physical properties of the particle. The v

6 contributions of this thesis are (1) development of a comprehensive biomass conversion model with a reduced number of assumptions, (2) assessment of various modeling approaches and assumptions of different sub-models by comparison with several experimental measurements, such as the different moisture evaporation models and assumptions regarding moisture diffusion and vapor re-condensation, the kinetic scheme and rate constants of biomass decomposition process and the amount of heat required in this process, (3) development of chemical kinetic mechanism for the release of alkali metals, potassium and sodium, from biomass during combustion in high CO 2 environment, (4) development of a multi pore structure for gasification of biomass char, and (5) development of a semi-empirical model for fixed-bed combustion. vi

7 List of Publications I. H. Fatehi, X.S. Bai A Comprehensive Mathematical Model for Biomass Combustion, Combustion Science and Technology, 2014, In Press. DOI: / II. H. Fatehi, X.S. Bai Gasification and Combustion of Biomass: Physical Description and Mathematical Modeling, Book chapter, HANDBOOK OF CLEAN ENERGY SYSTEMS, John Wiley Press, In Press. III. H. Fatehi, Y. He, Z. Wang, Z.S. Li, X.S. Bai, M. Alden, K.F. Cen LIBS Measurements and Numerical Studies of Potassium Release during Biomass Gasification, Accepted for Oral Presentation in Symposium of the Combustion Institute, 2014, San Francisco, USA. IV. H. Fatehi, X.S. Bai, A simplified model for char combustion at fixed bed environment, Manuscript to be submitted. V. H. Fatehi, X.S. Bai, Effect of pore size distribution on the gasification of biomass char, Manuscript to be submitted. VI. H. Fatehi, X.S. Bai, Modeling of coarse fuel (biomass) combustion in fixed bed grate system, 9th European Conference on Industrial Furnaces and Boilers, Portugal, April Other publications H. Fatehi, X.S. Bai, Thermodynamic data and heat of pyrolysis of biomass, 4th international Conference of Applied Energy, 5-8 July 2012, China X.S. Bai, H. Fatehi, CFD Analysis of a high moisture biomass combustion processes in a 10MW grate fired furnace, 9th European Conference on Industrial Furnaces and Boilers, Portugal, April vii

8 Contents Abstract Popular Science Description List of Publications Nomenclature INTRODUCTION 1 Motivations 1 Biomass: from source to bio-energy 1 Thermochemical conversion 2 Challenges and approaches 3 Thesis contribution 4 Thesis outline 5 FUNDAMENTALS OF BIOMASS THERMOCHEMICAL CONVERSION 7 PHYSICAL PROCESSES AND MATHEMATICAL MODELING 11 Modeling approach 11 Basic assumptions 15 Kinetic rates 18 Fluxes 20 Fick s law and Darcy flow 20 Porosity 22 Specific surface area 25 Multi pore structure 28 Ash layer around the particle 30 Physical and thermochemical properties of biomass 31 Boundary and initial conditions 32 Derivation of Energy Equation 34 RESULTS AND DISCUSSION 37 Model parameters 37 Drying, pyrolysis and char combustion 37 viii iii v vii x

9 Heat of pyrolysis 43 Release of alkali metals from biomass 46 Experimental measurements 47 Numerical analysis 50 Alkali metal release sub-model 51 Formation of ash around the particle 56 Gasification of biomass 58 A semi-empirical model for fixed-bed combustion 60 CONCLUSION AND FUTURE PERSPECTIVES 61 SUMMARY OF THE PAPERS 65 ACKNOWLEDGMENTS 69 REFERENCES 71 ix

10 Nomenclature A Pre exponential factor 1/s or m/sk Bi x Biot number c Heat capacity J/kgK Da Damköhler number J/kgK D Effective diffusion coefficient m /s d Characteristic particle length m E Activation energy kj/mole F Flux of detached char fragments kg/m s G(r) Surface area of the particle at radius r m h Convective heat transfer coefficient W/mK h Mass transfer coefficient m/s h Enthalpy of formation J/kg HHV Higher heating value J/kg J Diffusion flux of species i kg/m s k Kinetic rate constant 1/s K Permeability m k Effective thermal conductivity W/mK LHV Lower heating value J/kg M Molecular weight kg/mole Nu Nusselt number N Total convective flux kg/m s p Pressure Pa q Pore growth rate m/s R Moisture evaporation rate kg/m s R Volatile release rate kg/m s R Char formation rate kg/m s R Char reaction rate kg/m s R Source term in energy equation J/m s r Radial position m r Pore radius m R Particle radius m Sh Sherwood number

11 S Char specific surface area m /m T Temperature K u Gas velocity m/s v ash Ash volume fraction m/s Y X Gaseous species mass fraction Conversion h Heat of reaction J/kg δ ε η θ Ash thickness Void fraction Effectiveness factor Shrinkage factor μ Viscosity kg/ms ρ Density kg/m θ τ Stoichiometric coefficient Tortuosity Υ Ash melting rate 1/s Φ Thiele module ψ Structural parameter in random pore model Subscripts bm Biomass c Char cond Conduction conv Convection fc Formation of char fc t Char produced by tar secondary reaction g Gas g t Gas produced by tar secondary reaction i Gas species/pore size/numerical cell int Intrinsic j Species k Reaction index Kn Knudsen diffusion coefficient m Moisture n Reaction order rad Radiation t Tar True Non-porous density v Volatile 0 Initial state xi

12 Abbreviations ASA CELL DAF HCL LIG RPM Active Site Area Cellulose Dry Ash Free Hemicellulose Lignin Random Pore Model xii

13 INTRODUCTION Motivations Increased demands for energy sustainability and concerns over global warming have led to a growing interest in sustainable and renewable sources of energy. The EU s Renewable Energy Road Map proposes creating a new legislative framework to enhance the promotion and use of renewable energy and sets a mandatory target of 20% energy consumption through renewable energy and a mandatory minimum target of 10% for biofuels by 2020 [1]. As a sustainable and renewable source of energy, biomass has a great potential in generating heat and power. In 2005, the breakdown of renewable energy produced in the EU was as follows: 66.1% from biomass, 22.2% from hydropower, 5.5% from wind power, 5.5% from geothermal energy and 0.7% from solar power (thermal and photovoltaic) [1]. Sweden is the leading country in using renewable energy among EU countries, having achieved over a 40% share of renewable energy in their total energy consumption [2]. The energy from biomass can be utilized either in direct heat production or via conversion to biogas, bio-oil and char as storable and portable fuels. The current breakdown of Sweden s use of bioenergy indicates that approximately 88% of bioenergy sources are used to produce heat through combustion in combined heat and power (CHP) plants and industry units, while 8% is used to produce electricity in CHP plants and 4% is used in the transport sector [2]. This indicates the importance of biomass combustion in bioenergy production. Although anaerobic digestion in wastewater treatment plants, co-digestion plants, farm plants and landfill plants are the major approaches in Swedish biogas production, there is a trend towards using gasification for the production of biogas, bio-oil and char [3]. Biomass: from source to bio-energy Biomass is referred to plant-derived materials such as trees, forest residues, grasses, agricultural crops, agricultural residues and straw. In a broader view, 1

14 municipal solid waste (MSW) and waste from food factories can also be included in the biomass category. The composition and thermo-physical characteristics of biomass may vary significantly depending on the source. To produce bioenergy from biomass, whether through direct heat production or conversion to biofuels, two main routes can be considered: biochemical and thermochemical conversion. In biochemical conversion, biomass is converted to biofuels by using microorganism and catalysts. In thermochemical conversion which is the focus of this study, conversion may take place through torrefaction, pyrolysis, gasification or combustion. Torrefaction results in high-energy-density solid fuels. Bio-oil is the main product of pyrolysis and gasification results in syngas production. In the following Section, the thermochemical conversion of biomass is described briefly. Thermochemical conversion Thermochemical conversion of biomass is a complex, multi-phase and multi-scale process. This process is a result of strong interaction between the heat and mass transport phenomena and chemical reactions. Figure 1 shows a schematic diagram of thermochemical conversion of biomass. Depending on the heating rate and the ambient gas composition the desired path from biomass to the product can be selected. Figure 1. Thermochemical conversion of biomass The sub-processes involved in thermochemical conversion of a biomass particle start with transient heat conduction to the particle which results in elevated temperatures. Depending on the presence of moisture, the evaporation of moisture starts soon after the particle temperature reaches the evaporation 2

15 temperature. Free water convection by capillary flow through the void spaces may be important in low-temperature drying [4, 5]. The dried portion of the particle, typically initiating from the surface, starts a thermal decomposition process known as devolatilization. If this process takes place in the absence of oxygen, it is referred to as pyrolysis. The pyrolysis products are lumped in three groups: gas, tar (volatile) and char [6]. Low temperature pyrolysis is in favor of char production whereas the tar reaches its maximum yield in intermediate temperature and gas production increases with increasing temperature [7]. As a result of devolatilization, a layer of char is formed which may further react with oxygen (oxidation process), carbon monoxide and water vapor (gasification process). The volatile passing through the char layer can undergo secondary reactions namely tar cracking and polymerization, at high temperatures and high residence times due to the catalytic effect of char. Pressure induced by formation of volatiles results in a convective flow towards the particle surface. Temperature gradient and pressure buildup induced by formation of volatiles may result in formation of cracks in the particle structure or, in some cases, primary fragmentation [8]. Particle size can undergo shrinkage/swelling due to loss of mass and structural reordering. The transport of mass into the porous structure of char can play an important role during the char gasification and oxidation. Char structure can drastically change by the heterogeneous reactions which lead to peripheral fragmentation. Challenges and approaches Different aspects of biomass conversion systems have been extensively studied [9, 10]. However, there are still challenges such as efficiency, pollutant emissions, ash deposition and corrosion that need to be understood and improved. Regarding biomass conversion systems, two approaches can be considered. The first approach is to substitute the coal with biomass either partially (co-firing) or completely in the coal-fired plants in a way that the existing infra-structure can be used with minimum modifications. The challenge arises as thermo-chemical properties of biomass are different from those of coal. Biomass has a lower energy density, a higher volatile content and much higher oxygen content. Moreover, biomass contains higher amounts of alkali metals such as potassium and sodium and deposits from co-firing with biomass are highly adherent to surfaces [11] which can cause serious damage to a combustion system. The second approach is to originally design and develop the conversion system for biomass. A thorough understanding of the processes occurring during the thermochemical conversion of biomass is essential for the design and optimization of biomass conversion systems in both approaches. Mathematical modeling of thermochemical conversion of biomass can provide a useful tool for addressing 3

16 these challenges where the full scale experimentation is expensive and difficult due to the harsh environment of such systems. Important parameters on a large-scale reactor level, such as distribution of temperature and species and overall conversion rate of biomass, are strongly dependent on the processes involved in the particle scale. A detailed spatially and temporally resolved numerical simulation of a single particle is helpful for understanding the interaction between transport, chemical kinetic and morphological transformation at particle scale. Moreover, a detailed numerical model of particle conversion can further be used to develop simplified models for CFD calculations of the combustion systems where due to the computational limitations, global modeling are often desired. Thesis contribution In this thesis, a numerical approach is chosen to study the thermochemical process of biomass. To this end, biomass conversion is divided into several sub-processes where mathematical description and numerical modeling of each sub-process is explained. Because of the complexity of the problem and the limited computational resources, it is inevitable to make certain assumptions and simplifications or to neglect certain physical effects. Therefore this thesis is aimed at assessment and improvement of these assumptions and simplifications in order to gain a better understanding of the underlying physics of biomass thermochemical conversion. The main contribution of this thesis is listed below. Development of a comprehensive biomass conversion model A comprehensive biomass model is developed and evaluated. The model takes into account the main features of biomass conversion. Various modeling approaches and assumptions of different sub-models were assessed by comparison with several experimental measurements. Specifically, different moisture evaporation models and assumptions regarding moisture diffusion and vapor recondensation, the kinetic scheme and rate constants of devolatilization process and the heat of pyrolysis were studied. In addition, the effect of shrinkage on the particle porosity and biomass density was considered in the mathematical modeling. This model can be used for predictive studies of different controlling parameters on biomass conversion and also can work as a benchmark reference for evaluating model simplifications. Mechanism and kinetic rate constants of release of alkali metals from biomass during oxy-fuel combustion 4

17 A joint numerical and experimental study is carried out, in which, a mechanism for release of alkali metals from low chlorine biomass was proposed and corresponding kinetic rate constants during devolatilization and char reaction stages were obtained. Development of a multi pore structure for gasification of biomass char The heterogeneous char reactions occur at the surface of porous structure of the char. For regime II where the intra-particle mass transfer and chemical kinetic rates have comparable time scales, the evolution of porous structure proved to be important. A model based on the classical capillary tubes was developed taking into account different conversion rates for pores having different radius. This model improves the understanding of the contribution of each group of pores (micro, meso and macro-pores) in the conversion of biomass char. Development of a semi-empirical model for fixed-bed combustion Tradeoff was made in favor of computational time by simplifying the processes in a fixed bed combustor and a semi-empirical model was developed. This model can provide acceptable boundary conditions for the CFD simulation at a very low computational cost and is most suitable for the studies where the focus is on the gas phase combustion of free board in a fixed-bed combustor rather than the processes inside the bed. Thesis outline The thesis is organized as a collection of papers that have been written during this project. The introductory overview of the topic is presented in the thesis followed by the Papers. An overview of the physical processes involved in the thermochemical conversion of biomass is presented in Chapter 2. Mathematical modeling of each process is discussed in Chapter 3. In Chapter 4, summary of the results and discussions is presented. The detailed results are given in the Papers appended at the end of the thesis. Recommendation for the future work and main conclusions of this work are presented in Chapter 5, followed by summary of the Papers as presented in Chapter 6. 5

18 6

19 FUNDAMENTALS OF BIOMASS THERMOCHEMICAL CONVERSION Various processes involved in thermochemical conversion of biomass are described briefly in this Chapter. The conversion process starts with evaporation of moisture. Moisture evaporation process involves simultaneous heat and mass transfer. Heat is conducted and radiated towards the particle center. The moisture starts to evaporate when the temperature reaches the evaporation temperature at the surface of the particle. By diffusion and capillary movement, the bound and free water move towards the dried surface of the particle. The vapor moves out of the particle by convection and diffusion. Some parts of the vapor may diffuse to the colder inner part of the particle where condensation may occur. It can be shown that for initial moisture contents below the free water continuity point (around 45% on dry bases) and high-temperature drying, the free water movement is negligible [4]. After the moisture leaves the particle, the dry portion of the particle starts to decompose due to external heating. This process is referred to as pyrolysis. Different approaches for modeling the kinetics of pyrolysis have been adapted in the literature. Di Blasi [6] provided a comprehensive review on the kinetics of pyrolysis of biomass. In general, four categories of reaction mechanisms can be identified for biomass pyrolysis: global mechanisms with single-step chemistry, e.g., [12], global mechanisms with three parallel reactions [13]; Functional Group- Depolymerization, Vaporization Crosslinking (FG-DVC) models [14, 15] and multi-component semi-detailed mechanisms [16, 17]. Figure 2 shows a schematic view of the different pyrolysis mechanisms. The pyrolysis products are usually lumped in three main groups: gas, tar and char. Among others, a review by Neves et al. [7] suggests that depending on the heating rate and reactor temperature, the yield of pyrolysis products may vary significantly. Low temperature pyrolysis can lead to more fixed carbon residue (char), while a high heating rate promotes the formation of gas and tar. The compositions of tar, char and gas are a complex function of various factors, including the elemental mass fractions of carbon, hydrogen and oxygen in the fuel, the particle size, the reactor temperature and the structure of the fuel. The primary 7

20 gas species produced from pyrolysis are CO 2, CO, H 2, CH 4, C 2 H 4 and H 2 O. Tar is composed of several high molecular weight molecules such as aliphatic hydrocarbons and aromatic and phenolic compounds, which at room temperature are in liquid form. Figure 2. Schematic view of different pyrolysis schemes; (a) single-step, (b) three parallel with secondary reactions; (c) FG-DVC, and (d) multi-component semi-detailed The heterogeneous reactions of char are the rate controlling processes in the overall conversion of biomass. The overall reaction rates of char depend on the following chemical and physical processes: the transport of reaction agents through the external boundary of particle or film diffusion, diffusion through the ash layer, the transport of reaction agents inside the pore structure of char and intrinsic kinetic rates. The intrinsic kinetics rates are referred to those rates that are not influenced by transport processes. Depending on which of the above processes is controlling the reaction rates, three different burning regimes or zones are distinguished for char combustion and gasification. The Zone I reaction is a regime where the chemical reactions are limiting the reaction rate. In the Zone III reaction regime, the mass loss rate is limited by the rates of oxygen and gasification agents diffusion to the outer surfaces of the particle. In the Zone II regime, the characteristic rates for pore 8

21 diffusion and chemical reactions are of similar magnitudes, and both effects are important in determining the mass loss rate. The changes in the char structure, the distribution of active sites and the catalytic effect of minerals have strong influence on the char reactivity, which introduces great challenges in the modeling of char conversion. To deal with this problem, the reaction rates of char are usually expressed considering two separate effects: the effect of temperature and partial pressure of the gasifying/oxidizing agents and the effect of the specific surface area of the particle, R = R (T, p )f(x ), where X is the char conversion and f(x ) takes into account the changes in the particle specific surface area as the conversion proceeds. In the expression above, it is assumed that the intrinsic kinetic rates are independent of conversion, which holds only for homogeneous distribution of active sites and minerals inside the char structure. However, in the empirical formulation of f(x ) these effects are included as it is obtained by fitting the char reactivity to a certain experimental data. In addition to the chemical and transport phenomena, changes in the structure of biomass during conversion are very important. During its thermochemical conversion, a biomass particle may undergo fracturing and fragmentation at any stages of conversion, i.e., pyrolysis and gasification. Fragmentation process can be divided into two main stages; primary fragmentation and peripheral fragmentation. Primary fragmentation is referred to breakage of particle during the pyrolysis process by forming fractures and smaller particles. Peripheral fragmentation is disintegration of the char particle from its surface. As fragmentation results in producing smaller particles, it can change the overall combustion behavior of the solid particles. Furthermore, the emission of aerosol and flying ash increases due to fragmentation. This indicates the importance of understanding the fragmentation process. The extent and the source of fragmentation depend on several factors such as the particle size and porosity, the biomass type, the volatile content and the furnace temperature. Several studies, mainly focused on coal particles, have been carried out to clarify this process. Several mechanisms have been proposed regarding the source of the particle fracturing and fragmentation such as: Fractures caused by pressure build-up during the pyrolysis stage Thermal stress induced fractures Compressive stress due to particle shrinkage Attrition due to particle collision or packed bed height pressure Loss of structural integrity due to heterogeneous reactions One or several of the mechanisms mentioned above can cause fragmentation. The possibility of fragmentation, the time and location of the fractures and the distribution of the smaller particles produced from the fragmentation were subject 9

22 of many numerical and experimental studies. Here, a brief overview of different numerical approaches will be presented. As coal was one of the major energy sources for a long time, various studies have been carried out on coal primary fragmentation in different combustion systems such as pulverized coal furnaces [18, 19] and fluidized bed furnaces [20, 21]. For the early stage or primary fragmentation, thermal stress and pressure induced fragmentation are introduced as the main mechanisms. The thermal stress due to temperature gradient inside the particle was reported as the main mechanism for fragmentation [22-24]. The volatile pressure built up was also considered as the fragmentation source [21, 23]. During coal particle pyrolysis, an internal pressure as high as 150 bar can be built up inside the particle. For biomass particles, on the other hand, in spite of the high volatile content, the higher permeability of wood structure prevents high volatile pressure built up as compared with coal. As discussed by Sreekanth et al. [8], it is not clear whether this low volatile pressure (2-3 bar) for biomass can cause primary fragmentation. Consequently, Sreekanth et al. [8] proposed that the high shrinkage of biomass particle, up to 50%, can be one of the sources of primary fragmentation in the later stages of pyrolysis. Peripheral fragmentation, in contrast to primary fragmentation, is mainly due to the loss of structural integrity during the char heterogeneous reactions, i.e., char oxidation and gasification. When the particle porosity reaches a critical value, namely the percolation threshold, the remaining structure may loss its macroscopic connectivity and the particle can break into hundreds or even thousands of smaller particles at the surface and the particle radius shrinks. This phenomenon is called percolative fragmentation. Depending on the regime in which the char particle is reacting, the particle porosity and density may change. For a particle burning in the chemical kinetic control regime, the particle porosity reduces throughout the entire particle and uniform percolation may disintegrate the whole particle structure. When a particle is in the diffusion control regime, i.e., high temperature oxidation of char, the peripheral percolation up to the penetration depth of oxygen may occur [25]. The peripheral fragmentation is studied, among others, by Feng and Bhatia [26], Miccio [27], Hurt and Davis [28], Golfier et al. [29]. 10

23 PHYSICAL PROCESSES AND MATHEMATICAL MODELING Modeling approach Different types of models have been developed to study thermochemical conversion of biomass with certain simplifying assumptions depending on the model purpose. When modeling the biomass and char conversion, there are different issues that need to be considered. In general, models concerning thermal degradation of biomass assume a uniform porosity while distribution of the pore structure is one of the main concerns in char conversion modeling. The level of complexity of the model depends on the specific case of interest. Non-dimensional analysis can provide useful information on the assumptions and simplifications that can be made for modeling of a given case. In some special cases the heat and mass transfer or the chemical kinetic rate is the controlling phenomenon. In these cases some simplifications can be made regarding the calculation of certain conversion parameters. To define nondimensional numbers, first the characteristic rates of reactions as well as the rates of internal and external heat transfer should be defined. Note that here the term reaction is used liberally to include drying, pyrolysis and char reactions. The characteristic pyrolysis reaction rate can be expressed as a first order Arrhenius rate. For char on the other hand, the reaction rate may be expressed as per unit of the reacting surface. Then a total reacting surface per unit of mass can be used to express the reaction rate in per second unit. The external heat transfer depends on the convection heat transfer coefficient (h ), the particle density (ρ ), specific heat (c ) and the characteristic particle size l = V /A. The internal heat transfer depends on the effective thermal conductivity (k, ) of the particle. The following equations describe the characteristic rates of pyrolysis (R ), external heat transfer (R, ) and internal heat transfer (R, ), R = A exp E, (1) RT 11

24 R, = h ρ c l, (2) R, = k, ρ c l. (3) The corresponding dimensionless numbers for these specific rates are listed in Table 1 [10]. Thermal and mass Biot number, Bi, quantifies the relation between the internal and external heat and mass transfer. Damköhler number, Da, for solid particle compares the rate of reactions and the internal or external heat or mass transfer to the particle. For instance, a thermal Biot number that is smaller than unity indicates that the external heat transfer rate is smaller than the internal heat transfer rate. Hence the temperature can be assumed uniform in the entire particle and the particle heating is controlled by the external heat transfer rate. Similarly, a pyrolysis Damköhler number that is larger than unity implies that the pyrolysis process is limited by the heat transfer rate rather than the kinetic rate. Based on this non-dimensional analysis, expressions to estimate the time of drying and pyrolysis for asymptotic cases can be extracted. Gómez-Barea and Leckner [10] provided these expressions for a variety of asymptotic situations a particle may encounter. Classifications of combustion regimes in a packed bed, based on these non-dimensional numbers is provided in [30]. Table 1. Definition of dimensionless numbers Thermal Bi External/Internal heat transfer rate Bi =,, = Mass Bi External/Internal mass transfer rate Bi =,, = Thermal Da Reaction/Heat transfer rate; Internal Da, =, External Da, =, Mass Da Reaction/Mass transfer rate; Internal Da, = External Da, = As discussed earlier using the non-dimensional numbers provided above, some simplifications can be made when modeling the biomass conversion under certain conditions. The simplest model developed for conversion of biomass [31, 32] and char [33] is the unreacted shrinking core model as presented in Figure 3. This model assumes that pyrolysis or char conversion takes place at an 12

25 infinitesimally thin surface propagating from the surface of the particle towards the center. When the reaction front passes through any given point in the particle, the conversion is completed and a char layer, in the case of pyrolysis, or an ash layer, in the case of char conversion, will remain. This means that across the reaction front there is a sudden change in the properties of the solid medium. In the case of pyrolysis the properties of the solid medium change from biomass properties to char properties and in the case of char conversion, from char properties to ash properties. For pyrolysis modeling, both finite rate chemistry in the form of Arrhenius and infinite rate chemistry in the form of constant temperature for pyrolysis have been used within the framework of this model [34]. For char conversion, the rate of diffusion of oxygen through the gas film at the particle boundary, the rate of diffusion of oxygen through the ash layer and the chemical kinetic rates can be taken into account when calculating the conversion rate [35]. This model is also modified to include the moisture evaporation process by assuming constant-temperature evaporation [32, 34]. An estimation of the thickness of the reaction front as k, T/q is given in [31] where T is the temperature difference across the reaction front and q is the heating rate. For the assumption of infinitely thin reaction front to be valid, the heat flux to the particle should be sufficiently high. This assumption is more suitable for reaction of char with oxygen at high temperatures or pyrolysis under high heating rates. In case of a slow reaction (kinetic control regime), the model assumption is not valid. For particle with low ash content or in the case of a fluidized bed where the particle dynamics does not allow for any ash layer to be formed around the particle, an unreacted shrinking particle model may be used. The assumptions are similar to the ones for the unreacted shrinking core model, except that the external radius of the particle will shrink according to the mass loss. These models assume a non-porous structure for the particle. The particle density remains constant and the particle or the core radius decreases in proportion to the particle conversion to the power of 1/3 [36, 37]. A uniform conversion model is another approach valid for Bi 1, where the internal field of gas species and temperature are uniform throughout the entire particle. In this model the particle size remains constant and the density decreases in proportion to the particle conversion. 13

26 Figure 3. Schematic view of unreacted shrinking core model for char reaction and pyrolysis For slow reactions or a high porosity char where neither the infinitesimal reaction front nor uniform-temperature field assumption is valid, progressive models are developed. Depending on the behavior of the ash, two types of progressive models have been proposed, shrinking particle progressive model and shrinking core progressive model [10, 35]. The solution of the differential equations for species gas profile inside the particle is required for this type of models. The difference between these models and models with the assumption of infinitesimal reaction front is the presence of a broader reaction zone that can grow inwards from the surface to the center of the particle. In this situation, both the particle size and density will change with the particle conversion. The powerlaw relations [38] may be used to provide the relation between changes in the particle density and the particle size in case of char reaction. In addition to the simplified models mentioned above, a spatially resolved model for particle conversion has been developed for pyrolysis [39-43], for the entire biomass conversion [44-47] or char conversion [48]. The main features of these spatially resolved models suggested in the literature are the same. These models consider transient partial differential equations for the mass and the energy equations and use the momentum equation in the form of the Darcy law. Although they differ in some details, the basic assumption for these studies is the assumption of continuum. In continuum models the macroscopic properties such as temperature and species concentrations are smooth functions of time and space. These types of 14

27 models denote the classical engineering approach which employs the concept of volume averaging to describe a system which is characterized by several length scales. The volume in which the macroscopic properties are represented by average values should be small enough compared with the characteristic length scale of the system such as scalar gradients and large enough compared with the length scale of the porous structure, e.g., pore radius, so that the transport equations would be valid in this volume. The concept of volume averaging introduces some restrictions on employing the continuum models. For instance, in a system with large, disconnected voids or drastic changes in the connectivity of the media, the continuum models are not well suited. Another approach which is free from the limitations and assumptions of the continuum approach is to use discrete models. Discrete models are based on the percolation theory which deals with the concept of a spread of a hypothetical fluid through a random medium [49]. In these models, the particle is represented by two or three-dimensional random networks of large number of interconnected pore bodies and pore throats. The transport and reaction phenomena in this system can then be studied in a more realistic way by employing the ideas of statistical physics of disordered media. The main disadvantages of these models are the required previous knowledge of the microstructure of the porous system and their high computational costs. These models will not be elaborated here. Readers are referred to [49] for a comprehensive review on these models. Although due to convenience and the lack of practical alternatives, the continuum approach has been used widely in coal/biomass combustion modeling, one should bear in mind the limitations of this approach, i.e., the porous system should be connected while the pores are small enough to allow for a meaningful volume averaging. The model developed in this study lies within the category of the continuum approach and the governing equations in the continuum model framework will be presented. The particle is represented by a one-dimensional system consisting of a solid structure with void spaces filled with liquid and gas. Basic assumptions As discussed in previous Chapter, thermochemical conversion of biomass is a complex process. At the moment, a comprehensive mathematical description of the process still requires substantial simplifications for various processes involved. The model presented in this study is based on the following assumptions: The biomass particle consists of multiple components, including solid, gas and liquid All these components are in thermal equilibrium 15

28 The particle geometry can be spherical, cylindrical or slab (parallelepiped plate) and it is represented in a one dimensional frame work Gases inside the porous structure of the particle obey the ideal gas law The momentum transport in the particle is governed by the Darcy law Homogeneous gas reactions inside the particle are neglected For a system of solid, liquid and gas, under the model assumptions stated above, the conservation of mass for solid and liquid phases are governed by the following equations: ρ t = R, (4) ρ t = 1 r r D ρ + R, (5) where ρ is the density, D is the diffusion coefficient, and R is the mass release rate source term. This system of equations is valid for an arbitrarily shaped particle within a one dimensional framework. Therefore, r is the distance from the particle center and n is the Lamé coefficient accounting for a non-cartesian coordinates system: n = 2 for spherical, n = 1 for cylindrical and n = 0 for large slab plate shaped particles. The right-hand side terms are the source terms, which are described for each process in the following Section. The conservation of gaseous species and mass and energy are as follows: t ερ Y + 1 r (r N Y ) = 1 r. (r J ) + θ, R, (6) t ερ + 1 r (r N ) = θ, R, (7) ρ c + ρ c + ρ c + ερ Y c T t + N Y c T = 1 r r k T + R, (8) where N is the total convective flux, J is the diffusion flux of species i and ε is the void fraction. Y and θ, are the mass frction of species i and the stoichiometric coefficient of species i in reaction k, repectively. The energy equation as presented in Eq. (8), is derived from the enthalpy form of the energy 16

29 conservation and R is the heat source term appearing when writing the equation in the temperature form. The source term include the heat of reaction of different processes and the changes in sensible enthalpy as presented in the following equation, R = R h h + R h h + R h h + R h h + R h h + R c c dt + R h Y h + R h Y h + R c Y c dt + R c c dt + R c c dt + R c Y c dt + R c c dt + R c Y c dt, (9) where h is the enthalpy of formation and h is the heat of reaction. The derivation of energy equation is presented at the end of this Chapter. The enthalpy of formation of biomass, char and tar are related to their lower heating values 17

30 (LHV) upon their complete reaction with oxygen. If the C x H y O z represents either biomass, char or tar, the enthalpy of formation is as follows, C H O + x + y 4 z 2 O xco + y 2 H O, (10) h = LHV + 1 θ M M h,, (11) where M is the molecular weight and θ is stoichiometric coefficient, e.g., x for CO. If the lower heating values are not known, the correlation of Channiwala and Parikh [50] for biomass and char and correlation of Mason and Gandhi [51] for tar provides the heating values with good approximation, HHV / = 34.9Y Y Y 10.3Y 1.5Y 2.1Y MJ kg, (12) HHV = 34.1Y, Y, 12Y, MJ kg, (13) where HHV is the higher heating value. The reaction rates are presented in the following Section. Kinetic rates There are three stages of mass loss during the thermochemical conversion of biomass describing the reaction source terms in the governing equations, drying, pyrolysis and char reaction. The liquid mass conservation is governed by the Eq. (5) where the moisture diffusivity (D, ) inside the particle porous structure depends on the moisture mass fraction and temperature [5], D, = ex p T + 9.8Y. (14) The evaporation source term is based on the equilibrium model that is shown to perform well at different conditions [43]. The model can be mathematically presented as, R = h, S ρ, ρ,, (15) 18

31 where ρ, is the saturation vapor density calculated based on the Clausius Clapeyron relation for the saturation pressure and ρ, is the partial density of water vapor in the pores of the particle. S is the specific surface area of the particle available for moisture and vapor diffusion. h, is the mass transfer coefficient of vapor in the pores of the particle, which can be determined using the following equation [45], h, = 3.66D, d, (16) where d is avarage pore diameter. The dry solid particle (or biomass) mass degradation is governed by Eq. (4), where the source term is calculated according to a pyrolysis mechanism. In general, the pyrolysis reaction is presented as follows, ρ t = ρ k = R, (17) where k is the rate of reaction k in Arheneous form. If any intermediate solid species are considered by the mechanism, the mass conservation yields, ρ t = ρ k ρ k Two secondary reactions are assumed, converting tar to char and gas,. (18) ρ t = ρ k + k = R + R. (19) Three heterogeneous reactions are considered for char; C + θ O 21 θ CO + 2θ 1CO, (R ) C + CO 2CO, (R ) C + H O CO + H, (R ) where θ is the stoichiometric coefficient which depends on the reaction temperature and is evaluated by empirical correlations. The first reaction is the oxidation of char, whereas the second and third are the main reactions in char 19

32 gasification. The heterogeneous char reaction rates are presented in the following form, R = M θ, S k, C, (20) where k, is the intrinsic rate of conversion per unit area of the carbon in reaction k (R C1 R C3 ) and i = O, CO, H O. M is the carbon molecular weight, S is the specific internal surface area of the char pores structure, θ, is the stoichiometric coefficient of carbon in reaction k, C is concentration of species i and n is the reaction order. Fluxes The diffusion and convection fluxes play an important role during the heterogeneous char reactions. During devolatilization and drying, these fluxes have a direct impact on the pressure built up inside the particle which may or may not cause fractures and fragmentation. There are different approaches for modeling the transport inside a porous media [52]. In general, two main approaches used extensively in biomass conversion models are extension of the Fick s law and the Dusty Gas Model (DGM). In both of these models the Knudsen diffusion, molecular diffusion and effect of pressure gradient can be taken into account. DGM is based on the Stefan- Maxwell equation for diffusion and the convection due to pressure gradient. It is not possible to derive an explicit equation from DGM to compute fluxes and therefore a numerical solution is required to solve the nonlinear system of equations. On the other hand, the Fick s law model has the advantage that explicit analytical expressions can be derived for fluxes. In addition, the deviation between the two models in calculating the fluxes in porous structure have been reported to be negligible [53]. Fick s law and Darcy flow When there is a pressure gradient contributing to the total flux of species, a Darcy law equation can be added to the classical Fick s law resulting in an extension of the Fick s law model, 20

33 Y m = ρ D, r + Kρ Y P μ r, (21) where m is the total mass flux of species i in gas phase, K is the permeability, P is the pressure and μ is the gas viscosity. Permeability is a measure of mass transfer resistance inside the pores structure of biomass [54]. The diffusion and convection fluxes in Eqs. (6)-(8) can be expressed as follows, J = ερ D, Y r, (22) Using Eq. (23) the gas velocity can be obtained, N = Kερ P μ r. (23) u = K μ P r. (24) In Eq. (11), D, is the effective diffusion coefficient. D, takes into account the molecule-wall interaction in small pores (Knudsen diffusion) and molecular diffusion as well as the effect of porosity and tortuous structure of capillaries based on the following equations, D, = εd τ, 1 = 1 + 1, D D D, (25) where D is the Knudsen diffusion coefficient, D, is a the diffusion coefficient of species i into a mixture of other species m, and τ is tortuosity, which can be estimated by τ = 1/ε [55]. The Knudsen diffusion coefficient is dependent on the pore radius and is calculated based on the following expression, D = 2r 3 8RT πm, (26) where M is the molecular weight of species i. Note that Eq. (26) is independent of mixture composition as it depends on the molecular weight of each species 21

34 individually. D, depends on the gas composition and temperature which can be calculated based on the binary diffusion coefficient of each species in the mixture, D, = 1 Y x D (27) where N is the total number of gas species. The temperature effect is included in the binary diffusion coefficients, D. For simplification, the mixture diffusion coefficient can be expressed as a function of temperature which provides sufficient accuracy [56], D, = D = D, T. P T P (28) where D, is the diffusion coefficient at reference temperature K which is equal to 2 10 m /s. Figure 4 shows the diffusion coefficients for individual species in the mixture together with the mixture diffusion coefficient as a function of temperature obtained based on Eq. (28). Except for highly diffusive hydrogen, the correlation for the mixture diffusion coefficient is close to the diffusion coefficient of each individual species in the mixture. Figure 4 Mixture diffusion coefficient as a function of temperature Porosity Biomass conversion can be studied as a process in a porous medium in which the physicochemical interaction of a fluid with a porous structure leads to continuous evolution of the solid structure. This process involves several length scales which 22

35 depending on the shape, size and morphology of the structure can play a fundamental role in transport and reactions in the system. The simplest quantity that can characterize the state of a porous system is the volume fraction. The volume fraction can be defined for as many phases as exist in the system. For phase η, gas, liquid or solid, one can write, ε = V V, (29) where V is the total volume of the system. The total porosity in the framework of this study is the sum of the gas and liquid porosities, ε = ε + ε. (30) It worth noting that the moisture content of the particle is treated as liquid when existing in excess to the free water continuity point (around 45% on dry basis). The time evolution of porosity is represented by the following equation, ε t = 1 θ ρ, R (31) where ρ is the non-porous solid density. To obtain Eq. (31) it is assumed that the true solid density is constant. During drying and pyrolysis stages, as the remaining solid has a different structure from the original one, this assumption is not valid. The change in the solid structure may cause particle shrinkage. The effect of the particle shrinkage should be included in Eq. (31). As changes in the porosity and particle shrinkage are linked together, an empirical correlation is used for the particle shrinkage. Figure 5. The relation between mass loss and changes in volume during particle shrinkage Two different types of shrinkage are distinguished in this study: volumetric shrinkage and peripheral shrinkage. While the volumetric shrinkage is caused due 23

36 to the structural reordering of the solid, the peripheral shrinkage is a result of the structural disconnection from the outer radius of the particle. For the volumetric shrinkage, a linear change in the particle volume is assumed with respect to the mass loss during each process. The shrinkage factor θ relates the volume change to the release of the biomass components namely moisture, volatiles and char. The slopes of these linear functions are empirically determined, θ = 1 + (1 θ ) ρ ρ, 1 + θ (1 θ ) ρ ρ, 1 + θ θ (1 θ )( ρ ρ, 1), (32) where θ is defined as the ratio between the current volume of the particle to its initial volume. θ, θ and θ are empirical parameters with values between a minimum value and 1, representing the extent of shrinkage due to moisture evaporation, volatiles release and char reactions, respectively. A minimum value represents a constant density and a value of one represents a constant volume during thermochemical conversion. The minimum value for θ is the mass fraction of moisture in as received basis of the particle and minimum value for the volatile release is the mass fraction of final char. The ash mass fraction is the minimum value for θ. The subscript 0 denotes the initial state of the particle. The effect of shrinkage on the porosity is taken into account by adding a pseudo-convective term to Eq. (31), where the corresponding velocity for this convective term is the rate of particle volume change, ε t = ε θ θ t + 1 θ ρ, R. (33) Peripheral fragmentation occurs when the porosity of the particle exceeds a critical value, ε. The particle radius (R ) then decreases according to the following equation, = when ε > ε. (34) The flux of detached char fragments (F ) that leaves the particle is given by [57], 24

37 F = ρ. (35) where ρ is the apparent char density at the onset of fragmentation. Specific surface area Although porosity provides some information regarding the state of the system, it cannot provide any knowledge on the morphology of the system such as the amount of surface of the porous structure accessible by other phases or the distribution of the pores radius in the structure. For situations that transport or reaction phenomena take place at the surface of the system, such as in char conversion, the interfacial surface density is an important parameter which can be defined as [49]: S = Area of contact between phase η and σ. (36) Volume of the system V This variable is referred to as the specific surface area and as it only appears in the equations during the char conversion, the following discussion only concerns the char particle. There are different theories developed for representing a porous medium. Capillary models are by far the most commonly used models in this field [49]. One of the well-known theories in capillary models is the random pore model (RPM) developed by Bhatia and Perlmutter [58, 59] and Gavalas [60]. The RPM assumes cylindrical capillaries randomly located and oriented with any distribution of pore radii. The surface area of the structure may increase due to the surface reactions or decrease due to different capillaries overlapping. The specific surface area is related to the conversion degree of carbon structure based on the following equation, S = S, (1 X)1 ψ ln(1 X). (37) where S, is the initial specific surface area, X is the carbon conversion and ψ is the structural parameter depending on the carbon structure which can be determined using the original pore size distribution. However, it is usually estimated by fitting to measured experimental data. For coal, RPM predicts a maximum in the specific surface area as the conversion proceeds. The specific surface area first increases as the pore radius is enlarged due to the ongoing reactions and at a certain point, starts to decrease as the pores merge and coalesce. For biomass, on the other hand, some experimental 25

38 evidence shows a constant or a continuously increasing trend for the specific area during the conversion. Dasappa et al. [61] measured the specific surface area of different biomass chars at various conversion degrees, up to 80% conversion, using the BET technique. The specific surface area varied weakly as the conversion proceeded. In another experiment, Lussier et al. [62] reported a constant adsorbed hydrogen concentration between 0.5% and 40% conversion of Saran char. A continuous increase in the reactive surface area was reported in [63] for the beech wood char. The reason for these different behaviors can be the different origin of coal and the biomass char. Coal is formed under pressure and heat which cause the structural order of the original source to be destroyed during the process of formation. Biomass char usually has two orders of magnitude greater specific surface areas than coal. Hence wood chars present a remarkably regular structure similar to a beehive network, whereas the coal chars are characterized by highly random structures due to their origin [61]. Because of these differences, the theories originally developed for evolution of a pore structure of coal during the conversion may fail for biomass [64]. Because of the insufficient experimental data and the contradictory results, it is hard to draw a conclusion on the evolution of the surface area. In addition, other complexities such as formation of cracks and presence of ash bring more difficulties to the problem. Some simple models have been used in the literature, one of which assumes a constant specific surface area during the conversion [44, 45]. This model produces good results for high-temperature oxidation, where the boundary diffusion plays a more important role and the evolution of the internal surface area is of less importance. Mermoud et al. [64] developed a simple model based on the assumption that the surface area increases in a way that the conversion rate of a gasifying char remains constant. This model showed good agreement with the few experimental data that were obtained for up to 60% conversion, c.f. Figure 6. In another study, Dasappa et al. [61] argued that assuming a monotonous increasing pore size without any merging and coalescing is more relevant for the ordered structure of wood than using the RPM that was originally developed for coal. This assumption implies that the specific surface area increases continuously. This model is based on the classical capillary tubes model which relates the specific surface area to the porosity and pore radius r by the following equation, S = 2ε r. (38) By some mathematical manipulations and the assumption that the pores grow due to the radial expansion, one can relate the specific surface area to the carbon conversion, 26

39 S = S, 1 X (1 1 ε ). (39) In Figure 6 a comparison between the RPM, capillary model and constant reaction rate model is presented together with some experimental data. As discussed previously, there is not enough data to draw a firm conclusion about which model is the most suitable for biomass conversion. However, the RPM can be ruled out, as none of the experiments suggest a remarkable decrease in the specific surface area as the conversion proceeds. It should be noted that not the entire available surface area participates in the reactions. The concentration of active sites on the char surface is an important parameter in char conversion kinetics and it is expected that the reactions take place only on the surfaces where the oxygen-complexes are formed [35]. Despite many attempts at measuring the active site area (ASA), e.g., [65], there are difficulties in correlating the ASA with the carbon gasification reactivity [66]. As a result, the reason for increasing specific surface area versus conversion may be due to the increase in concentration of the active sites. The model based on the classical capillary tubes shows acceptable behavior compared with the experimental data and therefore it is used in this study. Based on this model, the volume fraction and pore radius increase with the rate of char conversion. Based on Eq. (38), the rate of increase in surface area yields, S t = R /r. (40) Figure 6. Evolution of specific surface area of some biomass driven char at different conversion degree 27

40 Multi pore structure There are three distinct groups of pores in the particle structure; Macro-pores (> 50 nm diameter) Meso-pores (2-50 nm diameter) Micro-pores (< 2 nm diameter) The reactions are taking place mainly in macro and meso pores. Although the contribution of micro-pores in the void fraction is high, their contribution to the conversion is limited due to the limited accessibility to the reactants. Therefore, the entire initial surface may not contribute to the reactions. On the other hand, the accessibility of the micro-pores may increase as the conversion proceeds due to the increase in the pore radius. Char can undergo different reactions with CO 2, H 2 O and O 2. It was observed that different reactants behave differently in micro-pores. Hurt et al. [67] reported that the char reaction with CO 2 primarily takes place outside the micro-porous network and according to Ballal and Zygourakis [68] micro-pores are not accessible to O 2. On the other hand, H 2 O can penetrate in small pores that are not accessible to oxygen [69]. The above discussion suggests that during the char conversion different pores may evolve with different rates and their accessibility to certain reactants may change during the conversion. To take into account the various pore sizes, a pore size distribution is assumed for the char particle. The pore sizes are divided into discrete bins of initial porosity ε (r ) corresponding to the pore radius, r, distribution. Similar to the work of Singer and Ghoniem [48, 70], a variable q is assumed which indicates the radial expansion of pore i at each time caused by the char conversion. By assuming that the reactions take place on the pore surface, the solid-phase conversion can be rewritten in terms of time variation of q, 28 dq dt = 1 ν ρ, R,,, (41) where R, is the intrinsic reaction rate of the char reactions. The pore radius at each time is given by, r = r, + q. (42) To identify the degree of participation of each pore in the reactions, the concept of pore scale Thiele module for individual capillary was developed by Gavalas [60]. The effectiveness factor based on the Thiele module is employed here. The effectiveness factor can verify whether or not a pore with a radius of r is in the kinetic control regime. For a first order reaction with an intrinsic rate

41 constant of k, and an effective diffusion of D,, the Thiele module Φ is given by, Φ, = l 2 2νk, r D,, (43) where l is the length of the pore i. Note that the Thiele module is different for each reaction and each pore radius. The effectiveness factor is given by, η, = tanh (Φ,) Φ,. (44) The effectiveness factor quantifies the participation of each reaction for each pore size, which can be used to modify the pore growth rate, dq dt = 1 η ρ, ν, R,,. (45) Equation (31) can now be employed for each porosity bin ε, dε dt = S ρ η, ν, R,,, (46) where S is the specific surface area corresponding to ε and r. The effective diffusion for each pore i can be obtained based on the combined Knudsen, Eq. (26), and bulk diffusion, Eq. (27). Contribution of each pore size in the total effective diffusion coefficient then should be considered. In a slightly different context, Wakao and Smith [71, 72] presented a relation for the diffusion coefficient in a porous catalyst. By assuming that the micro-pores are distributed in the solid structure between the macro and meso pores, the contribution of each group of pores in diffusion can be calculated based on the following equation [73], D = ε 1 D, + 2ε ε 1 D + 1, (47), D, where effective diffusion coefficient for each pore radius is, D, = (48) D, D, 29

42 Ash layer around the particle The ash liberated from the particle during the char conversion can form a layer around the particle. This layer affects the char conversion by restricting the oxygen diffusion to the particle. Moreover, the particle radius changes due to the ash layer which has an impact on the heat transfer rate to or from the particle. Two important parameter of ash should be considered, the ash thickness and ash porosity, the first one affects the heat balance and the second one affects the species concentration. It is assumed that the ash evaporation or reaction has no effect on the temperature of the particle. In other words, in the particle model the ash is assumed inert. Moreover, it is assumed that the ash is distributed uniformly inside the particle and is liberated as the char reaction proceeds. If ρ denotes the ash density inside the particle structure and ρ, represents the density of liberated ash in the particle surface, dρ, dt = ρ Y a,0 R, ρ, R,, (49) where Y, is the initial ash mass fraction inside the particle and R, is the rate of ash evaporation. The thickness of ash is defined as the char (core) radius minus particle radius and depends on the ash porosity. By assuming a uniform distribution of ash inside the particle, the ash thickness can be obtained according to the following equation [57], R v = R v + 3 r (1 ε )dr, (50) where v ash is the ash volume fraction and it depends on the initial ash mass fraction, initial particle apparent density (ρ, ) and true density of ash (ρ, ), 30 v = Y a,0ρ, ρ, θ,, (51) where θ c,0 is the extent of particle shrinkage at the start of char reactions, as it is assumed that ash is liberated by char reactions. In Eq. (50), R 0 is the radius of particle at the instant of char reaction initiation. Hurt et al. [74] assumed that before the ash layer reaches a critical thickness (100 μm), it is in the form of detached grains. They proposed a relation for the ash porosity of these detached grains. When the thickness of ash layer is above the

43 critical value, the porosity remains constant. The correlation is developed base on the observation for pulverized coal particles. The correlation is suitable for a particle in the pulverized regime. Here it is assumed that the ash is formed with the same porosity as char. The ash porosity, on the other hand, can change with time due to ash melting. With the assumption that the ash porosity is uniform in space, Eq. (50) yields, dr dt = R (1 ε ) v dr dt + 3R (1 ε ) R R dε dt, (52) where dε ash /dt represents the variations in the ash porosity due to the ash melting. By assuming that the ash porosity is linearly decreasing in time when the ash is exposed to a temperature higher than melting temperature, the following equation yields the ash porosity [48], dt ε = ε, Υ, (53) 1 + exp 10(T T ) where Υ is the ash melting rate and t is the time that the ash is exposed to temperature higher than the melting temperature (T melt ). The effective diffusion coefficient through the ash layer is related to the ash porosity according to the following equation [74], D, = D, ε.. (54) The effect of Knudsen diffusion through the ash layer is neglected. Physical and thermochemical properties of biomass Heat conductivity is a function of temperature, density, porosity, heat flow direction with respect to the grain and the moisture content [12]. The effective thermal conduction coefficient (k ) is a sum of the radiative (k ) and conductive (k ) heat transfer coefficients. The conductive heat transfer coefficient is a mass weighted sum of moisture (k ), char (k ), biomass (k ) and the gas (k ) coefficients, 31

44 k = k + k, k = εk + (1 ε)[(ρ k + ρ k + ρ k )/(ρ +ρ + ρ )], (55) k = 2εσT r. ε where σ is the Stefan-Boltzmann constant and ε is the emissivity. Permeability is a measure of the mass transfer resistance inside the pore structures of the biomass. The following expression is often used to predict the permeability [54], K = r ε 45(1 ε). (56) The specific heat of biomass, char [75] and tar [76] are calculated based on the following equations, c = exp( 0.004T), (57) c = T T, (58) 3.33 c = T T 10 T. (59) Boundary and initial conditions The boundary layer around the particle can affect the heat transfer because of the combustion of volatiles and the formation of flames; this can have a drastic effect during the pyrolysis if oxygen is present. However, due to the high importance of the radiation heat transfer to the particle surface, the boundary layer effect is often negligible. In the set of equations outlined in the previous sections, the boundary layer effect is considered through the heat and mass transfer coefficients by using a Stefan correction. At the particle s outer surface, the heat and mass transfer fluxes, along with ambient pressure, determine the boundary condition: 32

45 εd C r = h C, C, k T r = h (T T) + σεt T, (60) where T and T are the gas temperatures at the particle boundary and the wall temperature, respectively. The mass transfer coefficient, h, and the heat transfer coefficient, h, can be calculated according to the following equations: h = m c exp m, h = Nuk, c d 1 h (61) h = m /ρ exp m, h = ShD, d 1 ρ h (62) where d is a characteristic particle length, e.g., the diameter of a spherical particle, c is the specific heat of the gas, ρ is the gas density and k and D are the gas conductivity and diffusivity, respectively. The gas outflow from a particle can influence the convective heat transfer and mass transfer coefficients. The Stefan correlation [52] accounts for the outflow of volatiles and steam (m ). Nusselt (Nu) and the Sherwood (Sh) number, which are functions of Prandtl (Pr) and Reynolds (Re) numbers were obtained based on the Ranz-Marshall correlation [77], Nu = a + cre Pr, (63) where a, c, m and n are constants determined by the fluid and particle geometry. At the center of the particle, the gradient of pressure, temperature and all species are zero due to particle symmetry: p T = 0, r r = 0, u = 0, C r = 0. (64) For initial condition, one may assume that the particle is in equilibrium with surrounding gas at ambient pressure and temperature: p = p, T = T, u = 0, C = C,. (65) 33

46 Derivation of Energy Equation The energy equation in the form of the total enthalpy of the system consisting gas, liquid and solid can be written as follows, () ρ H + ρ H + ρ H + ερ Y H + ερ H + ερ H + G(r)u ερ Y H + ερ H + ερ H = G(r)k () + () G(r) ερ D, H. where H is the specific enthalpy and H = h, + c, dt. By simple mathematics, one can write, = h, + c, dt = c,dt = c,. = c, The energy equation can further be expanded. The first term becomes, ρ H + ρ H + ρ H + ερ Y H + ερ H + ερ H =, ρ + ρ + ρ + ερ + ερ + ερ + H + H + H + Y H + H + H The convection term in the energy equation can be expanded to () G(r)u ερ Y H + ερ H + ερ H =. ερ u + ερ + ερ + () Y H () () + H + H () Two terms in the following equation require further expansion ερ = ερ Y + ερ H 34

47 ερ u = ερ u Y + ερ u H. By substituting = c, and = c, and moving the term in the bracket to the right side of the energy conservation equation, the equation becomes, ρ c + ρ c + ρ c + ερ Y c + ερ c + ερ c ερ u Y c + ερ c + ερ c = () G(r)k + ερ u H H () H + H () + G(r) ερ D, H H Y H () + H + Y H () + H + H + H + () Here, the terms that appear after the conduction term are examined. By mathematical manipulation and substitution by the conservation of species, one can write,. () G(r) ερ D, H ερ H ερ u H θ, R Y θ, R + () G(r) ερ H = D,. H θ, R Y θ, R is equal to zero if no homogeneous or heterogeneous reaction occurs in the system, e.g., if the case only considers pyrolysis. Otherwise, this term represents the enthalpy changes by the changes in the gas composition due to heterogeneous reactions, e.g., k char reactions. The second term in the equation is due to the Dufour effect, but its contribution to the overall energy balance is negligible. The last two terms in the brackets can be simplified further by substituting the mass conservation equations into the source term on the right side of the energy equation. Hence, H + H + H + Y H Y () () H + H + H () H () + Y H () + () () + H () + H + H + () = H (). + + H + + H = H R R R + H R + R + H ( R ) + Y H R + R + H R R R + H (R ) + H Y θ, R H θ, R. The energy equation now can be written in the following form θ, R + 35

48 ρ c + ρ c + ρ c + ερ Y c + ερ c + ερ c + ερ u Y c + ερ c + ερ c = G(r)k () + R H Y H + R H H + R H H + R H Y H + R H H + R H H + R θ, H H 36

49 RESULTS AND DISCUSSION Model parameters The model parameters, including the physical and thermodynamic characteristics of the biomass, are taken from the reported values in each experimental set-up. If a value is not reported, then standard values are used. Particle density, thermal conductivity, diffusivity, emissivity, void fraction, pore size, permeability and specific heats for all components are required for each simulation. Apart from these parameters, the kinetic constants for pyrolysis and char reactions are also assigned based on existing literature. The apparent density of the biomass particle is usually reported and the values range between 330 to 1000 kgm. The thermal conductivity is calculated based on Eq. (55). The biomass thermal conductivity is set to 0.2 Wm K ; as for the moisture inside the biomass, k is set to 0.65 Wm K. The thermal conductivity of the char is a function of temperature and is calculated based on the expression from [78] as presented in Eq. (66): k = (1 ε )( T), (66) where ε is the char porosity. The ash conductivity is 1.2 Wm K and the gas conductivity is Wm K. The particle emissivity is given as a function of temperature. For temperatures below 450 K, the emissivity is set to 0.7, and for temperature above 550 K, the emissivity is set to For temperatures in between these values, a linear interpolation between the two values is used [79]. In this Chapter, the main results of the model discussed in the previous Chapter are presented followed by a summary of the papers and additional discussions. Drying, pyrolysis and char combustion In Paper I, a comprehensive model for the combustion of biomass was presented. The set of equations for the thermochemical conversion of a biomass particle described in the previous Chapter was used to study physical and chemical 37

50 processes that occur during biomass conversion. Various sub-processes involved in thermochemical conversion of biomass are described and mathematical modeling of these processes is presented. To validate the model, the experimental data of Lu et al. [45] was examined where temperature and mass loss were examined for a cylindrical particle with a diameter 9.5 mm and an aspect ratio (length to diameter) of 4 in a reactor filled with nitrogen with a gas temperature of 1050 K and an average wall temperature of 1276 K, c.f., Figures 1 and 2 in the Paper. To investigate the evaporation behavior, the case with a high moisture particle with 40% moisture based on the dry fuel mass was chosen. The reaction fronts of different chemical processes involved in the particle conversion for the case discussed in Paper I are illustrated in Figure 7. The changes in the particle radius due to shrinkage are noticeable in this figure. Drying starts immediately after that the particle is exposed to the hot furnace environment. As the drying front propagates inwards to the center of the particle, dry fuel is left behind. The pyrolysis front separates the dry fuel and the newly formed char layer. It takes 40 s for the moisture to evaporate and the devolatilization is completed after 70 s. The char oxidation front is very thin and it does not spread inside the particle as the process is controlled by mass transfer. Figure 7. Propagation of reaction front during biomass conversion Release of mass from the particle influences both the particle porosity and radius. The changes in porosity for a particle during the evaporation and devolatilization stages are presented in Figure 8 for four different time instants. The effect of particle shrinkage can be observed by comparing the solid lines and dashed lines in this figure. If the particle shrinkage was not included in the model, the final char porosity is calculated to be around 95%, whereas by including the particle shrinkage, this value is reduced to 82%, which is a reasonable value for this char particle as the final char yield is around 10% of its initial mass. The 38

51 distinct effects of drying and devolatilization can be distinguished in this figure as there is a time lag between the reaction fronts passing by each point along the particle radius. Biomass particles usually have high moisture content. In Paper I, a sensitivity study on different drying models is performed to evaluate the models for prediction of the drying behavior of biomass. Three different evaporation models were considered in this Paper: the heat flux model, the equilibrium model and the chemical kinetic model. The heat flux model is based on the assumption that evaporation takes place in an infinitely thin region and at normal boiling point of water. The equilibrium model is based on the assumption of equilibrium between water vapor and liquid water. In the chemical kinetic model, evaporation rate is expressed using the temperature dependent Arrhenius expression. These three models were compared with experimental measurements, c.f. Figure 5 in Paper I. The results of a detailed model developed by Di Blasi [5] was also included in the comparison. The Di Blasi model takes into account the water vapor convection and diffusion, capillary water convection due to pressure gradient and the bound water diffusion in the pores of the particle, c.f. Figures 6 and 7 in Paper I. The equilibrium model was found to be the best option for evaporation modeling both in comparison with the experimental data and the detailed model. In addition to reproducing the experimental measurements, the equilibrium model is numerically stable and capable of taking into account the moisture diffusion and recondensation process. Two assumptions associated with moisture evaporation modeling, moisture diffusion and vapor condensation were further studied in this Paper. It was shown that the diffusion of moisture towards the drying front speeds up the moisture evaporation and hence decreases the drying time. Vapor condensation increases the interior particle temperature and decreases the surface temperature, c.f., Figures 8 and 9 in Paper I. However, these effects are found to be relatively small especially at high-temperature drying conditions. In Paper II, fundamental processes involved in the gasification and combustion of biomass are discussed in details. The main focus of this Paper is on the devolatilization process. There are major concerns that should be addressed regarding modeling of the devolatilization process; chemical kinetic scheme, kinetic rate constants, secondary reactions and pyrolysis products. In this Paper these topics are discussed extensively. Three different kinetic schemes for pyrolysis of biomass, as presented in Figure 2, are discussed in this Paper. A single-step global pyrolysis model where biomass is converted to its products, volatiles and char through a one-step reaction, a three-step global scheme where conversion is achieved through three competitive reactions, and a multicomponent semi-detailed pyrolysis kinetic model were studied. The latter assumes that the cellulose, hemicellulos and lignin as the main components of biomass are independently pyrolyzed, forming intermediate products. The intermediates are further decomposed into gas, tar and 39

52 char. Each of these pyrolysis schemes has advantages and disadvantages, that was discussed in this Paper. Figure 8. Evolution of particle porosity and effect of particle shrinkage Different kinetic rate constants were employed when using the single-step global pyrolysis scheme to reproduce the experimental data of Park et al. [79], as presented in Figure 10 in Paper II. It was shown that none of the reported values for kinetic rates was able to predict the mass loss of the particle in this experimental set-up. It was possible to tune the kinetic rates to match the numerical results and the experimental measurements. However, these kinetic constants are considered case-dependent and cannot necessarily predict other experimental results. The assumption associated with the single-step pyrolysis scheme is that the particle has a uniform temperature. However, for large particles, there is a considerable temperature gradient inside the particle. Figure 9 shows the final char yield and the particle temperature along its radius at the end of pyrolysis. The temperature between the center and surface of the particle differs by 100 K and the final char yield production changes between 0.14 and 0.18 on dry mass basis. The reason for the different final char yield along the particle radius is that the pyrolysis is undergone at different temperatures and heating rates inside the particle. The single-step pyrolysis scheme is hence not suitable for large particles associated with large temperature gradient. 40

53 Figure 9. Effect of internal particle temperature on the final char yield The three-step global scheme for biomass pyrolysis is well established and proved to be a useful tool for prediction of particle mass loss and the major products of pyrolysis. There are plenty of kinetic rate constants associated with this pyrolysis scheme as reviewed in Ref. [6]. As presented in Figures 12 and 13 in Paper II, the kinetic rate constants of Di Blasi and Branca [80] can predict the mass loss behavior of the particle within acceptable tolerances, while the kinetic data of Thurner and Mann [81] under-predicts the mass loss and suggests a higher solid residue after pyrolysis. These results concur with the findings of other studies [82, 83]. The kinetic data of Thurner and Mann [81] for low-temperature pyrolysis and of Di Blasi and Branca [80] for high-temperature pyrolysis can provide accurate predictions of the process. The shortcoming of the three-step global scheme is that the pyrolysis products are presented as three main groups and detailed information about the gaseous species and tar composition is not provided by this scheme. The kinetic mechanism of Ranzi et al. [17], on the other hand, provides detailed information on gas and tar species produced during biomass pyrolysis. In this semi-detailed mechanism, 20 gaseous species have been carefully chosen to represent the gas and tar yields of the pyrolysis [16, 17]. According to this pyrolysis mechanism, cellulose, hemicellulose and lignin decompose independently. Although they may have interactions during decomposition and also extractives may play a role in the biomass pyrolysis, this is a state-of-art method of simplifying a very complex process. This pyrolysis scheme was considered in the particle model discussed previously and the results were compared with different experimental measurements. The results showed that the multi-component scheme performed well when compared with experiments at low heating rate conditions, c.f., Figures 12 to 14 in Paper II, while the results deviated at high heating rates, Figures 15 and 16 in the Paper. 41

54 In the multi-component scheme, to have a correct estimation of cellulose, hemicellulose and lignin contents of the particle, biochemical analysis of biomass is required. It is possible to estimate the biochemical analysis based on the ultimate analysis; however, this estimation brings more uncertainty into the calculations. Due to this reason, to have an estimation of biomass pyrolysis products, it is possible to link the three-step global scheme with empirical correlations. In Paper II, a method based on the empirical data of Neves et al. [7] and solving a set of algebraic equations to satisfy elemental mass and energy balance was investigated. It was shown that the prediction of the empirical correlations concurred with the measurements, Figure 17 in Paper II. The advantage of using empirical correlations over the detailed kinetic mechanism is simplicity. The results and comparison with different experimental setups show that there are a lot of uncertainties on the final products of the pyrolysis. The yield and composition of gas, char and tar produced in the pyrolysis process depend on several parameters, such as the heating rate, temperature, particle size, composition of biomass, residence time of the particle, sweeping gas flow rate and composition. A mechanism including all the details will be computationally heavy and will not be suitable for CFD applications of real furnaces with many particles. For a typical hardwood particle in millimeter size and in a temperature range of 500 K to 1000 K the model based on the detailed kinetic mechanism of Ranzi et al. [17] is used to predict the main products of the pyrolysis. The algebraic set of equations along with the empirical correlations proposed by Neves et al. [7] were also used to predict the yield of typical hardwood pyrolysis products in this temperature range. The comparison between the results of these two models is presented in Figure 10. Despite the simplicity of the empirical model, the results are in a good agreement. The model can predict the main products yield very well according to the detailed kinetic model. 42

55 Figure 10. Behavior of empirical model compared with detailed kinetic mechanism behavior for typical hardwood; symbols: detail kinetic model prediction, solid lines: empirical model prediction Heat of pyrolysis The heat of pyrolysis depends on pyrolysis products yield and composition. In a previous publication [84] the heat of biomass pyrolysis was investigated. Empirical correlations and thermodynamic data of biomass and the main pyrolysis products were provided to offer a systematic way for computing the heat of pyrolysis. The heat of pyrolysis can be defined as, Q = Y, h h = LHV bm Y j,f LHV j J kg, j (67) where Y, is the final yield of the pyrolysis products on a DAF basis of the original fuel (biomass) and the subscript j denotes char, tar and all of the gaseous species in the system. The correlation presented in Eqs. (12) and (13) can be used to compute the heat of pyrolysis. The pyrolysis products yield and composition are strongly dependent on the fuel source, temperature and heating rate to the particle. This results in a huge variation in the heat of pyrolysis of different biomass sources under different reactor conditions. Table 2 presents values and correlations for heat of pyrolysis of different biomass sources reported in different experimental setups. 43

56 Table 2. Heat of pyrolysis of different biomass sources at different experimental setups Sample Heat of pyrolysis KJ/Kg Ref Spruce 342 [85] Eucalyptus 400 [85] Poplar 207 [85] Sawdust 434 [85] Corn 267 [85] Sunflower 292 [85] Straw 375 [85] Sewage sludge 385 [85] Straw 500 [86] Cellulose Y, [87] Cellulose [88, 89] Pine 200 [90] Oak 110 [90] Spruce wood 3827Y, (1 Y, ) [91] Beech wood 3525Y, + 936(1 Y, ) [91] Maple, low heating rate 610 [92] Maple, High heating rate -105 to -395 [92] By means of Eq. (67) and correlations (12) and (13) and using a detailed mechanism [14], the effect of heating rate and maximum temperature on the heat of pyrolysis was investigated. The correlation between the heat of pyrolysis and the final char yield at different heating rates is presented in Figure 11. The heat of pyrolysis decreases monotonically as the final char yield increases; the slope of the change in the heat of pyrolysis depends on the heating rate. A very high final char yield indicates a very short residence time of the particle or a very low reactor temperature. In such cases, the solid residual is not primarily a carbon structure but is an intermediate composition between biomass and char. These data are discarded from Figure 11. Although there is a clear correlation between the final char yield and the heat of pyrolysis, the increase in the char yield does not cause the exothermicity of the pyrolysis process. This is because the char has higher carbon content and hence a higher heating value compared to the original biomass. Higher char production indicates higher heating value of the products, which results in a higher endothermicity of the pyrolysis process. The variations in the heat of pyrolysis are due to the quality of the produced char and tar at different heating rates or temperatures. Figure 12 shows the changes in the lower heating value of the main 44

57 pyrolysis products as a function of the heating rate to the particle. As the char yield increases, the resulting char has a lower LHV. In addition, for conditions that favor a higher char yield, less energetic tar and gas components are produced. The effect of heating rate on the LHV of the pyrolysis products can explain the behavior of the pyrolysis heat with respect to the final char yield. This behavior can be explained by considering the changes in the LHV of products at different heating rates. As the heating rate increases, the LHV of char and gas increases, which results in higher heating values of the pyrolysis products. The decrease in the LHV of tar as the heating rate increases is smaller than the increase in the LHV of the gas and char at high char yields. Figure 11. Relationship between the heat of pyrolysis and char yield (wt% of DAF) at different heating rates This study showed that higher heating values of the biomass, tar and char can be correlated with the elemental mass fractions of their components with an uncertainty of 20%. In spite of these uncertainties, the calculation of heat of pyrolysis based on this method can guaranty a strict elemental and energy conservation between the biomass and the pyrolysis products. The model results are consistent with the trends of the experimental data. 45

58 Figure 12. Changes in the lower heating values of char, tar and gas at different heating rates; the thick line in each group represents the highest heating rate. Release of alkali metals from biomass Release of alkali metals during biomass pyrolysis and gasification can cause serious damage to the combustion system due to fouling, slagging and agglomeration [93-96]. Previous investigations show that fouling is mainly due to the sticky deposition layer of alkali metals released from biomass combustion [97]. High-temperature corrosion in furnaces and in the turbine blades in combined cycle technology or reduced heat transfer in boilers due to deposits in convective surfaces are examples of these damages. The deposition of potassium and sodium on the active sites of vanadium and tungsten NO x reduction catalysts in power plants reduce the catalysts efficiencies and lead to more emissions [98]. To prevent damages caused by these species it is thus important to understand the release of alkali metals during thermochemical conversion of biomass. In Paper III, a joint numerical and experimental investigation of the release of potassium from biomass during gasification stage was presented. In this Section the method is further expanded to investigate the release of sodium and potassium during the devolatilization stage. 46

59 Experimental measurements A laser-induced breakdown spectroscopy (LIBS) was adopted to measure quantitatively the concentration of alkali metals from gasifying biomass. The effect of temperature on the potassium and sodium release was investigated in hot gas mixture of CO 2 and H 2 O with different concentrations of O 2. Experiments were carried out to measure the elemental mass of potassium and sodium of a biomass particle placed in a hot flue gas from combustion of methane/oxygen/carbon-dioxide under fuel-lean and fuel-rich conditions. The experimental results and the numerical results were used to delineate the various stages of potassium and sodium release. To measure the release of potassium and sodium from gasifying biomass, a cylindrical wood particle is placed on top of a laminar flame burner. The details of the experimental set-up and measurement technique can be found in Paper III. The experimental conditions are presented in Table 3. Two equivalence ratios and four different temperatures are studied. At an equivalence ratio of 0.9, there will be oxygen available in the flue gas at the particle boundary layer, while at an equivalence ratio of 1.2, biomass undergoes only gasification with water and carbon dioxide. Table 3. Temperature profile along radius of the burner exit plane for various test cases Equivalence ratio Temperature profile along radius [K] a Label r=0 r=10mm r=20mm r=30mm r=33mm F F F F4 a measured by an R-type thermocouple Figure 13 shows direct images from the experimental study providing qualitative information of different stages of conversion and correspondent release behavior of the particle. At the onset of this figure, the devolatilization stage is presented. The peak of the potassium release corresponds to the highest particle mass loss rate. Towards the end of the devolatilization stage, the potassium concentration decays to a very low level. When the char reactions start, the concentrations of potassium and sodium increase continuously and reach their peak values almost at the end of the char reaction stage. The concentrations monotonically decay during the ash-cooking stage. The yellow light emitted during the ash-cooking stage is due to the presence of sodium. 47

60 Figure 13. Direct images from experimental setup; a particle at three stages of conversion The measured concentrations of K and Na for four cases, F1-F4, are presented in Figure 14. The experimental data for each flame condition are marked with symbols and the solid lines show the filtered data. Each measurement point is an average of data from three different biomass samples. The measured data presented in Figure 14 are divided into two parts: the data obtained during the devolatilization stage and the data obtained during the char reaction and ashcooking stages as they show different behavior. Note that the process of heating, melting, shrinking and splitting of the ash particle in a high temperature environment is referred to as ash-cooking. The Na concentration is considerably lower during the char reaction and ashcooking stages. The Na concentration for the F3 case was below the sensitivity of the measurement technique. Moreover, for the F4 case, the signal to noise ratio of the sodium concentration is very high due to the low values of sodium release in this case. On the other hand, the sodium concentration during the devolatilization stage is higher than that of potassium. Considering the ultimate analysis of the particle, as presented in Table 1 in Paper III, the initial sodium mass fraction in the particle is lower than the potassium mass fraction, which indicates that a higher percentage of the sodium is released during the devolatilization stage. 48

61 Figure 14. Measured mass concentration of potassium and sodium at four different flame conditions From Figure 14, the concentrations of K and Na show two distinct peaks associated with the devolatilization and char combustion. After devolatilization finishes, the concentration of K and Na first decreases significantly and then starts to increase as the char reactions proceed. This phenomenon indicates that two different mechanisms are involved in the release of K and Na during the devolatilization and char reaction stages. During devolatilization, the release of K and Na is mostly due to evaporation of inorganic species inside the particle, whereas during char reaction the main source for K in the gas phase is the decomposition of carboxylic acid groups formed during the devolatilization and reaction of carbonate and sulfate groups. The mechanism for release of alkali metals is elaborated in Paper III. 49

62 Numerical analysis In case F4 the ambient gases around the particle are the products of a rich methane flame, meaning that no oxygen is available in the vicinity of the particle. For the other cases, F1-F3, the oxygen left from the lean combustion of methane can affect the burning behavior of the particle, which can be the reason for the longer burnout time in F4 compared with the other three cases. To understand the role of oxygen, possibility of air entertainment and accurate boundary conditions of the particle model, a CFD calculation of the flow and transport and reactions of chemical species around the particle was performed. Boundary conditions for the CFD simulation were extracted from the particle model described earlier. The temporal evolution of the mass fractions of the gaseous products as well as the velocity of the gaseous mixture at the surface of the particle was used as the boundary conditions for the CFD simulations. The burner out-flow condition was calculated by a CHEMKIN type solver and the GRI mechanism [99]. Figure 15 shows the heat release rate, temperature and mass fractions of O 2 and CO around the particle at an instance of time during the char combustion stage for the F3 case. The hollow region in the middle of the CFD domain represents the particle boundary. The CFD domain at the top is restricted by the presence of a cap. The reason for placing the cap was to stabilize the flame formed around the particle especially during the devolatilization stage. Figure 15. Heat release rate, temperature, oxygen and carbon monoxide concentration in gas phase around the particle 50