Chapter 14
Chapter 14
This example assumes that the feed specification if bio-oil is known and the objective is to
synthesise a flowsheet in tandem with the establishment of basic mass and energy balances.
Step 1
feed into syngas and fuel products (methanol). Bio-oil can be modelled using three chemical
constituents, i.e. acetic acid, acetol and guaiacol1-3. The methodology for selecting a set of
suitable representative components to model bio-oil can be found in chapter 10, Gibbs free
Table 1. List of components involved in the process modelling for bio-oil integrated
Category Component
H2
H2O
Intermediate Product
CO
(Syngas)
CO2
CH4
1
Biorefineries and Chemical Processes: Design, Integration and Sustainability Analysis, First Edition.
Jhuma Sadhukhan, Kok Siew Ng and Elias Martinez Hernandez.
© 2014 John Wiley & Sons, Ltd. Published 2014 by John Wiley & Sons, Ltd.
Gasifying medium O2
Physical property package - Since the modelling of integrated gasification system only
involves homogenous gaseous reactions and separation, the choice of suitable physical
property package is rather straightforward. Most of the common physical property package
such as Non-Random Two Liquids (NRTL) and Peng-Robinson can be applied. Most
importantly, the results obtained from the simulation should be compared and validated to
Gasifier
Select equilibrium reactor models, if there is no kinetic data available. There are four types of
reactor (REquil), Gibbs reactor (RGibbs) and Yield reactor (RYield). Which reactor should
be chosen for gasifier, water-gas shift reactor and synthesis reactor, depends on what
pressure are known. However, chemical reactions may not always be easy to define. Table
shows the reactor models adopted in Aspen Plus and the rationale for choosing the respective
models. It should be reminded that the choice of models is case-specific depending on the
desired level of details of the model and the need for subsequent analysis.
2
Biorefineries and Chemical Processes: Design, Integration and Sustainability Analysis, First Edition.
Jhuma Sadhukhan, Kok Siew Ng and Elias Martinez Hernandez.
© 2014 John Wiley & Sons, Ltd. Published 2014 by John Wiley & Sons, Ltd.
Reactor Aspen Rationale
Plus
model
Gasifier RGibbs Gasification reactions are complex (r.f. Chapter 9). REquil and
Water-gas REquil It can be assumed that only one reversible reaction occurs in the
shift reactor, i.e. water-gas shift. RGibbs can also be used when the
flow rate.
Methanol REquil For modelling methanol synthesis reaction, REquil is the most
synthesis suitable choice, since the synthesis reaction involves three main
reactor reversible reactions. RStoic and RYield are not appropriate in this
3
Biorefineries and Chemical Processes: Design, Integration and Sustainability Analysis, First Edition.
Jhuma Sadhukhan, Kok Siew Ng and Elias Martinez Hernandez.
© 2014 John Wiley & Sons, Ltd. Published 2014 by John Wiley & Sons, Ltd.
RGibbs may not be a right choice, because the product is generated
1. Gasification modelling
RGibbs model in Aspen Plus is used to model the bio-oil gasification process (GASIFIER),
shown in Figure . Two operating parameters are needed in RGibbs, i.e. temperature and
pressure. In this case, use a temperature of 1112°C and pressure of 30 bar. HE1 is used for
The bio-oil feed stream and oxygen stream are specified according to the values shown in
Stream BIO-OIL O2
Temperature (°C) 25 25
4
Biorefineries and Chemical Processes: Design, Integration and Sustainability Analysis, First Edition.
Jhuma Sadhukhan, Kok Siew Ng and Elias Martinez Hernandez.
© 2014 John Wiley & Sons, Ltd. Published 2014 by John Wiley & Sons, Ltd.
Acetic Acid 0.109
Acetol 0.109
Guaiacol 0.109
H2O 0.673
O2 1.0
The product composition is calculated based on Gibbs free energy minimisation method
(chapter 10). By default, RGibbs considers all the components as possible products. Simulate
the reaction. The product flow rate and composition is calculated, shown in Table .
Stream SYNGAS
Pressure (bar) 30
H2 0.307
H2O 0.251
CO 0.318
CO2 0.124
In a bio-oil integrated gasification system, the syngas should be free (or has negligible
5
Biorefineries and Chemical Processes: Design, Integration and Sustainability Analysis, First Edition.
Jhuma Sadhukhan, Kok Siew Ng and Elias Martinez Hernandez.
© 2014 John Wiley & Sons, Ltd. Published 2014 by John Wiley & Sons, Ltd.
reaction. This high level of heat can be recovered to generate steam and the steam can be
flown through to back pressure steam turbine to generate power. Water-gas shift (WGS)
reaction is important for adjusting the H2/CO molar ratio of the syngas, so that it is suitable to
be used in the synthesis reaction. The stoichiometric number for methanol synthesis, R
defined as (H2 − CO2) / (CO + CO2) should be manipulated to 2. Water-gas shift reaction is
also vital in concentrating CO2 in the outlet stream so that it can be removed by CO2 capture
unit (CO2SEP) more economically. Water removal unit (H2OREM) should be in place to
reduce the water content in the syngas, before the CO2 capture and synthesis reaction. This
unit can be modelled by using flash separators. The degree of water-gas shift, degree of water
removal as well as the degree of CO2 capture are three major factors in determining the
stoichiometric number, R and hence the yield of methanol from the synthesis reaction. After
capturing, CO2 is compressed to supercritical condition (~74 bar and 32°C), 80 bar and 35°C
are used in the modelling. This condition is needed to reduce the volume of CO2 to ease
transportation through pipeline and storage. The syngas conditioning system is shown in
Figure .
6
Biorefineries and Chemical Processes: Design, Integration and Sustainability Analysis, First Edition.
Jhuma Sadhukhan, Kok Siew Ng and Elias Martinez Hernandez.
© 2014 John Wiley & Sons, Ltd. Published 2014 by John Wiley & Sons, Ltd.
Table shows the model specification for each unit operation in the syngas conditioning
system.
cooler)
(CO2 capture)
7
Biorefineries and Chemical Processes: Design, Integration and Sustainability Analysis, First Edition.
Jhuma Sadhukhan, Kok Siew Ng and Elias Martinez Hernandez.
© 2014 John Wiley & Sons, Ltd. Published 2014 by John Wiley & Sons, Ltd.
3. Methanol Synthesis Reaction Modelling
Methanol reaction is carried out at 100 bar and 250°C. Due to the high pressure requirement
in the methanol reactor, the syngas stream (SYN3) needs to be compressed in SYNGCOMP
from 30 bar to 100 bar. After reaction, the product is a mixture of light gases (e.g. H2, CO)
and methanol. The pressure of the product stream (METPOUT) is lowered in SYNGEXP to
40 bar (minimum pressure to maintain the stream at gaseous phase since presence of liquid
has an adverse effect on the expander) and subsequently separated into gas and liquid streams
in METSEP. The liquid can be sent to distillation column for further methanol purification.
Gases can be combusted in gas turbine to generate power. Methanol synthesis reactor model
is shown in Figure .
Table shows the model specification for each unit operation in the methanol synthesis
process.
8
Biorefineries and Chemical Processes: Design, Integration and Sustainability Analysis, First Edition.
Jhuma Sadhukhan, Kok Siew Ng and Elias Martinez Hernandez.
© 2014 John Wiley & Sons, Ltd. Published 2014 by John Wiley & Sons, Ltd.
Table 6. Model specification for methanol synthesis process.
Table .
9
Biorefineries and Chemical Processes: Design, Integration and Sustainability Analysis, First Edition.
Jhuma Sadhukhan, Kok Siew Ng and Elias Martinez Hernandez.
© 2014 John Wiley & Sons, Ltd. Published 2014 by John Wiley & Sons, Ltd.
Final Flowsheet and Stream Results
10
Biorefineries and Chemical Processes: Design, Integration and Sustainability Analysis, First Edition.
Jhuma Sadhukhan, Kok Siew Ng and Elias Martinez Hernandez.
© 2014 John Wiley & Sons, Ltd. Published 2014 by John Wiley & Sons, Ltd.
Table 7. Stream results for integrated gasification and methanol synthesis model.
Stream Name
Component
Mole BIO-OIL O2 SYNGAS SYN1 SYN2 WATER SYN3 SYNGIN METPOUT METGAS METLIQ CO2
Fraction
Acetol 0.109
Guaiacol 0.109
H2O 0.673 0.253 0.099 0.004 0.912 0.005 0.005 0.007 0.020
CO2 0.122 0.276 0.298 0.082 0.060 0.060 0.102 0.138 0.037 1.000
O2 1.000
Molar flow 2.10 1.14 6.20 6.20 5.55 0.65 4.14 4.14 2.46 1.58 0.88 1.41
11
Biorefineries and Chemical Processes: Design, Integration and Sustainability Analysis, First Edition.
Jhuma Sadhukhan, Kok Siew Ng and Elias Martinez Hernandez.
© 2014 John Wiley & Sons, Ltd. Published 2014 by John Wiley & Sons, Ltd.
rate
(kmol s−1)
Mass flow
303.76 131.32 435.08 435.08 387.87 47.22 164.92 164.92 164.92 63.64 101.28 222.95
−1
rate (t h )
12
Biorefineries and Chemical Processes: Design, Integration and Sustainability Analysis, First Edition.
Jhuma Sadhukhan, Kok Siew Ng and Elias Martinez Hernandez.
© 2014 John Wiley & Sons, Ltd. Published 2014 by John Wiley & Sons, Ltd.
2. Description of bio-oil integrated gasification system for the production of FT liquids
Using the same analogy as in the methanol synthesis process, develop a simulation model of
an integrated gasification and Fischer-Tropsch (FT) synthesis system. The current case
Step 1
Chemical components –The current model assumes that FT product consists of C1-C30
desired level of detail of the model. The component list is updated and shown in Table .
Table 8. List of components involved in the process modelling for bio-oil integrated
Category Component
H2
H2O
Intermediate Product
CO
(Syngas)
CO2
CH4
Gasifying medium O2
13
Step 2 and Step 3
The fundamental modelling technique for gasification is same as in the methanol synthesis
process. Only the differences in terms of conditioning the syngas composition and FT
Figure shows the integrated gasification and FT synthesis flowsheet, consisting of a gasifier
unit, syngas cleaning and conditioning and FT synthesis reactor. The stream results are
shown in
Table . The model specification for the integrated gasification and FT synthesis system is
shown in Table .
14
Final Flowsheet and Stream Results
15
Table 9. Stream results for integrated gasification and FT synthesis model.
mole fraction O2 BIO-OIL SYNGAS WATER1 WATER2 SYN CO2 SYNGIN FTOUT FTLIQ RECYC FTGAS
Acetol 0.109
Guaiacol 0.109
H2O 0.673 0.251 0.986 0.864 0.001 0.002 0.161 0.866 0.003 0.003
CO2 0.124 0.008 0.128 0.004 1.000 0.007 0.009 0.006 0.010 0.010
O2 1.000
CH4
Molar flow
1.14 2.10 6.20 0.56 0.30 3.89 1.45 8.06 6.13 1.12 4.18 0.83
−1
rate(kmol s )
16
Table 10. Model specification for integrated gasification and FT synthesis system.
Type: Isentropic
Pressure = 25 bar
Pressure = 30 bar
Pressure = 30 bar
Pressure = 30 bar
Pressure = 30 bar
Pressure = 30 bar
Pressure = 80 bar
Pressure = 25 bar
17
Pressure = 25 bar
Pressure = 30 bar
Pressure = 25 bar
Pressure = 30 bar
Note: Process enhancement through recycling the FT offgas is considered in this case. In
order to result in an overall conversion of CO of 80%, discussed later, the split ratio to
The H2/CO molar ratio required in the FT reaction is approximately 2. WGS is used to adjust
the ratio. In this case, 2 water removal units are used. One of them is placed before the WGS
(H2OREM1) to manipulate the steam required by the water-gas shift reaction while the
second water removal unit (H2OREM2) is used to reduce the amount of water content in
It is assumed that only the straight-chain paraffins C1 to C30 are produced from FT synthesis
reactor. Stoichiometric reactor (RStoic) is selected for modelling FT reaction in Aspen Plus.
18
There are 30 reactions in total, inputs to the reaction list for the reactor model. The general
where n = 1, 2, 3,….., 30, is the carbon number. The fractional conversion of each reaction is
required in each of the reaction input. These are estimated based on the weight distribution of
αn-1 (1−α)2 n, where α the chain growth probability has a value between 0.7-0.9.
Song et al. shows an empirical model to relate α with temperature, as in Equation 11. This
creates a more flexible reactor model. α responds to different temperature and hence provides
a sensible prediction on the FT product distribution. α value in this case is 0.76, which
B 1 0.0039T 533
yCO
A (1)
y H 2 yCO
where
yCO and y H 2 denote the mole fraction of CO and H2, respectively at the inlet stream into FT
reactor.
Figure shows the modelling approach for estimating product distribution from FT synthesis.
19
Fractional conversion (30 reactions)
ASPEN Plus
Syngas
Stoichiometric Reactor
(CO, H2) Simulate Product
model (RStoic) for FT
reaction distribution
T
Match
Excel environment
Figure 6. Modelling approach for Fischer-Tropsch reactor using combined Aspen Plus
By applying ASF relation, the weight and molar distribution of each FT product are obtained.
An overall conversion of the syngas (assume CO) is assumed to be 0.8. The consumption of
CO for each reaction can be determined by solving Equation 2, where the total molar flow
rate of paraffin produced, Fp is determined using “Solver” function in Excel. Ultimately, the
fractional conversion based on CO for each reaction (CO consumption for a reaction divided
n 30
Total CO converted = y n 1
p ,n Fp n (2)
where
20
After specifying all the fractional conversion for each reaction into the reactor model (refer to
Table ), the simulation is thereby completed. The simulation results of the weight and molar
Table 11. Reaction specification in RStoic (Aspen Plus) for modelling Fischer-Tropsch
reaction.
21
17 0.0104 35 H2 + 17 CO → 17 H2O + C17H36
12
10
Weight distribution (wt %)
0
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30
Carbon number, n
(a)
22
25
20
Molar distribution (mol %)
15
10
0
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30
Carbon number, n
(b)
Figure 7. (a) Weight distribution of FT products predicted using ASF model. (b) Molar
distribution can be evaluated, as shown in Figure . These figures show that higher diesel
23
50
45
40
Weight distribution (wt%)
35
T = 160°C, α = 0.98
30
T = 200°C, α = 0.87
25 T = 240°C, α = 0.76
T = 300°C, α = 0.60
20
T = 400°C, α = 0.32
15
10
0
0 5 10 15 20 25 30 35
Carbon number, n
Figure 8. Variation of product distribution with temperature, with constant H2/CO molar
Research in Bio-oil
Aston University (UK) and Dynamotive (Canada) undertakes extensive research and
development of biomass pyrolysis technologies for the production of bio-oil. The Energy
research Centre of the Netherlands (ECN) focuses on the thermochemical production routes
from biomass to syngas and FT diesel production. Other specialist biomass gasification and
pyrolysis research centres include BTG. Gasification of bio-oil and bio-oil and char (bio-
slurry) has received considerable interest and research has been carried out extensively by
24
Test your skills with the following questions.
1. The offgas from methanol synthesis reactor in a bio-oil integrated gasification and
turbine2. Using the simulation model and specification in methanol synthesis process as the
basis, develop and simulate the models for power generation within the system using a
simulation package. Estimate the power generated from the expansion of the offgas through
gas turbine. Assume a discharge pressure of 2 bar and an isentropic efficiency of 90% for the
gas turbine. A combustor is required, operating at 1200°C and 14 bar. Include a heat recovery
steam generator (HRSG) to recover the heat from the exhaust gas of gas turbine, given that
the outlet temperature from the HRSG is 100°C. Air and an auxiliary fuel such as natural gas
can be used to help the combustion of offgas. Report results appropriately, including energy
efficiency.
and methanol synthesis system2 generated in Question 1 to improve the yield of methanol, by
implementing a recycle stream within the methanol synthesis reactor. However, this will
result in a reduction in power generation. Assume that 90% of the offgas from methanol
synthesis reactor can be recycled to obtain higher yield of methanol, while 10% of the offgas
can be sent to the gas turbine power generation system. Determine the amount of methanol
and power produced and compare the results with those from the methanol synthesis process
and Question 1. The same assumptions can be made for the power generation section as in
Question 1.
3. The amount of products generated from the integrated gasification and methanol
synthesis system varies with the type of feed used. Bio-oil can come from different types of
25
biomass feedstocks. The methanol synthesis process and Questions 1-2 used the bio-oil from
poplar wood as the feedstock. Estimate the yields of products (methanol and electricity) if a
different biomass, miscanthus is used. Carry out a simulation study to investigate this, using a
miscanthus are shown in Table . The bio-oil feed can be represented using acetic acid, acetol
and guaiacol. Refer to a study by Ng and Sadhukhan for the procedure of defining the feed
composition3.
value (wt%)
Carbon
processes such as methanol synthesis and power generation. Carry out a process simulation
using a different operating condition of the gasifier, i.e. 1400°C and 20 bar. Report results.
system in terms of removing CO2 from the system and also manipulating the stoichiometric
number for optimum yield of methanol. Comment, qualitatively, the impact on energy, yield
26
References
2. K.S. Ng and J. Sadhukhan, Process integration and economic analysis of bio-oil platform
for the production of methanol and combined heat and power, Biomass Bioenergy, 35, 1153-
1169 (2011).
Fischer-Tropsch and CHP synthesis platform, Biomass Bioenergy, 35, 3218-3234 (2011).
27