Contribution of Gas-phase Reaction on Heat Output from Smoldering Packed-bed of Dried Leaves INOUE, Yasunobu* 1,2, SUZUKI, Masataro 2 1

Similar documents
Department of Mechanical Engineering, University Visvesvaraya College of Engineering, Bangalore University, Bangalore, Karnataka, India

Numerical Examination of Two-Dimensional Smolder Structure in Polyurethane Foam

Combustion. Indian Institute of Science Bangalore

Transactions on Engineering Sciences vol 5, 1994 WIT Press, ISSN

Modeling of the pyrolysis of plywood exposed to heat fluxes under cone calorimeter

Experimental and Theoretical Study of the Ignition and Smoldering of Wood Including Convective Effects

Chapter 5 Test. Directions: Write the correct letter on the blank before each question.

11B, 11E Temperature and heat are related but not identical.

Measuring and Expressing Enthalpy Changes. Copyright Pearson Prentice Hall. Measuring and Expressing Enthalpy Changes. Calorimetry

Critical mass flux for flaming ignition of dead, dry wood as a function of external radiant heat flux

Chemical Energetics. First Law of thermodynamics: Energy can be neither created nor destroyed but It can be converted from one form to another.

SUPPLEMENTARY TOPIC 3 ENERGY AND CHEMICAL REACTIONS

CHAPTER 16 REVIEW. Reaction Energy. SHORT ANSWER Answer the following questions in the space provided.

Balancing chemical reaction equations (stoichiometry)

Section 1 - Thermochemistry

Elements and Their Oxides

1.4 Enthalpy. What is chemical energy?

Hydrogen addition to the Andrussow process for HCN synthesis

Chemical reactions. C2- Topic 5

Quantitative Study of Fingering Pattern Created by Smoldering Combustion

Personalised Learning Checklists AQA Chemistry Paper 2

Date: SCH 4U Name: ENTHALPY CHANGES

Chemistry Slide 1 of 33

Chapter 6 Thermochemistry 許富銀

17.2 Thermochemical Equations

Laminar Premixed Flames: Flame Structure

(g) 2NH 3. (g) ΔH = 92 kj mol 1

3.2 Alkanes. Refining crude oil. N Goalby chemrevise.org 40 C 110 C 180 C. 250 C fuel oil 300 C 340 C. Fractional Distillation: Industrially

Lecture Presentation. Chapter 6. Thermochemistry. Sherril Soman Grand Valley State University Pearson Education, Inc.

CONTROLLING MECHANISMS IN THE TRANSITION FROM SMOLDERING TO FLAMING OF FLEXIBLE POLYURETHANE FOAM

GLOWING AND FLAMING AUTOIGNITION OF WOOD

Burning a Hydrocarbon II

COMBUSTION OF FUEL 12:57:42

ONE-DIMENSIONAL MODEL OF PYROLYSIS AND IGNITION OF MEDIUM DENSITY FIBERBOARD SUBJECTED TO TRANSIENT IRRADIATION

Period 11: Chemical Energy and Fossil Fuels

Chemistry 101 Chapter 10 Energy

Careful observations led to the discovery of the conservation of mass.

Interactions between oxygen permeation and homogeneous-phase fuel conversion on the sweep side of an ion transport membrane

Most hand warmers work by using the heat released from the slow oxidation of iron: The amount your hand temperature rises depends on several factors:

Chemistry: The Central Science. Chapter 5: Thermochemistry

MARIYA INTERNATIONAL SCHOOL. Work sheet III. Term I. Level 8 Chemistry [MCQ] Name: CHEMICAL REACTIONS & SULFUR

BONDING FIRE RETARDANTS TO WOOD. PART I. THERMAL BEHAVIOR OF CHEMICAL BONDING AGENTS

AQA Chemistry Checklist

Reacting Gas Mixtures

Enthalpy Changes. Note: 1 calorie = 4.2 Joules

N L N G : C O M B U S T I O N

Gravity is a force which keeps us stuck to the earth. The Electrostatic force attracts electrons to protons in an atom.

Name Date Class THERMOCHEMISTRY

The Atmosphere Exam questions

Chapter Objectives. Chapter 9 Energy and Chemistry. Chapter Objectives. Energy Use and the World Economy. Energy Use and the World Economy

HKCEE Past Paper Questions: Part 9 Rate of Reactions Part A: Multiple Choices

Name Class Date. As you read Lesson 17.1, use the cause and effect chart below. Complete the chart with the terms system and surroundings.

Experimental investigation of Methane Partial Oxidation for Hydrogen Production

Methane contains atoms of two elements, combined chemically. Methane is a mixture of two different elements.

Onset of smoldering and transition to flaming fire

An acid made from sulphur, oxygen and hydrogen. A chemical with a ph value. of 8 or more.

THERMOCHEMISTRY & DEFINITIONS

Chemistry Chapter 16. Reaction Energy

Notes: Balancing Chemical Equations

Prediction of the Heat Release Rate of Douglas Fir

Lecture 0503 Calorimetry and Hess s Law

1º ESO UNIT 4: Chemical and physical changes. Susana Morales Bernal

Understanding Equations

Enthalpies of Reaction

Name /100. 1) Matter is defined as anything that is visible to the human eye. 1) 2) An amorphous solid has long range, repeating order.

Paper Reference. Sample Assessment Material Time: 2 hours

Thermochemistry. Energy. 1st Law of Thermodynamics. Enthalpy / Calorimetry. Enthalpy of Formation

A FLAMELET MODEL FOR DIFFUSION FLAMES IN POROUS MEDIA. Max Akira Endo Kokubun. Fernando Filho Fachini

AQA Chemistry (Combined Science) Specification Checklists. Name: Teacher:

Chapter 3. Stoichiometry: Calculations with Chemical Formulas and Equations

AQA TRILOGY Chemistry (8464) from 2016 Topics T5.1 Atomic structure and the periodic table (Paper 1) To pic. Student Checklist

Mr Chiasson Advanced Chemistry 12 / Chemistry 12 1 Unit B: Thermochemical Changes

2. Enthalpy changes. N Goalby chemrevise.org

OH, is an important feedstock for the chemical industry.

The Nature of Energy. Chapter Six: Kinetic vs. Potential Energy. Energy and Work. Temperature vs. Heat

A Numerical Study of Solid Fuel Pyrolysis under Time Dependent Radiant Heat Flux Conditions

First Law of Thermodynamics: energy cannot be created or destroyed.

8. Energetics I. N Goalby chemrevise.org 1

Chapter 5 Thermochemistry

What does rate of reaction mean?

Flame Spread and Extinction over Thermally Thick PMMA in Low Oxygen Concentration Flow

Well Stirred Reactor Stabilization of flames

SCIENCE 9 CONCEPT 4 CHEMICAL REACTIONS

Introductory Chemistry Fourth Edition Nivaldo J. Tro

5/14/14. How can you measure the amount of heat released when a match burns?

Thermochemistry: Energy Flow and Chemical Reactions

Enthalpy changes

Chapter 6. Energy Thermodynamics

DARS overview, IISc Bangalore 18/03/2014

Chemical Reaction (IGCSE Chemistry Syllabus )

Underground Coal Fire Extinction Model using Coupled Reactive Heat and Mass Transfer Model in Porous Media

Energetics. These processes involve energy exchanges between the reacting system and its surroundings.

ALE 26. Energy Changes ( E) and Enthalpy Changes ( H) in Chemical Reactions

s Traditionally, we use the calorie as a unit of energy. The nutritional Calorie, Cal = 1000 cal. Kinetic Energy and Potential Energy

Thermochemistry Ch. 8

Calorimetry. Chapter 5. Week 2 Unit 1. Calorimetry. Since we cannot know the of the reactants and products, we measure H through, the of.

Based on the kinetic molecular theory of gases, which one of the following statements is INCORRECT?

Microgravity Opposed Flow Flame Spread in Polyvinyl Chloride Tubes

ADVANCED DES SIMULATIONS OF OXY-GAS BURNER LOCATED INTO MODEL OF REAL MELTING CHAMBER

1.4 Energetics. N Goalby chemrevise.org 1. Standard Enthalpy Change of Formation. Standard Enthalpy Change of Combustion

Transcription:

Contribution of Gas-phase Reaction on Heat Output from Smoldering Packed-bed of Dried Leaves INOUE, Yasunobu* 1,2, SUZUKI, Masataro 2 1 Japan Tobacco Inc. Tobacco Science Research Center 2 Nagaoka University of Technology Abstract The heat generation from gas-phase reactions in vacant spaces inside a smoldering packed-bed is evaluated experimentally. The bed is a pile of cut leaves of dried plants, having a bulk density of 0.22 g/cm 3, in a thermallyinsulated vertical quartz tube of 17 mm inner diameter and 150 mm height. Its top is ignited so that the smoldering progresses vertically downward. Air is supplied from the bottom of the tube against the smoldering direction at a fixed flow rate of 30 mm/s. The temperature in the packed-bed is measured and concentrations of exhaust gas components are analyzed by means of gas chromatography, results of both which confirm that the smoldering process reaches a quasi-steady state. The smoldering is assumed to consist of three reactions: the pyrolysis of unburned leave generating char and pyrolysis gas, the surface combustion of the char generating exhaust gas and residue, and the gas-phase combustion of the pyrolysis gas also generating exhaust gas. The amount of heat generated from each of these reactions is evaluated from the balance of enthalpies between reactants and products in the reaction. Enthalpies of unburned leaves, the char, and the residue are measured by means of bomb calorimetry after those samples are collected from the corresponding regions in the tube. The exhaust gas is collected in a gas bag, and concentrations of oxygen, carbon dioxide, carbon monoxide and methane are measured by gas chromatography with thermal conductivity detection. Particulates in the exhaust gas are collected by a glass fiber filter, and its enthalpy is measured by means of bomb calorimetry. The pyrolysis gas, on the other hand, is simulated by a gas that is generated from unburned leaves while those are heated from room temperature to 400 ºC in air, and is collected and analyzed in the same way as the exhaust gas. The heat generated in the gas phase in vacant spaces in the bed is estimated at 8.1 J/s, which is about 40% of the total heat output of 20.1 J/s from the smoldering of the packed-bed. Introduction 1 Smoldering is a slow, low-temperature, flameless and self-sustaining combustion [1] ; it can occur through variety of materials such as wood [2][3], sawdust [4], flexible foam [5][6], and other biomass combustibles [7][8]. It includes a variety of chemical reactions and heat and mass transports in gas and solid phases, and it emits various combustible gases [7][8]. Smoldering of a biomass combustible generally involves two stages of chemical reactions: pyrolysis (thermal decomposition) and char oxidation [6]. In order for the smoldering to continue, sufficient heat must be transferred from the region of exothermic reactions to the adjacent unburned material so that the temperature becomes high enough to initiate the oxidation of the material [1]. Thus, oxidation reactions play an important role in the propagation of smoldering. One of typical self-propagating forms is the opposed downward smoldering, in which the reaction front propagates in the direction opposite to the oxidant flow. Under this condition, oxidant gas flows into the fresh part of the combustible material. The fresh material is heated and dried owing to conductive and radiative heat transfers from the downstream, and then is thermally decomposed under the presence of oxygen. Oxidation reactions follow depleting oxygen in the flow while it passes through the region of char and residue and heat is generated there [9] [10]. As it is noted in the review [11] of the previous studies, reactions have * Corresponding author: yasunobu.inoue@jt.com Proceedings of the European Combustion Meeting 2015 been regarded as two overall reactions are assumed representing two stages, namely: 1. Pyrolysis under oxidative atmosphere producing char and gas, and 2. Oxidation of the char. Ohlemiller and Lucca [12] have examined their reaction model for polymer and cellulosic material and have found that the first stage is predominant. This assumption can simplify it to a one-step reaction model. Solid Fuel + O 2 Solid Product + Gaseous Product + Heat (1) They concluded that the propagation rate of the smoldering is dependent on the heat transfer of both conduction and radiation, and that the intensity, or the rate of heat output, is dominated not only by oxygen supply but also by the reactivity of the fuel. Dasanjh et al. [13] also have adopted this one-step reaction scheme to formulate an opposed smoldering model. They pointed out that solid products are reactive and oxidation of those may vary depending on its pyrolysis reaction hysteresis. The rate of smoldering propagation is calculated by the model with consideration taken into the detail oxidation reaction of solid products. However, the intensity of smoldering measured experimentally is not always in good agreement with that calculated theoretically [8]. Further information about oxidation reactions will be required for improving the theoretical understanding of the biomass smoldering. While oxidation reactions on the solid surface have been examined considerably, those in the gas phase have been paid little attention even though it may occur in vacant spaces in the bulk of the biomass. There is a possibility that the discrepancy between

experiment and theory narrows if the gas-phase reactions are added into the theory. Thus, in this study, the contribution of gas-phase reactions on the heat generation from the smoldering of a packed bed of a biomass combustible is investigated. Reaction model and estimation of heat output The reaction model adopted in this work is schematically shown in Fig.1. It consists of three processes: the pyrolysis of unburnt solid, the char oxidation, and the partial combustion of the pyrolysis gas. The pyrolysis process generates gas species such as CO, CO 2, H 2, and CH 4 [14][15][16]. The generated gas from this process containing these species is referred to as a pyrolysis gas in this paper. Although pyrolysis reactions are commonly endothermic, this process may be exothermic since oxygen is present in the atmosphere and the process may include oxidation reactions. Subsequent to that, in the gas phase, the pyrolysis gas is oxidized partially into the combustion gas like CO 2 and H 2O before it is discharged. The char left in the solid phase also undergoes oxidation due to surface combustion. The total heat output H t is written as follows: H t = H 1 + H 2 + H 3 (2) where subscripts 1, 2, and 3 are for processes of the pyrolysis of unburnt solid, the char oxidation, and the partial combustion of the pyrolysis gas, respectively. The amount of heat generated from each reaction is evaluated on the basis of the balance of enthalpies of the reactant and the product of the reaction. Some of those products are sampled directly from the smoldering test, and others are done from a simulating test before enthalpies of those are measured using a bomb calorimeter (Shimadzu CA-4AJ) or evaluated on the basis of the ingredient analysis by GC-TCD (INFICON 3000 micro GC). Estimation of heat generated from solid and gas phase oxidation The heat generated from the pyrolysis process is calculated as follows: H 1 = H s (Y c H c +Y p H p ), (3) The heat generated from char oxidation is also calculated by using the equation. H 2 = Y c H c (Y r H r +0) (4) The heat generated from gas-phase oxidation is estimated by the following equation. H 3 = H t ( H 1 + H 2 ) (5) where H and Y are the enthalpy and the yield, respectively, and subscripts s, c, p, and r correspond to the unburnt biomass, the char, the pyrolysis gas, and the residue respectively. Each yield is the mass ratio of each product to the unburnt biomass. The enthalpy of each substance is regarded as the thermal energy generated from its complete combustion: Unburnt biomass +Oxygen Combustion gas +Ash + H s (6) Fig 1 Reactions of heat generation for smoldering biomass Char + Oxygen Combustion gas+ash + H c (7) Pyrolysis gas+oxygen Combustion gas + H p (8) Residue+Oxygen Combustion gas + Ash + H r (9) Experimental Sample Material The sample combustible adopted in this work is a packed bed of dried and cut leaves. An average size of a cut leaf is 1.0 mm in width, 10 mm in length and 0.1mm in thickness. Table 1 shows some properties of the sample; mass fractions of elements are measured by an elemental analyzer (SCAS Sumigraph NC-220F) and the apparent density of the sample by a porosimeter (Quanta chrome Pore Master GT60). Table 1 Properties of dried leaves Elemental analysis C (wt. %) 44.8 H (wt. %) 6.3 N (wt. %) 4.8 Apparent density 0.6 g/cm 3 Smoldering Test Smoldering experiments are conducted in a room with an air conditioning facility that controls the air to the standard ambient temperature and pressure (SATP). Figure 2 shows a schematic of the experimental setup. Cut leaves of dried plants are filled in a quartz tube of 17 mm inner diameter, which is covered with a heatinsulator composed of activated carbon fibers. Pieces of leaves are packed in the tube such that longitudinal directions of those pieces are parallel to the axial direction of the tube. The bed of those pieces has the bulk density at 0.22 g/cm 3, the void fraction at 0.8, and the height at 150 mm. The bed is ignited at the top by means of an electrically-heated nichrome wire. The combustion propagates towards vertically-downward in countercurrent contact with the air flow, the velocity of which is about 30 mm/s. This is the lowest value to prevent the combustion from extinguishing.

90 mm 100 mm 150 mm Temperatures on the radial center of the test tube are measured at the heights of 90, 95, and 100 mm from the bottom with thermocouples of 0.5 mm in diameter. Data are collected by a data logger at 1.0 s intervals. Simulating Test With regard to the pyrolysis gas, direct in-situ sampling is not possible because gas-phase reactions are supposed to occur in this work; the pyrolysis gas generated from the solid may be oxidized in the gas phase of vacant spaces inside the smoldering packedbed before it is collected. Therefore, the pyrolysis gas is simulated by a gas generated from 1 g of sample in the atmosphere of dry air while it is heated by an infrared image furnace (ULVAC image furnace RHL- E25) from ambient temperature to 400 ºC, below at which oxidation reactions are not supposed to react shown in fig.3. The pyrolysis test is conducted using the same tube as smoldering test and the collection regime is the same as the smoldering test. Sampling of Products Char and residue are collected after the smoldering bed is extinguished by supplying nitrogen gas into the tube. A typical temperature distribution and sampling regions are shown in Fig.3. The material in the pyrolysis region, where the temperature is about 400 ºC is collected as the char. The residue is collected from the region where the temperature once raised over 400 ºC and then went down below 400 ºC. The exhaust gas containing particulate matter is collected directly while the smoldering is taking place. The particulate matter is captured by a glass fiber filter and its enthalpy is measured by means of bomb calorimetry. The filtered gas, on the other hand, is collected into a gas bag and its ingredients are directly analyzed by means of gas chromatography with thermal conductivity detection. Measurement of Enthalpies The enthalpies of solid materials; sample (H s), char (H c), residue (H r) are measured using a bomb calorimeter (Shimadzu CA-4AJ). The enthalpies of pyrolysis gases (H p) are the sum of the enthalpies of gaseous materials and particulate matter. The enthalpies of gaseous materials are estimated as the sum of the enthalpy of their components (CO, CH 4). The enthalpy of particulate matter is collected by a glass filter pad and measured by a bomb calorimeter. The enthalpies of these materials per unit weight of original sample are calculated by multiplying the each enthalpy of materials and recovery coefficient of each material per unit weight of original sample. The recovery coefficient is calculated by the following equation. Y i = g i / g unburnt (i = s, c, and r) (11) Where, g i is a weight of each sampling material per unit length of the combustion tube after extinguish. g unburnt indicates a weight.of unburnt material per length of the tube before ignition. The subscripts s, c, Combustion Tube Φ = 17 mm L = 200 mm Thermocouples Sample packed-bed Heat Insulator Gas Flow GC-TCD Filter Igniter MFC Air Air cylinder Fig 2 Schematic of experimental apparatus (smoldering test) Fig 3 Schematic of experimental apparatus (simulating test) Fig 4 A typical result of temperature distribution in the tube and sampling region based on it and r mean the sample, the char, and the residue, respectively. Estimation of total heat output Total heat output ( H t [J/g-unburnt]) is calculated by the same method as a cone calorimeter as followings [17]. H t = 13100(C 0 C) 1.105 1.5C V airam o t (12) here C 0 [-] is the initial oxygen concentration of inlet gas, C [-] is the oxygen concentration of exhaust gas, V air [cm/s] is the air flow rate, A = 2.27 cm 2 is the cross-sectional area of the tube, M o = 1.43 10 3 g/cm 3 is the oxygen density, t = 236 s is the combustion duration time consumed until unburnt sample 1g burns out. Results Temperature Distribution The typical temperature distribution in the tube is shown in Fig.4. As shown, these three temperature curves at three different positions are similar to each other. That means the smoldering has progressed in

quasi-steady state. The smoldering propagation rate V s is estimated to be 9.3 10 3 cm/s based on this measurement. Gas Concentration in Exhaust Gas Figure 5 shows the variations of each gas concentration in exhaust gas. The concentration of each component becomes to be stable after 240 sec. from ignition. Due to the stability of each gas concentration, we obtain the quasi-state gas concentration for the smoldering bed. The exhaust gas is contained in a concentration of almost 0% oxygen during quasi-state smoldering. Heat of combustion The heats of combustion and yields of unburnt, char, residue and pyrolysis gas are shown in table 2. The heat output from each reaction per second ( h) is calculated by using following equation; h i = H i V s A ρ (13) V s is the smoldering propagation rate, A is the crosssectional area of the tube, ρ is the apparent density of unburnt in the bed. The heats of combustion of CO and CH 4 are obtained by reference to the literatures [18]. Table3 shows the heats output from the reactions calculated by using the equations as described previously. The total heat output from the smoldering bed, from pyrolysis reaction, from oxidation reaction on surface of char, and the heat generated in the gas phase are estimated at 20.1 J/s, 4.7J/s, 7.4J/s, and 8.0 J/s. The contribution of gasphase reaction on heat output from the smoldering packed-bed is about 40 % of the total heat output from the smoldering. Table 2 heat of combustion and yield for each material Material Yield Heat of combustion Component g/gunburnt J/g Unburnt 1 5251 Char 0.41 5720 Residue 0.11 5620 Pyrolysis compound CO 0.11 10100 CH 4 0.005 55640 Particle 0.089 4510 Table 3 Heat generation from reactions Reaction Heat output J/gsample Heat output J/s Total heat output 4741 20.1 Pyrolysis 1115 4.7 Surface combustion 1747 7.4 Gas-phase combustion 1879 8.0 Fig 5 Concentration profile of each component in exhaust gas analyzed by GC-TCD Fig 6 Conversion rate of methane for catalyst layer at various temperatures of catalyst layer and space velocities. Discussion The result, which is the heat generated in the gas phase is about 40% of the total heat output, is indicated that the oxidation of pyrolysis gases occurs in gasphase at comparatively low temperature to flame type processes which is famous for oxidation of flammable gas in gas-phase. Naganuma et al. [19] investigate the conversion rate of methane through catalyst layer under mixed gas of 3 % methane helium balance. The conversion rate of methane is increasing with increasing temperature of the catalyst layer or with decreasing space velocity shown in Figure 6. In our study, the length of region at over 600 ºC is 10mm based on temperature distribution. Therefore, the space velocity through the region over 600 ºC is estimated to be 10440 h -1 on our packed-bed. At the same as temperature and space velocity as our study (plotted in Figure 6), the conversion rate is estimated to be about 80%. The metal salts, which are comprised of biomass materials, sometimes show catalytic activities of oxidation reactions at smoldering process [20], this means, the oxidation of pyrolysis gas, such as carbon monoxide and hydrocarbons, which are the main byproducts of biomass smoldering combustion, may occur at the temperature in smoldering biomass process. This result and estimation imply that gas-phase reactions should be taken into account for an accurate prediction of heat output from a smoldering packedbed.

Conclusion Previous studies focus on the heat generated from surface oxidation in solid-phase for smoldering processes. In this study, we evaluate the heat generation from gas-phase reactions in vacant spaces inside a smoldering packed-bed experimentally. Dried plant leaves are filled in a quartz tube of 17 mm inner diameter with a bulk density of 0.2 g/cm 3. The air is supplied against the smoldering direction at a fixed flow rate of 30 mm/s. The heat of combustion of gasphase reactions is estimated from the balance of enthalpies of materials picked from the tube. The heat generated in the gas phase is estimated at 8.1 J/s, which is about 40% of the total heat output of 20.1 J/s from the smoldering. References [1] T.J. Ohlemiller, Modeling of Smoldering Combustion Propagation, Progress in Energy and Combustion Science 11 277-310 (1985) [2] K.N. Palmer, Smoldering combustion in dusts and fibrous materials, Combustion and Flame 1:129-54 (1957) [3] R. Bilbao, J.F. Mastral, M.E. Aldea, J. Ceamanos, M. Betran, Experimental and theoretical study of the ignition and smoldering of wood including convective effects. Combustion and Flame 126 1363-72 (2001) [4] W.C. Sun, M. Xie, Propagation and gas components of vertical smoldering with natural convection. Journal of Combustion Science and Technology 8 188-91 (2002) [5] J.H. Wang, C.Y.H. Chao, W.J. Kong Experimental study and asymptotic analysis of horizontally forced forward smoldering combustion. Combustion and Flame 135 405-19 (2003) [6] J. Porteiro, D. Patino, J. Collazo, E. Granada, J. Moran, J.L. Miguez Experimental analysis of ignition front propagation of several biomass fuels in a fixedbed combustor. Fuel 89(1) 26-35 (2010) [7] M.S. Saidi, M.R. Hajaligol, A. Mhaisekar, F. Rasouli Numerical simulation of burning cigarette during puffing. Journal of Analytical and Applied Pyrolysis 72 141-52 (2004) [8] A. Rostami, J. Murthy, M.R. Hajaligol Modeling of smoldering cigarette Journal of Analytical and Applied Pyrolysis 66 281-301 (2003) [9] J.L. Torero, A.C. Fernandez-Pello, M. Kitano Opposed forced flow smoldering of polyurethane foam Combustion Science and Technology 91 95-117 (1993) [10] A. Bar-llan, G. Rain, D.C. Walther, A.C. Fernandez-Pello, J.L. Torero, D.L. Urban The effect of buoyancy on opposed smoldering Combustion Science and Technology 176 2027-2055 (2004) [11] F.E. Roger, T.J. Ohlemiller, Smoldering Characteristics of flexible Polyurethane Foams, Journal of Fire and Flammability 11 32-44 (1980) [12] T.J. Ohlemiller, D.A. Lucca An experimental comparison of forward and reverse smolder propagation in permeable fuel beds. Combustion and flame 11 277-310 (1983) [13] S.S. Dosanjh, P.J. Pagni, A.C. Fernandez-Pello Forced concurrent smoldering combustion Combustion and Flame 68 131-142 (1987) [14] L. Bennamoun L, A. Belhamri Mathematical description of heat and mass transfer during deep bed drying: Effect of product shrinkage on bed porosity. Appl Therm Eng. 28(17-18):2236-44 (2008) [15] M. Asadullah, M.A. Rahman, M.M. Ali, M.S. Rahman, M.A. Motin, M.B. Sultan et al. Production of bio-oil from fixed bed pyrolysis of bagasse. Fuel 86(16):2514-20 (2007) [16] A. Zabaniotou, T. Damartzis, Modeling the intraparticle transport phenomena and chemical reactions of olive kernel fast pyrolysis. J Anal. Appl. Pyrol. 80(1):187-94 (2007) [17] C. Huggett, Estimation of rate of heat release by means of oxygen consumption measurements Fire and Materials 4(2) 61 65 (1980) [18] M. Himmelblau, Basic Principles and Calculations in Chemical Engineering, 4th Ed., Prentice-Hall, New Jersey (1982) [19] Y. Naganuma, Investigation of the Burning Conditions of Catalytic Combustion in a Low Concentration of Methane, Transactions of the Japan Society of Mechanical Engineers B52(473),479-484 (1986) [20] Vassilev, David Baxter, Christina G. Vassileva, An overview of the behaviour of biomass during combustion: Part I.-mineral transformations of organic and inorganic matter, Fuel, 112, 391-449 (2013)