Final Corrected Thesis
Final Corrected Thesis
Final Corrected Thesis
Equilibrium Modeling
by
Gopal Gautam
Auburn, Alabama
December 13, 2010
Approved by
Biomass has already emerged in the renewable energy arena as one of the promising
candidates for the future. Biomass has been a major source of fuel for human from the existence
of mankind. Rapid urbanization and widespread use of fossil fuels in the industrial world has
relegated it to the status of a minor source of energy. The rejuvenation, however, started with
increasing concerns over reducing carbon footprints and also due to strong causative connections
between non-renewable fossil fuels and global warming. Biomass gasification is a thermo-
chemical process of converting biomass into the producer gas or syngas (used interchangeably)
which can be subsequently used for heat, power and liquid fuels production through various
synthesis processes. The major objective of this study was to better understand the effect of
various parameters on the syngas composition from a stratified downdraft gasifier. The study is
primarily experimental but supplemented by the mathematical modeling that explains various
Chapter 1 provides basic literature review on the gasification process, various types of
gasifiers and elaborated discussion about the effect of various parameters on syngas composition
for different types of gasifiers. The effects were primarily discussed for updraft, downdraft and
fluidized bed gasifiers which currently cover more than 98% of the total biomass gasification
market.
Chapter 2 presents a thermodynamic model of the process inside the gasifier. Syngas
composition is predicted assuming thermodynamic equilibrium condition inside the gasifier. The
ii
thermodynamic equilibrium can be assumed because residence time is high in the fixed bed
gasification process. The effect of moisture content as well as temperature was determined. The
model was run for nearly 100 samples. Based upon the results of the simulations, using linear
regression analysis, equations were derived to predict the syngas composition of the biomass
moisture content and biomass flow rate on syngas composition in the downdraft gasifier.
Parameters studied are moisture content and biomass flow rate inside the gasifier. A mass,
energy and exergy analysis is also done to corroborate the experimental results as well as to
visualize the carbon, heat, and availability loss inside the gasifier in the process.
Chapter 4 discusses tar downdraft gasifier. Although the amount of tar from a downdraft
gasifier is always assumed to be small, it is more stable and might adversely affect when used for
power generation. Significant amount of toluene, o/p-xylene, naphthalene, phenol, styrene and
indene was observed. Tar concentration in the syngas from the gasifier was found to be 0.34-0.68
mg/Nm3.
iii
Acknowledgments
I would like to express my sincere gratitude to Dr. Sushil Adhikari and Dr. Sushil
Bhavnani for their constant support, guidance and help throughout this thesis work. Discussions
with them have formed the basic background for this work. Without their help, this thesis could
I am also very grateful to the Center for Bioenergy and Bioproducts and Alabama
Agricultural Experiment Station of Auburn University for providing funding for some portion of
this study. I would like to thank Mr. Christian Broadbeck for his immense help while conducting
the experiments with the gasifier. Without his help, much of this thesis would have remained
undone. Also Jonathan Griffith deserves special thanks for helping me with the gasification
experiments. I would also like to acknowledge the help from my colleague, Ms. Suchithra
Thangalazhy-Gopakumar for help in various experimental works conducted for this thesis.
Finally, I would like to thank my parents for being the source of motivation and support
for me throughout my entire academic career. They are always the source of my perseverance,
iv
Table of Contents
Abstract ........................................................................................................................................... ii
Acknowledgments.......................................................................................................................... iv
v
1.4.9 Bed Material ....................................................................................................27
vi
3.3.1 System Description .........................................................................................77
3.4.6 Carbon, Energy and Exergy Analyses with Commercial Wood Pellets .........95
vii
Appendix A. MATLAB Code for Syngas Equilibrium Modeling for Adiabatic Conditions ..122
A.3 Function File for Finding the Enthalpy Change in Gases .......127
Appendix B. Function File for Finding Syngas Composition at Constant Equivalence Ratio....128
Appendix C. Syngas Composition from MATLAB Simulation Used for General Formula
Derivation .... 132
Appendix D. Supplemental Data for Selected Figures........ 135
Appendix E. Sample Calculations ..........................................................146
E.1 Sample Calculations for Equilibrium Modeling....................................................................146
E.2 Carbon, Energy and Exergy Analyses ....148
viii
List of Tables
Table 1.1 Various characteristics, properties and difference between common types of
gasifier. .................................................................................................................................8
Table 1.2 Effect of moisture content upon major syngas constituents ..........................................11
Table 1.3 Optimal ER for some feedstocks in Downdraft and Fluidized gasifiers .......................14
Table 1.4 Ash content and its elemental composition for some common feedstocks.
(% dry basis) ......................................................................................................................18
Table 2.1 Coefficients of specific heat capacity for various gases ...............................................52
Table 2.2 Enthalpy of formation and coefficient for Eqn. (19) .....................................................53
Table 2.3 CO and H2 composition for most common feedstocks available in the U.S from
MATLAB model ................................................................................................................55
Table 3.6 Moisture content and biomass flow rate for different feedstocks .92
Table 3.8 Coefficients for specific heat capacity of various gases ................................................97
Table 3.9 Carbon, energy and exergy analyses of commercial wood pellets ................................99
ix
Table 4.2 Ultimate and proximate analysis of wood pellets ...110
Table C.1 Syngas composition from MATLAB simulation used for general formula
derivation 132
Table D.7 Snapshot of temperature of one typical run in the gasifier ........141
Table E.2 Syngas composition at different biomass flow rate for commercial wood pellets......148
x
List of Figures
Figure 2.2 Effect of moisture content on HHV of syngas under adiabatic condition....60
Figure 2.4 Effect of temperature on number of moles of syngas species (variable ).62
Figure 2.5 Effect of temperature on equivalence ratio in adiabatic condition (variable xg)..62
Figure 2.6 Effect of temperature on syngas species at fixed equivalence ratio of 0.396 ..63
Figure 2.7 Effect of temperature on syngas species at fixed equivalence ratio of 0.396...64
Figure 3.1 Schematic of the Auburn mobile downdraft gasifier designed by the CPC.....78
xi
Figure 3.2 Photograph of the Auburn gasifier along with dimensional sketch of thermocouples
position (Measurements are not to scale; A1 and T1 to T4 represent the position of
thermocouples in the gasifier) . ....79
Figure 3.3 The Auburn mobile downdraft gasifier parked outside the capitol building at
Montgomery, AL...80
Figure 3.8: Ash agglomeration in the grate of gasifier after the gasification of poultry litter...94
Figure 3.9 Biomass flow rate versus product gases for wood pellets .....100
Figure D.1 Temperature recorded by thermocouple at T1 from its start-up to steady state....142
Figure D.2 Temperature recorded by thermocouple at T2 from its start-up to steady state142
Figure D.3 Temperature recorded by thermocouple at T3 from its start-up to steady state ...143
Figure D.4 Temperature recorded by thermocouple at T4 from its start-up to steady state ...143
Figure D.5 Temperature recorded by thermocouple at T1 from its start-up to steady state ...144
xii
CHAPTER 1
LITERATURE REVIEW
In this chapter, various aspects of biomass gasification were reviewed. The most widely
used configurations of biomass gasifiers and the effect of various operating parameters on the
1.1 INTRODUCTION
Biomass has evolved as one of the most promising sources of fuel for the future. This has
spurred the growth of research and development efforts in both federal and private sectors [1].
This impetus is motivated by several factors; dwindling fossil fuels and thus an increase need of
energy security, environmental concerns and promotion of socioeconomic benefits to rural areas.
Another important fact is somewhat uniformly distributed nature of biomass worldwide which
means it is available locally and is helpful in reducing the dependence upon the fossil fuel [2].
The United States target has potential to produce 60 billion gallons of bioethanol per year if all
the available biomass is converted to produce bioethanol [2]. This requires one billion tons of dry
biomass per year which can be reasonably achieved. A recent study shows that the United States
has a potential of 1.3 billion tons of dry biomass per year from its forest and agricultural
resources [2-3]. Studies show the world-wide recoverable residues to be 31 exajoules per year
which is almost equivalent to 10% of the commercial energy use [4-5]. Figure 1.1 shows the
distribution of different renewable energy generation in the United States. The energy derived
1
from biomass is significant and accounts for more than half of all the renewable energy
generation among which, two-thirds is used for heat, power or combined heat and power (CHP)
[6].
The end-use of products from biomass conversion can be mainly in any one of heat and
power applications, transportation fuels (biodiesel, bioethanol) and chemicals for subsequent
processing [7]. Up to present, the only long-term solution for the carbon based fuels and
chemical is biomass and can be effectively converted into solid, liquid and gas [8-9]. Huber et al.
[10] opine that among all the renewable energy sources, biomass is the most optimal long-term
fuel for transportation. Biomass can be converted into biofuels using either thermal or chemical
processes. Among thermal conversion processes, gasification has received the most attention.
This is due to the higher efficiency compared to processes such as direct combustion, pyrolysis
and liquefaction [11-13]. Figure 1.2 shows the power generation from overall gasification
(including coal and biomass) from 1970 to 2004 [14]. This industrial rate of power generation
2
using gasification process can be expected to rise with advances in clean coal technologies and
Different forms of thermal treatment of biomass are distinguished from each other by the
amount of air supplied, residence time, temperature, and consequently the heat transfer rate in
the process. Supplying excess air results in combustion while treatment without air/oxygen
results in pyrolysis products [15]. Gasification is the conversion of biomass into the mixture of
combustible and non-combustible gases (referred as syngas hereafter) by partial oxidation at high
temperature around 800-900C in the presence of a gasifying medium such as air, oxygen or
steam. Syngas from biomass is a mixture of carbon monoxide (CO), carbon dioxide (CO2),
hydrogen (H2), water (H2O) and a small amount of methane (CH4). The use of syngas for power
Warnecke [17] has classified the gasifiers in four categories which are based on the fluid
3
i. Quasi non-moving or self-moving feedstock
a. Downdraft gasifier
b. Updraft gasifier
c. Cross-draft gasifier
c. Entrained-bed gasifier
b. Cyclone gasifier
Among those listed above, downdraft, updraft, BB and CFB gasifiers are the most
common as also shown by studies [18]. Figures 1.3-1.6 show schematics of various gasifiers that
are widely used in the commercial market. Commercially, about 75% of the gasifiers sold are
downdraft gasifiers, 20% fluidized bed, 2.5% updraft, and 2.5% of the other types [18].
The updraft gasifier is popular for application choice when the primary purpose of
gasification is heating only (below 10 MWt) due to its high thermal efficiency and ability to
handle feedstock with wide variation in size and moisture content as high as 50% [19].
Downdraft gasifiers are preferred for small scale power generation due to low amount of tar
content in the syngas. The problem with fixed-bed gasifiers is their inability to maintain uniform
radial temperature which results in local slag, bridging and clinkering problems. Lack of uniform
4
radial temperature is one of the reasons why this kind of gasifier cannot be scaled up rendering
Fluidized bed gasifier provides higher throughput than those with a fixed bed.
Fluidization enhances mass and heat transfer from the fuel thereby increasing heating value of
the output and higher efficiency rendering it excellent for low-rank coal and biomass
gasification. Entrained bed gasification is similar to fluidized bed gasification except for the
operation range temperature which is usually higher than 1900C. This type can have a even
higher throughput capacity but is limited to coal use only due to the particle size constraint on
5
Figure 1.5 Bubbling bed fluidized gasifier Figure 1.6 Circulating fluidized bed gasifier
The various types of gasifiers shown in Figure 1.3-1.6 are discussed in the sections
below. BB and CFG gasifiers are discussed in a single section as fluidized gasifier due to minor
The movement of the feedstock and the gasifying agent are in opposite directions in this
kind of gasifier (also called a counter-current gasifier). Since the syngas formed is not forced to
pass through the hot high temperature zone, the tar content is high in the syngas from this
gasifier. On the other hand, the temperature of syngas exiting from this gasifier is lower around
(200-300C) and hence the thermal efficiency of this kind of gasifier is high. Due to high tar
content in the syngas, a subsequent tar cleaning system is needed, which can become a major
6
1.2.2 DOWNDRAFT GASIFIER
In a downdraft gasifier, the feedstock and gasifying agent both move in the same
direction. The gases have to pass through the high-temperature so amount of tar is significantly
lower than that in an updraft gasifier. The particulate content is however higher for downdraft
gasifier and the thermal efficiency is lower since syngas draws an appreciable amount of energy
In fluidized bed gasifiers, feedstock is fluidized with some bed material like sand/silica
with a gasifying medium which can be air or steam. Fluidized bed gasifiers can further be
classified into two types: bubbling and circulating. Circulating fluidized bed adds one more
feature to bubbling bed such that solid material trapped in the gas phase is trapped and re-
circulated back to the gasification bed. This provides the significant advantages over the
bubbling bed gasifier in terms of mass conversion efficiency and reduces particulate content in
7
Table 1.1 Various characteristics, properties and difference between common types of gasifier ([17, 20-22])
Gasifier Type
Characteristics
Downdraft Updraft BBG CFBG
Less space required due to
enhanced heat transfer resulting
High space requirement for higher throughput due to
Gasifier size in much faster gasification and
modular design of the gasifier and high residence time
lower residence time inside the
gasifier
Not uniform temperature distribution in the radial Uniform temperature
Temperature Profile
distribution distribution inside the gasifier
Permissible particle size/ Size < 5 mm/more sensitive to
sensitivity < 50 mm/good feedstock size
Reaction zone temperature 800-1100C 800-1000C
Ability to handle fine paricles Limited Good
8
8
1.3 CHEMICAL REACTIONS IN THE GASIFICATION PROCESS
gasification sequence as drying and evaporating processes of biomass followed by pyrolysis, and
finally oxidation and reduction [23]. However, the overall process can be reasonably described
(1)
(2)
(3)
(4)
(5)
(6)
(7)
(8)
(9)
Among the reactions described above, the char-oxidation (Eq. 3) and partial-oxidation
(Eq. 4) reactions are slowest, and consequently the rate controlling factor in the overall
gasification process [24]. Pyrolysis also results in liquid which is resistant to the cracking due to
temperature increase though most of the pyrolyzed liquid does so at higher temperature. This
9
requires subsequent cleaning set-up for the tar, which can be a substantial investment in many
cases [23].
Syngas composition varies widely and mostly depends upon the gasifier type, feedstock,
feedstock pre-treatment, gasifying medium and operating parameters like temperature, pressure,
and nature of interaction between reactants in the gasification process [20, 26]. The effects of
major parameters affecting the quality of syngas are discussed in the sections below.
Biomass contains moisture in both ways: intrinsically by its nature, and extrinsically
wherein moisture is absorbed from the surrounding atmosphere [27-28]. Moisture content in the
biomass, during gasification, increases CO2 concentration by the water-shift reaction (Eq. 8)
which consumes CO and liberates H2 [27, 29-30]. While the equilibrium constant for water-shift
reaction varies little over a wide range of temperatures, the direction tends to reverse at higher
temperature. Since more heat is required for moisture evaporation than the small amount of heat
gained due to the exothermic behavior of the water-shift reaction, thermal energy inside the
gasifier reduces when gasifying biomass with higher moisture content [24]. Thus, the decrease in
temperature further exacerbates the scenario and forms more CO2 since the water-shift reaction
is improved at lower temperature. The overall effect is the reduction in calorific value of syngas
because, the small increase in H2 is not sufficient to compensate the loss of significant amount of
CO with increase in moisture content [27, 29-33]. However, the negative effect of moisture
content on the calorific value of syngas is lower at lower equivalence ratio (ER). The ER is the
ratio of actual air fuel ratio to the stoichiometric air fuel ratio which provides the basis for
10
evaluating the amount of air supplied for the gasification with respect to the amount of air
required for the complete combustion of the feedstock. Roy et al. [27] have observed that, in a
downdraft gasifier, when the moisture content is increased from 0 to 40%, heating value of
syngas decreases by 8.72% at ER of 0.45 while the decrease was of 4.7% when the ER used was
0.29. This result was reported from their equilibrium model and thus is applicable to any
gasification process. Table 1.2 summarizes the effect of moisture content in three common
gasifier types.
Maximum limit
Parameter Gasifier type CO CO2 H2 CH4
(% w.b)
Updraft -a +b + ~c <50 [19]
Moisture Content
Downdraft - + + ~ <40 [34]
(M.C)
Fluidized - + + ~ <10 [20]
a b c
decreases with increase in M.C, increases with increase in M.C, no significant change
literature beyond which self-sustaining gasification is not possible due to an enthalpy deficiency
for vaporization. In fact, supplemental fuel is required for most of the combustor when the
moisture content is greater than 50% on a wet basis [35-36]. Moisture content up to 30% (wet
basis) can be used for downdraft gasifier [21, 34]. When air is used as the gasification agent, the
amount of methane produced is small and stays almost constant with change in moisture content
[29, 37]. Thus the temperature decrease inside the gasifier due to moisture also results lower
mass conversion efficiency and increases tar content [30, 38-40]. Sheth et al. [41] report the
decrease in biomass consumption rate with increase in moisture content which is due to the
higher amount of heat necessary for drying those wood chips inside the reactor before they can
11
be pyrolyzed. However, some moisture content is always desirable since it enhances steam
reforming and helps to crack tar, and at higher temperature, also enhances other reactions such as
char gasification [42-43]. Steam injection is widely used in industrial applications to adjust
syngas composition in the gasification process but often, in the presence of higher temperature
Equivalence ratio (ER) is the most influential parameter in any gasification process and
often has significant impact on syngas composition. Increase in ER increases the temperature
inside the gasifier while ER decrease increases char formation inside the gasifier. As can be
seen from Figure 1.7, all combustible products reduce with an increase in ER with the formation
of higher amount of CO2 as well as total gas yield greatly diminishing the heating value of the
final syngas [45-48]. Zainal et al. [49] compared the best optimal value for the downdraft
gasifier with respect to equivalence ratio using furniture wood and wood chips as feedstock. The
effect of equivalence ratio for each syngas component was analyzed with the conclusion of an
optimal equivalence ratio of 0.38 for the gasifier performance for that particular feedstock. At
this equivalence ratio, CO, CH4 and calorific value each attain their maximum outputs while CO2
12
Figure 1.7 Effect of ER upon syngas composition in the reduction zone [39]
Heat conduction inside the gasifier is often limited by the thermal diffusivity of charcoal
particles- 0.7 10-7 m2/s which is lower than that of wood [50]. The optimum equivalence ratio is
necessary for accelerating pyrolysis and drying rate due to conduction and convection process
Both Skoulou et al. [51] and Sheth et al. [41] report an optimal equivalence ratio of 0.2
for downdraft gasification of olive kernels and olive tree cutting and furniture wood. The
optimum equivalence ratio varies for different biomass due to the amount of oxygen elementally
present in the biomass as well as the ash content. For example, coal requires far more oxygen
than common biomass materials for gasification due to its lower oxygen content [40]. The
existing literature shows that equivalence ratio should be around 0.3-0.4 for the successful
gasification. Table 1.3 shows the optimal equivalence ratio for selected feedstocks. Also from
13
Table 1.3, difference in optimal ER can also be observed for the feedstock with same elemental
composition (pine wood chips and saw dust) in fluidized bed. This is due to the difference in
gasification temperature which was lower for pine saw dust (780-830C) than pine wood chips
(>900C). Optimal equivalence ratio for an updraft gasifier is not shown in the table due to the
Table 1.3 Optimal ER for some feedstocks in downdraft and fluidized gasifiers
Gasifier
Feedstock Optimal ER References
type
Furniture wood + charcoal 0.38 [49]
Olive kernels and olive tree
0.2 [51]
Downdraft cutting
Hazelnut shells 0.28 [34]
Any compound with molecular weight greater than that of benzene is called tar [54]. Tar
concentration decreases with increase in ER. This is mainly due to two reasons: (a) higher
temperature as a result of higher ER increases reaction rates of the chemical products; and (b)
high ER supplies additional oxygen for cracking of tar into lower hydrocarbons, CO2 and H2O.
Thus, at some point between the applicable ranges (0.15-0.4), a shift between types of tar is also
reported. Light tar increases while heavy tar decreases [46, 55]. Corella et al. [56] suggests an
equivalence ratio above 0.36 for pine wood in a fluidized bed to reduce the tar content below 2
g/m3, a level below which coke formation does not exceeds its removal rate.
14
The effect of superficial velocity (SV) is worthy of discussion under the topic of
equivalence ratio due to its direct relation with the amount of gasification/fluidization medium
injected inside the gasifier. The SV is the ratio of volume flow rate of syngas to the cross-
sectional area of the gasifier and can be thought as one independent parameter unconstrained to a
particular gasifier size. Higher SV promotes burning as well as reaction rate and decreases the
residence time of biomass in the system [57]. Higher burning rate increases the temperature of
the gasifier. Yamazaki et.al [58] recommended SV greater than 0.4 m/s for syngas appropriate
for internal combustion engines. An overall increase in combustibles (except CH4 which shows
no appreciable change) is reported with increase in SV. An initial decrease reaching the
minimum level followed by an increase is reported with SV, the optimum SV being 0.4 m/s.
Increase tar after the increase beyond optimum SV is due to the short residence time of the tar
1.4.3 TEMPERATURE
Increase in temperature reduces the tar content as well as decreases char inside the
gasifier [51, 59]. Gas yield increases due to higher tar cracking. One of the means of increasing
temperature is by internal recirculation of syngas [60]. Tar cracking temperatures are often
reported to be around 1000-1100C with some dependency on gasifier design [34, 54]. Other
methods of tar cracking are also employed such as multi-stage gasifiers [61-63]. Multi-stage
gasifiers, as shown in Figure 1.8, have separated pyrolysis and gasification zones and make use
of partial oxidation of pyrolysis gas obtained in the pyrolysis zone for tar cracking and thus, tar
content can be reduced as low as 15mg/Nm3 (Nm3 volume at STP) [61]. CO content increases
with increase in temperature because endothermic reactions are more favored at higher
temperature [64]. Mass conversion efficiency decreases with decrease in temperature [60]. An
15
oxidation zone below a temperature of 725C gives significantly lower mass conversion
efficiency [65].
very important for efficient mixing in a fluidized bed. Generally, less than 100C difference in
total riser height is acceptable [43]. Heating value as well as syngas yield is found to increase
operating temperature driven by an external supply of heat in the gasifier for constant ER [43,
45, 47, 66-67]. However, this is different when the temperature is increased due to increase in
ER inside the reactor which actually reduces the combustibles [47]. Temperature control cannot
16
be independent in any gasification process and is an output variable, with the exception of small
lab-scale or pilot plants which can be heated with external heat. The temperature of the reactor is
dependent on various factors such as moisture content of the fuel, ER, heat losses from the
system, and amount of steam added [56, 68-69]. Thus, the temperature inside the gasifier should
represent an optimal compromise with ER. The best approach is the proper insulation of the
reactor and using waste heat. Higher temperature also reduces tar content significantly due to
thermal cracking [43, 45, 47, 67, 70-72]. In addition, Cao et al. [70] report higher reduction in tar
with same increase in top part of reactor than in bottom section. However, Drift et al. [55]
suggest that the tar that is cracked due to temperature is mostly the heavy tar while light tar is not
decomposed. Heavy tars are the product of pyrolysis process which has not gone through
cracking while light tars are the cracking products of heavy tar. In certain cases, light tar seems
to increase due to the subsequent breakdown of heavy tar into light tar and other compounds.
Typical temperatures suggested for biomass gasification in a fluidized bed are around 800-900C
by various studies [67, 72-73]. Although, high temperature increases carbon conversion
efficiency of the overall gasification system, consideration should be given to prevent the
formation of ash-melt, made not to form ash-melting, especially when used to gasify biomass
material with high ash content like rice husk [55]. Seggiani [74] has developed the empirical
relations that can be used to predict ash-fusion temperature of biomass-ash based on its elemental
composition under reducing conditions. Eq. (10) shows the general form of the relation.
(10)
In the above relation, to are the coefficients for calculating ash fusion temperature
17
1.4.4 BIOMASS TYPE
ratio and increases when these ratios increase, especially with an increase in Hydrogen/Carbon
ratio [35]. A higher oxygen concentration in biomass needs lower ER for gasification because of
its inherent oxygen that will also be available for gasification [40].
Table 1.4 Ash content and its elemental composition for some common feedstocks
(% dry basis)
Another important factor is the ash content of the feedstock. Table 1.4 provides the ash
content (% dry basis) and the elemental composition of various common biomass feedstocks [35,
75]. Although formation of clinkers can cause problems for the gasifier operation with biomass
having ash-content above 5%, successful gasification with ash-content up to 25% is reported
[22, 76]. Higher ash content causes slagging, and consequently ash agglomeration due to fusion,
the rate of which is dependent upon the ash content in biomass and ash composition [35, 77-78].
Thus high ash content biomass should be gasified at the temperature below the oxidation or
reducing temperature of the minerals constituents in the ash, often which is not possible if the
constituents have relatively low ash-fusion temperature [54, 62]. Common ash minerals in
18
biomass are silica, potassium, calcium, aluminum, magnesium, iron, sodium and chlorine. These
minerals present in biomass can exist as salts and vaporize during the gasification process
contaminating the syngas. Also, it is highly possible for these minerals to react with silicon in
presence of oxygen to create low-temperature melting silicates which can create a severe
deposition problem. Alkali metals such as potassium and calcium silicates have melting
temperatures even below 7000C [35]. One other way to tackle the problem is to resort to some
kind of removal process like leaching for alkali metal removal which has been reported to reduce
these minerals by more than 80%. Removal of these alkali metals will increase the ash fusion
The presence of ash in biomass requires careful control over the operating temperature.
Neither should it be high enough to fuse minerals in the ash forming a barrier to further
gasification by formation of clinkers, nor too low leading to unburnt carbon resulting in lower
Fixed bed gasifiers have lower biomass feedstock size restrictions compared to fluidized
bed gasifiers. Usually, feed size less than 51 mm and 6 mm is recommended for fixed bed and
fluidized bed, respectively [20]. Use of larger size feedstock has been tried and reported by
several authors [79-81]. Saravanakumar et al. [80] have successfully gasified long sticks with
particle size suggested for a conventional downdraft gasifier with throated design is one-eighth
of the reactor throat diameter [82]. The larger particles form bridges preventing the efficient flow
19
of biomass inside a gasifier while smaller particles interferes with the air/gasifying agent passage
creating high pressure drop and consequently can result in gasifier shut-down [22].
Sharma [83] reports increase in the temperature of oxidation and reduction zone with
decrease in particle size of the biomass feedstock in a downdraft gasifier. Decrease in particle
size reduces the heat loss due to radiation and enhances the thermal conductivity in the oxidation
and reduction zones. On the other hand, decrease in particle size increases pressure drop inside
the gasifier. Burning rate and thus the char oxidation period of fuel particles decrease with
increase in bulk density and particle size [33, 77]. Biomass consumption rate is inversely related
to particle size [57]. In other words, higher residence time is recommended for larger biomass
particle size. Decrease in CO with increase in CO2 concentration is observed. Ryu et al. [84]
report decrease from 18% to 13.5% CO when the size of wood cubes used in the experiments
were increased from 10 mm to 35 mm. Their model predicts a decrease in CH4 and an increase in
H2 with increase in size of biomass particles. Also, the temperature gradient decreases thus
increasing time taken for diffusion of heat. This will result in poor temperature distribution
which is also one of the reasons for the increase in CO2 concentration with increase in particle
size.
Carbon conversion efficiency is not strongly affected by particle size except the lower
biomass size increases tar concentration because of high entrainment susceptibility during
fluidization [85]. This is because particles can be easily transported to the upper part of the
reactor, leaving little time for tar cracking. Multi-staging can prevent this as demonstrated the
novel concept developed by Kersten et.al [86] using a gasifier design consisting of several cone
shaped structures welded together with the base of each connected to the next tubes of equal
diameters. The design, as shown in Figure 1.9, makes it possible to maintain numerous fluidized
20
sections in one reactor, and thus effectively control back-mixing of solids and gases. On the
other hand, the axial temperature drop increases significantly with decrease in size. This is due to
the easy passage of feed particles from the feed point and thus little or no reaction taking place
below the feed point. Thus, the homogeneity of the bed material cannot be maintained
throughout the reactor [55, 69]. Wiman and Almstedt [87] report increase in gas-particle
1.4.6 PRESSURE
High-pressure gasification reduces the size of the reactor for the same amount of
feedstock and can act to reduce the need for further compression when the gasification products
are intended for subsequent use in Fischer-Tropsch process or other chemical synthesis which
requires high pressure [88]. Pressure drop across the gasifier increases with smaller particle size
21
due to increased porosity [77]. The general recommendation of biomass particles size for various
Increase in pressure in a fluidized bed increases turbulence and thus increase in gas-
particle interaction is observed. Increase in pressure also results in bubble instability and bubble
splitting in fluidized bed. Wiman and Almstedt [87] have defined a parameter called bed
Their finding shows a significant increase in with increase in pressure. However, the
rate of increase drops with increase in pressure and levels off once the pressure reaches around 1
MPa. Among the two particles size of 0.45 mm and 0.70 mm selected for their experiment, is
Valin et al. [88] have studied the effect of pressure upon syngas composition with
pressure from 2 to 10 bar in fluidized bed with wood sawdust as a feedstock. With increase in
pressure, an increase in CO2, CH4 and H2 were observed, while CO decreased. In their lab-scale
reactor using steam and N2 as the gasification medium, with increase in pressure from 2 to 10
bars, an increase of 16%, 53% and 38% of H2, CO2 and CH4, respectively and decrease of CO by
33% was reported. Overall, increase in dry gas yield increase is reported to be 20% with increase
22
in pressure from 2 to 7 bars after which the gas yield remains constant. The increase in various
gases and total gas yield is due to the increase in char hold-up rate which increases catalytic
Air
Oxygen
Steam
Gasification with air results in syngas with low higher heating value (HHV) due to
inherent dilution with N2 present in the air. Conversely, gasification with oxygen yields syngas
with a heating value of 10-12 MJ/Nm3 and steam gasification results in syngas with heating
value even higher, 15-20 MJ/Nm3 [23]. Air gasification is widely used compared to oxygen and
Secondary air reduces the tar content in the syngas due to partial combustion of syngas.
This in turn establishes local high temperature zone due to exothermic reactions. Thermal
cracking of tar is thus due to these high temperature zones in the periphery of the secondary air.
Since heat is obtained from the combustion, HHV decreases with increase in secondary air [43,
70-71]. Pan et al. [71] suggest 20% secondary air (% of the primary air) for the minimum tar in
23
1.4.8 GASIFICATION OF WASTES AND BIOMASS CO-GASIFICATION
Bacaicoa et al. [90] have studied the co-gasification of a biomass and high density
decreased with increase in HDPE fraction in the mixture. Since HDPE vaporizes instantly at high
lignocellulose present in the mixture. A higher fraction of HDPE in the mixture consumes more
oxygen from the air supplied and thus leaves less air for lignocellulose to gasify consequently
reducing the biomass flow rate. They also report increase in CH4 and CO concentration with
increase in HDPE fraction, while H2 varies only slightly and CO2 decreases. This is attibuted in
other study to the strong affinity of HDPE towards a Boudouard reaction (Eq. 5) compared to a
Research has also been conducted with co-gasification biomass and coal [91-96]. Pan et
al. [97] have gasified poor quality coal (carbon content less than 38%) and biomass (pine wood
chips) in a CFB gasifier using air/steam as the gasifying agent. They have done experiments with
varying biomass/coal ratio from 0 to 1 at the interval of 0.25. Low temperature resulted when the
proportion of pine chips was increased in the blend due to increase in endothermic reaction rate
between wood charcoal of high reactivity and that of steam. The poor quality coal, when used
alone was not able to sustain gasification and only produced flue gas (gas consisting mostly CO2,
water vapors and nitrogen from the air). Increased combustibles, gas yield and increase in
heating value were reported with increase in the proportion of pine chips.
(PET) and brown coal was examined in a bench-scale FB with 10% O2 and 90% N2 as the
24
gasifying agent. The blend consists of 23% PET and 77% brown coal since difficulties in
gasification were reported when the proportion of PET was increased beyond 23%. The effect of
bed temperature and freeboard temperature was reported. Increase in CO and H2 was found
while CH4 showed slight decrease with increase in fluidized bed temperature. On the contrary,
fluidized bed temperature. It is due to the fact that all reactions will be completed in close
Xiao et al. [99] have studied plastic waste gasification in a FB gasifier with air as the
gasifying medium. Higher ER led to an increase in temperature in the reactor. The temperature
distribution inside the reactor showed gradual decrease of temperature from bed to freeboard.
Temperature difference reported by them is 80-100C. Gas yield also increased with ER, while
char and tar concentration decreased. The effect of bed height is also examined and an optimum
bed height suggested ensuring long residence time to facilitate the effective cracking of heavy tar
and char gasification. CO and H2 increased initially with an increase in ER due to efficient
thermal cracking at higher temperature obtained by higher ER and then subsequently decreased
with further increase in ER due to combustion of these products. As expected, HHV decreased
with increase in ER. The effect of bed height was also examined on syngas composition.
Velez et al. [100] report the co-gasification results with sub-bituminous coal and biomass
(sawdust, rice husk, coffee husk) in a fluidized bed gasifier with air/steam mixture with the
objective of finding the optimum proportion of biomass/coal yielding highest heating value.
Tests were run at 6%-15% of biomass beyond which proper fluidization was not possible due to
25
density difference in two fuel types resulting in fuel segregation upon gasification. Reactor
temperature decreases with increase in biomass concentration due to the lower heating value of
biomass compared to that of coal. Increase in H2 and CO was obtained with increase in
steam/mixture ratio. Ash agglomeration and sintering was also reported in their study.
Mclendon et al. [101] report lower gas yield from co-gasification of coal and biomass
with respect to biomass gasified alone. Another interesting research on waste gasification is the
gasification of waste tires. Similar effect of various parameters on syngas composition and
product yield, consistent with the above discussion were reported from the gasification studies
Recently, co-firing of biomass in coal-fired plants has also emerged as one of the biomass
utilization options. Though biomass has higher reactivity than coal and can provide some
improvements in overall coal gasification process, there are several problems associated with co-
firing of biomass in conventional coal-power plants [103]. The major issues associated with
High moisture content in biomass entails the need for prior drying before using into coal-
powered plant.
Biomass has low bulk density than coal and thus, might require additional handling
The ash in biomass has low melting point than those compared to coal. In addition,
biomass-ash is aggressive in nature and might corrode the gasifier and associated gas
supply system.
26
Biomass is hydrophilic and a non-friable which makes it very difficult in deal in
Impacts of these problems depend upon various factors such as coal/biomass ratio, actual
Proper consideration of bed material in a fluidized bed is important for achieving proper
homogenization of feed particles and efficient heat transfer so that minimum temperature
gradient is realized within the riser. In many cases, bed material can itself act as a catalyst
facilitating efficient tar cracking [46, 56]. Skoulou et al. [46] compared the effect of olivine over
silica sand, latter of which is reported to have adverse effect upon effective fluidization due to
agglomeration and tar formation when operating at the temperature below 8000C. Pfeifer et al.
[104] have studied in-bed catalysts in a dual bed fluidized bed reactor with Ni/olivine as the
catalyst and observed significant tar reduction. Use of catalyst for tar cracking is itself a vast
subject and further discussion is avoided here to remain within the scope of the review. The
Although biomass gasification is not a new concept by itself, current energy scenario and
significant interests in renewable energy has spurred the industrial and academic research in this
field. Various configurations of biomass gasifiers have been studied to achieve the maximum
efficiency from the process. An exhaustive amount of literature can be found in this field. The
27
Proper utilization of biomass through gasification can increases the energy security
Drying is usually cost-intensive. Utilization of waste heat to dry biomass can be very
helpful.
The equivalence ratio plays important role in determining the overall syngas quality.
While using air as the gasifying agent, high amount of sensible energy is lost in
heating the nitrogen from air. Although steam or oxygen gasification is possible, the
cost associated with the process make them economically unfeasible. Identifying and
Tar content has remained as one of the major issues in biomass gasification.
Although primary or secondary tar treatment can be done to reduce the level of tar
from the biomass gasifier, costs associated with the process might be considerable.
in turn, reduces quality of the syngas. Preventing heat losses from the gasifier by
proper insulation can reduce air need to maintain the sustainable gasification
temperature.
The type of biomass affects significantly in the overall syngas composition and
sometimes, also in the operational issues in the biomass gasification plants. High ash
28
content material is not desirable. However, many pre-treatment processes exists that
can be used to cure the biomass before feeding into the reactor.
High pressure gasification is very significant in decreasing the overall reactor size
and increasing the quality of syngas from the gasifier. However, costs and
Fluidized beds offer an excellent advantage over fixed bed gasifier in terms of
scalability. However, the constraints on particle size and moisture content often
make it unsuitable at some cases. On other hand, fixed beds are suited more for
small-scale application. Particle size constraints can impose enough restriction due to
though the focus of this new power plants may be more efficient utilization of coal
rather than biomass, current status of global warming and public awareness in
utilizing renewable energy, may create ample opportunities for biomass co-firing.
To develop a model that can predict syngas composition for wide variety of feedstocks
To study the selected operating parameters and their effect on syngas composition in a
29
To conduct detail study on tar concentration in the syngas from a stratified downdraft
gasifier
1.6 REFERENCES
[1] H.L. Chum, R.P. Overend, Biomass and renewable fuels, Fuel Processing Technology, 71
(2001) 187-195.
[2] R.D. Perlack, L.L. Wright, A.F. Turhollow, R.L. Graham, B.J. Stokes, D.C. Erbach, Biomass
as feedstock for a bioenergy and bioproducts industry: the technical feasibility of a billion-ton
[3] A. Demirbas, Biofuels sources, biofuel policy, biofuel economy and global biofuel
[4] S.P. Babu, Thermal gasification of biomass technology developments: End of task report for
[5] A. Demirbas, Biomass resource facilities and biomass conversion processing for fuels and
[6] D. Peterson, S. Haase, Market assessment of biomass gasification and combustion technology
for small- and medium-scale applications, in, National Renewable Energy Laboratory, Golden,
CO, 2009.
[7] K.B. Cantrell, T. Ducey, K.S. Ro, P.G. Hunt, Livestock waste-to-bioenergy generation
[8] H. Boerrigter, R. Rauch, Syngas production and utilization; Review of applications of gases
from biomass gasification, in: H.A.M. Knoef (Ed.) Handbook biomass gasification, Biomass
30
[9] New renewable energy resources : a guide to the future, Kogan Page, London, 1994.
[10] G.W. Huber, S. Iborra, A. Corma, Synthesis of Transportation Fuels from Biomass:
[11] A.V. Bridgwater, Renewable fuels and chemicals by thermal processing of biomass,
[12] A. Demirbas, Combustion characteristics of different biomass fuels, Progress in Energy and
[13] G.J. Stiegel, R.C. Maxwell, Gasification technologies: the path to clean, affordable energy
[14] D.R. Simbeck, Report on SFA Pacific gasification database and world market report, in:
[15] A. Demirbas, Biofuels securing the planet's future energy needs, Energy Conversion and
[16] A. Demirbas, Progress and recent trends in biofuels, Progress in Energy and Combustion
[17] R. Warnecke, Gasification of biomass: comparison of fixed bed and fluidized bed gasifier,
[18] H.A.M. Knoef, Inventory of Biomass Gasifier Manufacturers and Installations, Final Report
31
[20] P. Basu, Combustion and gasification in fluidized beds, Taylor and Francis Group, LLC,
[21] M. Dogru, A. Midilli, C.R. Howarth, Gasification of sewage sludge using a throated
downdraft gasifier and uncertainty analysis, Fuel Processing Technology, 75 (2002) 55-82.
[22] P. McKendry, Energy production from biomass (part 3): gasification technologies,
[23] A.V. Bridgewater, Renewable fuels and chemicals by thermal processing of biomass,
[24] S. Lee, Alternative fuels, 1st ed., Taylor & Francis, Washington, (1996). pp.134-44.
[25] R.C. Brown, Biorenewable resources-engineering new products from agriculture, 1st ed.,
[26] J.M. Prins, Thermodynamic analysis of biomass gasification and torrefaction, in, Eindhoven
[27] P.C. Roy, A. Datta, N. Chakraborty, Modelling of a downdraft biomass gasifier with finite
rate kinetics in the reduction zone, International Journal of Energy Research, 33 (2009) 833-851.
[28] P. McKendry, Energy production from biomass (part 1): overview of biomass, Bioresource
[29] Z.A. Zainal, R. Ali, C.H. Lean, K.N. Seetharamu, Prediction of performance of a downdraft
gasifier using equilibrium modeling for different biomass materials, Energy Conversion and
32
[31] A.K. Sharma, Equilibrium modeling of global reduction reactions for a downdraft (biomass)
[32] C.R. Altafini, P.R. Wander, R.M. Barretoa, Prediction of the working parameters of a wood
waste gasifier through an equilibrium model Energy Conversion and Management, 44 (2003)
2763-2777.
[33] D. Shin, S. Choi, The combustion of simulated waste particles in a fixed bed, Combustion
[34] M. Dogru, C.R. Howarth, G. Akay, B. Keskinler, A.A. Malik, Gasification of hazelnut
[35] B.M. Jenkins, L.L. Baxter, T.R.M. Jr., T.R. Miles, Combustion properties of biomass, Fuel
[36] T.B. Reed, A. Das, Handbook of biomass downdraft gasifier engine systems, SERI, Golden,
CO, 1988.
[38] T.H. Jayah, L. Aye, R.J. Fuller, D.F. Stewart, Computer simulation of a downdraft wood
[39] J.K. Ratnadhariya, S.A. Channiwala, Three zone equilibrium and kinetic free modeling of
[40] K.J. Ptasinski, M.J. Prins, A. Pierik, Exergetic evaluation of biomass gasification, Energy,
32 (2007) 568-574.
[41] P.N. Sheth, B.V. Babu, Experimental studies on producer gas generation from wood waste
33
[42] I. Narvaez, A. Orio, M.P. Aznar, J. Corella, Biomass gasification with air in an atmospheric
bubbling fluidized bed. effect of six operational variables on the quality of the produced raw gas,
[43] X.T. Li, J.R. Grace, C.J. Lim, A.P. Watkinson, H.P. Chen, J.R. Kim, Biomass gasification
[44] X. Li, J.R. Grace, A.P. Watkinson, C.J. Lim, A. Ergdenler, Equilibrium modeling of
gasification: a free energy minimization approach and its application to a circulating fluidized
Biomass Gasification in a Fluidized Bed Reactor, Energy & Fuels, 22 (2008) 3493-3498.
olive kernels in a 5-kW fluidized bed reactor for H2-rich producer gas, International journal of
[47] P.J. van den Enden, E.S. Lora, Design approach for a biomass fed fluidized bed gasifier
using the simulation software CSFB, Biomass and Bioenergy, 26 (2004) 281-287.
[48] A. van der Drift, J. van Doorn, J.W. Vermeulen, Ten residual biomass fuels for circulating
[49] Z.A. Zainal, A. Rifau, G.A. Quadir, K.N. Seetharamu, Experimental investigation of a
[50] J. Larfeldt, B. Leckner, M.C. Melaaen, Modelling and measurements of heat transfer in
charcoal from pyrolysis of large wood particles, Biomass and Bioenergy, 18 (2000) 507-514.
34
[51] V. Skoulou, A. Zabaniotou, G. Stavropoulos, G. Sakelaropoulos, Syngas production from
olive tree cuttings and olive kernels in a downdraft fixed-bed gasifier, International Journal of
[52] E. Natarajan, A. Nordin, A.N. Rao, Overview of combustion and gasification of rice husk in
[53] J. Gil, J. Corella, M.P. Aznar, M.A. Caballero, Biomass gasification in atmospheric and
bubbling fluidized bed: Effect of the type of gasifying agent on the product distribution, Biomass
[54] T.A. Milne, R.J. Evan, Biomass gasification "tars"; their nature, formation and conversion,
[55] A.v.d. Drift, J.v. Doorn, Effect of fuel size and process temperature on fuel gas quality from
CFB gasification of biomass, in: A.V. Bridgwater (Ed.) Progress in Thermochemical Biomass
[56] J. Corella, J.M. Toledo, G. Molina, Calculation of the conditions to get less than 2 g tar/mn3
[57] F.V. Tinaut, A. Melgar, J.F. Prez, A. Horrillo, Effect of biomass particle size and air
superficial velocity on tar from downdraft gasification of biomass, Energy & Fuels, 19 (2005)
1186-1191.
[59] M.M. Kk, A. Demirbas, Biomass conversion processes, Energy Conversion and
35
[60] P.R. Wander, C.R. Altafini, R.M. Barreto, Assessment of a small sawdust gasification unit,
[61] P. Brandt, E. Larsen, U. Henriksen, High tar Reduction in a two-stage gasifier, Energy &
(2008) 1383-1393.
[63] S. Monteiro Nunes, N. Paterson, D.R. Dugwell, R. Kandiyoti, Tar formation and destruction
in a simulated downdraft, fixed-bed gasifier: reactor design and initial results, Energy & Fuels,
21 (2007) 3028-3035.
[64] A. Zabaniotou, O. Ioannidou, V. Skoulou, Rapeseed residues utilization for energy and 2nd
[65] A. Rogel, J. Aguilln, The 2D eulerian approach of entrained flow and temperature in a
biomass stratified downdraft gasifier, American Journal of Applied Sciences, 3 (2006) 2068-
2075.
materials on fluidized bed steam gasification of wood biomass, Bioresource Technology, 100
(2009) 1419-1427.
[67] C. Wu, X. Yin, L. Ma, Z. Zhou, H. Chen, Design and operation of a 5.5 MWe biomass
integrated gasification combined cycle demonstration plant, Energy & Fuels, 22 (2008) 4259-
4264.
[68] J. Corella, A. Sanz, Modeling circulating fluidized bed biomass gasifiers. A pseudo-rigorous
36
[69] S.R.A. Kersten, W. Prins, A. van der Drift, W.P.M. van Swaaij, Experimental fact-finding
in CFB biomass gasification for ECN's 500 kWth pilot plant, Industrial & Engineering
[70] Y. Cao, Y. Wang, J.T. Riley, W.-P. Pan, A novel biomass air gasification process for
producing tar-free higher heating value fuel gas, Fuel Processing Technology, 87 (2006) 343-
353.
[71] Y.G. Pan, X. Roca, E. Velo, L. Puigjaner, Removal of tar by secondary air in fluidised bed
[72] Z. Wu, C. Wu, H. Huang, S. Zheng, X. Dai, Test results and operation performance analysis
of a 1-MW biomass gasification electric power generation system, Energy & Fuels, 17 (2003)
619-624.
[73] R.C. Saxena, D. Seal, S. Kumar, H.B. Goyal, Thermo-chemical routes for hydrogen rich gas
from biomass: A review, Renewable and Sustainable Energy Reviews, 12 (2008) 1909-1927.
[74] M. Seggiani, Empirical correlations of the ash fusion temperatures and temperature of
critical viscosity for coal and biomass ashes, Fuel, 78 (1999) 1121-1125.
[75] M.J. Fernandez Llorente, J.E. Carrasco Garca, Comparing methods for predicting the
[76] X.L. Yin, C.Z. Wu, S.P. Zheng, Y. Chen, Design and operation of a CFB gasification and
power generation system for rice husk, Biomass and Bioenergy, 23 (2002) 181-187.
[77] C. Ryu, Y.B. Yang, A. Khor, N.E. Yates, V.N. Sharifi, J. Swithenbank, Effect of fuel
properties on biomass combustion: Part I. Experiments--fuel type, equivalence ratio and particle
37
[78] W.R. Livingston, Biomass ash characteristics and behaviour in combustion, gasification and
[79] K.M. Bryden, K.W. Ragland, Numerical modeling of a deep, fixed bed combustor, Energy
[80] A. Saravanakumar, T.M. Haridasan, T.B. Reed, R.K. Bai, Experimental investigation and
modelling study of long stick wood gasification in a top lit updraft fixed bed gasifier, Fuel, 86
(2007) 2846-2856.
[81] R. Bilbao, A. Millera, M.B. Murillo, Temperature profiles and weight loss in the thermal
decomposition of large spherical wood particles, Industrial & Engineering Chemistry Research,
32 (1993) 1811-1817.
[82] D.M. Earp, The gasification of biomass in downdraft reactor, in, Aston University, UK,
1988.
[83] A.K. Sharma, Modeling fluid and heat transport in the reactive, porous bed of downdraft
(biomass) gasifier, International Journal of Heat and Fluid Flow, 28 (2007) 1518-1530.
[84] Y.B. Yang, V.N. Sharifi, J. Swithenbank, Effect of air flow rate and fuel moisture on the
burning behaviours of biomass and simulated municipal solid wastes in packed beds, Fuel, 83
(2004) 1553-1562.
[85] D.Y.C. Leung, C.L. Wang, Fluidized-bed gasification of waste tire powders, Fuel
[86] S.R.A. Kersten, W. Prins, B. van der Drift, W.P.M. van Swaaij, Principles of a novel
multistage circulating fluidized bed reactor for biomass gasification, Chemical Engineering
Science, 58 725-731.
38
[87] J. Wiman, A.E. Almstedt, Influence of pressure, fluidization velocity and particle size on the
2167-2176.
[88] S. Valin, S. Ravel, J. Guillaudeau, S. Thiery, Comprehensive study of the influence of total
pressure on products yields in fluidized bed gasification of wood sawdust, Fuel Processing
[89] A.V. Bridgwater, The technical and economic feasibility of biomass gasification for power
(2008) 5485-5491.
[91] F. Pinto, C. Franco, R.N. Andr, C. Tavares, M. Dias, I. Gulyurtlu, I. Cabrita, Effect of
experimental conditions on co-gasification of coal, biomass and plastics wastes with air/steam
[92] D.R. McIlveen-Wright, F. Pinto, L. Armesto, M.A. Caballero, M.P. Aznar, A. Cabanillas,
combustion and gasification power plant technologies for processing mixtures of coal, biomass
[93] M.L. Mastellone, L. Zaccariello, U. Arena, Co-gasification of coal, plastic waste and wood
[94] M.P. Aznar, M.A. Caballero, J.A. Sancho, E. Francs, Plastic waste elimination by co-
gasification with coal and biomass in fluidized bed with air in pilot plant, Fuel Processing
39
[95] J. Fermoso, B. Arias, M.G. Plaza, C. Pevida, F. Rubiera, J.J. Pis, F. Garca-Pea, P. Casero,
High-pressure co-gasification of coal with biomass and petroleum coke, Fuel Processing
Technology, 90 926-932.
[96] M. Lapuerta, J.J. Hernndez, A. Pazo, J. Lpez, Gasification and co-gasification of biomass
wastes: Effect of the biomass origin and the gasifier operating conditions, Fuel Processing
[97] Y.G. Pan, E. Velo, X. Roca, J.J. Many, L. Puigjaner, Fluidized-bed co-gasification of
residual biomass/poor coal blends for fuel gas production, Fuel, 79 (2000) 1317-1326.
Koutsk, K. Svoboda, Gasification of coal and PET in fluidized bed reactor, Fuel, 85 (2006)
2458-2468.
[99] R. Xiao, B. Jin, H. Zhou, Z. Zhong, M. Zhang, Air gasification of polypropylene plastic
waste in fluidized bed gasifier, Energy Conversion and Management, 48 (2007) 778-786.
[100] J.F. Vlez, F. Chejne, C.F. Valds, E.J. Emery, C.A. Londoo, Co-gasification of
Colombian coal and biomass in fluidized bed: An experimental study, Fuel, 88 (2009) 424-430.
[101] T.R. McLendon, A.P. Lui, R.L. Pineault, S.K. Beer, S.W. Richardson, High-pressure co-
gasification of coal and biomass in a fluidized bed, Biomass and Bioenergy, 26 (2004) 377-388.
[102] G. Xiao, M.-J. Ni, Y. Chi, K.-F. Cen, Low-temperature gasification of waste tire in a
[103] A. Maciejewska, H. Veringa, J. Sanders, S.D. Peteves, Co-firing of biomass with coal:
constraints and role of biomass pre-treatment, in, Institute for Energy, Netherlands, 2006.
[104] C. Pfeifer, R. Rauch, H. Hofbauer, In-Bed Catalytic Tar Reduction in a Dual Fluidized Bed
Biomass Steam Gasifier, Industrial & Engineering Chemistry Research, 43 (2004) 1634-1640.
40
[105] Z. Abu El-Rub, E.A. Bramer, G. Brem, Review of Catalysts for Tar Elimination in
Biomass Gasification Processes, Industrial & Engineering Chemistry Research, 43 (2004) 6911-
6919.
[106] L. Devi, K.J. Ptasinski, F.J.J.G. Janssen, A review of the primary measures for tar
41
CHAPTER 2
MODELING
2.1 ABSTRACT
Carbon monoxide (CO), carbon dioxide (CO2), hydrogen (H2) and methane (CH4) are the
major gases produced from biomass gasification. The composition of CO, CO2 and H2 in syngas
from the biomass gasification process was calculated via equilibrium modeling. Methane
concentration predicted by the equilibrium model was almost negligible (<0.15 vol. %) at
temperatures above 800C. Nearly one-hundred biomass samples were used to calculate synthesis
gas composition and the generalized equations were obtained by multiple regression analysis to
predict synthesis gas composition using elemental analysis of biomass. Equilibrium results were
compared with the experimental data. Effect of temperature and moisture content on synthesis
gas composition is also presented. Although perfect chemical equilibrium conditions cannot be
achieved in an actual gasification process, the derived formula generally predicts the syngas
42
2.2 INTRODUCTION
for meeting future energy demand. In addition, biomass is the only source of carbon-based
renewable fuels, and the proper and sustainable exploitation of this resource is essential to secure
the United States energy security. Among various biomass conversion technologies within
process for syngas production. Biomass gasification has received the highest interest among
various biomass conversion technologies because it is almost feedstock-agnostic and can be used
to produce electricity and liquid fuels such as green gasoline and diesel using the Fisher-
Tropsch process.
Further, biomass gasification shows a higher efficiency compared to other processes such
as direct combustion, pyrolysis, and liquefaction [1-2]. The product gas (also known as synthesis
gas or syngas; hereafter syngas) from the biomass gasification is a mixture of carbon dioxide
(CO2), carbon monoxide (CO), hydrogen (H2), methane (CH4), water (H2O) and nitrogen (N2) if
air is used as a gasifying agent. Syngas has been mostly accepted for power generation and is
considered to be more mature technology compared to other biomass conversion processes [3].
Fixed bed reactors are widely used for gasification of coal, biomass including municipal waste
utilization, because of their simplicity in design and efficiency [4].The choice of biomass for
gasification depends upon demographic factors. In the United States, Midwestern states have
abundant agricultural residues such as corn stover and wheat straw, whereas southern states have
more forest residues. These biomass feedstocks vary in their composition, which ultimately
affects the syngas composition. An exhaustive amount of literature is already present for biomass
gasification using various regionally appropriate feedstocks. Most of the time, syngas
43
composition is unknown until the gasification work is conducted. Experimental work is often
resource-intensive (time and money) and a mathematical model predicting syngas composition
(concentration of H2, CO, CH4, and CO2) using elemental analysis of biomass would be helpful.
semi-transient and transient that can be used to determine the syngas composition [5]. Among
these techniques, the thermodynamic equilibrium model is the simplest of all type and gives
syngas composition for various biomass types at selected gasification temperatures with
mechanical and chemical equilibrium. Chemical equilibrium is the state of minimum Gibbs free
energy and maximum system entropy. Mechanical equilibrium occurs when the system is not
prediction when the reaction temperature is sufficiently higher [6]. Equilibrium conditions are
difficult to achieve in practical operating conditions and results obtained from thermodynamic
equilibrium modeling can serve as the maximum limit on syngas composition. A few studies
have been conducted to determine syngas composition and heating value of syngas using
Watkinson et al. [9] have developed a thermodynamic equilibrium model and compared
their result with various types of gasifiers used for coal. The study found the best prediction for
entrained bed gasifier with a lower degree of accuracy in predicting syngas composition from
fluidized bed and moving bed gasifiers. Jarungthammachote and Dutta [7] and Melger et al. [8]
have predicted syngas composition from various biomass types using thermodynamic
temperature through an iterative process and the syngas composition at given equivalence ratio.
44
The objective of this chapter is to develop a mathematical expression to determine syngas
composition based on carbon, hydrogen and oxygen that can be applicable to any biomass type.
In addition, the analysis also includes the effect of moisture content in the biomass.
Thermodynamic results will be compared with the experimental data available for selected
biomass types.
2.3 METHODOLOGY
Empirical relations were developed for predicting the individual fraction of major
combustible species of the syngas, CO and H2. Although these equations can be used for any
type of gasifier, it is more accurate for a downdraft gasifier due to its low tar content. Although
there are several factors affecting syngas composition from biomass, it mainly depends on the
gasifier type, feedstock and feedstock pre-treatment, gasifying medium and operating parameters
such as temperature, pressure and equivalence ratio [10]. In this study, the formulation of the
- All carbon content in biomass is converted into gaseous form and the residence time
actual gasification process; however the degree of error introduced by this assumption
11]. The products taken into account are CO, CO2, H2, CH4, N2 and water. These are
the major gaseous compounds formed during the gasification process. Hydrocarbons
other than CH4 were assumed negligible in syngas and were not taken into account.
45
- Ash in the feedstock was assumed inert in all gasification reactions although it holds
true typically only for reaction temperatures less than 700C [12]. Herbaceous
biofuels contain silicon and potassium as the major mineral content which lowers ash
higher than 700C. Therefore, the relations derived in this study cannot be used
- All the gaseous products are assumed to behave as ideal gases. This will lead to
high temperature and low pressure. Also, the pressure drop inside the gasifier was
assumed to be negligible.
- The reaction was auto-thermal and no external source of heat was applied. The
process is completely adiabatic so that no heat losses occur from the gasifier. The
amount of air was varied to achieve the desired reaction temperature in the gasifier.
- The amount of tar in syngas was assumed to be negligible. This places the restriction
upon the use of this model for various configurations of gasifier design. For
downdraft gasifier, this assumption is valid since the tar concentration is significantly
low compared to other configurations [13-15]. For an updraft gasifier, where the
higher hydrocarbons produced are not cracked and hence contain high fraction of tar,
the results of this modeling cannot be applied [13]. Modifications such as subtracting
the amount of volatiles from the biomass and treating the rest as char for the
46
gasification can be done to improve the model [9]. However, this will lead to
increasing amount of error in terms of the final results. Output from the gasification
was assumed only to be permanent gases free of oxygen which is true because the
oxygen supplied is far less than that needed for combustion in a gasification process.
Sulfur and chlorine content in biomass were also neglected since they are less than
The chemical composition of biomass was taken to be in the form and the
(1)
(2)
The major reactions that occur inside the downdraft reactor are as follows:
(3)
47
(4)
The two reactions shown above can be combined into one single reaction (Eqn.5) known
(5)
The other reaction that is prominent in the gasification process is formation of methane as
shown below:
(6)
Eqns. (5) and (6) are the two major reactions that occur in the gasification process [6-8,
16].The equilibrium constant for these two above equations (5 and 6) as the function of their
(7)
(8)
48
Gibbs free energy is used in determining the value of K1 and K2 as presented in Eqn. (9).
For the given ideal gas, the Gibbs free energy is a strong function of the reaction temperature and
(9)
(10)
Eqns. (11-13) can be written by balancing carbon, hydrogen and oxygen moles,
(11)
(12)
(13)
Now, there are five equations (7, 8 and 11-13), and six unknowns ( ). The
final equation was obtained by an enthalpy balance inside the gasifier. Total enthalpy content in
49
any chemical species is the sum of its chemical enthalpy and sensible enthalpy and can be
written as follows:
(14)
Zainal et al. [6] have used HHV for predicting syngas composition from biomass.
However, the use of LHV for finding heat of formation is also common [7-8, 18]. In this study,
LHV is used for evaluating heat of formation of biomass. Heat of formation of biomass is
50
(15)
LHV is calculated in dry basis of biomass and was calculated using the following equation [18]:
(16)
The above equation (Eqn. 13) can be reduced to following form since
are zero at the reference temperature and pressure of 298 K and 1 atm.
(17)
51
Eqn. (17) acts as the constraint for the gasification process and forms the basis for
adjusting the amount of air to be supplied. The amount of air is adjusted in such a way that total
Cp can be determined using an empirical relation that holds for a wide range of temperature.
(18)
The sensible heat of each gas species was found by integrating Eqn. (18) from the
ambient temperature to gasification temperature. The value of c1-c4 is taken as reported by Reid
et.al. [19].
Species
52
Similarly, the change in Gibbs free energy for an individual gas is given by:
(19)
The values of a-g are taken from Probstein and Hicks [20] and are shown in Table 2.2
along with enthalpy of formation at standard reference state of 298 K and 1 atm pressure.
Species a b c d e f g
The model was run with an elemental composition of nearly 100 biomass feedstocks, which
were documented in a governmental database [21]. The feedstocks used in the model includes
pine, poplar, eucalyptus, corn stover, rice husk and many other common types of biomass. The
elemental compositions of these feedstocks were obtained from Syngas composition was
determined by solving six equations (Eqns. 7-8, 11-13 and 17) in MATLAB [22]. Newton-
Jacobi iteration was used for solving these equations. Complete MATLAB coding for these
53
overall equilibrium model are attached in Appendix A and B. Once the syngas composition was
determined from all feedstocks, a linear equation was developed to calculate the concentration of
each gas species. Syngas composition from all biomass feedstocks run in the model along with
their elemental and ash wt.% is reported in Appendix C. Multiple linear regression analysis was
performed to determine the coefficients for the linear equation using MS-EXCEL spreadsheet.
Table 2.3 shows the prediction of CO and H2 (combustible gases in syngas) for different
common types of biomass on dry basis using equilibrium model. These syngas composition are
computed from running the model at 800C. The equivalence ratio for all of these simulations is
automatically adjusted such that the pre-set temperature of 800C is achieved. Thus equivalence
ratio is not constant but a function of the elemental composition of biomass. Equivalence ratio
for the gasification of various biomasses as reported in Table 2.3 is in the range of 0.39-0.48. It
was found that increase in oxygen concentration in biomass reduces the equivalence ratio
because of inherent supply of oxidizing agent from biomass itself. The composition of syngas
predicted by the empirical formula is generally higher than the observed concentration at many
cases.
54
Table 2.3 CO and H2 composition for most common feedstocks available in the U.S from
MATLAB model.
Expressions for CO, H2 and CO2 were obtained in terms of three variables such as C, H
and O except for CO2 which is expressed as the function of two variables, C and O. The p-value,
which shows the significance of a parameter in regression analysis, was less than 0.00002 for
each of the independent variables asserting its influence in the individual syngas component.
55
Goodness of fit (R2) value obtained is higher than 0.98 showing significance of all independent
(20)
(21)
(22)
equilibrium modeling prediction was always less than 0.15% for biomass at a temperature of
800C and therefore, it is not presented here. Similar observations were reported in other
thermodynamic modeling studies [6-8]. Nonetheless, the methane concentration is in the range of
3-4 vol. % in actual gasification. The three relations shown above were derived assuming the
gasifier temperature to be 800C. On the other hand, the equivalence ratio was self-adjusted in
the model to maintain 800C so that the various gasification reactions inside the gasifier are self-
sustained. Temperature inside the gasifier is the optimal compromise between moisture content
56
very likely [23]. The Eqns. (20-22) gives the syngas composition for the temperature of 800C
which is adjusted for the particular biomass at the equivalence ratio enough to sustain
endothermic reaction and maintain the pre-set temperature of 800C. Moisture content can be
accounted in above correlation by the use of following values of C, H and O if the data for
(23)
(24)
(25)
Table 2.4 shows the comparison between the experimental results available from
literature and predicted values from equation derived from (20-21). Results are compared with
the corresponding references [16, 24-26]. H2 and CO composition as reported by Zainal et al.
[24], is the average of 57 test runs with the temperature around 700-900C for most of the
experimental duration. They used furniture wood as their feedstock with equivalence between
0.268-0.43. H2 and CO data from Bacaicoa et al. [27] is from experiment conducted in the
57
downdraft gasifier with capacity of 25-50 kg/hr and equivalence ratio of 0.247. The data taken
from Jayah et al. [25] is among one of their conducted experiments in the downdraft gasifier. The
syngas reported by Jayah et al. is between 18.4-22.1% of CO and 13-18.3% of H2 with the
temperature of the gasification zone in the range of 700-1000C. Comparison for both Bacaicoa
et al. [27] and Jayah et al. [25] was done with syngas composition at temperature close to this
model. CO and H2 composition from Altafini et al. [16] is the average of 10 test runs with the
reaction temperature around 832C and average air/sawdust ratio of 1.829. As can be seen from
Table 2.4, the predicted result is in good agreement with the experimental results, but the
composition of syngas predicted by the empirical relations is generally higher than the observed
The effect of various parameters are shown in Figure 2.1 to 2.8 for the particular biomass
with 50 wt.% carbon, 6 wt.% hydrogen and 44 wt.% oxygen which is the typical composition of
dry woody biomass. Figure 2.1 shows the effect of moisture content on the syngas composition.
The concentration of H2 increased from 16.9 vol.% to 17.8 vol.% with the change in moisture
58
content from 0 to 28 wt.% and started decreasing thereafter with further increase in moisture
content. The concentration of CO decreased monotonically with increase in moisture content and
the change in CO was more pronounced compared to the change in H2 with the same change in
moisture. CO decreased from 23.2 vol.% to 8.9 vol.% with an increase in moisture from 0 to 43
wt.%. As expected, CO2 concentration increases with increase in moisture content from 9.4
vol.% to 18.1 vol.% as moisture content increases from 0-43 vol.%. The methane concentration
is less than 0.15 vol.% for over the entire range of moisture content.
60 0.15
50 0.12
syngas species (% dry vol.)
H2
40 CO
CO2 0.09
N2
30 CH4
0.06
20
0.03
10
0.00
0 10 20 30 40 50
Moisture Content (% wet basis)
The effect of moisture content as shown in Figure 2.1 needs further elaboration at this
point. The model is based on the assumption that the process is completely adiabatic thus
additional air flow is required with an increase in moisture content to generate the heat required
to sustain the desired temperature. This equivalence ratio increases as seen in Figure-2. This can
59
be seen in the increase in concentration of N2 with increase in moisture content. In an actual
gasification process, if this air flow is not supplemented, decrease in gasifier temperature is
with increase in moisture content. The overall effect is the decrease in HHV of syngas with
increase in moisture content, which can be seen from Figure 2.2. HHV of syngas decreases from
5.1 MJ/m3 to 3.4 MJ/m3 with an increase in moisture content from 0 to 43 wt.%.
6 0.60
HHV
Eq. ratio 0.56
5
Higher heating value (MJ/m )
0.52
3
Equivalence ratio
0.48
3
0.44
2 0.40
0.36
1
0.32
0
0 10 20 30 40 50
Moisture Content (% wet basis)
Figure 2.2 Effect of moisture content on HHV of syngas under adiabatic condition
Figure 2.3 depicts the decrease in volumetric fraction of CO, H2 and CH4 (not appreciable
due to the inherently small concentration of CH4) with increase in temperature. This decrease is
due to the increase in dilution by N2 at higher temperature because the equivalence ratio adjusts
60
itself to meet the adiabatic condition set in Eqn. (17). The effect of temperature on equivalence
ratio can be seen graphically in Figure 2.5. Figure 2.4 depicts the actual number of moles of each
number of moles of H2 decreased monotonically. The number of moles of H2O and CO2
increased with the increase in temperature. Thus, effect of temperature presented in Figs. 2.3
and 2.4 is not the effect of increasing temperature alone but also the effect of increase in
75 1.0
60 0.8
Syngas species (% dry vol.)
H2
45 CO 0.6
CO2
N2
30 CH4 0.4
15 0.2
0 0.0
600 700 800 900 1000 1100 1200
o
Temperature ( C)
61
2.5 0.030
0.025
Number of moles of syngas species 2.0
H2 0.020
1.5 CO
CO2
H2O 0.015
N2
1.0
CH4
0.010
0.5
0.005
0.0 0.000
600 700 800 900 1000 1100 1200
o
Temperature ( C)
0.6
0.5
Equivalence ratio
0.4
0.3
0.2
0.1
0.0
600 700 800 900 1000 1100 1200
Temperature (oC)
Figure 2.5 Effect of temperature on equivalence ratio in adiabatic condition (variable xg)
62
The effect of temperature alone at fixed equivalence ratio (constant ) is shown by Figs.
2. 6 and 2.7. The model is run with the equivalence ratio of 0.396. This equivalence ratio is the
self-adjusted equivalence ratio for the particular biomass at 800C. Since the equivalence ratio
was fixed for developing Figs. 2.6 and 2.7, adiabatic condition is not valid. The increase in
temperature alone at fixed equivalence ratio results in an increase in the volumetric concentration
of CO as well as vapor in syngas, while concentration of CO2 decreases after reaching its
maximum value at around 850C. Concentration of H2 and CH4 decreases with the increase in
temperature and the CH4 concentration reach to negligible amount after 900C.
60 0.8
50
H2
0.6
CO
Syngas species (% dry vol.)
40 CO2
N2
CH4
30 0.4
20
0.2
10
0 0.0
600 700 800 900 1000 1100 1200
o
Temperature ( C)
Figure 2.6 Effect of temperature on syngas species at fixed equivalence ratio of 0.396
63
2.0 0.025
0.020
Number of mols of species
1.5
H2
CO
CO2 0.015
H2 O
1.0
N2
CH4 0.010
0.5
0.005
0.0 0.000
600 700 800 900 1000 1100 1200
o
Temperature ( C)
Figure 2.7 Effect of temperature on syngas species at fixed equivalence ratio of 0.396
64
Figure 2.8 compares the heating value of syngas for constant and variable fuel/air ratios.
For self-adjusted value of , HHV decreases rapidly with increase in temperature but increases
5
HHV (MJ/m3)
2
600 800 1000 1200
Temperature (oC)
significantly lower than those encountered in practical gasification tests. Typical CH4
concentration in downdraft gasifiers is 2-5 vol. % (dry basis). This is due to the fact that CH4
from the thermal cracking of tars and volatiles of biomass is not incorporated in the model. Since
CH4 is a very stable hydrocarbon, it is one of the main products of thermal breakdown of higher-
order hydrocarbons. Neglecting this phenomenon results in a lower CH4 prediction from the
65
2.5 CONCLUSIONS AND FINAL REMARKS
Empirical relations derived for CO, H2 and CO2 predict syngas composition with a
reasonable degree of accuracy. However, the relations become less accurate with increase in the
ash content in the biomass materials because a reaction with ash and heat absorbed by ash is
ignored in the model. Also, perfect adiabatic conditions are difficult to achieve in practical
gasifiers resulting in some discrepancy between predicted and experimental results. As the
temperature increases, the predicted values from this model and relation become more realistic.
Moisture content reduces CO fraction in syngas significantly and thus reducing HHV of the gas.
Thus, for the gasification process, it is essential to have the moisture content below a threshold
that meets the end-requirements. The amount of hydrocarbons produced during the gasification
process depends upon the temperature of the gasification and decreases rapidly with increase in
temperature. It is also seen that the concentration of CH4 above 900C is negligible as predicted
66
2.6 NOMENCLATURE
Carbon, hydrogen, oxygen, nitrogen and sulfur fraction in biomass (dry basis)
67
Number of moles of species i
Normalized coefficient of atomic hydrogen, oxygen and nitrogen for biomass molecule
68
2.7 REFERENCES
[1] A.V. Bridgewater, Renewable fuels and chemicals by thermal processing of biomass,
[2] A. Demirbas, Combustion characteristics of different biomass fuels, Progress in Energy and
[3] A. Demirbas, Progress and recent trends in biofuels, Progress in Energy and Combustion
[4] Y.B. Yang, V.N. Sharifi, J. Swithenbank, Effect of air flow rate and fuel moisture on the
burning behaviours of biomass and simulated municipal solid wastes in packed beds, Fuel, 83
(2004) 1553-1562.
[5] T.B. Reed, Principles and technology of biomass gasification, in: K.W. Boer, J.A. Duffie
(Eds.) Advances in solar energy, An annual review of research and development, Plenum press,
[6] Z.A. Zainal, R. Ali, C.H. Lean, K.N. Seetharamu, Prediction of performance of a downdraft
gasifier using equilibrium modeling for different biomass materials, Energy Conversion and
[9] A.P. Watkinson, J.P. Lucas, C.J. Lim, A prediction of performance of commercial coal
69
[10] J.M. Prins, Thermodynamic analysis of biomass gasification and torrefaction, in, Eindhoven
[12] B.M. Jenkins, L.L. Baxter, T.R.M. Jr., T.R. Miles, Combustion properties of biomass, Fuel
[13] R. Warnecke, Gasification of biomass: comparision of fixed bed and fluidized bed gasifier,
[14] S.M. Nunes, N. Paterson, D.R. Dugwell, R. Kandiyoti, Tar formation and destruction in a
simulated downdraft, fixed-bed gasifier: reactor design and initial results, Energy and fuels, 21
(2007) 3028-3035.
superficial velocity on tar from downdraft gasifier, Energy and fuels, 19 (2005) 1186-1191.
[16] C.R. Altafini, P.R. Wander, R.M. Barretoa, Prediction of the working parameters of a wood
waste gasifier through an equilibrium model Energy Conversion and Management, 44 (2003)
2763-2777.
[17] J.M. Smith, H.C.V. Ness, M.M. Abbott, Introduction to chemical engineering
[18] M.L.d. Souza-Santos, Solid fuels combustion and gasification Modeling, simulation and
[19] R.C. Reid, J.M. Prausnitz, B.E. Poling, The properties of gases and liquids, 4 ed., McGraw-
Hill, 1987.
70
[21] Biomass Property Database, in, U.S Department of Energy.
[23] J. Corella, M.J. Toledo, G. Molina, Calculation of the conditions to get less than 2 g tar/
Nm3 in a fluidized bed biomass gasifier, Fuel processing technology, 87 (2006) 841-846.
[24] Z.A. Zainal, A. Rifau, G.A. Quadir, K.N. Seetharamu, Experimental investigation of a
[25] T.H. Jayah, L. Aye, R.J. Fuller, D.F. Stewart, Computer simulation of a downdraft wood
[26] P.R. Wander, C.R. Altafini, R.M. Barreto, Assessment of a small sawdust gasification unit,
[27] G.P. Bacaicoa, R. Bilbao, J. Arauzo, M. L. Salvador, Scale-up of downdraft moving bed
gasifiers (25-300 kg/h) Designs, experimental aspects and results, Bioresource Technology, 48
(1994) 229-235.
71
CHAPTER 3
DOWNDRAFT GASIFIER
3.1 ABSTRACT
Auburn University and its collaborator, Community Power Corporation, have developed
a mobile 25 kWe downdraft gasifier. In this study, gasification of various biomass feedstocks
such as pine wood chips, commercial wood pellets, saw dust, peanut hulls and poultry litter (the
last four in pelletized form) were conducted. Ultimate and proximate analyses were carried out to
characterize the biomass feedstock used for gasification. The syngas obtained from various
feedstocks and different operating conditions were analyzed using the on-site gas analyzer.
Temperature distributions inside the gasifier for different feedstocks and operating conditions
were also examined. A minimum temperature difference across the reduction and combustion
zone was found in this gasifier. Gasification tests with commercial wood pellets were more
closely examined at various flow rates and carbon, enthalpy and exergy analyses were made.
syngas.
72
3.2 INTRODUCTION
combustible gases which can be subsequently used for energy application along with other
byproducts, such as volatiles, char and ashes. Under a broad classification, gasification systems
can be classified as moving bed and fluidized bed. Moving-bed is the oldest and simplest of all
gasification technologies and is generally more suitable for small-scale energy generation with
capacity less than 10-15 tons of biomass per hour [1]. When moving bed gasifiers have larger
capacities, there is non-uniformity of temperature distribution in the gasifier that results in low
quality synthesis gas (hereafter syngas). The gasification process is also more difficult to
optimize [2]. The most general forms of moving-bed gasifiers are updraft and downdraft, which
are only differentiated by the direction of flow of the gasifying (often called oxidizing agent)
agent with respect to fuel. The flow of the oxidizing agent is counter to fuel in the former case
The future prospects of gasification seem to be very promising in the United States. The
share of biomass in the total energy supply was 3.23% in 2007 and is expected to increase by an
annual average growth rate of 4.2% from 2007-2030, the highest growth rate amongst all other
energy sources [3]. The total biomass available in the United States is about 1.3 billion tons per
year with the southern states in the country accounting for about 423 million tons/year from
forest and agricultural residues [4]. Thus, the proper gasification of these residues will be
instrumental in reducing the nations dependency on fossil fuels, thereby increasing the energy
73
There are several parameters such as equivalence ratio, temperature, pressure, and
moisture content that influence the quality of synthesis gas produced from biomass gasification.
Equivalence ratio is one of the most studied parameters affecting syngas composition. It is
defined as the ratio of amount of air supplied to the biomass to the amount of stoichiometric
amount of air needed for complete combustion. Zainal et al. [5] observed the influence of
equivalence ratio upon the constituents of syngas as well as the calorific value and gas output
rate. The optimum equivalence ratio suggested was 0.38 for the downdraft gasification of
woodchips. Skoulou et al. [6] also investigated the effect of equivalence ratio and temperature
upon the quality of syngas in the downdraft gasifier for olive tree kernels and olive tree cuttings
and concluded an equivalence ratio of 0.42 optimal for the downdraft gasification of olive tree
cuttings and kernels. Based on previous studies, it can be assumed that the optimum equivalence
ratio is feedstock dependent [5-6]. The experiment conducted in the temperature range of 750-
950C also showed an increase in carbon monoxide (CO) and hydrogen (H2) with increase in
gasifier temperature; the overall effect was an increase in the heating value of syngas. Jayah et al.
[7] used a computer model calibrated by experimental results obtained from a typical medium-
scale downdraft gasifier for examining the effects of various parameters upon the biomass
conversion efficiency. A longer gasification zone, lower heat losses from the gasifier, lower
moisture content and higher air-inlet temperature all had a positive impact upon the conversion
efficiency. Longer gasification zone (hence increase in residence time of carbonaceous material)
inside the gasifier facilitates the conversion process but also increases the cost of building a new
plant. A gasification zone of 330 mm with a capacity of 80 kW was suggested for the downdraft
gasifier [7]. Bacaicoa et al. [8] studied the gasification of polyethylene and wood particle
mixtures in a downdraft gasifier with varying ratios of wood chips/polyethylene chips as well as
74
air flow rate. They found an increase in the calorific value of syngas with an increase in
in cold gas efficiency (ratio of chemical energy in the syngas to that of fuel) was also observed.
Lin et al. [9] gasified rice-husk in a downdraft gasifier in a baseline experiment for a pilot plant
design and concluded that the husk feed rate should be about 28 kg/hr for obtaining 10 kW of
power in a downdraft gasifier. Experiments by Wander et al. [10] with a downdraft, stratified,
open top gasifier showed independence between air/biomass ratio and mass conversion
efficiency in the temperature range between 500-900C while the mass conversion efficiency
depended upon the temperature and decreased below 800C. Sheth et al. [11] examined the effect
of various parameters on the performance of the gasifier and syngas composition. They found a
decrease in biomass consumption rate with an increase in moisture content and an increase with
increase in air-flow rate. Syngas composition, calorific value and gas output ratio with respect to
equivalence ratio were also reported. The optimal equivalence ratio was 0.20, beyond which a
decrease in calorific value was observed. Sharma [12] proposed an equilibrium model for a
downdraft gasifier which showed a decrease in CO with an increase in CO2 and H2 for feedstock
with higher moisture content. In another study, Zainal et al. [13] reported an increase in CO2 and
H2 for an increase in moisture content from 0% to 40% while noting a decrease in CO for the
Mass and energy analyses are very important since they serve as a validation of overall
gasification process. These are the applications of mass conservation and the first law of
thermodynamics. Carbon closure can be done to serve the purpose of the mass balance if the
latter is not possible due to experimental reasons such as inability to measure ash. The ratio of
input carbon from biomass to the sum of output carbon in various carbonaceous syngas
75
constituents is the measure of carbon closure. Total energy content in biomass can be measured
by the knowledge of higher heating value and total mass used in the experiment. Energy output
is the sum of chemical energy of the syngas output and the sensible energy gained in the
gasification process. The ratio of output energy to input energy gives the energy ratio of the
overall gasification process. While energy ratio can provide us the validity of the experiments, it
is not sufficient to measure the quality of energy that can be obtained from the gasification
process. Exergy is the amount of energy that can be used for useful work. Since exergy also
accounts for the losses due to irreversibilities of the process, exergy ratio is usually lower than
This chapter reports syngas composition and its heating value from wood chips and other
feedstocks such as pelleted wood, peanut hulls, sawdust, and poultry litter. Detail study was
conducted to calculate the carbon closure, energy and exergy ratio for the commercial wood
pellets due to their uniformity in size and moisture content. Low bulk density of many types of
agricultural residues entails the need for a densification process such as pelletizing which
increases the bulk density thereby improving the handling characteristics and significantly
reducing the space required to store and transport biomass [15]. Furthermore, the effect of
moisture content and biomass flow rate on composition and heating value of syngas obtained
To examine the syngas composition and heating value of syngas from selected feedstocks
76
To conduct mass, energy and exergy analyses of the overall downdraft gasification
process
schematic of an overall system is depicted in Figure 3.1 and a photograph of the gasifier is
presented in Figure 3.2. Figure 3.3 depicts the picture of the mobile gasifier. The interface
between the controls and the gasifier were created using LABVIEW. Biomass was loaded into
two bins (1 and 2) and passed through a sorting screen, which rejected biomass above 44.45 mm
(1.75 inch) and below 6.35 mm (0.25 inch). An auger was used to feed biomass into the gasifier
once it was sorted. The feeding rate was controlled based on a specified level of biomass inside
the gasifier. A level sensor detected the level of biomass in the gasifier and turned the feeder
(auger) once the biomass level fell below the set value. Air was used as an oxidizing agent for
biomass gasification. Primary air in the gasifier was obtained from the open top of the gasifier
(Figure 3.1). The gasifier had multiple secondary air injection nozzles where the secondary air
was fed with the 248.67 W (1/3 HP) air blower. The secondary air supplied was meant to
improve the combustion reaction and also to maintain uniformity in temperature along the
region. The grate was shaken at an adjustable regular interval via grate-shaker mechanism to
remove the ashes formed during the operation. The gasifier was also shaken at a regular interval
to facilitate the smooth flow of biomass inside the gasifier and prevent rat holes and bridging,
inside the gasifier. Charcoal left from the previous run (or fresh charcoal if the run was for the
77
first time) inside the gasifier was ignited with a cal rod (igniter/heater) before the fresh biomass
was fed inside the gasifier. The height and inside diameter of the gasifier reactor were 1200 mm
Figure 3.1 Schematic of the Auburn mobile downdraft gasifier designed by CPC
78
Figure 3.2 Photograph of the Auburn gasifier along with dimensional sketch of thermocouples
79
Figure 3.3 The Auburn mobile downdraft gasifier parked outside the capitol building in
Montgomery, AL
A certain protocol was followed prior to starting the gasifier for beginning of each
experiment, including conducting a pre-run check, involving checking leaks, various hose
connections, motors and electrical devices associated with the system, charcoal level, complete
flare raising etc. Heat exchangers and filters were preheated to 40C before igniting charcoal
inside the gasifier to avoid any condensation during the process. The heat exchangers and filters
were heated using electricity from the grid or using electricity generated from the generator (that
runs on propane gas). Once the filters and heat-exchanger were heated to 40C, igniter located
inside the gasifier was turned on. Igniter heats the charcoal inside the gasifier caused the
temperature inside the gasifier to rise. Temperatures across various locations inside the gasifier
80
seconds intervals. Among many thermocouples in the system, only the location of various
thermocouples of interest as measured from the grate is shown in Figure 3.2 schematically.
Secondary air supplies, via a char-air blower, were injected for each level, as the temperatures,
for the corresponding levels, reached 350C. Feeding of wood chips started only when
thermocouples measured above 800C in any three locations inside the gasifier.
Hot syngas coming out from the gasifier was cooled using the heat exchanger (HEX).
This results in the heating of air supplied by gas cooling blower and cooling of syngas. Char and
particulates were then removed by passing the cooled syngas through the filters. Even though the
cleaned syngas could be burned in an engine to produce power, the cleaned syngas was flared in
this study since the focus on the study was on effect of biomass feedstocks and operating
conditions on syngas composition. Syngas was sampled with online gas analysis system after
passing through the filters, the details of which are discussed in the following section.
The steady state of the system was indicated by steady temperature across the various
levels of the gasifier, most commonly 800C at any three locations among T1-T4. The time
required to attain steady-state generally varied from 30 minutes to one hour and was affected by
biomass feedstock type, and operating conditions. Once the gasifier reached steady-state,
parameters to be considered were the differential pressure in the HEX and the temperature of
syngas out of the HEX. An increase in higher pressure differential in the HEX suggests clogging,
which prevents smooth flow of syngas and reduces the effectiveness of the HEX. As the
effectiveness of HEX decreases, syngas temperature after passing through it is a concern as high
temperatures may ignite filter bags. Ideally, the amount of oxygen in the syngas has to be zero
but the system used in this experiment was bigger than that used in lab-scale experiments,
making this difficult. It was found that the oxygen level was between 0.5-2.0 (vol./vol.,% dry). If
81
the amount of oxygen in the syngas was more than 2 %, a leak inspection was carried out. A
typical temperature profile of the gasifier is attached in Appendix D (Table D. 7 and Figure D.1-
D.5).
After the experiment was completed, proper shut-down was carried out. During the
shutdown process, syngas flow rate increases because of an increase in the roots-blower speed.
This is pre-set to make sure that no gas remains inside after leaving the gasifier. More air is
drawn to obtain higher combustion so that fresh wood chips inside the gasifier can be burned
faster and are converted into charcoal for the next run. This also removes any smoke that may
occur. The feeder is automatically turned off and temperature increases to rise near the upper
level of the gasifier due to a higher combustion rate. The feed-gate and roots-blower will turn off
once the upper level of the gasifier attains a certain temperature that is sufficient to sustain slow
pyrolysis. This varies according to the operating conditions. The system is thus completely shut
down. Most of the experiments were carried out for 4 to 5 hours except in the case of poultry
litter where steady state conditions could not be supported for more than two hours.
This study was conducted to examine the effect of various feedstocks on syngas
composition. Five types of feedstock (pine wood chips, peanut hulls, poultry litter, saw dust
pellets and commercial wood pellets) shown in Figure 3.4 were tested. Commercial wood pellets
were obtained from American Wood Fiber, Columbia, MD. Furthermore, the effects of moisture
content and biomass-flow rate in gas composition and its calorific value were analyzed for pine
wood chips.
82
Carbon, energy and exergy analyses for the gasification tests were done only when
commercial wood pellets was used as a feedstock. Wood pellets were fed into the gasifier by an
external biomass feeder for an accurate measurement of the mass used in each experiment.
Instead of analyzing syngas after filters, syngas was sampled immediately after the gasifier
Experiments with commercial wood pellets were run for almost 4 hours for all tests.
Since the design of this gasifier makes it difficult to control the biomass feed rate
directly, an alternative approach is taken to achieve varying biomass feed rates. Since syngas
output rate depends upon the biomass feed rate in the system, the syngas output rate (which can
be easily computer controlled in this gasifier) was varied to control the biomass feed rate.
Syngas composition was measured using a gas analyzer (Nova 7905AQN4, Niagara Falls,
NY) on site and the gas composition was measured in a volumetric basis. The gas analyzer used
has the capacity of measuring oxygen (O2):0-25%, CO: 0-25%, CO2: 0-20%, methane (CH4) 0-
10%, 0-20% H2. It uses non-dispersive infrared (NDIR) detector for CO, CO2 and CH4 and
temperature controlled thermal cell for H2 and electrochemical sensor for O2. The accuracy of
this instrument is 1% of full scale. The gas analyzer was calibrated with air for oxygen. For
83
other gases, a known mixture of gases with the following composition was used: CO 25.16%,
CO2 20.05%, CH4 9.968%, H2 20.04% and the balance was nitrogen.
Once the gasifier reached a steady state temperature, syngas data (CO, CO2, H2, CH4 and
O2) were logged into a computer at 15 s intervals via the data logging software supplied with the
gas analyzer. The remaining volumetric proportion was assumed to contain only nitrogen.
Gasification temperature is not high enough to form nitrogen oxides (NOX) and the nitrogen
content in biomass is also fairly low except in the case of poultry litter. Therefore, NOx
measurement was not carried out in this study. Air-flow rate was calculated assuming that the
source of nitrogen was from air only and thus atmospheric mass proportion of nitrogen was
utilized for the calculation. The volumetric content of syngas constituents multiplied by their
individual higher heating value (HHV) gave the overall higher heating HHV of the syngas as
(1)
In above equation, HHVsyngas is the heating value of syngas while HHVi and yi are higher heating
value and volumetric fraction of syngas constituents (i=H2, CO, CH4). The HHV of H2, CO and
CH4 are 12.76 MJ/m3, 12.6 MJ/m3 and 39.8 MJ/m3, respectively [16].
Moisture content was measured following ASTM standard E871-82 where a representative
sample of biomass feedstock used for the experiment was heated for 16 hr at 103C to calculate
the mass difference and hence the moisture content [17]. Bulk density of biomass feedstock was
84
measured by determining a known quantity in a standard container with volume of 946.35 mL (1
quart). Ash content was measured according to the ASTM standard E 1755-01. This involves
heating of biomass sample (0.5 g-1 g) in a muffle furnace to 575 25C for three hours and
finding the amount that remains in the container [18]. Biomass samples were sent to Midwest
Microlab, LLC (Indianapolis, IN) for an ultimate analysis. Although the procedure for measuring
elemental composition (ultimate analysis) varies from instrument to instrument, the basic
principle for almost all is the combustion of small biomass sample in a pure oxygen environment
and subsequent measurement of C, H, N and S in the output stream which can easily be found.
HHV of the biomass was calculated by Dulong and Petits Formula given in Eqn. (2) using
results from ultimate analysis and also experimentally using an oxygen bomb calorimeter (IKA,
model C200, Wilmington, NC) with reference to ASTM D 2015-96 for verification [16, 19].
Moisture content in biomass samples was determined by calculating the weight loss of samples
by heating in an oven at 103 oC for 16 hours using ASTM E 871 Standard [20]. Ash fusion
temperatures were determined using ASTM D 1857 Standard in Alabama Power General Test
laboratory (Birmingham, AL) and Hazen Research Inc. (Golden, CO) [21]. Results of ultimate
and proximate analyses along with HHV are shown in Table 3.1.
(2)
where C, H, O and S are the carbon, hydrogen, oxygen and sulfur content of biomass in dry
basis.
85
Table 3.1: Characteristics of biomass feedstock used for gasification
Pellets
(% mass)
Ultimate Commercial Wood Chips
Poultry Peanut hulls Saw dust
Analysis, wood*
wt.%
C 22.1 47.8 45.2 47.7 45.2
H 4.3 5.5 5.8 6.0 5.6
N 3.1 0.8 0.3 0.0 0.1
S 0.6 - 0.0 0.0 0.0
O 31.3 43.1 46.3 45.8 47.7
Ash, wt.% 33.3 2.78 2.29 0.44 0.33
HHV,
11.21 (10.34) 18.67 (15.91) 18.07 (15.48) 18.34 (16.51) 18.82 (15.05)
MJ/kg
Bulk density,
680 790 725 750 210
kg/m3
M.C wt.% 8.5 5.1 4.7 2.5-5.3 17.6-25
*: Ultimate analysis done for dry pellets
ash free basis, by difference, -: Not detectable, values within the parenthesis are calculated
Syngas composition from selected feedstocks was evaluated at a constant syngas output
flow rate and the average values are reported in Table 3.2. Although equal moisture content for
all the feedstocks was not achieved as a proper moisture controlling set-up was not available, the
results show comparison between syngas from various feedstocks. All other feedstocks except
wood chips were tested as they were received. Wood chips were tested at the moisture content of
19.6 wt. %. Although syngas flow rate was set to 65 Nm3/hr, mass flow rate varied from 26.5
kg/hr for wood chips to 31 kg/hr for peanut hulls pellets under same experimental conditions. As
expected, pellets have slightly higher feed rate than wood chips due to better flow characteristics.
86
Syngas from peanut hulls pellets showed the largest fraction of CO as well as H2 and thus
had the highest HHV among all the feedstocks selected for this study. The carbon fraction of
peanut hulls, as can be seen from the ultimate analysis, had the highest fraction of carbon
amongst all the feedstocks under consideration for the current experiment. This could perhaps be
the reason why gasification of peanut hulls showed the highest heating value of syngas. Other
impacts could be due to the moisture variation and some difference in mass flow rate. Further,
research with accurate control on moisture and other conditions are already initiated and the
results will be reported in the future. The overall HHV as well as total volumetric combustibles
(CO, CH4 and H2) were found to be the lowest for poultry litter due to its high level of ash
vol. %
Feedstock HHV(MJ/m3 )
O2 CO CO2 CH4 H2
Peanut hulls (5.1) 0.50.1 22.80.7 8.90.5 2.70.3 20.10.3 6.10.2
Saw dust (4.7) 1.10.3 22.20.5 8.30.3 3.00.2 19.40.3 6.00.2
Poultry litter (8.5) 0.80.1 20.92.3 8.81.6 1.20.4 16.22.0 4.80.7
Commercial pellets (3.5) 0.50.2 22.10.8 10.40.7 1.90.2 16.61.0 6.10.2
Wood chips (19.6) 0.70.1 21.11.3 12.20.9 2.30.4 20.40.5 5.70.4
sign followed by numerical values are standard deviation, number given within the
To examine the effect of moisture content, syngas flow rate was set to 65 Nm3/hr thus
adjusting almost equal mass flow rate for all experiment which was about 26.5-27.5 kg/hr. Table
3.3 depicts syngas composition at various moisture content with pine wood chips.
87
Table 3.3: Effect of moisture content in syngas composition
Although biomass gasification is a complex process, the following reactions typically can
be used to represent the gasification process inside the gasifier [6, 22].
(3)
(4)
(5)
(6)
88
(7)
(8)
(9)
As can be seen from Eq. 7, an increase in moisture content decreases the amount of CO
and produces CO2 and H2. As expected, a decrease in CO was observed in the experiments with
an increase in moisture content but the increase/decrease in CO2, H2 and CH4 were not
significant. The gasifier used for current research was temperature-controlled which tries to
adjust its pre-set temperature (usually set at 800C) at various locations by increasing/decreasing
the amount of secondary air through proportional valve opening. However, H2 proportion in
syngas is also a strong function of temperature. As can be seen from Eqns., (4), (6) and (7),
reactions producing hydrogen are highly endothermic in nature so they demand high heat. The
temperature distribution for various moisture content discussed in a later section shows a similar
temperature profile. Due to this reason, there could be a small change in the hydrogen
concentration.
89
3.4.3 EFFECT OF BIOMASS FEED RATE IN SYNGAS COMPOSITION
The moisture content of wood chips used to analyze the effect of biomass feed rate varied
from 17.6 wt.% to 19.6 wt.%. Biomass feed rate was varied by setting the syngas flow rate which
automatically adjusts the biomass flow. Three syngas flow rates selected were 65, 55 and 45
Nm3/hr, which changed the biomass feed rate. The CO and H2 were found to increase slightly
with an increase in biomass feed while oxygen decreased. The effect of biomass feed rate upon
CH4 did not show any significant pattern. Syngas composition from different biomass feed rate is
Biomass
Moisture vol. %
flow rate HHV
content
( kg/hr, wet MJ/m3
(% wet basis)
basis) O2 CO CO2 CH4 H2
16.4 17.7 1.8 0.4 15.6 1.7 11.4 1.1 2.2 0.6 19.4 0.8 4.9 0.5
21.5 17.6 1.1 0.3 20.9 2.2 11.4 1.4 2.6 0.7 20.2 0.8 5.8 0.6
26.5 19.6 0.7 0.1 21.1 1.3 12.2 0.9 2.3 0.4 20.4 0.5 5.7 0.4
Moisture content reduces the reactor temperature of the gasifier due to heat absorption for
its evaporation. However as seen in Figure 3.5, deviation in temperature at various heights inside
the gasifier was found to be less than 30C even at a 5.4 wt % increase in moisture content. This
was probably due to automatic adjustments in the gasifier which tries to maintain the pre-set
temperature at different zones by increasing the amount of secondary air flow hence promoting
90
combustion around that region. Thus, Equivalence ratio was higher biomass with higher moisture
content. The equivalence ratio for various moisture contents is shown in Table 3.5.
1000
950
Temperature (oC)
900
850 25% MC
19.6% MC
800
23% MC
750
700
T1 T2 T3 T4 Grate
Location inside the gasifier (distance from the grate)
The temperature variation for various feedstocks in the gasifier is shown in Figure 3.6.
The moisture content and biomass flow rate corresponding to different biomass in Figure 3.6 is
shown in Table 3.6. While the temperature along the gasifier height was approximately equal for
saw dust and wood chips, a lower temperature distribution was found for the gasification of
91
peanut hulls. On the other hand, the temperature distribution for poultry litter gasification had a
different profile than the rest of the feedstocks with a sudden increase from T1 to T2. This
sudden increase in temperature was due to ash fusion around the vicinity of initial start-up
ignition and thus heat localization at one point which resulted in high temperatures. Since this
fused ash had high thermal resistance, no heat diffuses to the upper part of the gasifier.
Table 3.6 Moisture content and biomass flow rate for different feedstocks
1000.0
900.0
800.0
Temperature ( 0C)
300.0
T1 T2 T3 T4 Grate
Location inside the gasifier (distance from the grate)
92
1000
950
Temperature (0C)
900
700
T1 T2 T3 T4 Grate
Location inside the gasifier (distance from the grate)
Decreasing the biomass flow rate into the gasifier increases the grate temperature. The
moisture content of woodchips (% wet basis) is 19.2% to19.7% for the three biomass flow rate
shown in Figure 3.7. Lower feed rate increases the residence time of the biomass inside the
gasifier, and hence promotes efficient reduction reactions inside the gasifier. These reactions are
endothermic, hence the reduction in temperature with decrease in mass flow rate.
Due to good flow characteristics of pellets and a higher density than wood chips, the
residence time was increased from experience. Setting the same residence time resulted in
incomplete combustion of pellets and also a high pressure differential inside the gasifier in many
cases. The increase in residence time was achieved mainly by adjusting the interval of the grate-
shaker and of the gasifier vibrator itself. Normally when running the experiment with pellets, the
93
frequency of the grate shaker was decreased by a factor of four and gasifier vibration frequency
by 1.3 compared to the settings used for wood chips. Despite the modifications made, problems
were encountered with the gasification of poultry litter. It has been reported that agglomeration
could occur in the gasification of the feedstock with moisture greater than 5 wt.% of ash-content
[22]. Poultry litter test runs could not be sustained for more than two hours without significant
agglomeration. The low-melting mixture created blockage inside the gasifier and thus further
gasification was not possible. The snapshot shown in Figure 8 shows one of the clinkers removed
Figure 3.8 Ash agglomeration in the grate of gasifier after the gasification of poultry litter
Formation of ash clinker was assumed to be due to the low ash fusion temperature of the
minerals inherent in poultry litter. Abelha et al., reported the ash fusion temperature to be 660 C
for poultry litter [23]. Gasification of poultry litter in a downdraft gasifier may also be greatly
affected by the temperature inside the gasifier. The temperature should thus be such that it can
sustain the gasification but at the same time be lower than the ash-fusion temperature of the
poultry litter. Surprisingly, the results in the current study reported in Table 3.7, showed that the
ash fusion temperature for poultry litter was significantly higher than the numbers reported in the
literature [16]. Analyses were performed in two different laboratories to validate the results and
94
they were within 5% variation. Further, ash-fusion temperature of peanut hulls and poultry litters
do not differ significantly although, no problem was noticed while gasifying peanut hulls.
Reducing Atmosphere, oC
Softening Hemispherical
Feedstock Initial Temp. Fluid Temp.
Temp. Temp.
Results obtained from Alabama Power General Test laboratory, Numbers in the parenthesis are
PELLETS
Since the source of carbon input is only from pellets, amount of carbon fed into the
gasifier is calculated from pellets-flow rate and ultimate analysis. The amount of carbon output
was measured from the flow rate of individual carbonaceous syngas constituents.
The following assumptions were made in calculating the energy and exergy of biomass
95
-Ideal gas consideration for the syngas and its constituent gases
The gasifier used in the experiment is operated at atmospheric pressure and the pressure
drop across the gasifier is not significant. Gases can be treated as ideal gases at low pressure and
high temperature. Thus, the assumptions stated above are valid and introduce a negligible
The following formulas are used to calculate energy and exergy of individual gases in
syngas [24].
(10)
(11)
Ei and Exi are the energy and exergy of the gas in MJ/kg at the temperature T in kelvin
while E0i and Ex0i represents energy and exergy of the gases at the reference or dead state (Td)
taken to be at 25C or 298 K, respectively. The chemical energy (also known as the enthalpy of
formation) is taken from the corresponding reference [25]. The specific heat capacity, Cp is in
(12)
The values for the coefficients for Eqn. (12) are listed in the corresponding reference and
96
Table 3.8 Coefficients for the specific heat capacity of various gases
Gases c0 c1 c2 c3
Exergy of dry ash-free biomass without any sulfur is found using Eqn. (13) and (14) [14].
(13)
(14)
In Eqn. (13) and (14), is the ratio of chemical exergy of the biomass to the lower
heating value of the organic fraction of biomass (LHVorg). H, C, O, N denotes the hydrogen,
carbon, oxygen and nitrogen fraction by weight in biomass. Lower heating value of feedstock
(15)
The tests were run at different biomass flow rate dictated by the syngas flow rate set for
the experiment. Carbon, energy and exergy analyses were done for each test and are reported in
97
Table 3.9. Ideally, carbon closure is expected to be unity since input should be equal to output.
However, carbon closure was found greater than one for some experiments which might be due
to some instrumental errors as well as the residual biomass inside the gasifier, which could not
be measured due to operational difficulties. Also, the size of the gasifier contributed to these
discrepancies in carbon closures because of the higher probability that significant amount of
biomass that can remain in the gasifier after the completion of experiment. Wander et al. [10]
reported the similar carbon closure in their experiments with downdraft gasifier with the biomass
flow rate capacity of 12 kg/hr. Carbon closure obtained in all experiments were higher than 0.89,
98
Table 3.9 Carbon, energy and exergy analyses of commercial wood pellets
Wet
Moisture
biomass Carbon Eout Ein Exout Exin
content SF Fs Eout /Ein Exout /Exin
flow rate Closure (MJ/kg) (MJ/kg) (MJ/kg) (MJ/kg)
(% w.b)
(kg/hr)
17.6 3.4 45 45.0 1.08 311.5 310.4 1.00 224.0 337.0 0.66
18.0 3.8 45 45.8 1.04 318.7 316.6 1.01 223.5 343.8 0.65
18.6 4.5 45 45.0 0.98 306.9 325.4 0.94 215.1 353.3 0.61
18.7 4.1 45 44.3 0.99 299.9 327.6 0.92 210.9 355.7 0.59
19.0 2.7 45 47.4 1.02 335.8 337.4 1.00 239.0 366.4 0.65
19.8 3.4 55 55.0 1.15 382.1 349.1 1.09 274.2 379.0 0.72
99
20.6 5.3 55 54.9 1.15 390.3 356.6 1.09 279.5 387.2 0.72
23.1 3.8 55 50.0 0.91 358.4 405.4 0.88 258.7 440.1 0.59
24.6 2.7 65 59.9 0.98 424.4 437.9 0.97 299.3 475.4 0.63
24.9 5.3 65 65.0 1.13 465.9 430.7 1.08 332.2 467.7 0.71
26.5 3.7 65 57.6 0.89 403.0 465.4 0.87 285.5 505.3 0.57
27.0 3.8 65 62.6 0.96 442.7 475.6 0.93 317.5 516.3 0.61
28.8 3.5 65 65.0 0.90 445.0 508.9 0.87 319.1 552.5 0.58
SF: Syngas flow rate set for the experiment (Nm3/hr), Fs: Actual flow rate of syngas (Nm3/hr), : in wet basis,
Ein, Eout, Exin, Exout are input energy, output energy, input exergy and output exergy, respectively.
99
High energy ratios observed in all of the gasification tests reported in Table 3.9 shows
that the heat losses from the gasification system is minimal, and almost all of the energy present
in biomass is retained in the syngas. An Exergy ratio varies from 0.57 to 0.72 which is similar to
those reported in literature. Rao et al. [25] reported exergy ratios from 64% to 66% for different
The moisture content of the pellets was measured after each experiment and moisture
content of the pellets used in the experiment was found to be in the range of 2.5-5.3% (wet
basis). Figure 3.9 shows the volumetric fraction of individual gases with respect to biomass flow
rate.
30.0
Product gases (% vol /vol)
25.0
20.0
CO
15.0
CO2
10.0
CH4
5.0
0.0 H2
Figure 3.9 Biomass flow rate versus product gases for wood pellets
Syngas composition from wood pellets do not show any specific pattern change with
change in biomass feed rate inside the gasifier. The automatic adjustment of the gasifier, which
is difficult to control manually, tries to pre-set the temperature pre-set so that the effect of
100
biomass flow rate alone cannot be seen at constant syngas output. Thus, the temperature tries to
remain consistent as seen from the Figure 3.10 in spite of the change in biomass flow rate.
1000
900
Temperature (oC)
800 T1
T2
700
T3
600 T4
grate
500
15 18 21 24 27 30
Biomass flow rate (kg/hr)
Figure 3.11 shows the effect of biomass flow rate on HHV of the syngas. For all the
experiments conducted with wood pellets, HHV lies between 5.7-6.1 MJ/m3.
101
7
HHV ( MJ/m3)
5
3
15 18 21 24 27 30
Biomass flow rate (kg/hr)
Results based on biomass gasification using selected feedstocks in the Auburn mobile
commercial gasifier were presented along with the extended study of the effect of biomass flow
rate on syngas composition for commercial wood pellets. Syngas obtained from the gasifier was
found to have appreciable HHV for atmospheric air gasification. This particular gasifier was
found to maintain almost a constant pre-set temperature over a wide range of moisture contents
investigated in this study while decreasing temperatures were observed as the biomass feed rate
Difficulties encountered during the gasification of poultry litter warrants further research on
finding the suitable operating parameters as well as feedstock treatment. Carbon closures were
greater than 0.89 for all of the experiments conducted with commercial wood pellets suggesting
high carbon conversion efficiency of the gasifier. High energy ratios were obtained which
indicates negligible amount of heat losses from the gasifier. Exergy ratio of the gasifier was from
0.63 to 0.78 indicating significant amount of useful energy that can be recovered from syngas.
102
3.6 REFERENCES
[2] I. Olofsson, A. Nordin, U. Sderlind, Initial Review and Evaluation of Process Technologies
and Systems Suitable for Cost-Efficient Medium-Scale Gasification for Biomass to Liquid
[3] Annual energy outlook 2009 in: D.o.E. Energy information adminstration (Ed.), Washington,
DC, 2009.
[4] A. Milbrandt, A geographic perspective on the current biomass resource availability in the
United States, in, National renewable energy laboratory, Golden, Colorado, 2005.
[5] Z.A. Zainal, A. Rifau, G.A. Quadir, K.N. Seetharamu, Experimental investigation of a
olive tree cuttings and olive kernels in a downdraft fixed-bed gasifier, International Journal of
[7] T.H. Jayah, L. Aye, R.J. Fuller, D.F. Stewart, Computer simulation of a downdraft wood
(2008) 5485-5491.
[9] K.S. Lin, H.P. Wang, C.J. Lin, C.-I. Juch, A process development for gasification of rice
103
[10] P.R. Wander, C.R. Altafini, R.M. Barreto, Assessment of a small sawdust gasification unit,
[11] P.N. Sheth, B.V. Babu, Experimental studies on producer gas generation from wood waste
[12] A.K. Sharma, Equilibrium modeling of global reduction reactions for a downdraft (biomass)
[13] Z.A. Zainal, R. Ali, C.H. Lean, K.N. Seetharamu, Prediction of performance of a downdraft
gasifier using equilibrium modeling for different biomass materials, Energy Conversion and
[14] K.J. Ptasinski, M.J. Prins, A. Pierik, Exergetic evaluation of biomass gasification, Energy,
32 (2007) 568-574.
[16] P. Basu, Combustion and gasification in fluidized beds, Taylor and Francis Group, LLC,
[17] ASTM, ASTM Standards E 871-82-standard test method for moisture analysis of particulate
[18] ASTM, ASTM Standards E 1755-01.standard test method for ash in biomass, in, PA, USA,
2007.
[19] ASTM, ASTM Standards D 2015-96Standard test method for gross calorific value of coal
and coke by the adiabatic bomb calorimeter., in, PA, USA, 1998.
[20] ASTM, ASTM Standards E 871 - 82. Standard test method for ash in biomass, in, PA, USA,
2006.
104
[21] ASTM, ASTM Standards D 1857-04. Standard Test Method for Fusibility of Coal and Coke
[22] P. McKendry, Energy production from biomass (part 3): gasification technologies,
Leahy, Combustion of poultry litter in a fluidised bed combustor[small star, filled], Fuel, 82
(2003) 687-692.
[24] R.E. Sonntag, C. Borgnakke, G.J.V. Wylen, Fundamentals of thermodynamics, Sixth ed.,
[25] M.S. Rao, S.P. Singh, M.S. Sodha, A.K. Dubey, M. Shyam, Stoichiometric, mass, energy
[26] M. Dogru, C.R. Howarth, G. Akay, B. Keskinler, A.A. Malik, Gasification of hazelnut
[27] M. Dogru, A. Midilli, C.R. Howarth, Gasification of sewage sludge using a throated
downdraft gasifier and uncertainty analysis, Fuel Processing Technology, 75 (2002) 55-82.
105
CHAPTER 4
4.1 ABSTRACT
A study was conducted to see the different tar compounds and the effect of biomass flow
rate on tar concentration in a stratified downdraft gasifier. Tertiary condensed tar products such
as toluene, o/p-xylene, naphthalene, phenol, styrene and indene were observed in significant
amount. Tar concentration in the syngas was found to be in the range of 0.34-0.68 g/Nm3 lower
4.2 INTRODUCTION
Milne and Evans [1] defined tar from the gasification process as a material in the syngas
which condense inside the gasifier or in the equipment used for handling the product stream to its
end use. Tar compounds are largely aromatic in nature. They further classify the tar obtained due
106
Table 4.1 Classification of tar from thermal cracking of biomass
Among the type of products in the tar classified above, condensed tertiary products are
formed as a result of consecutive conversion of primary tar at high temperature, and thus these
two types, condensed tertiary and primary tar products are not usually found in the syngas at the
The maximum limit of tar concentration in syngas varies depending upon its end use. The
tolerable limit of tar concentration in syngas is 50-500 mg/Nm3, 50-100 mg/Nm3, less than 0.5
mg/Nm3 and less than 5 mg/Nm3 for compressors, internal combustion engines, methanol
synthesis and gas turbines, respectively [1]. Tar production in a downdraft gasifier is much
lower than in both updraft and fluidized bed gasifiers although it may not meet the requirements
needed to be used directly without prior treatment in power generation applications and liquid
fuel synthesis processes [2]. While liquid fuel syntheses from syngas requires purity in the
reacting gases and thus tar removal, the major problem with tar, when used in power generation,
is condensation at low temperature which creates blocking as well as fouling in power plant
equipment [3]. Hence, subsequent treatment is usually warranted depending upon the end use of
the syngas. Also the nature of tar from gasifier varies according to its design. Downdraft
gasifiers produce tertiary tar while tar from updraft gasifiers contain mostly primary tar due to
107
lower possibility of tar cracking inside the gasifier [1]. Syngas from fluidized bed gasifiers
contain tar which is the mixture of secondary and tertiary tar [1]. Tar content in a downdraft
gasifier is usually in the range of 0.01-6 g/Nm3 while updraft and fluidized bed gasifiers usually
have the tar content about 50 g/Nm3 and 6-12 g/Nm3 in average, respectively [1]. The residence
time and temperature in the gasification and reduction zones is the most important factor that
determining the level of tar in a downdraft gasifier [4]. With increase in temperature, tar content
in the syngas decreases due to thermal cracking [5]. Li et al., have reported that with increase in
temperature from about 700C to 820C, tar content decreases significantly from 15 to 0.54
g/Nm3 in a circulating fluidized bed [6]. Figure 4.1 shows the relationship between gasifier
reaction temperature and tar yield [7]. It can be observed from the figure that increase in
temperature significantly reduces the fraction of liquid, and hence the tar from the gasification
process. Increase in equivalence ratio also decreases tar content at the expense of higher level of
combustion inside the gasifier and results in a higher concentration of CO2, which is an
undesirable product [8]. Although tar concentration in syngas from a downdraft gasifier is
usually lower, these tars are also more stable and might be difficult to crack and remove
depending upon the end-need [9]. For use in an internal combustion engine, concentration of tar
should be less than 100 mg/Nm3 for successful long-term operation [1, 10].
108
Figure 4.1 Effect of maximum reactor temperature on tar production (Baker et.al [7])
Milne and Evans have discussed the tar reduction procedure as any one of physical,
thermal and catalytic techniques [1]. Han and Kim [5] have classified tar reduction methods into
cracking and plasma method. The mechanism method can effectively remove tar from 40-99% in
syngas but the useful energy that can be achieved from tar conversion is lost. In other methods,
tar is converted into other gases which increase the heating value of the syngas thus increasing
the energetic efficiency of the process. Devi et al. [3] suggest three methods for tar removal
which are adjustments of the operational parameters, addition of bed additives/catalysts and
gasifier modification. One-lump model, as shown in Figure 4.2 by Li and Suzuki [11], considers
all tar compounds lumped together as tar which disappears after simultaneous application of
various cracking and reforming processes and finally appears as secondary gases.
109
Figure 4.2 One-lump model for tar reduction (Li and Suzuki [11] )
This chapter quantifies different tar compounds present in tar from stratified downdraft
gasifier. Also reported is the tar concentration from the gasification of wood pellets as a function
Experiments to quantify the tar concentration in syngas were done with wood pellets
from the commercial source as the feedstock. The ultimate and proximate analysis of wood
110
Table 4.2 Ultimate and proximate analysis of wood pellets
Figure 4.3 shows the schematic diagram of the experimental set-up used to measure tar
concentration in a syngas sample from the gasification of wood pellets in a stratified downdraft
gasifier. Syngas is sampled from the port immediately after the downdraft gasifier and passed
through impinger bottles each containing 50 mL of isopropyl alcohol. The first impinger bottle is
kept at ambient conditions while the other two are kept in an ice-bath in order to maintain the
temperature around the freezing point of water. The tar present in the syngas condenses under
these conditions in the impinger bottles and later can be quantified. The water absorber after the
impinger bottles attracts all the moisture present in the syngas stream after condensation and the
syngas leaves dry after passing through the water-absorber. A flow-meter placed after the water-
absorber measures the syngas flow rate which is required to find the tar concentration per
standard volume.
111
The tar components were analyzed with an Agilent 7890 GC/5975MS using DB-1701
column (30 m; 0.25 mm i.d.; 0.25 mm film thickness). Thirty-one compounds were selected for
quantification and five data points were generated in such a way that concentration of tar
compounds fell within those five points. The tar which was already dissolved in isopropyl
alcohol was further diluted 5 times with dichloromethane. A dilute tar sample was injected into
the column and each sample was injected twice. Splitless injection was selected. The injector and
the GC/MS interface were kept at constant temperature of 280oC and 250oC, respectively. The
initial temperature of the column, 40oC, was maintained for 2 min and the temperature was
subsequently increased to 250oC at 5oC/min and the final temperature was held for 8 min.
Helium of ultra high purity (99.99%) supplied from Airgas Inc. (Charlotte, NC) was used as a
carrier gas and flowed at 1.25 mL/min. Compounds were ionized at 69.9 eV electron impact
conditions and analyzed over a mass per change (m/z) range of 50 550. Tar compounds were
identified by comparing the mass spectra with the NIST (National Institute of Standards and
Technology) mass spectral library and were reported as mg/Nm3 of syngas flow rate.
112
4.4 RESULTS AND DISCUSSIONS
Table 4.3 shows the various tar compounds along with the amount obtained from the test
runs from the gasifier. Individual concentration of tar compounds from each run of the
experiment is attached in Appendix F (Table F.1 and F.2). The major constituents observed in tar
are similar to those observed by other similar studies. Bari et.al [12] reported toluene,
ethylbenzene, styrene and p-xylene as a major tar constituents in the syngas obtained from the
gasification of feedstocks such as almond shells and oak in a downdraft gasifier using air as a
gasifying medium. Similar results was reported by Yamazaki et al. [13] on the experimental
investigation of the effect of superficial velocity on tar concentration in downdraft gasifier using
fir wood chips as a feedstock. As expected, the majority of tar compounds observed in higher
proportions are tertiary condensed tar products due to thermal cracking inside the gasifier. Figure
others Toluene
29% 28%
Styrene
9%
Phenol
Indene
7% Naphthalene
9%
18%
113
Table 4.3 Quantification of tar constituents in syngas from experiments
Toluene 76.8-198.3
o/p-Xylene 9.3-111.6
Naphthalene 62.3-126.1
Phenol 6.9-67.2
Styrene 21.0-65.1
Indene 15.7-55.8
Ethylbenzene 2.5-25.0
Phenol, 3-methyl- 1.3-25.4
Benzofuran 8.5-24.9
Biphenylene 7.1-22.2
Benzofuran, 2-methyl- 0-23.8
Benzene, 1-ethenyl-3-methyl-; (m-Methylstyrene) 6.6-18.8
Naphthalene, 2-methyl- 5.1-16.2
Naphthalene, 1-methyl- 5.9-14.6
Biphenyl 2.6-10.1
Phenol, 2-methyl- 0.5-8.9
Naphthalene, 2-ethenyl-; (2-Vinylnaphthalene) 0.4-6.7
Furfural 0-4.0
Naphthalene, 1,8-dimethyl- 0.6-3.6
Naphthalene, 1,5-dimethyl- 0-3.6
Dibenzofuran 0.4-3.4
.alpha.-Methylstyrene 1.5-3.1
Benzene, 1-ethyl-2-methyl-; (2-Ethyltoluene) 0.6-3.0
Benzene, 1,2,3-trimethyl- 1.4-2.4
Phenol, 2,4-dimethyl- 0-2.4
Acenaphthene 0.3-2.1
Phenol, 3,5-dimethyl- 0-1.9
Naphthalene, 2,3-dimethyl- 0-1.4
Phenol, 3-ethyl- 0-1.3
Phenol, 4-ethyl- 0-1.0
Naphthalene, 1,8-dimethyl- 0-0.8
Total 340-680
(Data for each experiment is attached in Appendix F)
114
Tar concentration in syngas from this stratified downdraft gasifier was found to be 0.34-
0.68 g/Nm3. Dogru et.al [14] and Phuphukrat et.al [15] reported tar concentration of 6.37-8.38
g/Nm3 for throated and throat-less downdraft gasifier respectively while using sewage sludge as
a feedstock. In another study conducted in the similar type of downdraft gasifier used for these
current experiments, Wei [16] reported the tar concentration of 0.054 mg/Nm3 when using wood
chips as a feedstock. This might be due to the difference in a bulk density of wood pellets and
wood chips. Since wood pellets are more than three times denser than wood chips, temperature at
the core of wood pellets might be lower than that in the surface and thus, producing higher tar
concentration.
Figure 4.5 shows the effect of biomass flow rate upon tar concentration in stratified
downdraft gasifier. Tar concentration shows the increase with increase in biomass flow rate from
17.6 kg/hr and it is observed to be highest at the biomass flow rate of 23.1 kg/hr. After an
increase in biomass flow rate from 17.6 kg/hr, tar concentration decreases with increase in
biomass flow and again increases after the biomass flow rate reaches to 26.5 kg/hr.
800.0
Tar concentration ( mg/m3)
600.0
400.0
200.0
0.0
15.0 18.0 21.0 24.0 27.0 30.0
Biomass flow rate ( kg/hr)
115
From Figure 4.5, it can be seen that biomass flow rate affects tar concentration in syngas
from a downdraft gasifier. However in this case, automatic temperature adjustment done by the
gasifier (already described in Chapter 3) injecting secondary air makes it difficult to analyze the
above tar concentration pattern due to the sole effect of biomass flow rate.
Tar concentration from the stratified downdraft gasifier was mostly condensed tertiary
mg/Nm3), and Indene (15.7-55.8 mg/Nm3). This shows that the primary and secondary tar
cracking is very efficient in the current configuration of the downdraft gasifier. Tar concentration
was also significantly lower than those reported by other in similar studies in conventional
downdraft gasifiers. Also, the tar concentration from the gasification of pellets was found to be
significantly higher than those with the gasification wood chips in similar type of gasifier due to
higher bulk density. The low tar concentration present in syngas from the downdraft gasifier
used for experiment makes it suitable for many synthesis process and power generation with
4.6 REFERENCES
[1] T.A. Milne, R.J. Evan, Biomass gasification "tars"; their nature, formation and conversion,
[2] R. Warnecke, Gasification of biomass: comparison of fixed bed and fluidized bed gasifier,
116
[3] L. Devi, K.J. Ptasinski, F.J.J.G. Janssen, A review of the primary measures for tar elimination
[4] S. Monteiro Nunes, N. Paterson, D.R. Dugwell, R. Kandiyoti, Tar formation and destruction
in a simulated downdraft, fixed-bed gasifier: reactor design and initial results, Energy & Fuels,
21 (2007) 3028-3035.
[5] J. Han, H. Kim, The reduction and control technology of tar during biomass
397-416.
[6] X.T. Li, J.R. Grace, C.J. Lim, A.P. Watkinson, H.P. Chen, J.R. Kim, Biomass gasification in
[7] E.G. Baker, M.D. Brown, D.C. Elliott, L.K. Mudge, Characterization and treatment of tars
from biomass gasifiers, in: Summer National Meeting, AIChE, Denver CO, 1988, pp. 111.
[8] P.M. Lv, Z.H. Xiong, J. Chang, C.Z. Wu, Y. Chen, J.X. Zhu, An experimental study on
[10] P. Hasler, T. Nussbaumer, Gas cleaning for IC engine applications from fixed bed biomass
[11] C. Li, K. Suzuki, Tar property, analysis, reforming mechanism and model for biomass
biomass in a downdraft fixed bed: a comparative study of the inorganic and organic products
117
[13] T. Yamazaki, H. Kozu, S. Yamagata, N. Murao, S. Ohta, S. Shiya, T. Ohba, Effect of
superficial velocity on tar from downdraft gasification of biomass, Energy & Fuels, 19 (2005)
1186-1191.
[14] M. Dogru, A. Midilli, C.R. Howarth, Gasification of sewage sludge using a throated
downdraft gasifier and uncertainty analysis, Fuel Processing Technology, 75 (2002) 55-82.
Characterization of tar content in the syngas produced in a downdraft type fixed bed gasification
[16] L. Wei, Experimental study on the effects of operational parameters of a downdraft gasifier,
Mississippi, 2005.
118
CHAPTER 5
predicting syngas composition depending upon the change in various operating parameters. The
linear equations derived to predict syngas composition based on knowledge of ultimate analysis
and moisture content of biomass is a significant achievement that can be applied to the
gasification process to find the upper limit of syngas production from an existing plant. The
effect of moisture content and temperature is also studied through the equilibrium model, which
gasifier. The results obtained, however are impacted significantly by the automatic factory pre-
set condition of the gasifier. Though clear effect of one parameter alone could not be seen from
the results, valuable information is provided about the syngas composition and temperature
distribution inside the gasifier under variable operating conditions. The final chapter which
discusses the tar content in syngas stream from the stratified downdraft gasifier is important
119
when it comes to the utilization of syngas for power generation and liquid-fuel synthesis from
syngas.
Future work can be done to develop a model that can accurately predict the syngas
composition depending on the configuration setting of the gasifier. Steady-state modeling of the
gasification process can be very useful for this purpose and can be used to predict the syngas
composition from the gasifier which cannot be predicted through the equilibrium model. With a
steady-state model, the effect of secondary air which is automatically injected inside the gasifier
can be more closely examined and validated. The steady state model can be also used to study
The experimental results obtained also shows some inconsistencies in temperature profile
of the gasifier. The temperature inside the gasifier, currently measured by the thermocouples
around the wall of the gasifier might not represent the true gasification temperature. Additional
temperature probes can be used to measure the temperature in the center of the reactor which will
give the better picture of temperature distribution inside the gasifier. Lack of temperature
uniformity has always remained one of the major problems of fixed bed gasification which can
be more closely studied with the addition of temperature measurement devices in the center of
the reactor.
One of the problems encountered during the operation of current gasifier used was the
limitations imposed by the automatic gasifier control system. Thus, effect of single parameter on
the downdraft gasification process could not be studied in detail. The modification of this current
configuration to accommodate the study of individual parameter or the design and fabrication of
120
new gasifier with the possibility of manual adjustments of various parameters will be immensely
helpful.
121
APPENDIX A
CONDITIONS
%THIS PROGRAM IS SET TO GIVE SYNGAS COMPOSITION IN DRY SYNGAS BASIS. IF FOR
%SOME REASON YOU WANT TO CHANGE, PLEASE MODIFY IN LINE 113 AND 114.
function[final_frac_comp]=eq_comp_model_gen(g_temp,ele_comp)
format short
tol=0.0001;
maxit=100;
%disp('elemental composition should be of the form [C, H, O, N, Ash]');
%ele_comp=input('Enter elemental composition of biomass: ');
%disp('Initial guess is of the form [H2 CO CO2 H2O CH4 3.76N2] ')
xx0=[0.1,0.1,0.1,0.1,0.1,0.1]';
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Heat of formation of different compounds at 25 C, kJ/kmol
H_f_H2O_g=-241818;H_f_H2O_l=-285830;H_f_CO2=-393509;H_f_CO=-110525;
H_f_CH4=-74520;H_f_H2=0;H_f_O2=0;H_f_N2=0;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Function for finding sensible heat for various gases
%constants
C_p_H2O=[32.24 0.1923e-2 1.055e-5 -3.595e-9];
C_p_H2=[29.11 -0.1916e-2 0.4003e-5 -0.8704e-9];
C_p_CO=[28.16 0.1675e-2 0.5372e-5 -2.222e-9];
C_p_CO2=[22.26 5.981e-2 -3.501e-5 -7.469e-9];
C_p_CH4=[19.89 5.204e-2 1.269e-5 -11.01e-9];
C_p_N2=[28.90 -0.1571e-2 0.8081e-5 -2.873e-9];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%finding general equations for calculating k1 and k2
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
G_CO=[3.376 0.557e-3 0 -0.031e5 -110525 -137169];
G_CO2=[5.457 1.045e-3 0 -1.157e5 -393509 -394359];
G_H2O=[3.470 1.450e-3 0 0.121e5 -241818 -228572];
G_H2=[3.249 0.422e-3 0 0.083e5 0 0];
G_C=[1.771 0.771e-3 0 -0.867e5 0 0];
G_CH4=[1.702 9.081e-3 -2.164e-6 0 -74520 -50460];
delta_ws_final=[];
delta_meth_final=[];
for iii=1:6
delta_ws=G_H2(iii)+G_CO2(iii)-G_CO(iii)-G_H2O(iii);
delta_meth=G_CH4(iii)-G_C(iii)-2*G_H2(iii);
delta_ws_final=[delta_ws_final delta_ws];
122
delta_meth_final=[delta_meth_final delta_meth];
end
T_0=298;
k1=exp(-((delta_meth_final(6)-
delta_meth_final(5))/(8.314*298.15)+(delta_meth_final(5)/(8.314*g_temp))...
+(int_eq_sp2(delta_meth_final,g_temp)/g_temp)-
int_eq_sp1(delta_meth_final,g_temp)));
k2=exp(-((delta_ws_final(6)-
delta_ws_final(5))/(8.314*298.15)+(delta_ws_final(5)/(8.314*g_temp))...
+(int_eq_sp2(delta_ws_final,g_temp)/g_temp)-
int_eq_sp1(delta_ws_final,g_temp)));
%%%%%%%%%%%%%%%%%%%%%%%%%%%%Continued below%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%function for calculating int_eq_sp
function int_for_gibbs_difff1= int_eq_sp1(difff,g_temp)
tau=g_temp/298.15;
int_for_gibbs_difff1=difff(1).*log(tau)+((difff(2).*T_0+...
((difff(3)*T_0^2+(difff(4)/(tau^2.*T_0^2)))*((tau+1)/2)))*(tau-
1));
end
function int_for_gibbs_difff2= int_eq_sp2(var_sp,g_temp)
tau=g_temp/298.15;
int_for_gibbs_difff2=var_sp(1).*T_0*(tau-1)+...
var_sp(2)*0.5*T_0^2*(tau^2-1)+var_sp(3)*T_0^3*(tau^3-1)/3+...
var_sp(4)*(tau-1)/(tau*T_0);
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%finding lambda and gamma for below calculation%%%%%%%%%%%%%%%%%%%%%%%%%
norm_1_C=ele_comp(1)/(12);
norm_1_H=ele_comp(2)/(1.008);
norm_1_O=ele_comp(3)/(16);
norm_1_N=ele_comp(4)/(14.007);
lambda=norm_1_H/norm_1_C;
gamma=norm_1_O/norm_1_C;
beta=norm_1_N/norm_1_C;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
x0=xx0;
iter=1;
iter_m=1;
sol_final=[];
%w=M_fs*Moisture_Content/(18*(1-Moisture_Content));
w=linspace(0,1,15);
%w=0;
Moisture_Content=[];
M_fs=12+lambda*1.008+gamma*16;
for N=1:length(w)
Moisture_Content=[Moisture_Content 18*100*w(N)/(M_fs+18*w(N))];
end
Moisture_Content
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Main Loop for solving the equations of interests
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for iter_m=1:length(w)
while(iter<=maxit)
y=-df1(x0)\f1(x0);
123
xn=x0+y;
err=max(abs(xn-x0));
if(err<=tol)
x=xn;
else
x0=xn;
end
iter=iter+1;
end
iter=1;
sol_temp=x;
%%%%%%%%%%%%%%%%%%%Continued below%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
sol_final=[sol_final sol_temp];
iter_m=iter_m+1;
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Multiplying m with 3.76 to get correct N2 mols
p=length(w);
frac_N2=[];
for l=1:p
frac_N2=[frac_N2 sol_final(6,l)*3.76];
end
final_comp=[sol_final(1:5,1:p);frac_N2]
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%fin_rep=input('Do you want to find syngas composition in dry syngas
basis(y/n): ','s');
fin_rep='y';
if fin_rep=='n';
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%finding total amount of product gas for each moisture content
total_frac_m=[];
for n=1:p
total_frac_m=[total_frac_m sum(final_comp(1:6,n))];
end
total_frac_m; %sum of all product gases
else
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%finding total amount of product gas on dry basis for each moisture
%content
dry_final_comp=final_comp;
dry_final_comp(4,:)=[];
total_frac_m=[];
for n=1:p
total_frac_m=[total_frac_m sum(dry_final_comp(1:5,n))];
end
total_frac_m;
final_comp=dry_final_comp;
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%expressing all the components in molar fraction or volumetric fraction
final_frac_comp=[];
for MM=1:length(total_frac_m)
124
final_frac_m=[];
if fin_rep=='y'
l_in=length(xx0)-1;
else
l_in=length(xx0);
end
for NN=1:l_in
final_frac_m=[final_frac_m;final_comp(NN,MM)/total_frac_m(MM)];
end
final_frac_comp=[final_frac_comp final_frac_m];
end
%%%%%%%%%%%%%%%%%%%Continued below%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function f=f1(X)
x_1=X(1); x_2=X(2); x_3=X(3);x_4=X(4);x_5=X(5);m=X(6);
val_1=x_2+x_3+x_5-1;
val_2=x_1+x_4+2*x_5-w(iter_m)-(lambda/2);
val_3=x_2+2*x_3+x_4-2*m-gamma-w(iter_m);
val_4=-k1*x_1^2+(x_5*(x_1+x_2+x_3+x_4+x_5+3.76*m));
val_5=x_2*x_4*k2-x_1*x_3;
val_6=x_1*t_en_gas(H_f_H2, C_p_H2, g_temp)+...
x_2*t_en_gas(H_f_CO, C_p_CO, g_temp)+...
x_3*t_en_gas(H_f_CO2, C_p_CO2,g_temp)+...
x_4*t_en_gas(H_f_H2O_g, C_p_H2O, g_temp)+...
x_5*t_en_gas(H_f_CH4, C_p_CH4, g_temp)+...
3.76*m*t_en_gas(H_f_N2, C_p_N2, g_temp)-...
heat_bio(ele_comp)-w(iter_m)*(H_f_H2O_l+1000);
f=[val_1; val_2;val_3;val_4;val_5;val_6];
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function df=df1(X);
x_1=X(1); x_2=X(2); x_3=X(3);x_4=X(4);x_5=X(5);m=X(6);
df=[0,1,1,0,1,0;1 0 0 1 2 0; 0 1 2 1 0 -2;-
2*x_1*k1+x_5,x_5,x_5,x_5,2*x_5+(x_1+x_2+x_3+x_4+3.76*m),3.76*x_5; -x_3
k2*x_4...
-x_1 k2*x_2 0 0; t_en_gas(H_f_H2, C_p_H2, g_temp)...
t_en_gas(H_f_CO, C_p_CO, g_temp) t_en_gas(H_f_CO2, C_p_CO2, g_temp)...
t_en_gas(H_f_H2O_g, C_p_H2O, g_temp) t_en_gas(H_f_CH4, C_p_CH4,
g_temp)...
3.76*t_en_gas(H_f_N2, C_p_N2, g_temp)];
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function dh_comp=t_en_gas(H_for, sp_heat, temp)
heat_coeff=sp_heat;
dh_comp=H_for+quad(@sensible,298,temp);
function sens_heat=sensible(t)
sens_heat=heat_coeff(1)+heat_coeff(2).*t+heat_coeff(3).*t.^2+heat_coeff(4).*t
.^3;
end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Standard heat of formation of various biomass
%H_f_XX is the heat of formation of XX compound, units in kJ/kmol
%LHV is lower heating value of biomass, kJ/kg
125
%LHV_mol is the lower heating value of biomass, kJ/kmol
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function H_f_bio= heat_bio(comp)
LHV=4.187*(81*comp(1)+300*comp(2)-26*comp(3)-54*comp(2));
LHV_mol=LHV*(12+lambda*1.008+gamma*16);
H_f_bio=(lambda/2)*H_f_H2O_l+H_f_CO2+LHV_mol;
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end
126
var_sp(2)*0.5*T_0^2*(tau^2-1)+var_sp(3)*T_0^3*(tau^3-1)/3+...
var_sp(4)*(tau-1)/(tau*T_0);
end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
127
APPENDIX B
EQUIVALENCE RATIO
The functions used for calculating the enthalpy change and equilibrium constants are
same as that for adiabatic condition which is already mentioned in appendix A. However, the
%THIS PROGRAM IS SET TO GIVE SYNGAS COMPOSITION IN DRY SYNGAS BASIS. IF FOR
%SOME REASON YOU WANT TO CHANGE, PLEASE MODIFY IN LINE 113 AND 114.
function[final_comp]=eq_model_const(g_temp,ele_comp,m)
format short
tol=0.00001;
maxit=100;
%disp('elemental composition should be of the form [C, H, O, N, Ash]');
%ele_comp=input('Enter elemental composition of biomass: ');
%disp('Initial guess is of the form [H2 CO CO2 H2O CH4] ')
xx0=[0.1,0.1,0.1,0.1,0.1]';
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Heat of formation of different compounds at 25 C, kJ/kmol
H_f_H2O_g=-241818;H_f_H2O_l=-285830;H_f_CO2=-393509;H_f_CO=-110525;
H_f_CH4=-74520;H_f_H2=0;H_f_O2=0;H_f_N2=0;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Function for finding sensible heat for various gases
%constants
C_p_H2O=[32.24 0.1923e-2 1.055e-5 -3.595e-9];
C_p_H2=[29.11 -0.1916e-2 0.4003e-5 -0.8704e-9];
C_p_CO=[28.16 0.1675e-2 0.5372e-5 -2.222e-9];
C_p_CO2=[22.26 5.981e-2 -3.501e-5 -7.469e-9];
C_p_CH4=[19.89 5.204e-2 1.269e-5 -11.01e-9];
C_p_N2=[28.90 -0.1571e-2 0.8081e-5 -2.873e-9];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%finding general equations for calculating k1 and k2
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
G_CO=[3.376 0.557e-3 0 -0.031e5 -110525 -137169];
G_CO2=[5.457 1.045e-3 0 -1.157e5 -393509 -394359];
G_H2O=[3.470 1.450e-3 0 0.121e5 -241818 -228572];
G_H2=[3.249 0.422e-3 0 0.083e5 0 0];
G_C=[1.771 0.771e-3 0 -0.867e5 0 0];
G_CH4=[1.702 9.081e-3 -2.164e-6 0 -74520 -50460];
delta_ws_final=[];
128
delta_meth_final=[];
for iii=1:6
delta_ws=G_H2(iii)+G_CO2(iii)-G_CO(iii)-G_H2O(iii);
delta_meth=G_CH4(iii)-G_C(iii)-2*G_H2(iii);
delta_ws_final=[delta_ws_final delta_ws];
delta_meth_final=[delta_meth_final delta_meth];
end
T_0=298.15;
%%%%%%%%%%%%%%%%%%%%%%%Continued below%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
k1=exp(-((delta_meth_final(6)-
delta_meth_final(5))/(8.314*298.15)+(delta_meth_final(5)/(8.314*g_temp))...
+(int_eq_sp2(delta_meth_final,g_temp)/g_temp)-
int_eq_sp1(delta_meth_final,g_temp)));
k2=exp(-((delta_ws_final(6)-
delta_ws_final(5))/(8.314*298.15)+(delta_ws_final(5)/(8.314*g_temp))...
+(int_eq_sp2(delta_ws_final,g_temp)/g_temp)-
int_eq_sp1(delta_ws_final,g_temp)));
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%function for calculating int_eq_sp
function int_for_gibbs_difff1= int_eq_sp1(difff,g_temp)
tau=g_temp/298.15;
int_for_gibbs_difff1=difff(1).*log(tau)+((difff(2).*T_0+...
((difff(3)*T_0^2+(difff(4)/(tau^2.*T_0^2)))*((tau+1)/2)))*(tau-
1));
end
function int_for_gibbs_difff2= int_eq_sp2(var_sp,g_temp)
tau=g_temp/298.15;
int_for_gibbs_difff2=var_sp(1).*T_0*(tau-1)+...
var_sp(2)*0.5*T_0^2*(tau^2-1)+var_sp(3)*T_0^3*(tau^3-1)/3+...
var_sp(4)*(tau-1)/(tau*T_0);
end
%k1=9.72e-02;k2=1.4561;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%finding lambda and gamma for below calculation%%%%%%%%%%%%%%%%%%%%%%%%%
norm_1_C=ele_comp(1)/(12);
norm_1_H=ele_comp(2)/(1.008);
norm_1_O=ele_comp(3)/(16);
norm_1_N=ele_comp(4)/(14.007);
lambda=norm_1_H/norm_1_C;
gamma=norm_1_O/norm_1_C;
beta=norm_1_N/norm_1_C;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
x0=xx0;
iter=1;
iter_m=1;
sol_final=[];
w=0;
Moisture_Content=[];
for N=1:length(w)
Moisture_Content=[Moisture_Content 18*100*w(N)/(24+18*w(N))];
end
Moisture_Content;
129
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Main Loop for solving the equations of interests
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for iter_m=1:length(w)
while(iter<=maxit)
y=-df1(x0)\f1(x0);
xn=x0+y;
err=max(abs(xn-x0));
if(err<=tol)
x=xn;
else
%%%%%%%%%%%%%%%%%%%%%%%Continued below%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
x0=xn;
end
iter=iter+1;
end
iter=1;
sol_temp=x;
sol_final=[sol_final sol_temp];
iter_m=iter_m+1;
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Multiplying m with 3.76 to get correct N2 mols
p=length(w);
frac_N2=[];
for l=1:p
frac_N2=[frac_N2 m*3.76];
end
final_comp=[sol_final(1:5,1:p);frac_N2];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%fin_rep=input('Do you want to find syngas composition in dry syngas
basis(y/n): ','s');
fin_rep='n';
if fin_rep=='n';
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%finding total amount of product gas for each moisture content
total_frac_m=[];
for n=1:p
total_frac_m=[total_frac_m sum(final_comp(1:6,n))];
end
total_frac_m; %sum of all product gases
else
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%finding total amount of product gas on dry basis for each moisture
%content
dry_final_comp=final_comp;
dry_final_comp(4,:)=[];
total_frac_m=[];
for n=1:p
total_frac_m=[total_frac_m sum(dry_final_comp(1:5,n))];
end
total_frac_m;
final_comp=dry_final_comp;
130
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%expressing all the components in molar fraction or volumetric fraction
final_frac_comp=[];
for MM=1:length(total_frac_m)
final_frac_m=[];
if fin_rep=='y'
l_in=length(xx0);
else
l_in=length(xx0)+1;
end
%%%%%%%%%%%%%%%%%%%%%%%Continued below%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for NN=1:l_in
final_frac_m=[final_frac_m;final_comp(NN,MM)/total_frac_m(MM)];
end
final_frac_comp=[final_frac_comp final_frac_m];
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function f=f1(X)
x_1=X(1); x_2=X(2); x_3=X(3);x_4=X(4);x_5=X(5);
val_1=x_2+x_3+x_5-1;
val_2=x_1+x_4+2*x_5-w(iter_m)-(lambda/2);
val_3=x_2+2*x_3+x_4-2*m-gamma-w(iter_m);
val_4=-k1*x_1^2+(x_5*(x_1+x_2+x_3+x_4+x_5+3.76*m));
val_5=x_2*x_4*k2-x_1*x_3;
f=[val_1; val_2;val_3;val_4;val_5];
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function df=df1(X)
x_1=X(1); x_2=X(2); x_3=X(3);x_4=X(4);x_5=X(5);
df=[0,1,1,0,1;1 0 0 1 2; 0 1 2 1 0;-
2*x_1*k1+x_5,x_5,x_5,x_5,2*x_5+(x_1+x_2+x_3+x_4+3.76*m); -x_3 k2*x_4...
-x_1 k2*x_2 0];
end
end
131
APPENDIX C
FORMULA DERIVATION
132
Moisture-free elemental composition Dry syngas composition
C H O N Ash CO CO2 CH4 H2
49.5 6.0 40.6 0.5 3.5 18.7 7.6 0.2 18.7
49.7 5.8 41.5 0.7 2.3 18.9 7.4 0.2 18.9
50.6 6.1 41.6 0.5 1.3 19.7 6.9 0.2 19.7
50.0 6.2 41.1 0.5 2.2 19.7 7.2 0.2 19.7
49.9 5.7 42.3 0.1 2.1 18.7 7.5 0.2 18.7
49.9 5.9 43.5 0.1 0.7 19.7 7.3 0.2 19.7
49.5 6.3 42.0 0.5 1.8 20.1 7.3 0.2 20.1
51.7 4.5 35.1 0.2 8.5 13.4 8.0 0.1 13.4
49.4 5.8 42.3 0.2 2.4 18.8 7.7 0.2 18.8
49.9 6.0 41.2 0.2 2.7 19.1 7.4 0.2 19.1
50.0 6.2 39.6 0.2 4.1 18.7 7.5 0.2 18.7
50.7 6.4 41.8 0.3 1.0 20.4 6.8 0.2 20.4
49.5 6.2 41.7 0.2 2.4 19.7 7.5 0.2 19.7
49.7 6.2 43.8 0.3 0.1 20.6 7.1 0.2 20.6
49.8 5.7 39.8 0.3 4.4 17.7 7.8 0.1 17.7
49.4 6.1 43.0 0.2 1.3 20.0 7.4 0.2 20.0
48.4 6.0 41.6 0.2 3.8 18.9 8.1 0.2 18.9
49.0 5.5 39.2 0.2 6.2 16.7 8.4 0.1 16.7
48.5 5.8 41.2 0.2 4.3 18.3 8.2 0.1 18.3
48.2 5.7 41.6 0.2 4.3 18.0 8.4 0.1 18.0
46.0 5.9 41.4 0.9 5.9 18.1 9.2 0.1 18.1
47.0 5.5 41.1 0.7 5.7 17.2 9.0 0.1 17.2
46.5 5.8 40.4 0.6 6.7 17.6 9.1 0.1 17.6
46.0 5.4 39.2 0.6 8.7 16.1 9.7 0.1 16.1
46.7 5.5 40.6 0.6 6.5 17.0 9.2 0.1 17.0
44.8 5.5 37.7 0.7 11.3 15.5 10.2 0.1 15.5
47.0 5.7 40.7 0.6 6.0 17.6 8.9 0.1 17.6
46.5 6.1 40.1 0.7 6.5 18.3 8.9 0.1 18.3
46.3 5.6 41.0 0.7 6.5 17.2 9.3 0.1 17.2
47.1 5.8 37.5 0.7 8.9 16.5 9.1 0.1 16.5
48.0 5.8 37.0 0.7 8.5 16.6 8.7 0.1 16.6
46.8 5.5 38.4 0.7 8.7 16.1 9.4 0.1 16.1
47.0 5.7 41.4 0.7 5.3 17.8 8.9 0.1 17.8
48.8 5.5 42.3 1.0 2.4 18.4 7.8 0.2 18.4
49.4 5.2 39.5 1.1 4.8 16.6 8.0 0.1 16.6
49.4 5.8 39.6 1.3 4.0 18.2 7.6 0.1 18.2
46.5 5.6 41.9 1.2 4.9 17.9 9.0 0.1 17.9
49.5 5.6 37.4 1.1 6.5 16.7 8.0 0.1 16.7
49.3 5.9 42.8 0.7 1.3 19.6 7.4 0.2 19.6
Continued in the next page
133
Moisture-free elemental composition Dry syngas composition
C H O N Ash CO CO2 CH4 H2
50.3 5.6 40.9 0.7 2.6 18.1 7.3 0.1 18.1
49.7 5.6 42.6 0.6 1.6 18.7 7.4 0.2 18.7
50.7 5.5 35.4 0.8 7.7 15.8 7.9 0.1 15.8
50.4 5.9 38.1 0.9 4.7 17.8 7.4 0.1 17.8
49.5 5.7 36.0 0.8 8.0 16.2 8.2 0.1 16.2
49.6 5.5 42.3 0.7 2.0 18.4 7.5 0.2 18.4
47.3 5.3 41.6 0.5 5.3 17.0 9.0 0.1 17.0
47.3 5.6 41.1 0.7 5.3 17.6 8.8 0.1 17.6
47.6 5.6 40.2 0.6 6.0 17.3 8.7 0.1 17.3
47.8 5.6 39.2 0.7 6.7 17.0 8.7 0.1 17.0
48.0 5.7 40.0 0.7 5.6 17.6 8.5 0.1 17.6
48.0 5.6 39.0 0.5 6.9 16.8 8.7 0.1 16.8
48.5 5.5 38.2 0.6 7.1 16.5 8.6 0.1 16.5
46.7 5.6 41.5 0.4 5.8 17.6 9.1 0.1 17.6
46.7 5.7 42.1 0.6 4.9 18.1 8.9 0.1 18.1
46.9 5.5 42.0 0.6 5.0 17.7 9.0 0.1 17.7
46.6 5.6 41.2 0.6 6.0 17.5 9.1 0.1 17.5
47.0 5.4 41.1 0.6 5.9 17.0 9.1 0.1 17.0
46.7 5.6 41.0 0.5 6.3 17.3 9.2 0.1 17.3
46.6 5.7 41.5 0.6 5.7 17.8 9.1 0.1 17.8
47.6 5.6 41.4 0.2 5.3 17.7 8.8 0.1 17.7
49.7 5.9 41.9 0.1 2.5 18.9 7.6 0.2 18.9
50.3 6.0 42.1 0.0 1.6 19.4 7.2 0.2 19.4
43.9 5.3 38.8 0.6 11.5 15.1 10.7 0.1 15.1
45.4 5.4 31.0 1.0 15.9 12.9 10.5 0.1 12.9
35.0 4.4 21.3 2.8 35.4 7.2 14.2 0.0 7.2
45.4 5.9 35.9 0.9 11.4 15.9 9.8 0.1 15.9
39.7 5.8 27.2 0.8 26.1 11.1 12.3 0.1 11.1
49.8 5.5 42.4 0.5 1.8 18.5 7.5 0.2 18.5
50.4 5.7 40.6 0.5 2.8 18.3 7.3 0.2 18.3
134
APPENDIX D
The data used for figures in different chapters is reported in this appendix. Data in each
table corresponds to the figure mentioned alongside.
Moisture
H2 CO CO2 CH4 N2 HHV Eq.
Content
(% vol.) (% vol.) (% vol.) (%vol.) (%vol.)((MJ/m3) ratio
(% wet basis)
0.0 16.9 23.2 9.4 0.1 50.2 5.1 0.40
5.1 17.3 21.5 10.6 0.1 50.4 5.0 0.41
9.7 17.6 19.9 11.6 0.1 50.7 4.8 0.41
13.8 17.8 18.5 12.5 0.1 51.0 4.7 0.42
17.6 17.9 17.2 13.3 0.1 51.4 4.5 0.43
21.1 17.9 16.1 14.0 0.1 51.8 4.4 0.44
24.3 17.9 15.0 14.7 0.1 52.3 4.2 0.45
27.3 17.8 14.0 15.3 0.1 52.8 4.1 0.47
30.0 17.7 13.1 15.8 0.1 53.3 4.0 0.48
32.5 17.5 12.2 16.3 0.1 53.8 3.8 0.49
34.9 17.3 11.5 16.7 0.1 54.4 3.7 0.50
37.1 17.1 10.7 17.1 0.1 54.9 3.6 0.51
39.1 16.8 10.1 17.5 0.1 55.5 3.5 0.52
41.1 16.5 9.4 17.8 0.1 56.1 3.3 0.53
42.9 16.2 8.9 18.1 0.1 56.7 3.2 0.54
135
Table D.2 Data for Figure 2.3
H2 CO CO2 CH4 N2
Temperature HHV(MJ/m3)
(% vol.) (%vol.) (% vol.) (% vol.) ( % vol.)
136
Table D.3 Data for Figure 2.4
137
Table D.4 Data for Figure 2.5
138
Table D.5 Data for Figure 2.6
H2 CO CO2 CH4 N2
Temperature HHV(MJ/m3)
(% vol.) (% vol.) (% vol.) (% vol.) (% vol.)
139
Table D.6 Data for Figure 2.7
140
Table D.7 Snapshot of temperature of one typical run in the gasifier
Figure D.1 to D.5 represents the temperature variations inside the gasifier from the initial
start-up to the steady state at which temperature almost remains constant.
141
1200
1000
Temperature ( oC)
800
600
400
200
0
6:43 9:07 11:31
Time (A.M)
Figure D.1 Temperature recorded by thermocouple at T1 from its start-up to steady state
1000
800
Temperature (C)
600
400
200
0
6:43 9:07 11:31
Time (A.M)
Figure D.2 Temperature recorded by thermocouple at T2 from its start-up to steady state
1000
Temperature (C)
800
600
400
200
0
6:43 9:07 11:31
Time (A.M)
Figure D.3 Temperature recorded by thermocouple at T3 from its start-up to steady state
142
1000
Temperature (C)
800
600
400
200
0
6:43 9:07 11:31
Time (A.M)
Figure D.4 Temperature recorded by thermocouple at T4 from its start-up to steady state
800
Temperature (C)
600
400
200
0
6:43 9:07 11:31
Time (A.M)
Figure D.5 Temperature recorded by thermocouple at T1 from its start-up to steady state
Moisture Content
T1(C) T2(C) T3(C) T4(C) Grate(C)
(% wet basis)
19.6 89746 86739 83934 83710 8118
23 91739 86727 84429 84414 7928
25 88149 86020 83619 82919 8102
143
Table D.9 Data for Figure 3.6
144
Table D.12 Data for Figure 3.10
145
APPENDIX E
SAMPLE CALCULATIONS
The input for the program is elemental composition of the feedstock and reaction
temperature. Program automatically calculates the molecular formula for the biomass in the form
of CHxOyNz.
The system of equations mentioned in Chapter 2 is solved using Newtons Jacobi method
using these stoichiometric numbers and accessing various function files for finding equilibrium
constants and other thermodynamic properties of various gases involved in the gasification
process.
Solving these equations, the number of moles of H2O, CO, CO2, CH4 and N2 required for
the gasification process is obtained. The fraction of gases can be either expressed as wet syngas
content. Table below shows the typical data when the program was run with elemental
146
Table E.1 Calculation of syngas composition from MATLAB
The calculation procedure for other tables in Chapter 2 are very similar to that explained above
147
E.2 CARBON, ENERGY AND EXERGY ANALYSES
Table E.2 shows the syngas composition at different biomass flow rates along with the
moisture content, grate temperature and actual syngas flow rate upon the gasification of
commercial wood pellets. This is also the supplemental data from the experiments with
Table E.2 Syngas composition at different biomass flow rate for commercial wood pellets
The experiment done with a biomass flow rate of 28.8 kg/hr is selected for sample calculation
purpose (highlighted above). Sample calculations of carbon closure, energy ratio and exergy
148
E.2.1 CARBON CLOSURE
Carbon closure for the experiment taken for sample calculation is 0.90 from Table 3.9 for the
The dry biomass flow rate (mdry) is calculated by subtracting the amount of moisture present
which is:
where, are mass of dry biomass flow rate, wet biomass flow rate and amount of
The amount of carbon present in the biomass can be found by multiplying the dry biomass flow
rate with its carbon content which is 47.7% (Reported in Chapter 3-Table 3.1).
(1)
molar density of ideal gas, molecular weight of carbon in kg/mol, and syngas flow rate (m 3/hr),
149
respectively. For ideal gas condition, molar density at standard temperature and pressure (STP) is
44.615 mol/m3.
Substituting the respective values in Eqn. (1) gives the following result.
Energy ratio for the experiment taken for sample calculation is 0.87 is taken from Table 3.9 for
E.3.
150
Table E.3 Properties of syngas constituents
The energy content in biomass can be found using Eqn. (2). Higher heating value (HHV) of
(2)
The volumetric flow rate of CO, CH4 and H2 is expressed in the form of mass flow rate as
151
is obtained by difference as follows.
As mentioned in Chapter 3-Eqn. (10), the total energy of a gas is given by Eqn. (3).
(3)
where Ei and E0i are the total energy and chemical energy at dead
state temperature (Td) taken as 25C, respectively. Cp is the specific heat capacity (kJ/kg-K) of
syngas while T is the syngas temperature (taken as grate temperature) in Kelvin as highlighted
Table E.2.
The specific heat capacity for given temperature can be calculated from Eqn. (4) where
(4)
In above equation, a, b, c and d are the coefficients of specific heat capacity which is reported in
152
Integrating Eqn. (3) gives the specific energy of an individual gas. Then, multiplying the specific
energy of biomass at given temperature with biomass flow rate will give the total energy of
syngas as follows.
By adding the enthalpy of individual gases, we can find Eout, the total energy of the output gas.
Exergy ratio for the experiment taken for sample calculation is 0.58 from Table 3.9.
The chemical exergy associated with biomass can be found from Eqn. (5).
153
(5)
(6)
Exin is the chemical exergy of biomass. LHV and are the lower heating value (MJ/kg)
and ash free fraction of biomass, respectively. H, C, O and N represent the fraction of hydrogen,
carbon, oxygen and nitrogen present in the biomass respectively. H/C, N/C and O/C are
calculated from the ultimate analysis of the feedstock which is 47.7% C, 6.0% H, 45.8% O and
LHV was calculated according to the Eqn. (7) where HHV should be expressed as (kJ/kg).
(7)
is the fraction of ash in the feedstock which was found to be 0.44% in the commercial
wood pellets, as can be seen in Table 3.1. Now, Exin can be calculated from above relations.
154
(8)
where Ex0i and Exi ( are the chemical exergy and total exergy of
the individual gases. The chemical exergy of selected gases is given in Table E.2.
Multiplying the specific exergy of biomass at given temperature with biomass flow rate will give
(9)
The following exergy of individual gases can be obtained by substituting the value of specific
The addition of above exergy of individual gases gives the total exergy output per hour
155
Exergy ration is the ratio of Exin to Exout can be now calculated.
156
APPENDIX F
Table F.1 and F.2 shows the concentration of various tar constituents in syngas from the
gasification of commercial wood pellets in a downdraft gasifier for eleven experiments. The
experiments were conducted at different biomass flow rate which reported along with its
157
Table F.1 Concentration of tar constituents in syngas (Supplemental data-A)
158
Table F.2 Concentration of tar constituents in syngas (Supplemental data-B)
159
APPENDIX G
UNCERTAINTY ANALYSIS
(1)
A Taylor series expansion can be used to find the infinitesimal change in due to the
expression.
(2)
associated with can be expressed as the root-sum-square which is shown in the following
equation.
160
(3)
The above expression is used for calculating the uncertainty associated with finding mass, energy
The following equation was implemented to calculate the mass closure in the gasifier which
can be obtained by combining cin and cout from carbon closure sample calculation.
Following are the uncertainties associated with each parameter in above equation.
161
= 1 m3/hr of a scale reading, = 10-2 kg of a scale
162
Thus the absolute uncertainty associated with the mass closure is given by the following relation.
Substituting the values found above, the uncertainty associated with mass closure was found to
be 0.009.
HHV is the higher heating value of biomass (commercial wood pellets) which is expressed in the
units of MJ/kg.
This specific energy can be converted into hourly flow rate as following.
163
The following is the sum of total energy in the syngas from the gasification process.
above equation which is already defined in Chapter 3- Eqn. (12). The uncertainty associated with
temperature ( is .
164
Thus, it can be seen from above discussion that:
Finally, the uncertainty in enthalpy balance can be calculated using the following relation:
165
G.3 UNCERTAINTY ASSOCIATED WITH EXERGY RATIO
Exergy of the input mass in a gasifier can be found by following relation assuming zero
, and with integration and some simplification, the following expression can be
166
The ratio of output exergy to input exergy is .
Finally, the uncertainty in exergy balance can be calculated using the following relation:
167
Using those values and the above equation, uncertainty in exergy balance was determined to be
0.0246.
REFERENCES
[1] E.O. Doebelin, Measurement Systems: Application and Design, 5th ed., McGraw-Hill, 2003.
[2] K.J. Ptasinski, M.J. Prins, A. Pierik, Exergetic evaluation of biomass gasification, Energy, 32 (2007)
568-574.
168