This work focuses on thermodynamic analysis of the autothermal reforming of palm empty fruit bunch (PEFB) bio-oil for the production of hydrogen and syngas. PEFB bio-oil composition was simulated using bio-oil surrogates generated from a mixture of acetic acid, phenol, levoglucosan, palmitic acid and furfural. A sensitivity analysis revealed that the hydrogen and syngas yields were not sensitive to actual bio-oil composition, but were determined by a good match of molar elemental composition between real bio-oil and surrogate mixture. The maximum hydrogen yield obtained under constant reaction enthalpy and pressure was about 12 wt% at S/C = 1 and increased to about 18 wt% at S/C = 4; both yields occurring at equivalence ratio Φ of 0.31. The possibility of generating syngas with varying H2 and CO content using autothermal reforming was analysed and application of this process to fuel cells and Fischer-Tropsch synthesis is discussed. Using a novel simple modelling methodology, reaction mechanisms were proposed which were able to account for equilibrium product distribution. It was evident that different combinations of reactions could be used to obtain the same equilibrium product concentrations. One proposed reaction mechanism, referred to as the ‘partial oxidation based mechanism’ involved the partial oxidation reaction of the bio-oil to produce hydrogen, with the extent of steam reforming and water gas shift reactions varying depending on the amount of oxygen used. Another proposed mechanism, referred to as the ‘complete oxidation based mechanism’ was represented by thermal decomposition of about 30% of bio-oil and hydrogen production obtained by decomposition, steam reforming, water gas shift and carbon gasification reactions. The importance of these mechanisms in assisting in the eventual choice of catalyst to be used in a real ATR of PEFB bio-oil process was discussed.
Citation: Lifita N. Tande, Valerie Dupont. Autothermal reforming of palm empty fruit bunch bio-oil: thermodynamic modelling[J]. AIMS Energy, 2016, 4(1): 68-92. doi: 10.3934/energy.2016.1.68
Related Papers:
[1]
Hafizah Abdul Halim Yun, Valerie Dupont .
Thermodynamic analysis of methanation of palm empty fruit bunch (PEFB) pyrolysis oil with and without in situ CO2 sorption. AIMS Energy, 2015, 3(4): 774-797.
doi: 10.3934/energy.2015.4.774
[2]
Lihao Chen, Hu Wu, Kunio Yoshikawa .
Research on upgrading of pyrolysis oil from Japanese cedar by blending with biodiesel. AIMS Energy, 2015, 3(4): 869-883.
doi: 10.3934/energy.2015.4.869
[3]
Sunbong Lee, Shaku Tei, Kunio Yoshikawa .
Properties of chicken manure pyrolysis bio-oil blended with diesel and its combustion characteristics in RCEM, Rapid Compression and Expansion Machine. AIMS Energy, 2014, 2(3): 210-218.
doi: 10.3934/energy.2014.3.210
[4]
Hasanudin Hasanudin, Wan Ryan Asri, Utari Permatahati, Widia Purwaningrum, Fitri Hadiah, Roni Maryana, Muhammad Al Muttaqii, Muhammad Hendri .
Conversion of crude palm oil to biofuels via catalytic hydrocracking over NiN-supported natural bentonite. AIMS Energy, 2023, 11(2): 197-212.
doi: 10.3934/energy.2023011
[5]
Wasinton Simanjuntak, Kamisah Delilawati Pandiangan, Tika Dwi Febriyanti, Aryani Putri Islami, Sutopo Hadi, Ilim Ilim .
Catalytic upgrading of palm oil derived bio-crude oil for bio-hydrocarbon enrichment using protonated zeolite-Y as catalyst. AIMS Energy, 2024, 12(3): 600-616.
doi: 10.3934/energy.2024028
[6]
Xianhui Zhao, Lin Wei, James Julson .
First stage of bio-jet fuel production: non-food sunflower oil extraction using cold press method. AIMS Energy, 2014, 2(2): 193-209.
doi: 10.3934/energy.2014.2.193
[7]
Vincent E. Efeovbokhan, Augustine O. Ayeni, Osuvwe P. Eduvie, James A. Omoleye, Oladotun P. Bolade, Ajibola T. Ogunbiyi, Victoria N. Anyakora .
Classification and characterization of bio-oil obtained from catalytic and non-catalytic pyrolysis of desludging sewage sample. AIMS Energy, 2020, 8(6): 1088-1107.
doi: 10.3934/energy.2020.6.1088
[8]
Lihao Chen, Kunio Yoshikawa .
Bio-oil upgrading by cracking in two-stage heated reactors. AIMS Energy, 2018, 6(1): 203-215.
doi: 10.3934/energy.2018.1.203
[9]
Dheeraj Rathore, Anoop Singh, Divakar Dahiya, Poonam Singh Nigam .
Sustainability of biohydrogen as fuel: Present scenario and future perspective. AIMS Energy, 2019, 7(1): 1-19.
doi: 10.3934/energy.2019.1.1
[10]
Shouyun Cheng, Lin Wei, Xianhui Zhao, Yinbin Huang, Douglas Raynie, Changling Qiu, John Kiratu, Yong Yu .
Directly catalytic upgrading bio-oil vapor produced by prairie cordgrass pyrolysis over Ni/HZSM-5 using a two stage reactor. AIMS Energy, 2015, 3(2): 227-240.
doi: 10.3934/energy.2015.2.227
Abstract
This work focuses on thermodynamic analysis of the autothermal reforming of palm empty fruit bunch (PEFB) bio-oil for the production of hydrogen and syngas. PEFB bio-oil composition was simulated using bio-oil surrogates generated from a mixture of acetic acid, phenol, levoglucosan, palmitic acid and furfural. A sensitivity analysis revealed that the hydrogen and syngas yields were not sensitive to actual bio-oil composition, but were determined by a good match of molar elemental composition between real bio-oil and surrogate mixture. The maximum hydrogen yield obtained under constant reaction enthalpy and pressure was about 12 wt% at S/C = 1 and increased to about 18 wt% at S/C = 4; both yields occurring at equivalence ratio Φ of 0.31. The possibility of generating syngas with varying H2 and CO content using autothermal reforming was analysed and application of this process to fuel cells and Fischer-Tropsch synthesis is discussed. Using a novel simple modelling methodology, reaction mechanisms were proposed which were able to account for equilibrium product distribution. It was evident that different combinations of reactions could be used to obtain the same equilibrium product concentrations. One proposed reaction mechanism, referred to as the ‘partial oxidation based mechanism’ involved the partial oxidation reaction of the bio-oil to produce hydrogen, with the extent of steam reforming and water gas shift reactions varying depending on the amount of oxygen used. Another proposed mechanism, referred to as the ‘complete oxidation based mechanism’ was represented by thermal decomposition of about 30% of bio-oil and hydrogen production obtained by decomposition, steam reforming, water gas shift and carbon gasification reactions. The importance of these mechanisms in assisting in the eventual choice of catalyst to be used in a real ATR of PEFB bio-oil process was discussed.
1. Introduction
The need to find alternative sources of energy to substitute fossil fuel has gained a lot of attention with increased awareness in the role greenhouse gases play in global warming [1]. These greenhouse gases, together with other pollutant gases, are produced during the combustion of fossil fuel and include CO2, SO2 and NOx[2]. Biomass, defined as biological material obtained from a living or recently living organism, is being considered as a substitute to fossil fuel since it is abundant and can be readily accessed and processed sustainably. Though the combustion of biomass/biofuel generates CO2, it is considered a near-neutral carbon process since the CO2 released during combustion is in part the same CO2 absorbed by plants when synthesizing carbohydrates during photosynthesis [3]. However, fossil CO2 emissions may occur during soil conditioning by use of synthetic fertilizers, and also during transportation and storage when using conventional vehicles and equipment for biomass production and processing. One other major drawback with the use of fossil fuel is that it is a finite source of energy and existing world reserves are known to be depleting fast as a result of increase in world demand [2,3,4,5]. Due to its complex organic nature, biomass can serve as feedstock, either directly or indirectly after processing, for the production of multiple bio-based chemical products [6,7]. Biomass can be valorised by thermochemical conversion into bio-oil which is a dark brown, polar, high-density and viscous organic liquid containing a complex mixture of oxygenated compounds such as sugars, carboxylic acids, phenols, esters, ketones, aldehydes and benzenoids [8,9,10]. Palm empty fruit bunch (PEFB) is obtained after oil extraction at palm oil mills. It is estimated that for every tonne (t) of palm oil produced from a fresh fruit bunch, approximately 1 t of PEFB is produced [11].
Bio-oil can be produced by fast/flash pyrolysis, slow pyrolysis and solvolysis of solid biomass feedstock [12]. Pyrolysis is the thermal decomposition of biomass in the absence of oxygen (or in the presence of a limited amount of oxygen) at temperatures ranging from 400 ºC to 600 ºC [13]. Many studies have been carried out to understand the thermochemical conversion of biomass into bio-oil (also called pyrolysis oil) with authors focusing on process parameters such as pre-treatment of biomass, particle size of feedstock, reaction temperature, reactor type, choice of catalyst and reaction time. Depending on process conditions and the presence or absence of catalyst, bio-oil yield from PEFB can range from 30-70 wt% of the solid feedstock, with non-condensable gases and solid char formed as co-products [14,15,16,17]. In a study carried out by Mantilla et al. [9], PEFB pyrolysis experiments were conducted in a fixed-bed reactor at temperatures 460-600 ºC, gas residence time 16-80 s and particle size <0.5 mm, as well as 0.5-1.4 mm. The maximum yield of bio-oil of 48.4 wt% was obtained at a temperature 540 ºC, gas residence time of 31 s and particle size <0.5 mm. Their study concluded that temperature was the most significant parameter, of the three considered, in determining the bio-oil yield. Abdullah et al. [18], studied the fast pyrolysis of PEFB using a fluidised bed system for which the reactor temperature was varied from 400 to 600 ºC, the residence time between 0.79-1.32 s, and particle size diameters (with corresponding ash content) were <150μm (8.49%), 150-250 μm (7.46%), 250-300 μm (6.70%) and 355-500 μm (4.83%). They obtained optimun bio-oil yields at a pyrolysis temperature of 450 ºC, retention time of 1.02 s and for the particle size 355-500 μm with the lowest ash content. The somewhat high ash content of PEFB (rich in potassium) particularly reduces the yield in bio-oil resulting in the need for pre-treatment options before pyrolysis. Abdullah & Gerhauser [14] were able to demonstrate that by washing PEFB feedstock with distilled water before pyrolysis, the yield in bio-oil (organic phase) obtained rose to 61.3 wt% which was significantly higher than that obtained for the unwashed PEFB of 34.7wt%. Solvolysis, or hydrothermal liquefaction, is the use of highly pressurised solvent (pressure > critical pressure of solvent) such as water, methanol, ethanol or mixture of water and organic solvent to cause the decomposition of biomass to yield bio-oil at mild temperatures, usually < 400 ºC [11,19]. A few studies have been carried out to determine the influence of different solvent and catalyst on the production of bio-oil by solvolysis of PEFB. Akhtar et al. [20], investigated the effect of different catalysts such as NaOH, KOH and K2CO3 on PEFB solvolysis with water as solvent. The highest bio-oil yield of about 67 wt% was obtained with K2CO3 catalyst with a concentration of 1 M. Table 1 presents a summary of some selected PEFB bio-oil composition published by various authors.
Table 1. Summary of PEFB bio-oil composition obtained in literature.
The quality of bio-oil obtained depends on the type of biomass used and the severity of process parameters such as temperature, gas residence time and heating rate [24]. Bio-oil can be directly used as fuel in boilers or gas turbines or subsequently upgraded to produce automobile fuels and bulk chemicals using several methods such as zeolite catalytic cracking, hydrogenation, and aqueous phase processing [25].
Synthesis gas (syngas) can be produced from hydrocarbons by either reforming or partial oxidation. The most widely used process is the steam reforming of methane carried out in tubular reactors. Steam reforming (SR) is an endothermic process causing it to be very energy intensive. Reformers are designed to optimize heat exchange and recovery leading to huge capital investments[26,27,28]. The SR of methane for example, is carried out at pressures of 1.4-4.0 MPa and temperatures ranging between 750-900 ºC [29]. A number of thermodynamic and experimental studies have been performed on the SR of oxygenated organic compounds (oxygenates) in general and bio-oil in particular [5,30,31,32,33,34]. The general equation for the SR of oxygenates is given by reaction 1. This is an endothermic reaction and the amount of heat required depends on the value of the coefficients and . More hydrogen is produced subsequently via the water gas shift (WGS) reaction as CO reacts with water (reaction 2).
CnHmOk+(n−k)H2O→nCO+(2n+m−2k)/2H2
(1)
CO+H2OCO2+H2ΔHθ298=−41kJ/molofCO
(2)
It is noteworthy that, given its mild exothermicity, the WGS reaction is favoured at temperatures well below those used for SR and so, if a single reactor is to be used under isothermal conditions, an optimum temperature has to be selected for which both reactions combine to give maximum hydrogen yield. Otherwise, maximum yields are obtained using separate reactors operated at high and low temperatures respectively [35,36]. Partial oxidation (POX) involves the use of sub-stoichiometric combustion oxygen to produced synthesis gas from hydrocarbons. It can be operated either thermally or catalytically but the overall hydrogen yield is lower than that obtained for SR due to the absence of steam as an extra source of hydrogen and because some hydrogen oxidizes to form water [37]. POX is exothermic and can appear to produce syngas at a lower cost since SR requires energy to sustain the endothermic SR reaction and generate steam [38]. The POX of oxygenates is given in reaction 3. The amount of energy released by this reaction depends on the particular oxygenated feedstock. The amount of oxygen used is very critical during POX as the use of stoichiometric oxygen will lead to complete oxidation (COX-reaction 4).
CnHmOk+(n/2−k/2)O2→(m/2)H2+nCO
(3)
CnHmOk+(n+m/4−k/2)O2→(m/2)H2O+nCO2
(4)
Autothermal reforming (ATR) is a process that uses steam and oxygen to produce synthesis gas. ATR reactors are designed to couple endothermic SR with exothermic POX so as to obtain a thermoneutral (adiabatic) or slightly exothermic process [39,40,41]. This coupling helps to reduce the overall cost of the synthesis gas produced compared to SR and also achieves a higher hydrogen yield compared to pure POX process [42,43]. ATR systems in which both partial oxidation and steam reforming reactions occur in a single catalytic bed are mostly useful for fuel cell applications. Those which have a separate burner for complete oxidation (combustion zone) followed by a catalytic bed for steam reforming are ideal for gas to liquid applications [44]. The overall reaction occurring during the ATR of bio-oil (or any oxygenated fuel) can be written as shown in reaction 5 [45].
CnHmOk+aO2+bH2O→cH2+dCO+eCO2+fCH4
(5)
The stoichiometric coefficients c, d, e, and f depend on the amount of oxygen and steam (a and b) and also on temperature and extent of side reactions [31]. Some important side reactions which may occur together with SR, POX, WGS and COX include:
Decomposition:
CnHmOk→kCO+m/2H2+(n−k)C
Methanation of carbon:
C+2H2CH4ΔHθ298=−75kJ/molofC
Boudouard reaction:
2COCO2+CΔHθ298=−86kJ/molofCO
Gasification:
C+H2OCO+H2ΔHθ298=−131kJ/molofC
C+1/2O2COΔHθ298=−110kJ/molofC
Methane dry reforming:
CH4+CO22CO+2H2ΔHθ298=247kJ/molofCH4
Unlike methane and other fossil fuels, the use of bio-oils as feedstock to generate hydrogen (or synthesis gas) presents some tough challenges because of their very heterogeneous composition and thermal instability [10,31,46]. The main challenge is the formation of coke due to dehydration and polymerisation reactions which can be mitigated by using excess steam, bio-oil blending and appropriate catalyst choice [38,47]. Czernik and French [37] performed ATR of oak, poplar and pine bio-oils using 0.5% Pt/Al2O3 catalyst. They were able to produce 9-11 g of H2 per 100 g of bio-oil (corresponding to 70 to 83% of stoichiometric potential) while varying temperature 800-850 °C, steam to carbon ratio 2.8-4.0 and oxygen to carbon ratio 0.9-1.1. Problems with volatility lead to the 11%-30% of bio-oil carbon forming deposits in the evaporator.
In order to optimise the ATR of bio-oil, it is important that the contributions of participating reactions in determining equilibrium product concentrations are known, bearing in mind the additional complexity stemming from using bio-oil as a feedstock. To the authors’ knowledge, no work has been published on this aspect of bio-oil ATR reforming and this work is an attempt to contribute to this knowledge gap.
2. Materials and Method
2.1. Bio-oil composition generation
Due to the enormous variability in chemical composition that exists among bio-oils produced from different biomass sources [48], most thermodynamic equilibrium simulations have been carried out with the use of model compounds to simulate bio-oil feedstock in SR, POX, and ATR studies. Among these model compounds, acetic acid has received the most attention [49,50,51]. Other compounds such as cresol, acetone and ethylene glycol have also been used [31,47,52]. Even when the same biomass feedstock is used, variations in process parameters lead to different bio-oil compositions. Sukiran et al. [21] studied the effects of pyrolysis temperatures, particle sizes and heating rates on the yield of bio-oil from PEFB. They found significant differences in the bio-oil obtained when these process parameters were varied in the following ranges: temperature 300-700 °C, heating rate 10-100 °C min−1, particle size <90, 91-106, 107-125 and 126-250 μm.
Using information from existing literature, the moisture free (mf) elemental compositions of some PEFB bio-oils were determined. The results are summarized in Table 2. Mean elemental compositions were determined excluding the results from [22] due to the low H/C ratio and high O/C ratio reported.
The mean over five mf elemental compositions obtained for PEFB bio-oil was C0.3724H0.5304O0.0886 for which the nitrogen content was neglected. Acetic acid, phenol, levoglucosan, palmitic acid and furfural were selected as representative compounds found in PEFB bio-oil, since their presence in PEFB bio-oil has been repeatedly detected in significant amounts via GC-MS semi-quantitative analyses [9,53,54]. Other authors have performed thermodynamic analysis of complex mixtures by using mixtures of simpler compounds with the same molar elemental compositions. Zin et al. [55] used mixtures of acetic acid, levoglucosan, vanillin and furanone to perform the thermodynamic analysis of pine bio-oil aqueous fraction steam reforming. Hanika et al.[56] used a mixture of glucose, vanillin, n-butyl-stearate, methionine and tri-ethyl-phosphate as representative compounds to simulate the partial oxidation of rape meal. Table 3 gives the physical properties of acetic acid, phenol, levoglucosan, palmitic acid and furfural.
Table 3. Physical properties of model compounds found in PEFB bio-oil.
Table 4. Model PEFB bio-oils obtained using different combinations of acetic acid, phenol, levoglucosan, palmitic acid and furfural. Target composition: C0.3724H0.5304O0.0886.
BOS1
BOS2
BOS3
BOS4
BOS5
Name
mol (%)
mol (%)
mol (%)
mol (%)
mol (%)
Acetic acid
12.2
23.2
0.0
0.3
0.0
Phenol
35.1
45.2
36.0
65.3
12.2
Levoglucosan
2.2
0.3
8.6
16.4
0.0
Palmitic acid
20.0
16.3
21.6
17.5
25.5
Furfural
30.5
15.0
33.8
0.5
62.2
TOTAL
100.0
100.0
100.0
100.0
100.0
C
0.3724
0.3724
0.3724
0.3724
0.3740
H
0.5390*
0.5390*
0.5372
0.5390*
0.5374*
O
0.0886
0.0886
0.0904*
0.0886
0.0886
Maximum relative error* (%)
1.6
1.6
2
1.3
1.3
* For a given surrogate, only the maximum elemental relative error, when compared to the average value, is shown. For BOS1, 2, 4 and 5, error shown for H, for BOS3, error shown for O.
Chemical Equilibrium with Applications (CEA) software developed by National Aeronautics and Space Administration (NASA) was used to perform thermodynamic equilibrium simulations. The software determines the equilibrium properties of a reaction mixture by using the Gibbs free-energy-minimisation method based on a known pool of reactant and species, and user defined initial composition, temperature and pressure [59]. This method takes into consideration the fact that the total Gibbs free energy of a reacting system reaches a minimum at equilibrium when varying the mixture composition at constant pressure and temperature.
For a given mixture with a number K of species, the Gibbs free energy can be written as:
G=K∑i=1μ ini
(6)
where is the Gibbs free energy, is the chemical potential of species and the number of moles of species . The condition of equilibrium is the minimisation of . In order to find the that minimize the value of , it is necessary that the values of satisfy certain constraints, one of which is the elemental mass balance given by
K∑i=1ajini−boj=0
(7)
or
bj−boj=0(j=1,…,M)
(8)
where aji are the number of gram atoms of element j per gram mole of species i and boj is the number of gram atoms of element j in the reaction mixture.
Using Lagrangian multipliers, G can be written as
G=g+M∑j=1λj(bj−boj)
(9)
where λi are Lagrangian multipliers and g the Gibbs free energy per gram of reaction mixture. Based on these equations, the condition for equilibrium can be expressed as:
G=K∑i=1niΔ G0i+RTK∑i=1nilnyi+RTK∑i=1nilnP
(10)
where ΔG0i is the standard Gibbs free energy of formation of species i, R is the universal gas constant, T is the tmperature in Kelvin, yi is the mole fraction of species i, and P is the total pressure.
The thermodynamic state for which the equilibrium composition is determined has to be specified by two intensive properties which, in principle, can be any combination of: temperature (T), pressure (P), specific enthalpy (h), specific entropy (s) and specific volume (v). In CEA, the “tp” setting is used for constant temperature and pressure processes, and the “hp” setting is used for constant pressure and enthalpy (adiabatic) processes. To solve equation 10, an iteration procedure is used with the Newton-Raphson method applied to solve for corrections to the initial estimates for composition, ni, Lagrangian multipliers, moles of gaseous species and (when required) temperature, T [59].
All inputs into the CEA software were based on a bio-oil feedstock carbon number of 1500 (user-chosen arbitrarily). This carbon number was also used to calculate the amount of water and oxygen needed for ATR based on the desired steam to carbon ratio (S/C) and amount of oxygen expressed as the equivalence ratio, Φ, as used to describe oxy-combustion processes (actual O2 to fuel molar ratio divided by stoichiometric combustion O2 to fuel molar ratio). All temperatures were entered in Kelvin (K) and pressures in atmosphere (atm).
Thermodynamic simulation executed in CEA generates an output file containing all relevant thermodynamic properties and an equilibrium composition in mole fractions. To obtain the molar yields of equilibrium products, the mole fractions were converted to moles using a carbon balance. The total number of moles at equilibrium was determined using equation 11.
total moles=nc∑ji=1cixi
(11)
where nc is the initial (input) moles of carbon in the feed which in this case is equal to the chosen carbon number of 1500, j is the number of carbon containing species at equilibrium, and ci and xi are the carbon number and mole fraction respectively of the species considered. Once the total moles of equilibrium species was determined, the equilibrium yields ni of each species present were calculated using equation 12.
ni=xi×total moles
(12)
For each set of process condition considered, the overall performance was evaluated using critical factors such as hydrogen yield (YH2, in wt% of the mf feed), and the percentage selectivity (Xi) to CH4, CO and CO2 (see equations 13-16).
YH2=mass of H2 at equilibriuminitial mass of bio−oil×100
(13)
XCH4=moles of CH4total moles of all carbon products×100
(14)
XCO=moles of COtotal moles of all carbon products×100
(15)
XCO2=molesofCO2totalmolesofallcarbonproducts×100
(16)
2.3. Modelling the global reactions of bio-oil ATR
It has been established in existing literature that several reactions are involved when biofuels (or organic fuels in general) undergo reforming to produce H2 or syngas [10,60,61,62]. In order to eventually optimise the ATR of bio-oil, it is helpful to devise a tool that can be used to predict the contribution of the different participating reactions. The CEA software determines equilibrium composition by applying numerical techniques which are independent of the actual reaction mechanisms taking place. To determine possible reaction mechanisms, different sets of reactions where proposed and then tested to see how well they could fit the equilibrium yields obtained using CEA. Mechanism testing was performed algebraically using the Solver function in MS Excel. A mechanism was accepted as correct if predicted equilibrium concentrations with the proposed mechanism were close to actual equilibrium concentrations for all species with a relative percent error less than or equal to 1%.
As an example, suppose we propose a mechanism for ATR of bio-oil for which only three reactions, POX, SR, and WGS are assumed to occur. The following equations can be written for BOS1:
C0.3724H0.5390O0.0886+0.1419O2→0.3724CO+0.2695H2
C0.3724H0.5390O0.0886+0.2838H2O→0.3724CO+0.5533H2
CO+H2→CO2+H2
Let the moles of BOS1 (nBOS1) consumed by POX and SR reactions be nPOX and nSR respectively; and nWGS the moles of carbon monoxide (CO) consumed by the WGS reaction. The general material balance equation for a particular species can be written as:
Input+Generation−Consumption=Output+Accumulation
(17)
Given that the calculations are performed for a system at equilibrium there is no accumulation term and the other terms will depend on the particular chemical species considered. Applying equation 17 to our example, we can write the following equations to predict equilibrium concentration (in moles) of the chemical compounds involved:
nCO=0.3724nPOX+0.3724nSR−nWGS
(18)
nCO2=nWGS
(19)
nH2=0.2695nPOX+0.5533nSR+nWGS
(20)
nBOS1=nPOX+nSR
(21)
The number of predicted chemical species has to be equal to the number of equations in order for system to produce a unique solution. This system consisting of 4 equations and 4 unknowns is entered in Excel such that
By substituting random values for nPOX, nSR and nWGS it is possible to generate estimated equilibrium yields of CO, CO2 and H2 corresponding to an input moles of bio-oil (i.e. nCO,nCO2,nH2 and nBOS1). These calculated values are then compared to the actual (desired) equilibrium concentration values and an error between both sets of values calculated. A solution is accepted if the errors for each species are below 1%, and the combined sum of errors is also below 5%.
Using the ATR equilibrium data generated by the CEA software, different combinations of reactions shown in Table 5 where tested using the methodology described above. Any number of the reactions given in Table 5 can occur during ATR of bio-oil. An acceptable mechanism should be able to account for all chemical species present at equilibrium and must contain equations which account for the following processes: bio-oil degradation, oxygen consumption, steam consumption, carbon formation, methane formation, carbon removal and methane removal.
Table 5. List of all reactions considered during bio-oil ATR mechanism modelling.
Name
Abbreviation
Reaction
1
Partial oxidation
POX
CnHmOk+(n-k)/2O2→nCO+m/2H2
2
Complete oxidation
COX
CnHmOk+(n+m/4–k/2)O2→nCO2+m/2H2O
3
Steam reforming
SR
CnHmOk+(n-k)H2O→nCO+(2n+m-2k)/2H2
4
Decomposition
DEC
CnHmOk→kCO+m/2H2 +(n-k)C
5
Water gas shift*
WGS
CO+H2O↔CO2+H2
6
Boudouard reaction
BO-RX
2CO→CO2+C
7
Methanation of C(s)
MEN
C(s)+2H2→CH4
8
Carbon gasification 1
C-GS1
C(s)+H2O→CO+H2
9
Carbon gasification 2
C-GS2
C(s) +0.5O2→CO
10
Methane steam reforming
ME-SR
CH4+H2O→CO+3H2
11
Carbon monoxide oxidation
CO-OX
CO+0.5O2→CO2
12
Hydrogen oxidation
H-OX
H2+0.5O2→H2O
* The reverse of the water gas shift reaction (R-WGS) was used in some cases
3.1. Sensitivity analysis and product distribution
The influence of bio-oil chemical composition on hydrogen yield was investigated. ATR was performed on the five bio-oil surrogates considered in this study by varying S/C ratio and amount of oxygen. The amount of oxygen used during ATR was expressed in terms of the equivalence ratio, Φ, which in this case was defined as the ratio of actual moles of oxygen/moles of carbon present in the reaction mixture to the stoichiometric moles of oxygen/moles of carbon needed for complete oxidation (COX) of the bio-oil feedstock. This is written mathematically as
The equivalence ratio was preferred over the more traditional O2/C ratio because it highlights the relative amount of oxygen in the system and indicates how far off the system is from complete oxidation (combustion). By dividing the moles of oxygen needed for the stoichiometric partial oxidation of the bio-oil with that needed for its stoichiometric complete oxidation (see reaction 3 and4) we can define the special value of equivalence ratio for stoichiometric partial oxidation, (equation 23).
ΦPOX was equal to 0.31 for all five bio-oil surrogates. Other equivalence ratios used in this study were obtained by considering 50%, 150%, and 200% of this value. Therefore for a given S/C ratio, the equivalence ratio considered were Φ = 0.15, 0.31 (ΦPOX), 46, and 0.61 (corresponding to molar O2/C of 0.19, 0.38, 0.57, and 0.76).
For all five bio-oil surrogates, equilibrium hydrogen yield and product concentrations were similar for all ATR equilibrium conditions examined. The maximum standard error obtained when comparing mean hydrogen yields from all five bio-oil surrogates was 0.324, corresponding to a percent error of 3.1%. This was obtained for S/C ratio = 4 and Φ=0.15 (Figure 1). This implies that equilibrium product distribution is insensitive to exact bio-oil composition. The minimal variations in mean hydrogen yield observed were due to the slight difference in elemental composition among the bio-oil surrogate mixtures.
Figure 1. Mean hydrogen yield (BOS1-5) and standard error for Φ (0.15−0.61) and S/C (1−4).
The bio-oil surrogate mixtures all undergo similar reactions when subjected to the same oxidizing conditions. Zin et al. [55] also found that chemical equilibrium products from the SR of different mixtures of simulated aqueous fraction of pine bio-oil were the same. Table 6 shows mean equilibrium temperatures and their standard deviations obtained during the ATR of the five bio-oil surrogate mixtures. The equilibrium temperatures are almost equal for similar conditions of steam and oxygen with the maximum percent error of 1.6% obtained for S/C=1 and Φ=0.30. This provided further evidence that ATR proceeds with a similar mechanism for all bio-oil surrogates and the equilibrium product composition depends on the final equilibrium (exit) temperature.
Table 6. Mean temperatures (K) and standard deviations obtained during ATR of the five bio-oil surrogates (BOS1-5) considered in this study.
At constant S/C ratio, increasing the amount of oxygen (equivalence ratio) causes the exothermic oxidation reaction to become more favourable leading to an overall increase in the temperature of the system. On the other hand, increasing S/C ratio at constant Φ reduces the equilibrium temperature due to the high heat capacity of water which absorbs some of the surrounding heat without causing a temperature increase.
Figure 2a shows the influence of the S/C ratio on the equilibrium hydrogen yield for the ATR of BOS2. BOS2 was used to discuss all remaining results because its composition is more realistic when compared to that of bio-oils found in published literature. As expected, the amount of equilibrium hydrogen increased with increase in S/C ratio. The maximum hydrogen yield obtained was about 12 wt% at S/C = 1 and increased to about 18wt% at S/C = 4. This was due primarily to increase in water gas shift reaction which shifts to the right (towards forming more products) as more steam is introduced in the system.
Figure 2. Influence of S/C ratio and the equivalence ratio on the amount of hydrogen produced during the ATR of BOS2. (a) Φ range 0.15-0.61, (b) refined Φ scale (0.21-0.39).
For all S/C ratios, the maximum hydrogen yield was obtained at values of Φ close to ΦPOX (0.31) that is, when the amount of oxygen in the system was close to that needed for stoichiometric partial oxidation. A closer look (Figure 2b) reveals that at low S/C ratios (1 and 2), the maximum hydrogen yields occurred at equivalence ratio lower than ΦPOX (Φ≈ 0.27) but attained this value at S/C = 4. ΦPOX is therefore an important parameter that can be used to determine the amount of oxygen to use during ATR in order to achieve maximum hydrogen yield.
3.2. Selectivity to carbon containing products
The influence of S/C ratio and amount of oxygen on the selectivity to carbon containing products is shown in Figure 3. Overall, as more steam is added, the product gas becomes increasingly rich in CO2. In Figure 3a, solid carbon (in the form of graphite), CH4, CO2 and CO are all present at the lowest equivalence ratio considered in this study. The presence of carbon and methane is an indication of possible bio-oil thermal decomposition and Boudouard reaction. Both reactions are known to occur under oxidant deficient conditions [47]. As more oxygen is added, carbon and CH4 essentially become negligible and CO remains as the major product.
Figure 3. Influence of S/C ratio and O2/C ratio on the selectivity of carbon and carbon containing products during the ATR of BOS2 at 1 atm. a) S/C = 1 b) S/C = 2 c)S/C = 3 d) S/C = 4.
The decrease in carbon and CH4 concentrations observed between Φ=0.15 and 0.30 (Figure 3a) is due to carbon gasification and methane steam reforming respectively (reaction 8 and 10 on Table 5). Both reactions directly contribute in increasing the H2 concentration and CO selectivity. The high CO content indicates that virtually no water gas shift reaction takes place under this process condition. The decrease in CO2 concentration that occurs between Φ=0.3 and 0.45 is due to the reverse water gas shift reaction which become favorable at high temperatures (1204-1963 K). The slight increase in CO2 observed Φ>0.45 is due to the bio-oil undergoing combustion (complete oxidation) producing CO2 and H2O. The same trends as just explained hold for Figure 3b-d. The main difference being that as the S/C ratio is increased to 2, 3, and 4 respectively, the water gas shift reaction becomes increasingly prominent, converting most of the CO in the system to CO2.
3.3. Synthesis gas composition
Hydrogen or synthesis gas can be used as a primary feedstock in fuel cells or as feed for downstream chemical synthesis. Depending on the end use of the synthesis gas produced, the ATR process can be operated by choosing appropriate values for S/C ratio and amount of oxygen (Φ) to give a desired synthesis gas composition [63].
3.3.1. Fuel cell feed
Fuel cells are electrochemical devices which convert the chemical energy of a chemical reaction directly into electrical and thermal energy [64]. Generally speaking, there are five types of hydrogen fuel cells divided into two main categories: low-temperature fuel cells and high-temperature fuel cells. The low-temperature fuel cells such as polymer electrolyte membrane fuel cell (PEMFC), alkaline fuel cell (AFC), and phosphoric acid fuel cell (PAFC) operate at temperatures ranging between 370-473 K and use hydrogen as their only fuel source with very little tolerance for CO (<20ppm for PEMFC) [43,65]. For such fuel cells, the output gas from an autothermal reformer will have to be purified and all CO2, CO and unreacted feedstock removed to give an essentially pure hydrogen stream. To reduce the cost of the downstream purification, the ATR process will have to be operated under conditions of maximum hydrogen yield for a given S/C ratio and amount of oxygen (see Figure 2). Purification can then be achieved by using CO2 absorbent, pressure swing adsorption systems and catalytic preferential oxidation [43,66].
The high-temperature fuel cells include molten carbonate fuel cell (MCFC) and sulphur oxide fuel cell (SOFC). These fuel cells operate at much higher temperatures ranging between 923-1273 K and show more flexibility in feedstock and catalyst requirements [64]. For these fuel cells, hydrogen competes with CO and even CH4 as fuel source making the combined H2+CO from ATR an important parameter. Figure 4 shows the influence of Φ and S/C ratio on total H2+CO yield.
Figure 4. Influence of S/C ratio and Φ (amount of oxygen) on the total H2 + CO obtained during ATR of BOS2 at 1 atm.
The maximum total H2+CO reduces as the S/C ratio is increased from 1 to 4. This happens because the formation of CH4 becomes favorable at high S/C ratios and Φ<ΦPOX (Figure 3). In the absence of methanation, the total H2+CO remains the same (Φ>0.33) due the equal mole to mole ratio between H2 and CO as one mole of CO converted to CO2 via the water gas shift reaction gives one mole of H2. Irrespective of the S/C ratio chosen, optimal yield for H2+CO is obtained at equivalence ratio approximately equal to ΦPOX.
3.3.2. Chemical synthesis feed
Synthesis gas is an important intermediate in the production of several important chemicals such as methanol, dimethyl ether (DME), ammonia and liquid fuels (via gas-to-liquid ‘GTL’ processes). These processes rely either on direct combination of reactants or Fischer-Tropsch (FT) chemistry and have different requirements in the amounts of H2, CO and CO2 in synthesis gas. A key parameter for such processes is the H2/CO ratio whose variation for the ATR of PEFB is shown in Figure 5.
Figure 5. Plot of S/C ratio versus H2/CO ratio at different values of Φ during ATR of BOS2. Equilibrium simulation carried out at 1 atm. The maximum H2/CO ratio shown is 10.
Gas-to-liquid (GTL) Fischer-Tropsch processes for which only H2 and CO are reactants require a H2/CO ≈ 2 [36,67]. The synthesis gas in this case can be produced from an autothermal reformer operating with a S/C ratio between 1 and 2 and Φ<ΦPOX (Figure 5). Synthesis of higher alcohols require H2/CO = 1 [35]. In this case ATR can be performed at a low S/C ratio (S/C ≤ 1) and equivalence ratio slightly larger than the value required for stoichiometric partial oxidation, that is, Φ>ΦPOX so as to avoid any carbon formation (Figure 3a). H2, CO and CO2 are all reactants in methanol, dimethyl ether and high temperature Fischer-Tropsch synthesis [68]. For such processes, the synthesis gas is made to have the same stoichiometry as the final product with its composition expressed as shown in equation 24.
M=H2−CO2CO+CO2
(24)
M is called the module and is equal to 2 for methanol and dimethyl ether synthesis [35,68]. For ATR, this value of M can only be obtained either by addition of H2 or removal of CO2. For BOS2 used in this study, the maximum value of M was 0.7 obtained at Φ=ΦPOX for all S/C ratios examined. The explanations given in this section are simplified and meant to serve as a guide only. The eventual choice in process parameters will depend on other important factors like process scale and amount of product recycle [35,69].
3.4. Reaction mechanism
The advantage of using the Gibbs minimization energy is that a very large pool of chemical species is used to determine eventual equilibrium composition. The alternative will be to assume the prevailing reactions under the given process conditions and then use their equilibrium constants to determine equilibrium concentrations. The limitation of the latter is that fewer numbers of reactions and potential products are considered, compared to the number of species used by the Gibbs minimization method. Two main types of mechanisms were successful in accounting for equilibrium species obtained by the CEA software in the ATR of PEFB bio-oil. The results presented are those obtained from using BOS2 as feedstock. Similar results were obtained for all five bio-oil surrogates considered in this study. All reactions in this section are identified using the reaction nomenclature introduced in Table 5.
3.4.1. Partial oxidation (POX) based mechanism
The reactions considered for this mechanism were: bio-oil partial oxidation (POX), bio-oil steam reforming (SR), water gas shift (WGS), Boudouard reaction (BO-RX), methanation of carbon (MEN), carbon gasification (C-GS1 and C-GS2), methane steam reforming (ME-SR), hydrogen oxidation (H-OX) and carbon monoxide oxidation (CO-OX). These reactions were used to fit equilibrium results obtained at S/C = 1-4 and Φ=0.31 and 0.46. Not all reactions were involved at the same time for a given process condition. With a relative error less than 0.1% on individual molar production rates, it was impossible to distinguish between the mechanism-predicted and the equilibrium (‘actual’) values of H2 yield as shown in Figure 6.
Figure 6. Comparing actual equilibrium hydrogen yield with predicted yield using POX based mechanism.
The contribution to hydrogen production by participating reactions for the POX based mechanism is shown in Figure 7. By way of this mechanism, hydrogen production was primarily from POX, WGS and carbon gasification (C-GS1) at Φ=0.31 and by POX and SR at Φ=0.46. The implication of this mechanism is that under low oxygen content (Φ≤ΦPOX), the bio-oil completely undergoes POX and that more hydrogen is produced by the water gas shift reaction with no significant contribution from SR. Further production of hydrogen is achieved by the gasification of all solid carbon formed. Meanwhile at the intermediate Φ=0.46 (>ΦPOX), POX and SR were the only hydrogen producing reactions.
Figure 7. Percentage contribution to hydrogen production by participating reactions for the POX mechanism at (autothermal temperatures given in the figure for each S/C).
For this mechanism, no significant hydrogen consumption occurred as methanation of carbon (MEN) was virtually zero. In the case of S/C = 1 and Φ = 0.46, H-OX and reverse WGS (R-WGS) where responsible for H2 consumption. Figure 8shows how the amount of steam and oxygen influenced the bio-oil consuming reactions.
Figure 8. Influence of S/C ratio and oxygen on Bio-oil consuming reactions.
Table 7 gives a summary of the different reactions that dominate in the POX based mechanism at various temperature ranges. This table provides valuable information for the eventual choice of catalyst. For ATR carried out at temperature below 1100 K, the catalyst should be very selective to POX and WGS. High temperatures typical of low S/C ratios and high oxygen content (Φ≫ΦPOX) should be avoided since they can lead to catalyst degradation.
Table 7. Summary of POX based mechanism based on temperature range. Only reactions which contribute to equilibrium products are included.
Based on this mechanism, it can be said that at Φ≤ΦPOX, SR reactions are minimal and the choice of catalyst should be based on the prevailing POX, WGS and C-GS1 reactions.
3.4.2. Complete oxidation (COX) based mechanism
Another mechanism was validated for which complete oxidation (COX) was the dominant oxygen-consuming reaction. There was a near perfect agreement between the predicted hydrogen concentration and actual equilibrium hydrogen yield (Figure 9). The reactions considered for this mechanism were bio-oil thermal decomposition (DEC), bio-oil complete oxidation (COX), bio-oil steam reforming (SR), WGS, reverse water gas shift (R-WGS), methane steam reforming (ME-SR), methanation of carbon (MEN), carbon gasification (C-GS1 and C-GS2), and hydrogen oxidation (H-OX), as shown in Figure 10.
Figure 9. Plots of predicted hydrogen and actual equilibrium hydrogen showing near match with maximum relative error of 0.1% for the COX based mechanism.
Figure 10. Percentage contribution to hydrogen production by participating reactions for the COX based mechanism (autothermal temperatures given in the figure for each S/C).
At Φ=0.31, hydrogen production stemmed from DEC, SR, WGS and C-CG1. As expected more hydrogen was produced from WGS as the S/C ratio was increased from 1 to 4. A similar trend was observed at Φ=0.46 except for the fact that there was no C-GS1 reaction. Carbon formed at the higher equivalence ratio is removed via reaction with oxygen (C-GS2).
For this mechanism, bio-oil was consumed almost in the same proportion among COX, SR and DEC independently of the amount of oxygen, as illustrated in Figure 11.
Figure 11. Influence of S/C ratio and oxygen on Bio-oil consuming reactions.
This mechanism relies on gas phase decomposition and combustion which occur significantly at all ATR conditions examined. Table 8 gives a summary of the different reactions that dominate in the COX based mechanism at various temperature ranges. A suitable catalyst for this mechanism will have to be very selective for COX, DEC, SR and C-GS1, as well as offering high thermal stability, although at high Φ, all reactions are expected to become homogeneous (non-catalytic) due to high autothermal temperatures.
Table 8. Summary of different COX mechanism based on temperature range. Only reactions which contribute to equilibrium products are included.
The POX and COX based mechanisms discussed above highlight the fact that there may be several routes leading to the formation of the desired H2 and CO products autothermally. Such schemes are typical for systems where several reaction equilibria occur simultaneously [70]. For the POX based mechanism to be realistic it has to be completely catalytic. An appropriate choice of catalyst can lead to the suppression of undesirable side reactions and products like carbon (coke). The COX based mechanism relies on homogenous oxidation and decomposition. Thermal decomposition accounts for about 30% of the bio-oil consumption and this can prove challenging to manage due to excessive carbon (coke) formation on reactor walls and catalyst and can be difficult to completely eliminate by gasification (depending on type of carbon formed). In addition, it might be necessary to operateunder conditions of Φ>ΦPOX in order to compensate for heat loss and feed preheating (sensible heat) [71].
4. Conclusion
This work demonstrates that ATR is, in theory, a viable process for the production of hydrogen-rich syngas from bio-oil. Using bio-oil surrogate mixtures with different compositions, it was established that hydrogen yield and concentration of other equilibrium products were insensitive to actual chemical composition. The molar elemental composition proved to be the determining factor for equilibrium hydrogen and syngas yield. The possibility of generating syngas with different H2 and CO compositions by varying the S/C ratio and the equivalence ratio makes ATR of bio-oil a feasible option for applications like fuel cells and chemical synthesis. Mechanisms were proposed to account for equilibrium product yields. A POX based mechanism was proposed in which H2 was produced from POX, SR and WGS. For this POX based mechanism, ATR can be viewed as partial oxidation combined with WGS instead of the more traditional notion of exothermic oxidation coupled with endothermic SR. Another mechanism validated was the COX based mechanism in which thermal decomposition accounted for about 30% of bio-oil consumption with hydrogen production assured by decomposition to carbon, steam reforming, water gas shift and carbon gasification reactions.
The equilibrium calculations performed in this study do not take into consideration the kinetic aspects of the reactions involved and could prove unrealistic in real ATR reactors [72]. The proposed mechanisms can only occur when equilibrium is attained for example working at low space velocities and high catalytic activity. Future work should therefore focus on kinetic studies and the influence of other process parameters like pressure and space velocity.
Conflict of interest
All authors declare no conflicts of interest in this paper.
References
[1]
Naik SN, Goud VV, Rout PK, et al. (2010) Production of first and second generation biofuels: A comprehensive review. Renew sust energ rev 14: 578-597. doi: 10.1016/j.rser.2009.10.003
[2]
Ni M, Leung DYC, Leung MKH, et al. (2006) An overview of hydrogen production from biomass. Fuel process technol 87: 461-472.
[3]
Ni M, Leung DYC, Leung MKH, et al. (2006) An overview of hydrogen production from biomass. Fuel process technol 87: 461-472.
[4]
Wu C, Huang Q, Sui M, et al. (2008) Hydrogen production via catalytic steam reforming of fast pyrolysis bio-oil in a two-stage fixed bed reactor system. Fuel process technol 89: 1306-1316.
[5]
Chattanathan SA, Adhikari S, Abdoulmoumine N (2012) A review on current status of hydrogen production from bio-oil. Renew sust energ rev 16: 2366-2372.
[6]
Czernik S, Bridgwater AV (2004) Overview of Applications of Biomass Fast Pyrolysis Oil. Energ fuel 18: 590-598. doi: 10.1021/ef034067u
[7]
Kamm B, Kamm M, Gruber PR, et al. (2005) Biorefinery Systems—An Overview. In: Kamm B, R GP, Kamm M, editors. Biorefineries-Industrial Processes and Products: Status Quo and Future Directions. Weinheim: Wiley-VCH Verlag GmbH, pp. 1-40.
[8]
Jacobson K, Maheria KC, Dalai AK (2013) Bio-oil valorization: A review. Renew sust energ rev 23: 91-106. doi: 10.1016/j.rser.2013.02.036
[9]
Mantilla SV, Gauthier-Maradei P, Gil PÁ, et al. (2014) Comparative study of bio-oil production from sugarcane bagasse and and palm empty fruit bunch: Yield optimization and bio-oil characterization. J anal appl pyrol 108: 284-294. doi: 10.1016/j.jaap.2014.04.003
[10]
Czernik S, Evans R, French R (2007) Hydrogen from biomass-production by steam reforming of biomass pyrolysis oil. Cataly today 129: 265-268. doi: 10.1016/j.cattod.2006.08.071
[11]
Chang SH (2014) An overview of empty fruit bunch from oil palm as feedstock for bio-oil production. Biomass bioenerg 62: 174-181. doi: 10.1016/j.biombioe.2014.01.002
[12]
Isahak WNRW, Hisham MWM, Yarmo MA, et al. (2012) A review on bio-oil production from biomass by using pyrolysis method. Renew sust energ rev 16: 5910-5923. doi: 10.1016/j.rser.2012.05.039
[13]
Bridgwater AV (2012) Review of fast pyrolysis of biomass and product upgrading. Biomass bioenerg 38: 68-94.
[14]
Abdullah N, Gerhauser H (2008) Bio-oil derived from empty fruit bunches. Fuel 87: 2606-2613. doi: 10.1016/j.fuel.2008.02.011
[15]
Sulaiman F, Abdullah N (2011) Optimum conditions for maximising pyrolysis liquids of oil palm empty fruit bunches. Energy 36: 2352-2359. doi: 10.1016/j.energy.2010.12.067
[16]
Kim SW, Koo B, Ryu J, et al. (2013) Bio-oil from the pyrolysis of palm and Jatropha wastes in a fluidized bed. Fuel process technol 108: 118-124.
[17]
Autaa M, Erna LM, Hameeda BH (2014) Fixed-bed catalytic and non-catalytic empty fruit bunch biomass pyrolysis. J anal appl pyrol 107: 67-72. doi: 10.1016/j.jaap.2014.02.004
[18]
Abdullah N, Gerhauser H, Sulaiman F (2010) Fast pyrolysis of empty fruit bunches. Fuel 89: 2166-2169.
[19]
Akhtar J, Amin NAS (2011) A review on process conditions for optimum bio-oil yield in hydrothermal liquefaction of biomass. Renew sust energ rev 15: 1615-1624. doi: 10.1016/j.rser.2010.11.054
[20]
Akhtar J, Kuang SK, Amin NS (2010) Liquefaction of empty palm fruit bunch (EPFB) in alkaline hot compressed water. Renewable energy 35: 1220-1227.
[21]
Sukiran MA, Chin CM, Bakar N (2009) Bio-oils from Pyrolysis of Oil Palm Empty Fruit Bunches. Am j appl sci 6: 869-875. doi: 10.3844/ajassp.2009.869.875
[22]
Pimenidou P, Dupont V (2012) Characterisation of palm empty fruit bunch (PEFB) and pinewood bio-oils and kinetics of their thermal degradation. Bioresource technol 109: 198-205. doi: 10.1016/j.biortech.2012.01.020
[23]
Khor KH, Lim KO, Zainal ZA (2009) Characterization of bio-oil: a by-product from slow pyrolysis of oil palm empty fruit bunches. Am j appl sci 6: 1647-1652.
[24]
Ayabe S, Omoto H, Utaka T, et al. (2003) Catalytic autothermal reforming of methane and propane over supported metal catalysts. Appl Catal a-gen 241: 261-269.
[25]
Rennard DC, Dauenhauer PJ, Tupy SA, et al. (2008) Autothermal Catalytic Partial Oxidation of Bio-Oil Functional Groups: Esters and Acids. Energ fuel 22: 1318-1327. doi: 10.1021/ef700571a
[26]
Dybkjaer I (1995) Tubular reforming and autothermal reforming of natural gas - an overview of available processes. Fuel process technol 42: 85-107.
[27]
Rabenstein G, Hacker V (2008) Hydrogen for fuel cells from ethanol by steam-reforming, partial-oxidation and combined auto-thermal reforming: A thermodynamic analysis. J Power sources 185: 1293-1304. doi: 10.1016/j.jpowsour.2008.08.010
Jonga Md, Reindersa AHME, Kok JBW, et al. (2009) Optimizing a steam-methane reformer for hydrogen production. Int j hydrogen energ 34: 285-292. doi: 10.1016/j.ijhydene.2008.09.084
[30]
Adhikari S, Fernando S, Gwaltney SR, et al. (2007) Athermodynamic analysis of hydrogen production by steam reforming of glycerol. Int j hydrogen energ 32: 2875-2880. doi: 10.1016/j.ijhydene.2007.03.023
[31]
Vagia EC, Lemonidou AA (2008) Thermodynamic analysis of hydrogen production via autothermal steam reforming of selected components of aqueous bio-oil fraction. Int j hydrogen energ 33: 2489-2500. doi: 10.1016/j.ijhydene.2008.02.057
[32]
Xie J, Su D, Yin X, et al. (2011) Thermodynamic analysis of aqueous phase reforming of three model compounds in bio-oil for hydrogen production. Int j hydrogen energ 36: 15561-15572. doi: 10.1016/j.ijhydene.2011.08.103
[33]
Garcia L, French R, Czernik S, et al. (2000) Catalytic steam reforming of bio-oils for the production of hydrogen: effects of catalyst composition. Appl catal a-gen 201: 225-239.
[34]
Trane R, Dahl S, Skjøth-Rasmussen MS, et al. (2012) Catalytic steam reforming of bio-oil. Int j hydrogen energ 37: 6447-6472.
[35]
Rostrup-Nielsen JR, Sehested J, Nørskov JK (2002) Hydrogen and synthesis gas by steam- and C02 reforming. Advance catalysis 47: 65-139.
[36]
Aasberg-Petersen K, Dybkjær I, Ovesen CV, et al. (2011) Natural gas to synthesis gas – Catalysts and catalytic processes. J nat gas sci eng 3: 423-459. doi: 10.1016/j.jngse.2011.03.004
[37]
Czernik S, French R (2014) Distributed production of hydrogen by auto-thermal reforming of fast pyrolysis bio-oil. Int j hydrogen energ 39: 744-750. doi: 10.1016/j.ijhydene.2013.10.134
[38]
Rennard D, French R, Czernik S, et al. (2010) Production of synthesis gas by partial oxidation and steam reforming of biomass pyrolysis oils. Int j hydrogen energ 35: 4048-4059. doi: 10.1016/j.ijhydene.2010.01.143
[39]
Kolios G, Frauhammer J, Eigenberger G (2000) Autothermal fxed-bed reactor concepts. Chem eng sci 55: 5945-5967. doi: 10.1016/S0009-2509(00)00183-4
[40]
Aasberg-Petersen K, Christensen TS, Stub Nielsen C, et al. (2003) Recent developments in autothermal reforming and pre-reforming for synthesis gas production in GTL applications. Fuel Process technol 83: 253-261.
[41]
Martin S, Wörner A (2011) On-board reforming of biodiesel and bioethanol for high temperature PEM fuel cells: Comparison of autothermal reforming and steam reforming. J power sources 196: 3163-3171. doi: 10.1016/j.jpowsour.2010.11.100
[42]
Ruiz JAC, Passos FB, Bueno JMC, et al. (2008) Syngas production by autothermal reforming of methane on supported platinum catalysts. Appl Catal a-gen 334: 259-267.
[43]
Semelsberger T, Brown F, Borup RL, et al. (2004) Equilibrium products from autothermal processes for generating hydrogen-rich fuel-cell feeds. Int j hydrogen energ 29: 1047-1064. doi: 10.1016/S0360-3199(03)00214-3
[44]
Zahedi Nezhada M, Rowshanzamira S, Eikanic MH (2009) Autothermal reforming of methane to synthesis gas: Modeling and simulation. Int j hydrogen energ 34: 1292-1300.
[45]
Haynes DJ, Shekhawat D (2011) Chapter 6 Oxidative Steam Reforming. In: Shekhawat D, Spivey JJ, Berry DA, editors. Fuel Cells: Technologies for Fuel Processing. Amsterdam: Elsevier: 129-190.
[46]
Garcia La, French R, Czernik S, et al. (2000) Catalytic steam reforming of bio-oils for the production of hydrogen: effects of catalyst composition. Appl catal a-gen 201: 225-239.
[47]
Wu C, Liu R (2010) Carbon deposition behavior in steam reforming of bio-oil model compound for hydrogen production. Int j hydrogen energ 35: 7386-7398. doi: 10.1016/j.ijhydene.2010.04.166
[48]
Hou T, Yuan L, Ye T, et al. (2009) Hydrogen production by low-temperature reforming of organic compounds in bio-oil over a CNT-promoting Ni catalyst. Int j hydrogen energ 34: 9095-9107. doi: 10.1016/j.ijhydene.2009.09.012
[49]
Resende KA, Ávila-Neto CN, Rabelo-Neto RC, et al. (2015) Thermodynamic analysis and reaction routes of steam reforming of bio-oil aqueous fraction. Renewable energy 80: 166-176. doi: 10.1016/j.renene.2015.01.057
[50]
Wang S, Cai Q, Zhang F, et al. (2014) Hydrogen production via catalytic reforming of the bio-oil model compounds: Acetic acid, phenol and hydroxyacetone. Int j hydrogen energ 39: 18675-18687. doi: 10.1016/j.ijhydene.2014.01.142
[51]
Latifi M, Berruti F, Briens C (2014) Non-catalytic and catalytic steam reforming of a bio-oil model compound in a novel “Jiggle Bed” Reactor. Fuel 129: 278-291. doi: 10.1016/j.fuel.2014.03.053
[52]
García-García I, Acha E, Bizkarra K, et al. (2015) Hydrogen production by steam reforming of m-cresol, a bio-oil model compound, using catalysts supported on conventional and unconventional supports. Int j hydrogen energ 40: 14445–14455.
[53]
Zin RM, Lea-Langton A, Dupont V, et al. (2012) High hydrogen yield and purity from palm empty fruit bunch and pine pyrolysis oils. Int j hydrogen energ 37: 10627-10638. doi: 10.1016/j.ijhydene.2012.04.064
[54]
Sembiring KC, Rinaldi N, Simanungkalit SP (2015) Bio-oil from Fast Pyrolysis of Empty Fruit Bunch at Various Temperature. Energy procedia 65: 162-169. doi: 10.1016/j.egypro.2015.01.052
[55]
Zin RM, Ross AB, Jones JM, et al. (2015) Hydrogen from ethanol reforming with aqueous fraction of pine pyrolysis oil with and without chemical looping. Bioresour technol 176: 257-266. doi: 10.1016/j.biortech.2014.11.034
[56]
Hanika J, Lederer J, Tukac V, et al. (2011) Hydrogen production via synthetic gas by biomass/oil partial oxidation. Chem eng j 176-177: 286-290. doi: 10.1016/j.cej.2011.06.050
[57]
Lima da Silva A, Malfatti CdF, Müller IL (2009) Thermodynamic analysis of ethanol steam reforming using Gibbs energy minimization method: A detailed study of the conditions of carbon deposition. Int j hydrogen energ 34: 4321-4330. doi: 10.1016/j.ijhydene.2009.03.029
[58]
Lwin Y (2000) Chemical Equilibrium by Gibbs Energy Minimization on Spreadsheets. Int j eng edu 16: 335-339.
[59]
Gordon S, McBride BJ (1994) Computer program for calculation of complex chemical equilibrium compositions and applications: National Aeronautics and Space Administration.
[60]
Haynes DJ, Shekhawat D (2011) Oxidative Steam Reforming. Fuel cells: technologies for fuel processing: 129-190.
[61]
Smith MW, Shekhawat D (2011) Catalytic Partial Oxidation. Fuel cells: technologies for fuel processing: 73-128.
[62]
Nahar GA (2010) Hydrogen rich gas production by the autothermal reforming of biodiesel (FAME) for utilization in the solid-oxide fuel cells: A thermodynamic analysis. Int j hydrogen energ 35: 8891-8911. doi: 10.1016/j.ijhydene.2010.05.042
[63]
Enger BC, Lødeng R, Holmen A (2008) A review of catalytic partial oxidation of methane to synthesis gas with emphasis on reaction mechanisms over transition metal catalysts. Appl catal a-gen 346: 1-27.
[64]
Williams MC (2011) Fuel Cells. Fuel cells: technologies for fuel processing: 11-27.
[65]
Pant K, Gupta RB (2008) Fundamentals and Use of Hydrogen as a Fuel. In: Gupta RB, editor. Hydrogen Fuel Production, Transport, and Storage: CRC Press.
[66]
Dejong M, Reinders A, Kok J, et al. (2009) Optimizing a steam-methane reformer for hydrogen production. Int j hydrogen energ 34: 285-292.
[67]
Wilhelm DJ, Simbeck DR, Karp AD, et al. (2001) Syngas production for gas-to-liquids applications technologies, issues and outlook. Fuel process technol 71: 139-148.
[68]
Aasberg-Petersen K, Bak Hansen JH, Christensen TS, et al. (2001) Technologies for large-scale gas conversion. Appl catal a-gen 221: 379-387.
[69]
Rostrup-Nielsen JR (2000) New aspects of syngas production and use. Cataly today 63: 159-164. doi: 10.1016/S0920-5861(00)00455-7
[70]
Enger BC, Lødeng R, A H (2008) A review of catalytic partial oxidation of methane to synthesis gas with emphasis on reaction mechanisms over transition metal catalysts. Appl catal a-gen 346: 1-27.
[71]
Kumar R, Ahmed S, Krumpelt M (1996) The low-temperature partial oxidation reforming of fuels for transportation fuel cell systems. Fuel cell seminar. Kissimmee, FL (United States), 17-20 Nov 1996: Argonne National Laboratory, Argonne, IL.
[72]
Vagia E, Lemonidou A (2007) Thermodynamic analysis of hydrogen production via steam reforming of selected components of aqueous bio-oil fraction. Int j hydrogen energ 32: 212-223. doi: 10.1016/j.ijhydene.2006.08.021
This article has been cited by:
1.
Walid Nabgan, Tuan Amran Tuan Abdullah, Ramli Mat, Bahador Nabgan, Yahya Gambo, Sugeng Triwahyono,
Influence of Ni to Co ratio supported on ZrO2 catalysts in phenol steam reforming for hydrogen production,
2016,
41,
03603199,
22922,
10.1016/j.ijhydene.2016.10.055
2.
J. Spragg, T. Mahmud, V. Dupont,
Hydrogen production from bio-oil: A thermodynamic analysis of sorption-enhanced chemical looping steam reforming,
2018,
43,
03603199,
22032,
10.1016/j.ijhydene.2018.10.068
3.
Lifita N. Tande, Erik Resendiz-Mora, Valerie Dupont,
Bioh2, Heat and Power from Palm Empty Fruit Bunch via Pyrolysis-Autothermal Reforming: Plant Simulation, Experiments, and CO2 Mitigation,
2021,
14,
1996-1073,
4767,
10.3390/en14164767
4.
Pedro Cerqueira, M.A. Soria, Luis M. Madeira,
Valorization of olive mill wastewater via autothermal reforming for hydrogen production,
2023,
219,
09601481,
119502,
10.1016/j.renene.2023.119502
Table 4. Model PEFB bio-oils obtained using different combinations of acetic acid, phenol, levoglucosan, palmitic acid and furfural. Target composition: C0.3724H0.5304O0.0886.
BOS1
BOS2
BOS3
BOS4
BOS5
Name
mol (%)
mol (%)
mol (%)
mol (%)
mol (%)
Acetic acid
12.2
23.2
0.0
0.3
0.0
Phenol
35.1
45.2
36.0
65.3
12.2
Levoglucosan
2.2
0.3
8.6
16.4
0.0
Palmitic acid
20.0
16.3
21.6
17.5
25.5
Furfural
30.5
15.0
33.8
0.5
62.2
TOTAL
100.0
100.0
100.0
100.0
100.0
C
0.3724
0.3724
0.3724
0.3724
0.3740
H
0.5390*
0.5390*
0.5372
0.5390*
0.5374*
O
0.0886
0.0886
0.0904*
0.0886
0.0886
Maximum relative error* (%)
1.6
1.6
2
1.3
1.3
* For a given surrogate, only the maximum elemental relative error, when compared to the average value, is shown. For BOS1, 2, 4 and 5, error shown for H, for BOS3, error shown for O.
* For a given surrogate, only the maximum elemental relative error, when compared to the average value, is shown. For BOS1, 2, 4 and 5, error shown for H, for BOS3, error shown for O.
Name
Abbreviation
Reaction
1
Partial oxidation
POX
CnHmOk+(n-k)/2O2→nCO+m/2H2
2
Complete oxidation
COX
CnHmOk+(n+m/4–k/2)O2→nCO2+m/2H2O
3
Steam reforming
SR
CnHmOk+(n-k)H2O→nCO+(2n+m-2k)/2H2
4
Decomposition
DEC
CnHmOk→kCO+m/2H2 +(n-k)C
5
Water gas shift*
WGS
CO+H2O↔CO2+H2
6
Boudouard reaction
BO-RX
2CO→CO2+C
7
Methanation of C(s)
MEN
C(s)+2H2→CH4
8
Carbon gasification 1
C-GS1
C(s)+H2O→CO+H2
9
Carbon gasification 2
C-GS2
C(s) +0.5O2→CO
10
Methane steam reforming
ME-SR
CH4+H2O→CO+3H2
11
Carbon monoxide oxidation
CO-OX
CO+0.5O2→CO2
12
Hydrogen oxidation
H-OX
H2+0.5O2→H2O
* The reverse of the water gas shift reaction (R-WGS) was used in some cases
Equivalence ratio, Φ
S/C =1
S/C = 2
S/C = 3
S/C = 4
Φ1 ≈ 0.15
874 ± 4
804 ± 4
755 ± 4
715 ± 4
Φ2 ≈ 0.31
1204 ± 19
1043 ± 13
942 ± 10
871 ± 7
Φ3 ≈ 0.46
1963 ± 19
1587 ± 14
1370 ± 11
1225 ± 9
Φ4 ≈ 0.61
2521 ± 12
2079 ± 13
1773 ± 11
1566 ± 9
T < 1100 K
1200 < T < 1600 K
T > 1800 K
Steam content
S/C > 2
2 < S/C < 4
S/C < 2
Oxygen content
Φ<ΦPOX
Φ>ΦPOX
Φ≫ΦPOX
Reaction condition
Catalytic
Catalytic and homogenous
Homogenous
Reactions
POX, WGS, BO-RX, MEN, C-GS1
POX, SR, C-GS2 CO-OX
POX, SR, R-WGS, H-OX, CO-OX
T < 1100 K
1200 < T < 1600 K
T > 1800 K
Steam content
S/C > 2
2 < S/C < 3
S/C < 2
Oxygen content
Φ<ΦPOX
Φ>ΦPOX
Φ≫ΦPOX
Reaction condition condition
Catalytic and homogenous
Catalytic and homogenous
Homogenous
Reactions
COX, DEC, SR, WGS, C-GS1,
COX, DEC, SR, C-GS1
COX, DEC, SR, R-WGS, C-GS2, H-OX
Figure 1. Mean hydrogen yield (BOS1-5) and standard error for Φ (0.15−0.61) and S/C (1−4)
Figure 2. Influence of S/C ratio and the equivalence ratio on the amount of hydrogen produced during the ATR of BOS2. (a) Φ range 0.15-0.61, (b) refined Φ scale (0.21-0.39)
Figure 3. Influence of S/C ratio and O2/C ratio on the selectivity of carbon and carbon containing products during the ATR of BOS2 at 1 atm. a) S/C = 1 b) S/C = 2 c)S/C = 3 d) S/C = 4
Figure 4. Influence of S/C ratio and Φ (amount of oxygen) on the total H2 + CO obtained during ATR of BOS2 at 1 atm
Figure 5. Plot of S/C ratio versus H2/CO ratio at different values of Φ during ATR of BOS2. Equilibrium simulation carried out at 1 atm. The maximum H2/CO ratio shown is 10
Figure 6. Comparing actual equilibrium hydrogen yield with predicted yield using POX based mechanism
Figure 7. Percentage contribution to hydrogen production by participating reactions for the POX mechanism at (autothermal temperatures given in the figure for each S/C)
Figure 8. Influence of S/C ratio and oxygen on Bio-oil consuming reactions
Figure 9. Plots of predicted hydrogen and actual equilibrium hydrogen showing near match with maximum relative error of 0.1% for the COX based mechanism
Figure 10. Percentage contribution to hydrogen production by participating reactions for the COX based mechanism (autothermal temperatures given in the figure for each S/C)
Figure 11. Influence of S/C ratio and oxygen on Bio-oil consuming reactions