M. Mozahar Hossain January 27, 2025: Hydrogen Production and Storage
M. Mozahar Hossain January 27, 2025: Hydrogen Production and Storage
M. Mozahar Hossain
January 27, 2025
When used in a PEMFC, the product gas from the water – gas shift must be further
processed to reduce CO to < 10 ppm.
H2 Production – SMR Mass/Energy/Economy
CHE-300
CHE-404
H2 Production – Steam Methane Reforming
CHE-300
CHE-404
Table of contents
Steam Reforming of Natural Gas
Process Descriptions
Unit Operations
Chemistry
Reactor design
Non-hydrocarbon gases
“wet’= which condense on compression at T ambient
Associated gas is invariably ‘wet’ than ‘dry’ natural gas
H2 Production – SMR Major steps
CHE-300
CHE-404
H2 Production – SMR Process
CHE-300
CHE-404
H2 Production – SMR Process Description
CHE-300
CHE-404
H2 Production – SMR Process Description
CHE-300
CHE-404
H2 Production – SMR Process Description
CHE-300
CHE-404
H2 Production – SMR Major Units Operation
CHE-300
CHE-404
Major units
Hydrotreating/Cracking/pre-
reformer
Primary Reformer
Secondary Reformer
It is accomplished in the
desulfurization unit (DSU).
A small split stream of the H2 In some cases, traces of halides (Cl, Br, F, I)
product is used as a are also present in the NG & these are
hydrogenating agent in this step. removed by an Al2O3 guard bed.
HDS (hydrodesulfurization) –
Unit operation and Chemistry CHE-300
CHE-404
Co–Mo catalyst
290–370°C In some cases, traces of halides (Cl, Br, F, I) are also
present in the NG & these are removed by an Al2O3
H2 Production – SMR Major Units Operation
CHE-300
CHE-404
Major units - Desulfurization unit (DSU)
HDS (hydrodesulfurization) –
Unit operation and Chemistry CHE-300
CHE-404
After the DSU, NG is fed to a catalytic
reforming Unit 2 where it reacts with
steam to produce synthesis gas.
Reactions -
H2 Production – SMR Major Units Operation
CHE-300
CHE-404
Methane Steam Reformer (SMR)
Objective of the steam reforming process
is to produce syngas.
𝐶𝐻4 + 𝐻2 𝑂 ↔ 𝐶𝑂 + 3𝐻2
(2) CO disproportionation
𝐶𝑂 + 𝐶𝑂 ↔ 𝐶 + 2𝐶𝑂2
An addition of the supplemental
steam shifts the reforming reaction
equilibrium away from “C” formation.
To supply heat for the endothermic SMR
reaction, the catalyst is loaded into a
bundle of reactor tubes (~ 15 m long & 12
cm ID) made out of heat-resistant Ni alloy.
The estimated that for the hydrogen plant with the Scholz (another group)- estimated the energy
capacity of 1.5 × 106 Nm3/day, the energy efficiency efficiency of SMR plant at 81.2%.
was 89% (on high heating value [HHV] basis).
Process Descriptions
Unit Operations
Chemistry
Reactor design
Non-hydrocarbon gases
“wet’= which condense on compression at T ambient
Associated gas is invariably ‘wet’ than ‘dry’ natural gas
H2 Production – SMR reactions
CHE-300
CHE-404
H2 Production – SMR reactions
CHE-300
CHE-404
H2 Production – SMR reactions
CHE-300
CHE-404
H2 Production – SMR reactions
CHE-300
CHE-404
H2 Production – SMR reactions
CHE-300
CHE-404
H2 Production – SMR reactions
CHE-300
CHE-404
H2 Production – SMR reactions
CHE-300
CHE-404
H2 Production – SMR reactions
CHE-300
CHE-404
H2 Production – SMR reactions
CHE-300
CHE-404
H2 Production – SMR equilibrium
CHE-300
CHE-404
Composition
2
K1 = = 𝑃 0.6
PCH4 )×(PH2O yCH4 )×(yH2O 0.5 H2 (1 bar)
PCO2 ) × (PH2 yCO2 ) × (yH2 0.4 CH4 (1 bar)
K2 = = 0.3 H2 ( 30 bar)
PCO ) × (PH2O yCO ) × (yH2O 0.2
CH4 ( 30 bar)
0.1
Methane conversion increases with
0.0
temperature. 500 700 900 1100 1300 1500 1700
Temperature (K)
The methane conversion increases with 100
Methane % Conversion
higher S/C ratios (S/C varies from 1 to 5) and 80
decreases with increasing pressures (1 – 20
60
bar pressures were studied).
40
A complete conversion of methane could be 20
achieved around 700 ° C at 1 bar pressure and 0
the S/C ratio of above 2.5, while a 200 400 600 800 1000 1200
temperature of above 900 ° C would be Temperature (℃)
required to achieve the complete methane For Steam/Methane Ratio = 1
conversion at 20 bar pressure. For Steam/Methane Ratio = 2.5
For Steam/Methane Ratio = 5
H2 Production – SMR equilibrium
CHE-300
CHE-404
Composition
0.4
The carbon dioxide content goes
through a maximum. This can be 0.3
explained as follows:
0.2
carbon dioxide is formed in
exothermic reactions only, while it is a 0.1
reactant in the endothermic reactions.
0.0
Hence, at low temperature carbon 500 1000 1500 2000
dioxide is formed, while with Temperature (K)
increasing temperature the
endothermic reactions in which
carbon dioxide is converted become
more important.
[Importance of Thermodynamics and
Kinetics]
H2 Production – SMR equilibrium
CHE-300
CHE-404
H2 Production – SMR reaction equilibrium
CHE-300
CHE-404
Effects of Pressure, Steam and
Temperature – on methane conversion
100 kmol/h of a stream containing 20 mol% CO, 20 mol% H2O, and 60 mol% H2 is combined
with additional steam and fed to a water–gas shift reactor to produce more hydrogen
according to the reaction:
The equilibrium constant of the water–gas shift reaction is 200, and the mole fraction of
hydrogen at the outlet is 0.75. Write and balance the chemical equation. How much
additional steam is fed to the water–gas shift reactor? (Answer: 5.1 kmol/h)
H2 Production – SMR Mass/Mole and Energy Balance
CHE-300
CHE-404
[Link]
H2 Production – SMR Mass/Mole and Energy Balance
CHE-300
CHE-404
Problem 2
A reactor is designed to take in a feed of 200 mol/min of methane and 50% excess steam, and
convert them to CO and H2 with a 65% conversion rate.
Determine the mole flow rates methane, hydrogen and CO in the reactor product stream.
Calculate the Energy Efficiency based on LHV
(For your verification - Answer: methane flow rate at the outlet is 70 mol/min)
[Link]
H2 Production – SMR Mass/Mole and Energy Balance
CHE-300
CHE-404
H2 Production – SMR Mass/Mole and Energy Balance
CHE-300
CHE-404
H2 Production – SMR Mole Balance in the Furnace
CHE-300
CHE-404
H2 Production – SMR Mole Balance in the Furnace
CHE-300
CHE-404
H2 Production – SMR Mole Balance in the Furnace
CHE-300
CHE-404
H2 Production – Example on mole Balance in the Furnace
CHE-300
CHE-404
H2 Production – Example on Mole balance with multiple reactions
CHE-300
CHE-404
H2 Production – Example on mole Balance in the Furnace
CHE-300
CHE-404
H2 Production – Example on mole Balance in the Furnace
CHE-300
CHE-404
H2 Production – SMR equilibrium
CHE-300
CHE-404
H2 Production – Example on Absorber/separation
CHE-300
CHE-404
H2 Production – Example on Absorber/separation
CHE-300
CHE-404
H2 Production – Example on Absorber/separation
CHE-300
CHE-404
H2 Production – Example on WGS reactor and
Absorber/separation mole balance CHE-300
CHE-404
H2 Production – Example on WGS reactor and
Absorber/separation mole balance CHE-300
CHE-404
H2 Production – Example on WGS reactor and
Absorber/separation mole balance CHE-300
CHE-404
H2 Production – Example on reactions and seperation
CHE-300
CHE-404
H2 Production – SMR reactions
CHE-300
CHE-404
H2 Production – SMR reactions
CHE-300
CHE-404
The CO2 Balance in the Hydrogen Plant
CHE-300
CHE-404
Membranes performances in
terms of removal efficiency and
purity are expected to be lower.
The CO2 Removal Technologies from Syngas
CHE-300
CHE-404
When a traditional plant is retrofitted with the addition of
such removal unit, the SMR burners have to be checked for
the new tail gas composition
which, being poorer in CO2, has a higher tendency to NOx
formation, with the consequence that low-NOx burners
Space for a future installation of a de-NOx catalytic system in
the convection section of the heater, could be considered at
the design stage.
As a consequence of the retrofit the PSA will have some
capacity margin due to the reduced inlet flow rate after CO2
removal.
Licensors active in these technologies include: BASF,
Giammarco-Vetrocoke, Linde, Lurgi, UOP, etc.
The CO2 Removal Technologies from Flue gas
CHE-300
CHE-404
CO2 in the SMR flue gas has a very low partial pressure (~0.2
bar abs) and the flue gas actual volumetric flow is 20 times or
more the raw syngas volumetric flow (m3/h).
These two factors & presence of O2, SOx, NOx, dust, makes the
CO2 removal from this stream complicated & expensive.
Nevertheless its quantity is very high & its capture is an
environmental & economical challenge. Main technologies -
Amines (MEA, MDEA); Chilled ammonia
Membranes.
The CO2 Removal Technologies from Flue gas
CHE-300
CHE-404
The basic principles of CO2 removal from flue gas are the
same as from raw H2, but the experience is so far limited.
The environmental concern related to the GHG emissions is
determining a growing interest for this application in many
industries, for example in the power sector.
In practical terms a CO2 removal of 85 - 90% is an acceptable
target, and the CO2 is also in this case obtained at high purity
and low pressure.
No pre-investment is needed in case of retrofit of a traditional
plant, however additional ducting and ID fan addition shall be
considered when the unit is added to the SMR plant.
Licensors active in these technologies are: Axens, Fluor, GTC,
MHI, MTR etc.
The CO2 Capture – Case Study
CHE-300
CHE-404
The CO2 removal efficiency from raw H2 and from flue gas is considered
100% (rounded off) and 90% respectively.
The pure CO2 stream is then compressed to 150 barg for sequestration,
EOR or other utilisation (out of the scope of the study).
Hydrotreating unit
Process Descriptions
Unit Operations
Chemistry
Reactor design
where η is the energy efficiency, EH2 the energy in H2 product, Esteam, 4.8 MPa the 4.8
Mpa steam energy (exported), ENG the NG energy, and Esteam, 2.6 MPa the 2.6 MPa
steam (required).
The estimated that for the hydrogen plant with the capacity of 1.5 × 106 Nm3/day,
the energy efficiency was 89% (on high heating value [HHV] basis).
If the steam were not included in the equation (e.g., 4.8 MPa steam could not be used
by another source), the energy efficiency would decrease to 79.2%. The hydrogen
plant efficiency drops to 69.1% if
If 2.6 MPa steam is considered as an energy input (i.e., it is not produced internally).
Scholz (another group)- estimated the energy efficiency of SMR plant at 81.2%.
H2 Production – Reactor Design
CHE-300
CHE-404
Hydrotreating unit
Product (B)
Boundary layer
Bulk fluid phase
Active sites
Product (B)
Reactant (A)121
Steps in heterogeneous catalytic reactions
CHE-404
Catalyst
Reactant (A) Product (B)
Active site
A section of a porous catalyst 122
4. Design of Isothermal Packed-Bed Reactor
Substituting (3) into (2) gives:
(4a)
(4b)
Remarks:
▪ While catalyst deactivation is unavoidable, as we have seen in Chapter VI, there are
many applications in which it is very slow (say, 2-3 years), so that catalyst activity is
essentially constant
▪ Accordingly, PBR design calculations should take into account variation of pressure
along the reactor, as described by the Ergun equation. The design problem then
becomes more complicated and in most situations we need to combine Equation
(4a) or (4b) with the Ergun equation and solve a set of differential equations
numerically
123
H2 Production – Reactor Design
CHE-300
CHE-404
The SMR reaction results in the production of H2 and CO2 and is a strongly
endothermic reaction and favored at high temperatures.
Heat transfer is critical to ensure fast kinetics as the reaction moves forward.
Integration of SMR reaction and products (H2 or CO2 ) separation steps to shift
the equilibrium and enhance the reaction rate in the forward direction.
This also lowers the reforming temperature to achieve the same conversion
level.
125
4. Design of Isothermal Packed-Bed Reactor
❑In the case of solid catalytic reactions, the amount of catalyst is what
is important for product formation
❑Therefore it is more suitable to define the bulk rate of the catalytic
reaction on the basis of catalyst mass (W) as follows :
❑The bulk rate of the reaction per mass of catalyst can be related to
reaction rate per volume of reactor, (-rA), using the bulk bed density,
ρB [mass of catalyst/volume of reactor]. Therefore,
126
4. Design of Isothermal Packed-Bed Reactor
▪ We need to distinguish between the density of the particle, ρp, and
bulk density of the bed ρB. The two densities are related as follows:
127
4. Design of Isothermal Packed-Bed Reactor
Assumptions:
❑Plug flow of the fluid phase; i.e. the velocity profile is flat across the
bed.
❑Even distribution of flow over the bed cross-section ( no radial
velocity profiles)
❑The heterogeneous nature of the system may cause appreciable
concentration differences between the bulk fluid and internal and/or
external surface of the catalyst particles. These differences are
accounted for using the effectiveness factor concept
❑Pressure drop across the catalyst bed is negligible
❑Catalyst deactivation is insignificant
128
4. Design of Isothermal Packed-Bed
Reactor (cont.)
We can write a steady-state mole balance for component A, over the differential
element ΔW, as follows:
in – out + generation – consumption = accumulation
(1)
where FA is the molar flow rate of component A, and is the effectiveness factor
129
4. Design of Isothermal Packed-Bed Reactor
Taking the limit of ΔW → 0 re-arrang; we obtain the differential form of the
performance equation:
(2)
we can integrate Equation (2) to obtain the integral form of the performance
equation:
dFA dX A
= 0 − FAo (3)
dW dW
130
4. Design of Isothermal Packed-Bed Reactor
Substituting (3) into (2) gives:
(4a)
(4b)
Remarks:
▪ While catalyst deactivation is unavoidable, as we have seen in Chapter VI, there are
many applications in which it is very slow (say, 2-3 years), so that catalyst activity is
essentially constant
▪ Accordingly, PBR design calculations should take into account variation of pressure
along the reactor, as described by the Ergun equation. The design problem then
becomes more complicated and in most situations we need to combine Equation
(4a) or (4b) with the Ergun equation and solve a set of differential equations
numerically
131
H2 Production – Reactor Design
CHE-300
CHE-404
Consider the above catalytic reaction. Using the rate equation found for
this reaction determine the amount of catalyst needed in a packed bed
reactor(assume plug flow) for 35% conversion of A to R for a feed of
2000 mol/hr of pure A at 3.2 atm and 117°C.
Example - 2
4. Design of Isothermal Packed-Bed Reactor
Example-1:
143
4. Design of Isothermal Packed-Bed
Reactor (cont.)
We can write a steady-state mole balance for component A, over the differential
element ΔW, as follows:
in – out + generation – consumption = accumulation
(1)
where FA is the molar flow rate of component A, and is the effectiveness factor
144
4. Design of Isothermal Packed-Bed Reactor
Taking the limit of ΔW → 0 re-arrang; we obtain the differential form of the
performance equation:
(2)
145
4. Design of Isothermal Packed-Bed Reactor
Example-1:
149
4. Design of Isothermal Packed-Bed Reactor
Example-1:
150
4. Design of Isothermal Packed-Bed Reactor
Example-1:
151
4. Design of Isothermal Packed-Bed Reactor
Example-1:
H2 Production – SMR equilibrium
CHE-300
CHE-404
H2 Production – SMR Kinetics/rate
CHE-300
CHE-404