ELSEVIER
PH: S0016-2361(96)00030-0
Fuel Vol. 75 No. 8, pp. 987-998, 1996 Copyright © 1996ElsevierScienceLtd Printed in Great Britain. All rights reserved 0016-2361/96 $I 5.00+ 0.00
Pyrolysis characteristics of biomass and biomass components K. Raveendran, Anuradda Ganesh and Kartic C. Khilar* Energy Systems Engineering, Department of Mechanica/ Engineering, *Department of Chemica/ Engineering, Indian Institute of Technology, Bombay--400 076, India e-mai/."
[email protected]/n (Received 22 May 1995; revised 2 January 1996) Biomass pyrolysis studies were conducted using both a thermogravimetric analyser and a packed-bed pyrolyser. Each kind of biomass has a characteristic pyrolysis behaviour which is explained based on its individual component characteristics. Studies on isolated biomass components as well as synthetic biomass show that the interactions among the components are not of as much significance as the composition of the biomass. Direct summative correlations based on biomass component pyrolysis adequately explain both the pyrolysis characteristics and product distribution of biomass. It is inferred that there is no detectable interaction among the components during pyrolysis in either the thermogravimetric analyser or the packedbed pyrolyser. However, ash present in biomass seems to have a strong influence on both the pyrolysis characteristics and the product distribution. Copyright © 1996 Elsevier Science Ltd. (Keywords: biomass; biomass components; pyrolysis)
For an agriculture-based economy like that of India, the prospect of being able to convert widely available biomass materials into various forms of fuel is most attractive. In every thermochemical conversion route, pyrolysis plays a vital role. Hence there is renewed interest in understanding the complex pyrolysis process. In the past few decades much work has been done to study the influence of the operating parameters, such as temperature, heating rate, pressure and residence time, on pyrolysis 1. However, there still remains a need to study the effect of the feedstock properties on the process. Different kinds of biomass, although consisting of the same major constituents, have different compositions. The influence of biomass composition and ash content and composition on pyrolysis characteristics is the object of the present study. This paper discusses the pyrolysis characteristics and product distribution of biomass and their correlation with those of the individual components. Most of the pyrolysis characteristics reported in literature are for woody materials. A few attempts have been made at correlating the pyrolysis characteristics of biomass with those of its constituents, viz. cellulose, hemicellulose, lignin and extractives 1. Thermal analysis curves for wood often exhibit three peaks and have led researchers to believe that the mechanism of wood pyrolysis is a superposition of the mechanisms of its components 2-7. Shafizadeh and McGinnis 8 argued that qualitatively the thermal behaviour of wood reflected the behaviour of its components. Beall 9 showed that the same major components isolated from wood by different methods behave
differently; and Roberts 1° showed that the structural properties of the components influence the pyrolysis characteristics, so he did not attempt to correlate the pyrolysis behaviour of biomass with that of its components. Other workers/1-13 attempted to correlate biomass pyrolysis characteristics with those of its components using simple overall pyrolysis kinetic models. Though good agreement was seen at a particular heating rate, specific for each biomass, the agreement was not satisfactory at other heating rates. This discrepancy was attributed to the influence of heating rate on char formation / . Antal 1, in his exhaustive review on pyrolysis of biomass, addressed the problem of 'whether the pyrolysis of lignocellulose can be represented as a simple superposition of the behaviour' of its components'. He proposed that a mathematical superposition of the components' t.g.a, curves should explain their interaction adequately. It is appropriate to mention here some of the work carried out on copyrolysis of various biomass and coal mixtures which give insight into the interaction of biomass during pyrolysis. Klose and Stuke 14 reported no interaction between coal and biomass during copyrolysis in t.g.a. However, Nikkhah et al. 15, in their detailed copyrolysis studies of various biomass and coal mixtures in a batch reactor, reported increased gas yields, as well as increased heating value and hydrocarbon content of the pyrolysis gases. McGee et al) 6 reported copyrolysis studies of the mixtures of poly(vinyl chloride) (PVC) and wood/straw, to simulate municipal solid waste pyrolysis char. They found that the interaction between PVC and wood/straw increased the char
Fuel 1996 Volume 75 Number 8
987
Pyrolysis characteristics of biomass and biomass components: K. Raveendran et al. was determined using standard geological rock analysis methods2°'2l; the results are presented in an earlier paper 22. All biomass samples taken for both t.g.a, and PBP studies were ground <250 #m. Isolated biomass components and synthetic biomass also were studied. Holocellulose, cellulose, hemicellulose, lignin and extractives were isolated from wood by methods described in the literature ~9. Synthetic biomass was prepared to simulate the results of original whole biomass (with and without ash).
yield but reduced the char reactivity. Copyrolysis studies conducted by Khan et al. ~7 on mixtures of coals and heavy petroleum residues and by Saxby and Sato ~8 on Australian oil shale and lignite, all in a packed-bed pyrolyser (PBP), revealed the prevalence of synergetic effects; they also showed that the initial composition of the feedstock mixture had a direct bearing on the product distribution and properties. On critical analysis of these studies, the interesting fact emerges that most of the studies conducted in a t.g.a. show no interaction among biomass components or between biomass and coal, whereas the studies conducted in packed-bed batch reactors show interaction in both cases. Given the range of knowledge reported above, the following objectives were set for the present study:
Synthetic biomass preparation. Synthetic biomass samples were prepared by mixing each of the individual biomass constituents (cellulose, xylan, lignin, extractives and ash) proportionately. The proportions of individual constituents were obtained from summative analysis. Standard cellulose, lignin and xylan samples were obtained commercially for this purpose. Extractives were isolated from each kind of biomass according to TAPPI standard Tllm. The ash was obtained by combusting the corresponding biomass in a muffle furnace. To simulate deashed biomass, ash was excluded from the mixture of other constituents.
(1) to investigate comprehensively the pyrolysis characteristics and product distribution of various kinds of biomass and their components in both t.g.a. and PBP, and to compare the results to attempt to reconcile the discrepant information and inferences reported in the literature; (2) to develop a mathematical model to predict the pyrolysis characteristics and product distribution of biomass from the characteristics of the individual components; (3) to study the influence of the composition of the ash as a parameter within (1) and (2).
Pyrolysis To fulfil the objectives of the study, pyrolysis experiments were conducted in both a t.g.a, and a PBP. In the t.g.a, experiments a small sample (5-10mg) was taken, thus reducing the effect of secondary reactions and heat and mass transfer on the product yield, whereas in PBP experiments a much larger sample was used (10-25 g), so that secondary reactions and other heat and mass transfer effects played greater roles.
EXPERIMENTAL Feedstock Biomass. Fourteen commonly available kinds of biomass in the region were selected for the present study. These included a variety of potential biomass fuels available from energy crops (wood) as well as residues of agricultural (husks, straw etc.) and food processing activities (bagasse, nut shells etc.). The chemical composition of these materials was determined by methods described in the literature ~9. The composition of ash GAS .EATER7 / "'" "J""
I~.'.~ ......
[ REACt0R ~ I!'"']
"~
'|;J. . . . . . . . ~ ' ~ I'.
Jel
v
I~ t / , r
-,, -~ ~. ' v " ,~* :
t
~1:1::'.~1
,,~'.: "..b-. o -
,0'EMPE"ATURE = * . , . 0 " E . It TEMPERATURE INDICATOR
, CERA.,C woo,
.0.,
..SU,AT.O. s .eatlnSEUEMe"T
,, WATER C,ReO,AT,O.P~MP
s CERAMIC TUBE
16 GAS FLOW METER
~','
I
//
J
- II
\\
'
ATMOSP.ER~,
~""-~_~_~
® Figure 1
988
.EDWATERTA.,
.... •
•1'
,--" ...........
, RO*AME,,R 3 PYROLYSISREACTOR
,ELECTR,CALFU..ACE' ,, CO.OE.SE" TURE
I'::~'111"J~
l : 1 ; ;':;i
~,-¢ . . . .
T THt=RMOCOUP!.E 8 CLEANING'PORT 1 PRESSURE REGULATOR 9 PRODUCT VAPOUR
w°",t*l')~l~ I; "b2~LI
I::,:~:,,I~I.~F;I-';.::.I ~ ~, ~. .~ tl'~" ~"
I F "H00K [ l J-*,-i
T
| e ~
I'~,~lff"l."i~:-.:~ -I:: ,il l
T.g.a. studies. A DuPont 9900 thermal analyser was used for all the pyrolysis studies. Dynamic t.g.a, studies were carried out at linear heating rate of 50Kmin -1, covering the temperature range from atmospheric to 1273K. All the studies were carried out in an inert atmosphere of flowing nitrogen (50 cm 3 min-1).
Schematic of the packed-bed pyrolyser test set-up
Fuel 1996 Volume 75 Number 8
@
Pyrolysis characteristics of biomass and biomass components." K. Raveendran et al. Table I Pyrolysischaracteristicsof biomass componentsin t.g.a.~
Yield (wt% daf)
Whatman cellulose Wood cellulose Alkali lignin Acid lignin Hemicellulose Xylan Extractives
Volatiles
Char
Max. rate (wt% K -l)
97.5 86.0 59.4 52.9 68.0 70.0 73.0
2.5 14.0 40.6 47.1 32.0 30.0 26.9
1.7 1.7 0.3 0.4 0.5 0.9 0.4
Temp. at max. rate (K)
Initial decomp. temp. (K)
Temp. at involution point (K)
682 667 695 693 573 584 668
573 573 413 473 448 463 393
713 697 773 700 550 500 575
a Heating rate 50 K min-1
PBP studies. Experiments were conducted in a packed-bed pyrolyser designed for the purpose, with provision for collecting the pyrolysis products. Figure 1 shows a diagram of the pyrolysis reactor set-up, consisting of an electrically heated stainless steel tubular pyrolysis reactor and another electrically heated stainless steel tubular gas heater. Nitrogen from the cylinder was first heated to the operating temperature in the gas heater and entered the pyrolysis reactor tangentially at the top, flushed the pyrolysis vapour and left through the bottom of the reactor through a double-walled glass condenser tube to a train of wash-bottles immersed in an ice bath. Cold water was circulated through the condenser. A programmable PID controller was used to control the temperatures of the reactor and heater. Thermocouples were incorporated at three places in the reactor and one in the sample boat to monitor the sample temperature. Five representative kinds of biomass were chosen for investigations in this reactor. Isothermal experiments 773K in an inert atmosphere of flowing nitrogen (500cm 3 min -1) were conducted for each of these five biomass samples, isolated biomass components and mixtures of the components. A known weight of ~ 1 0 - 2 5 g of powdered samples was taken in the stainless steel mesh 'boat' and placed in the furnace, which was already maintained at the predetermined temperature of 773 K with the nitrogen flow. The volatiles evolved were collected and quenched, and the non-condensable gases were passed through a flow meter. The experiments were continued until the evolution of gases ceased. The furnace was switched off and the flow of nitrogen was maintained until the sample attained room temperature. The char remaining in the furnace was carefully removed and weighed. The yield of liquids was obtained by weighing the wash-bottles before and after the experiment, and the yield of gas was obtained by difference. Repeatability was ensured by repeating the experiments until the same product yield fractions were obtained with a precision of 4-2%. RESULTS A N D DISCUSSION
Pyrolysis characteristics of biomass components in t.g.a. Figure 2 presents the integral and differential thermograms of the basic components. The pyrolysis data are
r
1,00
Untreated biomass, isolated biomass components and synthetic biomass were all subjected to these analyses.
r
If • • •
0,75
r
Acid l i g n i n AIk=tl o + llgnin 0 Extractivcs w h a t r n a n ce(lulos¢ Wood c e l l u l o s e Herni cellulose Xylon
g u ~ 0,50
o, 0.25
.==.
I=,
~
t
0.00
1.50
+ O o A v • •
DTG
1.25
1.00
Acid tignin Alkali lignin Extmctlves Whatmon cellutos~ Wood cellulose Hemi cellulose Xylan
E m
0.75
>
O,SO
Q
0'25
0.00 150
250
350
450 550 Temperature ( C )
650
750
850
Figure 2 T.g.a.and d.t.g, curvesof biomass components
presented in Table 1. Several observations can at once be highlighted: (1) cellulose decomposes within a narrow temperature range, 573-703 K; (2) cellulose decomposition rate is the highest and char yield the lowest; (3) Whatman cellulose has a higher rate as well as a narrower temperature range for decomposition than cellulose extracted from wood; (4) lignin on the other hand decomposes over a wider temperature range, 523-823 K; (5) the char yield from lignin is the highest, in the region of 45-50 wt%; (6) hemicellulose and xylan are thermally the most unstable and start decomposing at a much lower temperature than every other component, with ~30 wt% char yield;
Fuel 1996 V o l u m e 75 Number 8
989
Pyrolysis characteristics of biomass and biomass components: K. Raveendran et al. Table 2 Pyrolysis characteristics of whole biomass in t.g.a,a
Volatiles
Char
Max. rate (wt% K-l)
79.7 81.1 69.8 70.7 56.8 73.5 70.9 80.6 68.7 70.1 70.0 74.7 76.3 72.8
20.3 18.9 30.2 29.3 43.2 26.5 29.1 19.4 31.3 29.9 30.0 25.3 23.7 27.2
0.9 0.6 0.8 0.8 0.6 1.1 1.1 1.3 0.7 0.9 0.8 1.0 0.9 0.9
Yield (wt% daf)
Bagasse Cashewnut shell Coconut coir Coconut shell Coir pith Corn cob Corn stalks Cotton gin waste Groundnut shell Millet husk Rice husk Rice straw Subabul wood Wheat straw
Temp, at max. rate (K)
Initial decomp, temp. (K)
Temp. at involution point (K)
677 638 672 615 622 603 634 679 662 653 666 651 683 604
483 473 513 518 483 533 498 523 493 523 518 518 498 493
688 673 673 678 663 653 653 688 683 663 663 673 663 663
Heating rate 50 K min-1
Table 3 Distribution of volatiles released during biomass pyrolysis in t.g.a,a
Distribution of volatiles (temp. ranges)
Bagasse Cashewnut shell Coconut coir Coconut shell Coir pith Corn cob Corn stalks Cotton gin waste Groundnut shell Millet straw Rice husk Rice straw Subabul wood Wheat straw
Moisture Zone I <373K
Zone II 373-523 K
Zone III 523-623 K
Zone IV 623-773 K
Zone V >773 K
Total volatiles
5.6 1.3 8.3 7.4 14.6 5.5 5.4 7.4 8.0 8.3 7.7 9.7 6.2 6.5
3.2 4.7 2.2 2.8 1.0 2.5 2.3 1.7 4.4 4.2 6.2 5.1 3.6 3.2
27.2 37.7 25.0 30.3 22.0 28.6 27.2 27.6 20.9 30.6 22.7 24.1 24.2 36.0
41.8 27.2 35.1 33.1 30.6 38.9 30.4 42.8 36.7 30.4 37.5 40.3 44.5 24.3
4.5 11.5 7.6 4.4 3.2 3.5 11.0 8.6 6.8 4.9 8.2 5.2 4.0 6.3
79.7 81.1 69.9 70.7 56.8 73.5 70.9 80.7 68.7 70.1 70.0 74.7 76.3 72.8
a Moisture, wt% ash-free biomass; volatiles, wt% daf biomass; heating rate 50 K min-l
(7) extractives decompose in a similar way to lignin b u t at a slightly higher rate a n d at a slightly lower temperature. F r o m the above observations, a ' z o n a t i o n ' scheme can be envisaged, relating to the d e c o m p o s i t i o n status of each c o m p o n e n t . The 'zones' are postulated as: zone I: <373 K m a i n l y m o i s t u r e evolution zone II: 373-523 K extractives start d e c o m p o s i n g zone III: 523-623 K predominantly hemicellulose decomposition zone IV: 6 2 3 - 7 7 3 K m a i n l y cellulose a n d lignin decomposition zone V: >773 K m a i n l y lignin d e c o m p o s i t i o n
990
Fuel 1996 Volume 75 Number 8
Pyrolysis characteristics of biomass in t.g.a. The fourteen potential b i o m a s s fuels selected for this study were subjected to pyrolysis in the t.g.a. The results are presented in Table 2. I n accordance with the z o n a t i o n hypothesis, the volatiles released o n pyrolysis were classified as shown in Table 3. T h a t each k i n d of b i o m a s s has its o w n typical decomposition characteristics is evident from these tables. F o r example coir pith releases only 5 6 w t % as volatiles, whereas bagasse a n d c o t t o n gin waste each release ~ 8 0 w t % . The m a x i m u m rate o f decomposition of c o t t o n gin waste is m o r e t h a n twice that of coir pith. Even t h o u g h it is evident from Table 3 that m o s t of the volatiles are released in zone IV, the
Pyrolysis characteristics of biomass and biomass components. K. Raveendran et al. t.O0
J
O
'~ 0.75
.
?,
-
\ll I
~'.'N"~. I \ x .
Synthetic .
/~l I ,,}'\I
.
.
.
Untreated
..
COR.CO
.
\\,~
"
....
~i
Untreated
\ i~
,
, II
100 J~"
Synthetic
-Untreated
U
0.75 o m
RICE HUSK
I
0.50
WOOD
,ti',
o.0.50
.C Oa
'~'
0.25
0.25 > >
~___
--
0.00
I
, . 1.00
i
• A
8 o.75
--
- Untrzotzd
-
~\/,l
-
150
Synthetic
-
350 A
o
550
o
750
Temperature ( C )
Untreated
035
u
RICE STRAW
COfR PITH
O
u. 0,50
"~,
/
0.50
0.zs
0.25
___¢' ,"
0,00 150
3 SO
550 Temperature (C)
.~
>o
,
750~150
350 Temperature
SSO ( C )
0.00 SO
Figure 3 T.g.a. and d.t.g, curves of synthetic and untreated biomass
numerical values range from 28wt% for coir pith to almost 81 wt% for corn stalk. The above observations bring about the specificity of individual biomass types. This specificity is held to be attributable to the differences in composition of the biomass. To investigate whether the effects of the individual components of a biomass are simply additive, synthetic biomass prepared by mixing the relevant components was studied. Figure 3 presents the weight loss and derivative weight loss curves for synthetic biomass samples (the curve for fresh biomass is also shown for comparison). It is appropriate to mention here that the curves of the synthetic biomass matched well (with minor deviations explained later) the curves for fresh biomass only when wood cellulose was used for preparing the synthetic samples. The deviation of the fresh biomass curve was much greater for synthetic biomass containing Whatman cellulose. To explain this, X-ray diffraction (XRD) analyses were carried out to determine the crystallinity of celluloses isolated from various kinds of biomass. The crystallinity index (Cr-I) was calculated by the method suggested by Segal and co-workers 23'24. Table 4 gives the Cr-I values obtained for different types of cellulose. It may be noted that crystallinity of cellulose extracted from biomass lies between the Cr-I values of Whatman cellulose and chromatograph cellulose. It is known 25-28 that crystallinity influences the pyrolysis characteristics. Hence, to synthesize biomass it is necessary to use cellulose having a crystallinity similar to that in the natural biomass. Hence wood cellulose was used. The above results indicate the interesting fact that the components themselves play individually significant roles in determining the pyrolysis characteristics of biomass. Also, the basic structure or degree of polymerization of the biomass is less significant than its composition. In other words, the way in which components are bound (chemically) is not as important as the actual amounts of individual components present in a particular biomass. However, the chemical structure of the individual components, for example the type of
Table 4 Crystallinity index of cellulose Cellulose
Cr-I
Whatman Chromatograph Wood Coir Corn cob Rice husk Rice straw G r o u n d n u t shell
34.4 91.5 68,9 49,6 68.9 56.5 52.9 54.7
cellulose (defined in terms of Cr-I, which incidentally is of the same order of magnitude for all types of biomass), has an effect on the pyrolysis characteristics of biomass.
Mathematical correlation for pyrolysis in t.g.a. An attempt was made to represent mathematically the correlation between the pyrolysis characteristics of biomass and those of its components. This simple additive correlation assumes that the overall pyrolysis behaviour ofa biomass is the weighted sum of the partial contributions of its components, the relative proportions in the total composition defining the respective weight losses of the components: [AWbIT]C = AWc]Z ° Xc + AWhIT" Xh + AW]IT" Xl + AWeIT'Xe
(1)
where [AWbtr]c is the correlated weight loss for any given biomass and X¢, Xh, Xl, Xe are the initial fractions of cellulose, hemicellulose, lignin and extractives, respectively, present in the corresponding biomass. The values of local weight loss at a temperature T for biomass and for all components, designated as AWblT, Awcl T, Awh[T, AWl[T and AwelT respectively, are obtained from individual t.g.a, curves. The correlated data for biomass are obtained through the following steps based on Equation (1): (1) The t.g.a, weight loss curve for each component is accurately read in a fixed temperature range to define
Fuel 1996 Volume 75 Number 8
991
Pyrolysis characteristics of biomass and biomass components: K. Raveendran et al. 1O 0
~,.-.~
Im
,
"~N~. /t
075-.
~
, ' [,
w
.
--"
~
~\
$1r1= 0.179 SE2=0.0098
SE,--" 0.130 $£2= 0"0097
\',~
~ '
S[1 ': 0"11'! $E2=0"0086
I
1,00 075 "~"
•
. . . .
/., 000
I
_..
,
~
\,, o
: ~0.S0-
\
ooo.---
150
Figure 4
550 Temperature ( C )
-:.
fi ~
750~150
~-.
--.
,
350 SSO Temperature( C)
_ Experiment=t ---Correlatedl Eqns.(l) ~,. (2) . . . . Correlated, Eqns.(3)~,. (4) $E1- Standard error, Eqn$. (1) ~, (2) SE2-Stondarderror , Eq~ (3). (4)
"
750
T.g.a.and d.t.g, curvesof untreatedbiomass,and correlateddata
the weight loss at that temperature; these weight losses are used as weightings. (2) The weighted sum of the weight loss of biomass over the specified temperature range is computed by multiplying the weight losses of the individual components obtained in step (1) by the corresponding initial fractions of the components present in the biomass and summing. (3) The resulting weighted sum is taken as a measure of weight loss at the specified temperature for the total biomass, which is summed over the temperature range 423-1023 K to obtain the complete weight loss curve.
Figure 4 compares the calculated weight loss with typical t.g.a, curves for few kinds of biomass, as well as the derivative curves for these cases. The equation for the derivative of weight loss curve is as follows (as may be deduced from Equation 1): dw
dw
q- (~)1T'XI-}- (~tt)e T"Xe
Fuel 1996 Volume 75 Number 8
To calculate the weight loss data using Equations (1) and (2), wood cellulose t.g.a, data were used (Figure 2). It is important to mention here that the cellulose chosen for study should have similar crystallinity, as discussed above. When Whatman cellulose data were used in Equations (1) and (2) the calculated weight losses did not match the experimental data. Figure 5 shows the experimental and calculated (using Whatman cellulose) t.g.a, data. As can be seen, the deviation obtained is quite significant, which is mainly attributed to the difference in the cellulose crystallinity. To put the above inferences on a firmer basis, correlation studies were also conducted with synthetic biomass samples prepared with and without ash; the results are shown in Figures 6 and 7. It is interesting and encouraging to find that the correlated data match very well with the calculated data in the case of synthetic biomass without ash. The deviation seen in the case of synthetic biomass with ash is obviously due to the influence of ash. To account for the influence of ash, a multiplicative correction factor was introduced in the above correlations:
(2)
where [(dw/dt)b[7-]c is the correlated derivative weight loss for any given kind of biomass and the values of the derivative weight loss at temperature T for the biomass and for all its components at a temperature are designated a s (dw/dt)blT, (dw/dt)clT, (dw/dt)llT, (dw/dt)hlr, and (dw/dt)elr, and are directly obtained from the individual d.t.g, curves. As can be observed, the correlated data are comparable with the experimental data, with a standard error of 0.13-0.18. The deviation in the correlated data is mainly in the region where the devolatilization rate is high and where a major fraction of the devolatilization occurs. This deviation may be attributed to the influence of ash elements, as shown below, besides some unaccountable factors such as density and porosity.
992
_ 0.50 =~n "~
-
750"
r.,',,p,~a,,,r,(C)
3.75 ~ _
550
350
~
Rice straw
f~
,
ISO
/ /\[1
~
//i~: ~ 350
= - J
0 ,It49
st.; -- o.oo-9,..
:l ' t] ; i~!:!
.. :
SEI =
":~ !I
~, ~ / ~ ] [ ~ i
"
~:4"~x. ~'
s~2 = o.oo9,.
/il
~02s-
,o
$EI = 0.137
,....)
[eXWblr]E = [aWb[r]C x
A-
~
XdT: xexp
(3)
[(~-~)b TJE= [(-~)b TIC
X (A- [()(a~nlxnz] Slj (4) where [AWbIT]E and [(dw/dt)blr]E are the experimental weight loss data and derivative weight loss data for any given biomass, Xa, Xsi and X1 are the initial silica-free
Pyrolysis characteristics of biomass and biomass components." K. Raveendran et al. 1.00
1.00 i
!
$E -- 0 o2241
0,75
~\
,,
f I
i 0,50 -
^
coRN
cob
_\/
~l
5E=0,2253
N
l
-
-
/
~l
~
.ICE
0-75 u
_
"~I~,
Husk
\\'
I
SE = 0"2656
-
-
wooo
, ~
1~, l,
_
I ?,~1
0.S0
-~
O.2S
.~
~ 0.25
-
-
i I
'.
~ . . . _ . m _ = _ = ~ _ = ~__ ~
0,00
-
~'X
i ~
"
,
I
"-~'x.
/~
" ~
i
,so
~so
,, - 0-00
~ - ~ -
sso
7s
Temperature (C) c 0,75
0,75 u
o
Experimental) Eqns. (1) ~. (2) Correlated
"= o.s,
0.50 .~
SE-Standard
(Whatmo error
Cellulose)
II
.c
o.2s -~
T, 0.2
0,00
0,00
150
350
SSO
Temperature
750,150
(C)
350 SSO Temperature (C)
750
5 T.g.a. and d.t.g, curves of untreated biomass,and correlated data
Figure
1.0
g o.7s
u~ --
5E = 0,0683
5E = 0.0732 Synthetic rice husk without ash
Synthetic corn cob with out ash
~ ~ ~ "
....
.,
5E "~ 0.7 o u
I
1 !SE
0.00
" 0.75
: 0"0 473
°
: 0.0732
Synthetic coir pith
~so "-¢'~
Synthetic rice straw" 0.?S ;" I without ash
" '
I
without ash
==
- 050
O.S,
350Temperature(cSSO )
7SO
Experimental .... Correlated, E q n s . ( t l & ( 2 ) SE - Standard error
.=
E en 0.25
0.25 \
0.00F=iS0
0.25 o.oo
o
Figure
Synthtti¢ wood wdhout ash
I
t,oo
O.SO
0-25
~
5E=0.0142
-~
350 550 Temperature(C)
_
/• "/5OptSO
r~ 350 SSO Temperature ( C )
~'0.00 750
6 T.g.a. and d.t.g, curves of syntheticbiomasswithout ash, and correlated data
ash, silica and lignin fractional contents in dry biomass on a weight basis, and A = 0.5, n 1 = 8.5 and n2 = 7.0. Constants A, n 1 and n2 were obtained by non-linear optimization, using least-squares criteria for convergence, using the Box complex algorithm 29. As can be seen from Figures 4 and 7, the weight loss data predicted using Equations (3) and (4) match the experimental results well. Equations (3) and (4) show the combined effect of ash elements and biomass constituents, particularly the lignin component. This may be explained in a similar manner to that in the earlier paper on the influence of ash 22, that lignin present in biomass forms more char during pyrolysis and this char is gasified by catalytic
action of certain ash elements in the presence of pyrolysis products such as water and carbon dioxide. The inorganic ash elements are known to play a catalytic role in gasifying char. Silica in ash does not play any catalytic role; however, it alters the thermal properties and the pore structure of char and hence its reactivity.
Mathematical correlation for pyrolys& in PBP: prediction of product distribution To investigate further, studies were carried out in the packed-bed pyrolyser with five representative types of biomass and also the standard biomass components such as cellulose, xylan and lignin. This was intended to provide insight into the interaction, if any, among the
Fuel 1996 Volume 75 Number 8
993
Pyrolysis characteristics of biomass and biomass components: K. Raveendran et al. 1,00
-~
~
~'~x
n
,
-- "~"~
SEI =0.2Or.
;o,,-
,,,=o.oo,,
~tl.~, 1, Synthetic corn ¢ob k ~1 ~'l w,,h a.h
t
o.,o-
_
.'. i
o~s-
I
°°°L /
~"
/ 0150I-
~ ~
~ 0.25[ ;
~ / ~
Synth¢tic ricehusk j~ll w.h ash
t 0.00~ 150
//I
-
~
I
//i:
I~
- 0.50 "~ :
i
750,150
I
350
.... ~
550
Temperature(C)
~ ,
~
" 550
- J
z
~ I
' " ' ' ~
350
/-i
~
* " " °'"
--
Temperature(C)
o oo,,
'~ ! Synth¢tic e,\.~L .,,h a.h
Synthcticrice straw 3.75
~
/
..#..-5' k
1.00
_o,,7
woo¢l
--oso
o.,.
Synthetic coir pith
i
SEt = 0,981
ti\4 /..:,,
A~'~ w,l, a.~ t ~ _
r
"~.
o oo,.
,_- o.
.9
~
\.,
......
"-~
I
5E! = 0.932
"
"...
). i
I
'"%,
~ - ~
"~"~------ooo ,,o
o
,o
~
Expcrimentol E,.s..~= C2~ ...... Correlettetl, Eqll,, (3)& (/~-) SE1- StandarderrorpEqBs(1)&(?.) SE2-Standard error, EqBS (3'&'&) ---¢orre,a,ed.
._.
0.,, 13 13n 750
Figure 7 T.g.a. and d.t.g, curves of synthetic biomass with ash, and correlated data
Table 5 Pyrolysis product yields (wt% daf) in PBP experiments
Biomass Coir pith Corn cob Groundnut shell Rice husk Rice straw Subabul wood Components Cellulose Lignin Xylan
Volatiles
Char
Liquids
Gas
70.5 79.9 72.9 82.7 78.8 80.7
29.5 20.1 27.1 17.3 21.2 19.3
29.5 37.4 40.5 41.2 47.0 22.6
41.0 42.5 32.5 41.5 31.8 58.1
88.9 58.3 79.3
11.1 41.7 20.7
46.8 26.8 40.5
43.0 30.5 38.8
components during pyrolysis. The information on product distribution would be useful in design and development of any pyrolysis process.
Biomass and components pyrolysis product distribution. Table 5 presents the product distribution obtained for biomass on pyrolysis in the PBP. The char yields from biomass vary from ~17wt% for rice husk to 30wt% for coir pith. The liquid yields vary from ,,~22wt% for wood to almost double that figure for groundnut shell and rice husk. The gas yields vary from ~30wt% for groundnut shell to ~58wt% for wood. As seen earlier in the t.g.a, studies, here also biomass showed specificity in terms of product distribution. (Note: for the sixth set of data, for rice straw, see below.) To gain more information to explain the specificity of biomass or the variation in product distribution, experiments with biomass components were conducted with the individual isolated components. The product distributions obtained for the biomass components studied are also presented in Table 5. As can be seen, wood cellulose gives the lowest char yield and the highest liquid and gas yields, whereas lignin gives the highest char yield and the lowest liquid and gas yields. The product distribution from xylan lies between that of cellulose and lignin. It is also interesting to note that the
994
Fuel 1996 Volume 75 Number 8
product distributions obtained from most of the biomass types are intermediate between those obtained with cellulose and lignin. In other words, the product yields obtained with cellulose and with lignin hold as the upper and lower limits respectively for the product yields that can be obtained with most types of biomass, for a given operating condition. This is an important criterion to be considered for any biomass pyrolysis process design.
Additive correlation for product distributhgn. Similar to the correlation discussed above for t.g.a, data, an additive model is suggested for predicting the product distribution obtained in the PBP. On similar lines as explained earlier, this model also assumes that the additive overall product distribution of biomass is the weighted sum of the partial contributions of its components, with the relative proportions in the total composition defining the respective weights of the components:
(rb, i)c =rc,
i.xc+
rx, i . ( x h + x o ) +
rl,;.(xL)
(5)
where Y is the product yield and subscripts b, c, x, h, and 1 denote the yields from biomass, cellulose, xylan, hemicellulose and lignin as obtained by additive analysis, on the dry and ash-free basis; subscript i denotes char, liquid or gas accordingly. The extracts and hemicellulose are seen to behave similarly, especially in terms of char yield (cf. Figure 2). Hence extractives can be combined with hemicellulose (Xh + Xe) for calculating the product yield. Table 6 presents the theoretical product distribution calculated by Equation (5) and the deviation from the experimental values. As can be seen, the deviation is greatest for wood and least for groundnut shell. To explain the deviation between the experimental and calculated product yields, the possible influence of two factors was considered: (1) interaction among the components; and (2) influence of ash elements. Each factor was investigated separately.
(1) Interaction among components. Since cellulose and lignin together form ,,~70wt% of the biomass, these two
Pyrolysis characteristics of biomass and biomass components." K. Raveendran et al. Table 6 Calculated product yieldsa Char
Liquid
Gas
Calc.
Diff.
Dev.
Calc.
Diff.
Dev.
Calc.
Diff.
Dev.
Coir pith
24.4
-5.2
- 17.5
46.8
17.3
Corn cob
21.1
1.0
5.2
42.4
4.9
58.9
36.5
-4.5
- 10.9
13.1
40.7
- 1.8
Groundnut shell
24.0
-3.1
- 11.4
38.8
-4.3
- 1.7
-4.5
38.0
5.5
Rice husk
21.2
3.9
22.4
17.1
41.5
0.3
0.7
39.9
-1.7
-4.0
Rice straw
21.1
-0.1
Wood
21.9
2.6
-0.5
41.8
-5.2
-11.0
39.0
7.2
22.7
13.3
39.2
16.6
73.2
38.0
-20.0
-34.5
a Calc., calculated product yield (wt% daf) by Equation (5); Diff., difference between experimental and calculated values (wt% daf); Dev., relative deviation from experimental data (%)
Table 7 Pyrolysis product yields (wt%) of component mixtures Cellulose: lignin ratio
Volatiles
Char
Liquids
Gas
81.9 70.2 63.9
18.1 29.8 36.1
40.5 33.1 31.5
41.3 37.0 32.8
3:1 1: 1 1:3
60
I
I
I
I
I
t
(0:1)
(1:3)
(1:1)
(3;1)
(1:0)
',.9 LIJ
~= 4O Q ,..I
_w ~ 20 U
o t, - GAS Y I E L D
O.
~ ,,,!
....
- 20
0
i ....
i ....
20
40
CE L L U L O S E ( ' / , I
.
120
l
l
,
i
. . . .
100
I
80
)
i
,
i
I
60
i ....
I ....
60
. . . .
80
i ....
100
120
WEIGHT) I
40
,
i
i
,
I
20
. . . .
(Yb, i)E ----=(Yb,i)C x [A XI nl j[-~2 Xsn3]
(6)
where ( Yb,i)E is the product yield obtained experimentally and (Yb,i)C is defined by Equation (5). Figure 9 shows a plot of %,i vs. [A XInl Xan2 Xs~3], where 'ffb,i (Yb, i)E/(Yb, i)C and the values of nl, n2, n3 and A are presented in Table 8 along with standard error and R 2 values. These values of nl, n2, n3 and A are difterent for the char, liquid and gas yields. The above correlation was developed using data obtained for the five biomass materials used in the PBP experiments. To validate the correlation, experiments were conducted for one more kind of biomass, namely rice straw; yields of char, liquids and gases matched well with those obtained from the correlation. It can be inferred from the results of both t.g.a, and PBP experiments that the secondary reactions or any other effect due to the increased sample mass mainly influence the product distribution and that the interaction of the individual components remains unaffected. This is clear from the fact that the t.g.a, and PBP results can be used to develop correlations to obtain pyrolysis characteristics. The differences lie in the correlation coefficients and not in the observed trend. =
(C-" L ) - C E L L U L O S E °LIGNIN RATIO
..,/.
ash and silica contents of the biomass, is in accord with previous studies22. The following relation describes this function:
i
0
. . . .
i
-20
LIGNIN 1% WEIGHT)
Figure 8 Product yields of cellulose-lignin mixtures vs. cellulose and lignin contents
major components were mixed in different proportions. The experimental product distributions obtained are presented in Table 7. As seen for biomass product distribution above, here also the product yields obtained lie between the upper and lower limits obtained from cellulose and lignin. It can also be readily seen from Table 7 that as the cellulose content of the mixture decreases, the char yield increases and the liquid and gas yields decrease. Figure 8 presents the product distribution of mixtures as a function of their cellulose or lignin content. It is clear that the product yields are direct functions of the initial composition of the mixture. Hence it may be said that there is no detectable effect of interaction among the components on pyrolysis product distribution.
(2) Influence of ash. To explain the influence of ash, a similar approach to that used for the correlation in the t.g.a, study was used here. The influence of ash was found to be a function of the lignin, silica-free
Simplified overall correlation. To be able to predict the biomass pyrolysis product yields using Equation (6), it is necessary to know not only the chemical composition of the biomass but also the pyrolysis product yields from the individual components at a specified temperature. To simplify this further, another correlation was attempted, keeping in view the influences of ash and of cellulose and lignin, which contribute the major fraction by weight. This is expressed as:
where the constants A, nl, n2 and n 3 were obtained by regression analysis and are presented in Table 8 along with the standard error and R 2 values obtained. The correlations were developed with the experimental data obtained for five kinds of biomass and were verified with the data obtained for rice straw. Figure 10 shows the product yields as a function of composition, and as can be seen, the calculated yields are in good agreement with the experimental values.
Fuel 1996 Volume 75 Number 8
995
Pyrolysis characteristics of biomass and biomass components: K. Raveendran et al. 1'60
WD
1,4,0 ~ 1.20 Gs
GS ~ A ~ C P
~,~ 1.00
.
.
.
.
....
.
"80 ,.2 -°''u
l o.s
I 0"8
', I 1-o
L 1-2
I 1.4
0.4.
5S o R H
1o,
...............
.
.
.
cP
/2~- - -6 /
.
I i
I
I
0'6
0'8
.
.
.
.
.
.
.0"8 8
"CB
I
1"0
I
I
1"2
1,4-
1"6
.
oc/cop e
.g
~. 0,92 0'56
o
1'300
I WDI 0'6
aWD
CP I
i 0'8
!
If 1'0
I
I
I
I
0'6 0'6 0'8 RH-RICE HUSK RS-RICE STRAW
0"2
CP- COIR PITH CP CB-CORN COB -GSGROUNDNUT SHELL f o G S
I
1'0
1"2 CP o j
-wo -WOOD
1"4.
s
1.10 J¢ u
.
.
.
.
0.90 o I RH I 0"70 i ~ i I I I I I 0,75 -85 -95 1"05 n 1 n2 n3 X t .Xa "Xsi C
.
.
.
.
WD /
3; f-oG8
.
I -
R Ho,,,~/
Jl I
I
1.15
I
I
0"8
I
"9 -X•
I II
1"0
-Xsi
]
I
1-1
I
I
1'2
1.5
.C
Figure 9 Ratio of product yields (experimental and correlated) vs, cellulose, lignin, ash and silica contents
Table 8 Regression coefficients for Equations (6) and (7) Correlation coefficients A
n1
n2
n3
Standard error
R2
0.79
Additive correlation (Eqn 6) Char
2.584
0.378
0.343
-0.048
0.13
Liquid
1.888
0.694
0.273
0.030
0.23
0.79
Gas
0.349
-0.562
-0.363
0.001
0.15
0.86
0.99
Overall correlation (Eqn 7) Char
49.90
-0.0258
0.258
-0.044
0.03
102.86
0.7028
0.335
0.033
0.16
0.90
Gas
13.93
-0.4168
-0.337
-0.002
0.15
0.87
Char (t.g.a.)
56.43
-0.2664
0.079
0.048
0.12
0.91
Liquid
The above correlations give a good fit for char yields obtained in the t.g.a, studies also. The correlation coefficients obtained for t.g.a, data are also presented in Table 8. As can be seen from Figure 10, the t.g.a, data also follow a similar trend to that of the PBP data. CONCLUSIONS Thermogravimetric studies show that each kind of biomass has unique pyrolysis characteristics, by virtue of the specific proportions of the components present in it. A simple additive correlation to define this phenomenon has been developed, which suggests that there is no detectable interaction among the components during pyrolysis. These results are supported by
996
Fuel 1996 Volume 75 Number 8
synthetic biomass studies. The influence of ash is reflected in the correlation developed. It is also shown that the individual components taken for study should be similar to those present in biomass. In this respect the crystallinity index is found to be a useful measure for comparing various celluloses. The specificity of biomass is also seen in the product distribution obtained in the packed-bed pyrolyser, which is further explained by the product distribution from the individual components. In this case also, a simple additive correlation incorporating the influence of ash has been developed which describes the product distribution well. These results are supported by studies on mixtures of components. The results re-emphasize the fact that the feedstock composition and properties play a
Pyrolysis characteristics of biomass and biomass components. • K. Raveendran et al. 60 WD
S
~" so
RH
40 >.
30
20 !o SO
i
i
,'o
i
s'o '
GS
,o
!
30~
'
4 ~0
:
510
i
60
RHH~"
o.. 4(] i/I ' ~ 313 >2(
I
;tl5
I
3iS
!
415
1
20
I
30
30
40
5b
oGS
m 2S 13.
"E
~ 20 oRH'~ 15
45
,¢
15
I
* oWO
,119
I
I
I
23
217
CP- COIR PITH 40 • CB-CORN COB G5- GROUNDNUT SHELL
I
I
I
I
19
1S
I
I
23
I
RH-RICE HUSK RS-RICE STRAW WD-WOOD
CP o
I
27
30 CP o
p. io
~>- 30
20
20
J
2,5
i
3
~
i
i 35
x l. x 2. x , in3 . c
I
w 37.5 20
I
25
t
/
30
I
!
35
I
37.5
rxq "1 x2.2 x, 3. c L~?J "
Figure 10 Product yields as a function of cellulose, lignin, ash and silica contents
significant role in determining pyrolysis characteristics a n d p r o d u c t distribution, a n d these are generically specific for each k i n d o f b i o m a s s , u n d e r given o p e r a t i n g conditions. Results o f b o t h t.g.a, a n d P B P experiments show t h a t s e c o n d a r y reactions o r a n y o t h e r effect due to an increase in s a m p l e m a s s m a i n l y influence the p r o d u c t distribution; the i n t e r a c t i o n o f the i n d i v i d u a l c o m p o n e n t s r e m a i n s unaffected. ACKNOWLEDGEMENTS T h e a u t h o r s wish to t h a n k P r o f e s s o r W o l f g a n g Klose, F G T h e r m o d y n a m i k , K a s s e l University, G e r m a n y , for useful discussions a n d v a l u a b l e suggestions. T h e y also t h a n k D r K. I. G a s n a s e k a r a n a n d D r A. S u n d a r e s a n , R e s e a r c h Scholars in the C h e m i s t r y D e p a r t m e n t , I.I.T., B o m b a y , for their help in o b t a i n i n g the X R D patterns.
REFERENCES Antal, M. J., Jr. In 'Advances in Solar Energy', Vol. 2 (Eds K. W. Boer and J. A. Duffle), American Solar Energy Society, New York, 1983, pp. 175-239
2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
Akita, K. Rep. Fire Res. lnst. Japan 1956, 9, 10 Stamm, A. J. Ind. Eng. Chem. 1956, 48, 413 Arseneau, D. F. Can J. Chem. 1961, 39, 1915 Radislav, D. and Frantisek, R. Holz Roh- Werkstoff 1962, 20, 473 Tang, W. K. Forest Service Research Paper FPL 71, US Dept. of Agriculture, 1967 Shafizadeh, F. and DeGroot, W. F. In 'Fuels and Energy from Renewable Resources' (Eds F. Shafizadeh, K. Sarkanen and D. Tillman), Academic Press, New York, 1977 Shafizadeh, F. and McGinnis, G. D. Carbohyd. Res. 1971, 16, 273 Beall, F. C. Wood Sci. Technol. 1971, 5, 156 Roberts, A. F. In 'Thirteenth Symposium (International) on Combustion', The Combustion Institute, Pittsburgh, 1971, p. 893 Mellotte, H. and Richard, J. R. In 'Energy from Biomass', EUR 8245, Commission of the European Communities, 1983, pp. 523-529 Vovelle, E. and Mellotte, H. In 'Energy from Biomass', EUR 8245, Commission of the European Communities, 1983, pp. 925-930 Ward, S. M. and Braslaw, J. Combust. Flame 1985 21,261 Klose, W. and Stuke, V. Fuel Process. TechnoL 1993, 36, 283 Nikkhah, K., Bakhshi, N. N. and MacDonald, D• G. In 'Energy from Biomass and Wastes XVI' (Ed. D. L. Klass), Institute of Gas Technology, Chicago, 1993 McGee, B., Norton, F., Snape, C. E. and Hall, P. J. Fuel 1995, 74, 28
Fuel 1996 Volume 75 Number 8
997
Pyrolysis characteristics of biomass and biomass components: K. Raveendran et al. 17 18 19 20 21 22
998
Khan, M. R., Heshieh, F. Y. and Heaclky, L. Am. Chem. Soc. Div. Fuel Chem. Preprints 1989, 34, 1167 Saxby, J. D. and Sato, S. Fuel 1990, 69, 1109
Browning, B. L. 'Methods of Wood Chemistry', Vols I & II, Wiley, New York, 1970 Shapino, S. and Brannock, W. W. 'Rapid Analysis of Silicate, Carbonate and Phosphate Rocks', US Geological Survey Bulletin 1144-A, Washington, 1962 Boar, P. L. and Ingram, L. K. Analyst 1970, 95, 124 Raveendran, K., Ganesh, A., Khilar, K. C. Fuel 1995, 74, 1813
Fuel 1996 Volume 75 Number 8
23
Segal, L., Creely, J. J., Martin, A. E., Jr and Conrad, C. M. Tex. Res. J. 1959, 21, 786
24 25 26 27 28 29
Loeb, L. and Segal, L. Tex. Res. J. 1955, 21, 516 Chatterjee, P. K. and Conrad, C. M. Tex. Res. J. 1966, 36, 487 Kilzer, F. J. and Broido, A. Pyrodynamics 1965, 2, 151 Broido, A. and Weinstein, M. Combust. Sci. Technol. 1970, 1, 279 Weinstein, M. and Broido, A. Combust. Sci. Technol. 1970, 1, 287 Box, M. J. Computer J. 1965, 8, 42