Ahmadreza Rahbari1, Mahinder Ramdin1, Leo J P van den Broeke1, Thijs J H Vlugt1. 1. Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands.
Abstract
Syngas is an important intermediate in the chemical process industry. It is used for the production of hydrocarbons, acetic acid, oxo-alcohols, and other chemicals. Depending on the target product and stoichiometry of the reaction, an optimum (molar) ratio between hydrogen and carbon monoxide (H2:CO) in the syngas is required. Different technologies are available to control the H2:CO molar ratio in the syngas. The combination of steam reforming of methane (SRM) and the water-gas shift (WGS) reaction is the most established approach for syngas production. In this work, to adjust the H2:CO ratio, we have considered formic acid (FA) as a source for both hydrogen and carbon monoxide. Using thermochemical equilibrium calculations, we show that the syngas composition can be controlled by cofeeding formic acid into the SRM process. The H2:CO molar ratio can be adjusted to a value between one and three by adjusting the concentration of FA in the reaction feed. At steam reforming conditions, typically above 900 K, FA can decompose to water and carbon monoxide and/or to hydrogen and carbon dioxide. Our results show that cofeeding FA into the SRM process can adjust the H2:CO molar ratio in a single step. This can potentially be an alternative to the WGS process.
Syngas is an important intermediate in the chemical process industry. It is used for the production of hydrocarbons, acetic acid, oxo-alcohols, and other chemicals. Depending on the target product and stoichiometry of the reaction, an optimum (molar) ratio between hydrogen and carbon monoxide (H2:CO) in the syngas is required. Different technologies are available to control the H2:CO molar ratio in the syngas. The combination of steam reforming of methane (SRM) and the water-gas shift (WGS) reaction is the most established approach for syngas production. In this work, to adjust the H2:CO ratio, we have considered formic acid (FA) as a source for both hydrogen and carbon monoxide. Using thermochemical equilibrium calculations, we show that the syngas composition can be controlled by cofeeding formic acid into the SRM process. The H2:CO molar ratio can be adjusted to a value between one and three by adjusting the concentration of FA in the reaction feed. At steam reforming conditions, typically above 900 K, FA can decompose to water and carbon monoxide and/or to hydrogen and carbon dioxide. Our results show that cofeeding FA into the SRM process can adjust the H2:CO molar ratio in a single step. This can potentially be an alternative to the WGS process.
One
of the consequences of the energy transition is that fossil fuel based
production of chemicals will be replaced with renewable energy based
processes.[1−3] The current infrastructure for producing chemicals
is predominantly based on hydrogen and carbon. This means that to
support the energy transition, a widely available and sustainable
C1 source is required. Therefore, the reuse of carbon dioxide
will be an essential part of future chemical production processes.[4−7] A range of efforts are underway to use carbon dioxide as a sustainable
and economical source of C1 to produce value-added chemicals.[5−8] There are basically two pathways for the conversion of carbon dioxide:
either by conventional hydrogenation or by electrochemical conversion.Formic acid (FA) is one of the simplest products that can be made
from carbon dioxide.[8] Recently, FA production
by electrochemical reduction of CO2 has gained significant
interest.[7,9−13] In this process, the overall reaction in the electrochemical
cell is the conversion of carbon dioxide with water to FA according
toThe main advantage of the electrochemical conversion
of carbon dioxide is that in the reaction water can be used as the
hydrogen source. The cathodic half-cell reduction of carbon dioxide
is described by the following reaction:[11]The formation of FA is a two electron reaction, and
the electric power to convert 1 kg of carbon dioxide to FA follows
from[14]where P is the power input in kWh per kg carbon dioxide, I (A) is the electric current, U is the
electrical potential which is on the order of 2.2–2.5 (V),
λ is the number of electrons, λ = 2, F is the Faraday coefficient which is equal to 96485 C molelectron–1, Q (C) is the total electric charge provided to
the reactor, t (s) is the time, and MCO (g mol–1) is the molecular
mass of carbon dioxide. For an overall energy efficiency, ϵ,
of around 70%, the energy required to convert 1 kg of carbon dioxide
into 1 kg of FA is ca. 4 kWh.[15] A simple
gross profit analysis using $5 per kg carbon dioxide and an electricity
price of $5 per kWh leads to a cost price of around $25 per kg FA.[16]For the hydrogenation reactions, the most
sustainable approach to produce the required hydrogen is by water
electrolysis, while traditional hydrogen production methods are based
on fossil fuels.[5,17] The conventional catalytic hydrogenation
of carbon dioxide to FA proceeds according to[4,18]FA is the simplest C1 carboxylic
acid, it is a nontoxic liquid between 281.55 and 373.15 K, and it
can be safely stored in aqueous solutions.[19,20] In addition, hydrogenation of biomass derived feedstocks has been
suggested as potential sustainable pathways to formate/formic acid
production.[6,21−26] Alternatively, value-added chemicals such as methanol, dimethyl
ether, and formate/formic acid can be produced by hydrogenation of
carbon dioxide.[4,18] To date, FA is mainly considered
as a hydrogen storage material via its decomposition to hydrogen and
carbon dioxide.[4,27−44] One of the key observations is that FA can be considered as a carbon
monoxide carrier as well via its decomposition to water and carbon
monoxide.[41,42] Basically, by combining the two main decomposition
pathways toward hydrogen and carbon monoxide, and additional products
such as water and carbon dioxide, FA can therefore be considered as
a source for syngas. Yoshida et al. have reported the presence of
FA as an intermediate in the water-gas shift reaction (WGS) reaction
reaction:[20,45,46]On a molecular weight basis, FA contains 4.3
wt % hydrogen and 60.9 wt % carbon monoxide. Using a FA density of
1.22 kg L–1 at standard conditions leads to 53 g
H2 per liter FA and 744 g of carbon monoxide per liter
FA. Based on the amount of 4.3 wt % or 53 g of hydrogen, FA is identified
as one of the most promising candidates for hydrogen storage.[38,47−49] Considering the high carbon monoxide fraction in
FA, it is interesting to explore the potential of FA as carbon monoxide
carrier.Typically, the WGS reaction is used together with steam
reforming of methane (SRM) to adjust the composition of the synthesis
gas (syngas). This is one of the most common and oldest methods for
syngas production.[5,50−56,58,59] The reaction pathways for the SRM and WGS areComparing reactions and R6 shows that by cofeeding FA to the SRM process, the WGS and the SRM
reactions can be performed in a single step.In this work, we
show that by using thermochemical equilibrium calculations, the syngas
composition (the H2:CO molar ratio) can be adjusted to
any required value between one and three by cofeeding FA to the SRM
reaction. FA in the reactant feed decomposes to water and carbon monoxide
and/or to hydrogen and carbon dioxide which are all involved in the
WGS reaction at high temperatures. This can potentially change the
conventional SRM and WGS reactions (R5 and R6) from a two-step process into a single-step process.This paper is organized as follows. In section , thermodynamic modeling of reactions and R5 is explained in detail. The Gibbs free energies of each component
is calculated at standard pressure and temperatures between 400 and
1400 K based on the partition function of isolated molecules. The
Gibbs minimization method is used to calculate the composition of
the product syngas at chemical equilibrium. Our results are summarized
in section . It is
shown that the H2:CO molar ratio can be adjusted to any
value between one to three based on the initial concentration of the
FA in the feed. Our conclusions are summarized in section .
Applications
of Formic Acid
Formic Acid Decomposition
The decomposition of FA can proceed according to two different
pathways: decarbonylation (or dehydration) into carbon monoxide and
water or decarboxylation (dehydrogenation) into hydrogen and carbon
dioxide:The selectivity
toward FA dehydration or dehydrogenation depends on the temperature,
pressure, and the type of catalyst. For the heterogeneous FA decomposition,
the dehydration/dehydrogenation selectivity of different solid catalysts
has been studied.[39,60] Metals and zinc oxide are predominantly
active for reaction , while other oxides are predominantly active for reaction .[60] Lopez et al. reported the results for different catalysts used for
the heterogeneous FA decomposition reactions in the temperature range
of T = 573–673 K.[39]Blake and Hinshelwood investigated the homogeneous decomposition
of FA acid in the gas phase for temperatures between T = 709 and 805 K and concluded that catalytic effects become negligible
at temperatures above T = 773 K.[41] Therefore, reactions and R8 are assumed to be in
equilibrium at high temperatures, which is a reasonable assumption
since kinetics are fast and of minor importance.[41] In the temperature range of T = 709 and
805 K, it was observed that reaction is of first-order while reaction is of second-order. The reaction rates for
packed and unpacked reactors were essentially the same for reactions and R8. In the beginning of the 1970s, Blake et al. extended
the experiment to the temperature range of T = 820–1053
K.[42] In this temperature range, reaction was also observed
to be a minor process, with typical CO:CO2 = 10:1 molar
ratios. Reaction is
of second-order for temperatures below T = 943 K
and has an order of 1.5 for higher temperatures. The difference in
yield of CO and CO2 was attributed to the water-gas shift
reaction.
Synthesis of Formic Acid
Current
industrial synthesis of FA is mainly based on fossil feedstocks using
methanol carbonylation/methyl formate hydrolysis and naphtha partial
oxidation.[47] On a large scale, FA is produced
in a two-step process of methanol carbonylation followed by methyl
formate hydrolysis. In 2014, this two-step process was used to produce
81% of FA acid worldwide.[61] In the first
step, carbon monoxide reacts with methanol at pressures around P = 4 MPa and temperatures around T = 353
K to produce methyl formate. FA and methanol are produced in the second
step by methyl formate hydrolysis. The produced methanol is recycled
back to the first step:[25,61]The sum of reactions and R10 reduces to the direct reaction of
carbon monoxide with water. FA synthesis based on methanol is a reliable
and established process; however, it has some drawbacks.[25,61] The process uses carbon monoxide produced at high temperatures which
is very energy intensive.[25] Also, a large
excess of water is required to decompose methyl formate to FA (R10).[25] The main application
of FA is for the production of preservatives and as antibacterial
agent;[62] it is also used for dyeing in
the leather industry. FA has received increased attention as a suitable
material for controlled hydrogen storage and release.[19,47,62−65] A relatively new application
is the use of FA in direct formic acid fuel cells (DFACF).[66,68,69] It has also been proposed to
use FA for storage and transportation of carbon monoxide[70] or carbon dioxide.[62,71]
Established Syngas Technologies
Syngas
refers to gas phase mixtures of hydrogen and carbon monoxide with
various H2:CO ratios.[5,52,53] Syngas can be produced by reforming almost any hydrocarbon source,
such as naphtha, heavy oil, natural gas, biomass, or coal.[52,56] Currently, steam reforming of light hydrocarbons (e.g., methane,
ethane, methanol, and ethanol) is the most commonly used method for
syngas production.[5,50−56,58,59] An alternative source for syngas production are coal reserves; however,
the investment costs associated with a coal-based syngas plant are
approximately 3 times higher as compared to a natural gas-based plant.[56] Therefore, natural gas remains the major source
for syngas production.[5,56] Syngas is an intermediate in
many industrial applications, and depending on the downstream process,
the optimal H2:CO molar ratio required in the syngas typically
lies between one and three.[51,72] The most common syngas
applications in the chemical process industry are methanol synthesis
(H2:CO = 2:1),[51,73] Fischer–Tropsch
(FT) synthesis (H2:CO = 2:1),[74−76] oxo-synthesis
or hydroformylation (H2:CO = 1:1),[77−81] and acetic acid synthesis (H2:CO = 1:1).[82] As an illustrative example, Figure shows different reaction pathways
leading to various syngas compositions by partial oxidation, steam
reforming, carbon dioxide reforming, and the combined FA and steam
option, as outlined in this work.
Figure 1
Different reaction pathways to reduce
methane to syngas using oxygen (R12), steam
(R5), carbon dioxide (R11), and an aqueous mixture of FA (R4) (proposed
in this work). Syntheses of different products require favorable syngas
H2:CO ratios.[145] Synthesis of
liquid hydrocarbons using the FT reaction (H2:CO = 2:1),[51] metal carbonyls, oxo-alcohols (H2:CO = 1:1),[81] acetic acid (H2:CO = 1:1),[82] methanol synthesis (H2:CO = 2:1),[73] and phosgene (H2:CO = 0:1).[146]
Different reaction pathways to reduce
methane to syngas using oxygen (R12), steam
(R5), carbon dioxide (R11), and an aqueous mixture of FA (R4) (proposed
in this work). Syntheses of different products require favorable syngas
H2:CO ratios.[145] Synthesis of
liquid hydrocarbons using the FT reaction (H2:CO = 2:1),[51] metal carbonyls, oxo-alcohols (H2:CO = 1:1),[81] acetic acid (H2:CO = 1:1),[82] methanol synthesis (H2:CO = 2:1),[73] and phosgene (H2:CO = 0:1).[146]To produce syngas from methane, various technologies have
been developed, such as SRM[5,58] and WGS,[5,58] carbon dioxide reforming of methane (CRM),[5,58] catalytic
partial oxidation of methane (POM),[74,75,81,83] combined partial oxidation
and carbon dioxide reforming of methane or autothermal reforming of
methane (ARM),[84,85] combined steam reforming, and
carbon dioxide reforming of methane (CSRCRM).[53,86]The first industrial SRM plant was commissioned in the early
1930s.[87,88] Methane is a very stable molecule, and only
at relatively high temperatures a high conversion rate to syngas is
obtained.[56,89] Syngas production from methane is divided
into two steps. First, at high temperatures ranging from T = 1073 to 1273 K and pressures ranging from P =
20 to 40 bar, the SRM reaction takes place. Second, the WGS is performed
after the SRM reaction to adjust the H2:CO molar ratio.[51,58,90] SRM is typically performed using
Ni-based catalysts.[5] This is related to
the low cost and favorable activity of the Ni-based catalysts as compared
to noble metals.[59,91] Although noble metals are more
coke resistant,[56] the high cost and the
limited availability make Ni catalysts a more practical choice in
commercial applications.[92] SRM has two
major drawbacks. In particular, the Ni-based processes suffer from
coke formation which leads to deactivation of the catalyst. To avoid
coke formation on the catalyst surface, excess steam is added which
results in H2 enriched syngas,[93] and this will lead to a syngas composition with a H2:CO
molar ratio larger than three.[54,55,84,94] The syngas compositions with
high H2:CO molar ratios do not meet the requirements for
many downstream petrochemical processes, e.g., FT synthesis,[56,74,75,88] acetic acid synthesis,[88] or methanol
synthesis.[72,94−97] The other disadvantage is that
the SRM reaction is highly endothermic and subsequently highly energy
intensive.[51,89,94,98,99]Alternatively,
in CRM (dry reforming), steam is replaced by carbon dioxide:[58]CRM is a CO2-consuming
reaction at temperatures between T = 1073 and 1273
K, resulting in syngas with H2:CO = 1:1 molar ratio.[5,56,100−102] This syngas composition is more suitable for some downstream processes.[56,74,75,100,101] To lower the H2:CO
molar ratio of the syngas, CRM is widely used as a secondary reforming
reaction after the SRM reaction.[103] CRM
synthesis using Ni-based catalysts, Co-based catalysts, and noble-metal-based
catalysts are reported extensively in the literature.[83,100,102] The main drawback of the CRM
reaction is the rapid coke deposition, which can be explained by the
Boudouard reaction[56,104] (2CO → C + CO2) taking place on the catalyst surface. Another disadvantage is the
direct decomposition of methane[56,104] (CH4 →
C + 2 H2) at high concentrations of CO2 in the
feed.[5,83,100,103,105]Catalytic partial
oxidation of methane (POM), also known as oxy-reforming, was introduced
as an alternative to obtain syngas with a H2:CO = 2:1 molar
ratio, suitable for producing long chain hydrocarbons,[5,74,75] and as a feed for methanol synthesis.[94−96]POM is favorable
for a wide range of temperatures allowing close to 100% methane conversion
to syngas.[89,106] The advantages include a short
residence time and mild exothermicity.[53,56,76,88,94,107] The main drawback is the presence
of hot spots as a result of the high conversion rates of methane.[53,76,92,93,108,109] Removing
the heat produced in the reactor is difficult for large-scale operations,
making the process difficult to control. From experiments by Prettre
et al. it was shown that the catalytic oxidation of methane, with
reactant feed composition CH4:O2 = 2:1, is not
accurately represented by (R12).[106,110] It seems that the POM reaction proceeds in two steps. The first
step is exothermic which involves deep oxidation (combustion) of a
part of the methane (approximately 25% of the starting moles) to carbon
dioxide and steam. All oxygen is consumed during this process. In
the second step, the residual methane reduces steam and carbon dioxide
to syngas.[89,94,106,107,110] This is an endothermic process. The POM reaction mechanism can be
described by the following three reactions:[106,110]The overall sequence of reactions
(R13–R15) using
a Ni/Al2O3 catalyst results in syngas with a
H2:CO = 2:1 ratio as reported by Dissanayake et al.[106] Yamamoto et al. have proposed the same reaction
mechanism for partial oxidation of C6+ hydrocarbons using supported Ni catalysts.[106,111] Different combinations of feedstock and catalysts can provide a
specific H2:CO molar ratio.[58,109]Autothermal
reforming of methane (ARM) is a combination of the POM and SRM-CRM
process.[112−115] ARM is performed either in one or two separate reactors to reduce
the energy consumption.[5,58] The combination of the exothermic
POM and endothermic SRM is energetically favorable.[85] ARM was originally designed for syngas production in ammonia
and methanol plants in the 1950s.[57] The
oxygen-steam flow is mixed with methane typically at around T = 2200 K,[56] and methane is
oxidized in a substoichiometric flame. Combustion products enter the
catalyst bed reactor with high thermal stability and with the temperature
in the range of T = 1200–1400 K:[56,57]Adding steam is crucial for the ARM process
as it prevents explosion hazards and suppresses coke formation.[5,85] Equilibration of the syngas is further governed by the SRM and WGS
reactions.[57] The H2:CO molar
ratio in the syngas can be precisely controlled by adjusting the H2O:CH4 and O2:CH4 molar ratios
in the feed.[56]Combined steam and
carbon reforming of methane (CSCRM) was proposed as an alternative
to directly control the syngas composition.[53,86] In this process, the H2:CO molar ratio is adjusted by
partially cofeeding carbon dioxide and steam with the reaction feed.
Adding steam to CRM process drastically reduces coke deposition on
the catalyst.[83,116] By changing the H2O:CO2:CH4 ratio in the reaction feed, a H2:CO ratio in the syngas between 1.5 and 2.5 is obtained.[53,57,76,83,93,94,116−118]
Modeling
and Methodology
For a single chemical reaction, the composition
of the reaction product at chemical equilibrium is calculated from
the method of equilibrium constants.[119−122] In this approach, mole fractions
are expressed as functions of a single variable called the reaction
coordinate (ϵ). The equilibrium constant is related to the individual
mole fractions of the components and the stoichiometric coefficients.
Therefore, ϵ is calculated for a single reaction.[119,122] The method of equilibrium constants becomes numerically more difficult
as the number of chemical species and reactions increases.[120,122,123]A necessary condition
for chemical equilibrium is that the total Gibbs energy of the mixture
reaches a minimum value at a given temperature and pressure. Based
on this principle, the Gibbs minimization method[119,122,124] is used as a robust method to
compute the composition of the reaction product at chemical equilibrium
for multicomponent systems with simultaneous reactions.[119,122,124,125] The solution obtained based on this method is less sensitive to
the initial guess as compared to other methods.[119,122,124] The composition of the reaction
product at chemical equilibrium is obtained by changing the initial
composition such that the Gibbs energy of the mixture is minimized.
The total number of atoms of each type should remain constant during
this minimization process. The Gibbs free energy, or the chemical
potential, of each component at the standard reference pressure, P° = 1 bar, can be evaluated from the isolated molecule
partition function:[126−129]with q(V,T)/Vo the temperature-dependent
part of the ideal gas partition function, kB is the Boltzmann constant, P° is the standard
reference pressure (1 bar), T is the temperature,
and the volume V° = kBT/P°. Details regarding the
calculation of the reference chemical potential from eq are provided in the Supporting Information.[129] The total Gibbs energy of a multicomponent mixture equals[122,130,131]where Gt is the total Gibbs energy of the mixture, S is
the number of components in the mixture, n is the number of moles of component i, μ is the chemical potential
of component i in the mixture, and S is the total number of components in the mixture. Considering the
standard state as an ideal gas, the chemical potential at any temperature
and pressure is obtained from[129,132]where R is the universal gas constant, y is the mole fraction of component i, and φ is the fugacity
coefficient of component i. The fugacity coefficient
can be obtained from experimental volumetric data or an equation of
state.[133] Combining eqs and 5 yieldsAt chemical equilibrium, the function Gt reaches a minimum. In a closed system, the minimization of eq is subject to the constraints
of the material balance.[119,122,124] In other words, the number of moles of each atom type remains constant
during the reaction. For k types of atoms in the
mixture, k independent mass balance equations are
applied as constraints:where A is the number of atoms of type k and α is the number
of atoms of type k present in molecule type i. Therefore, calculating the mixture composition at chemical
equilibrium is reduced to minimizing eq subject to the constraint of eq . The objective function, eq , is minimized using the function fmincon implemented in the MATLAB Optimization Toolbox.[134] In every iteration, the Peng–Robinson
equation of state (PR-EoS)[133,135,136] is used to evaluate the fugacity coefficients φ in eq . The mixture parameters are based on pure component parameters and
van der Waals mixing rules.[137,138] The effects of the
binary interaction parameters (BIPs) are negligible for gaseous mixtures
at high temperatures.[132] Therefore, the
BIPs are set to zero in this work. Further details of the PR-EoS modeling
and pure component parameters are provided in the Supporting Information. The standard Gibbs energies of reactions –R8 at P° are obtained based
on the computed chemical potentials of individual components, eq , and the corresponding
stoichiometric coefficients of the reactionwhere ν is the stoichiometric
coefficient of component i. The standard reaction
enthalpy ΔHr° is directly computed using the Gibbs–Helmholtz
equation:[126]
Results and Discussion
Main Reactions
In Figure , the values
for ΔGr° and ΔHr° (eqs and 9) are plotted as a function of the temperature. The data in Figure are obtained for reactions –R8. For more details on the computing of ΔGr° and ΔHr°, the reader is referred to the Supporting Information. The SRM reaction (R5) is endergonic, ΔGr° > 0, at
temperatures below T = 880 K,[89] and exergonic, ΔGr° < 0, at temperatures
above T = 880 K. This indicates that the syngas production
in the SRM reaction is favorable at high temperatures. The FA decomposition
reactions (R7 and R8)
are also endergonic for the temperature range of T = 400–1400 K. Therefore, thermodynamic equilibrium favors
high conversion of FA to water, hydrogen, carbon dioxide, and carbon
monoxide at high temperatures.[41,42] The WGS reaction is
endergonic at temperatures above T = 1100 K. At high
enough temperatures, higher conversion of carbon dioxide and hydrogen
to carbon monoxide and water is favored.[139,140] The reaction enthalpies are calculated directly from the Gibbs–Helmholtz
equation (eq ). From
the reaction enthalpies, ΔH, it is clear that reactions and R7 are endothermic and reactions and R8 are exothermic.
Figure 2
(a) Standard
Gibbs energies of reaction and (b) reaction enthalpies for reactions and R8 (per mole of FA), reaction (per mole of methane), and reaction (per mole of water) as a
function of temperature at P° = 1 bar. The equilibrium
constant is related to the Gibbs free energy change of the reaction.[122,126] The symbols indicate SRM (downward-pointing triangles), WGS (circles),
dehydration of FA (squares), and dehydrogenation of FA (upward-pointing
triangles). A dashed line is used as a reference line at zero. Standard
Gibbs energies of carbon monoxide, water, carbon dioxide, hydrogen,
and FA are provided in Table S1 of the Supporting Information.
(a) Standard
Gibbs energies of reaction and (b) reaction enthalpies for reactions and R8 (per mole of FA), reaction (per mole of methane), and reaction (per mole of water) as a
function of temperature at P° = 1 bar. The equilibrium
constant is related to the Gibbs free energy change of the reaction.[122,126] The symbols indicate SRM (downward-pointing triangles), WGS (circles),
dehydration of FA (squares), and dehydrogenation of FA (upward-pointing
triangles). A dashed line is used as a reference line at zero. Standard
Gibbs energies of carbon monoxide, water, carbon dioxide, hydrogen,
and FA are provided in Table S1 of the Supporting Information.The Gibbs minimization
method is used to compute the syngas equilibrium composition for the
SRM and WGS reactions (R5 and R6). The reaction is studied with an equimolar feed mixture
of water and methane, H2O:CH4 = 1:1, in the
temperature range of T = 800–1500 K at P = 1 and 25 bar. The results are shown in Figure . As expected, the H2:CO molar ratios in the syngas are larger than three for the two
pressures. It follows from Figure that full conversion of methane is achieved at T = 1200 K at P = 1 bar, while nearly full
conversion of methane at P = 25 bar is not achieved
until temperatures above T = 1500 K. For both pressures,
low concentrations of carbon dioxide are observed in the syngas mixture
at high temperatures. This is because the WGS equilibrium shifts toward
carbon monoxide and water at high temperatures.[51,139−141]
Figure 3
Equilibrium composition of syngas as a function
of temperature computed using the Gibbs minimization method (reactions and R6): (a) at pressure of 1 bar and H2O:CH4 = 1:1 and (b) at pressure of 25 bar and H2O:CH4 = 1:1. In both panels: mole fractions of hydrogen (squares),
mole fractions of carbon monoxide (circles), mole fractions of methane
(downward-pointing triangles), mole fractions of water (diamonds),
and mole fractions of carbon dioxide (upward-pointing triangles).
Equilibrium composition of syngas as a function
of temperature computed using the Gibbs minimization method (reactions and R6): (a) at pressure of 1 bar and H2O:CH4 = 1:1 and (b) at pressure of 25 bar and H2O:CH4 = 1:1. In both panels: mole fractions of hydrogen (squares),
mole fractions of carbon monoxide (circles), mole fractions of methane
(downward-pointing triangles), mole fractions of water (diamonds),
and mole fractions of carbon dioxide (upward-pointing triangles).
FA Combined
with the SRM Process: SRM-FA
To reduce the carbon footprint
of hydrogen and syngas production, alternative process schemes need
to be developed. In Figure , we propose a process scheme in which FA is combined with
the SMR process to provide a wide range of H2 to CO ratios.
In this way both the CH4:H2O and the HCOOH:H2O molar ratios can be varied. By using essentially both CH4 and CO2 as the C1 feedstock, the overall
consumption of methane will be reduced.
Figure 4
Comparison between a
typical SRM layout and the layout for the proposed combined SRM-FA
process. In the exiting SRM process, steam reforming is followed by
the WGS process to adjust the H2:CO ratio. In the alternative
process, first FA is synthesized, and second the FA is added to the
SMR to adjust the H2:CO ratio. FA can be synthesized either
by electrochemical conversion of CO2[9−13] or by conventional catalytic hydrogenation of CO2.[4,18]
Comparison between a
typical SRM layout and the layout for the proposed combined SRM-FA
process. In the exiting SRM process, steam reforming is followed by
the WGS process to adjust the H2:CO ratio. In the alternative
process, first FA is synthesized, and second the FA is added to the
SMR to adjust the H2:CO ratio. FA can be synthesized either
by electrochemical conversion of CO2[9−13] or by conventional catalytic hydrogenation of CO2.[4,18]For existing hydrogen and syngas production processes, there
are two sources of carbon dioxide. To obtain the required product
specifications for the hydrogen or the syngas, pressurized carbon
dioxide is removed from the SRM and the WGS processes. Additionally,
carbon dioxide is produced during heat generation and is present in
the flue gas stream. The pressurized carbon dioxide stream from the
existing hydrogen or syngas production units can be used as feedstock
for the synthesis of FA, for both the electrochemical conversion and
the hydrogenation of carbon dioxide. It should be noted that large
scale conversion of carbon dioxide to FA is not yet available. The
advantage of the electrochemical route is that the product will be
an aqueous FA stream. Various aqueous FA solutions, with different
FA wt %, can be fed to the SRM-FA process, where the final syngas
composition can be adjusted by the operating conditions for the pressure
and temperature.From the SRM process, syngas with a molar ratio
of H2:CO = 3:1 is generally obtained. However, for most
applications a lower H2:CO molar ratio is required (see Figure ). To assess the
potential of FA as a carbon monoxide carrier, the thermodynamic equilibrium
of combining the FA decomposition reactions and the SRM reaction was
evaluated. The composition of the feed mixture was defined by the
molar ratio between water and methane, H2O:CH4, and varying the molar ratio between FA and water, HCOOH:H2O. Two cases for the H2O:CH4 molar ratio are
considered: H2O:CH4 = 1:1 and H2O:CH4 = 2:1. For the FA, a HCOOH:H2O molar ratio in
the range from 0.49 to 5.66 has been used (see Table ). The equilibrium composition of the syngas
is calculated using the Gibbs minimization method based on reactions and R5.
Table 1
Different Molar Ratios
of FA in FA–Water Mixtures Used in the Reactant Feeda
FA:H2O
FA:(FA + H2O) (%)
FA (wt %)
0.11
10
22
0.49
33
56
1.00
50
72
1.50
60
79
5.66
85
94
The corresponding mole percentage and weight percentage
of FA (wt %) in the mixture is calculated based on the molar ratio
between FA and water. The molar ratios between water and methane used
in the simulations are H2O:CH4 = 1:1 and H2O:CH4 = 2:1.
The corresponding mole percentage and weight percentage
of FA (wt %) in the mixture is calculated based on the molar ratio
between FA and water. The molar ratios between water and methane used
in the simulations are H2O:CH4 = 1:1 and H2O:CH4 = 2:1.The results for the equilibrium syngas composition for the temperature
range of T = 900–1500 K at P = 1 bar are shown in Figure , and the results for P = 25 bar are shown
in Figure . At P = 1 bar, full conversion of methane is achieved at temperatures
up to T = 1100 K. By increasing the temperature further,
the equilibrium favors conversion of hydrogen and carbon dioxide to
water and carbon monoxide. This is in agreement with the equilibrium
of the WGS reaction at high temperatures.[36,140,142] In addition, thermodynamic equilibrium
favors complete FA decomposition (R4) in this
temperature range. This leads to an increase in the mole fractions
of water, carbon dioxide, and carbon monoxide compared to the SRM-WGS
process. Because the mole fraction of hydrogen is decreasing with
the increase in temperature, contrary to the mole fraction of carbon
monoxide, different H2:CO molar ratios are obtained at
different temperatures.
Figure 5
Equilibrium composition of syngas as a function
of temperature obtained by cofeeding FA to the SRM reaction at 1 bar
and H2O:CH4 = 1. The Gibbs minimization method
is used to obtain the syngas equilibrium composition using eqs and R5. Initial mole fraction of FA relative to the mole fraction
of water: (a) HCOOH:H2O = 0.49, (b) HCOOH:H2O = 1.00, (c) HCOOH:H2O = 1.50, and (d) HCOOH:H2O = 5.66. In all panels: mole fractions of hydrogen (squares), mole
fractions of carbon monoxide (circles), mole fractions of methane
(downward-pointing triangles), mole fractions of water (diamonds),
mole fractions of carbon dioxide (upward-pointing triangles), and
mole fractions of FA (crosses).
Figure 6
Equilibrium composition of syngas as a function of temperature obtained
by cofeeding FA to the SRM reaction at 25 bar and H2O:CH4 = 1. The Gibbs minimization method is used to obtain the
syngas equilibrium composition based on reactions and R5. Initial mole
fraction of FA relative to reactions mole fraction of water: (a) HCOOH:H2O = 0.49, (b) HCOOH:H2O = 1.00, (c) HCOOH:H2O = 1.50, and (d) HCOOH:H2O = 5.66. In all panels:
mole fractions of hydrogen (squares), mole fractions of carbon monoxide
(circles), mole fractions of methane (downward-pointing triangles),
mole fractions of water (diamonds), mole fractions of carbon dioxide
(upward-pointing triangles), and mole fractions of FA (crosses).
Equilibrium composition of syngas as a function
of temperature obtained by cofeeding FA to the SRM reaction at 1 bar
and H2O:CH4 = 1. The Gibbs minimization method
is used to obtain the syngas equilibrium composition using eqs and R5. Initial mole fraction of FA relative to the mole fraction
of water: (a) HCOOH:H2O = 0.49, (b) HCOOH:H2O = 1.00, (c) HCOOH:H2O = 1.50, and (d) HCOOH:H2O = 5.66. In all panels: mole fractions of hydrogen (squares), mole
fractions of carbon monoxide (circles), mole fractions of methane
(downward-pointing triangles), mole fractions of water (diamonds),
mole fractions of carbon dioxide (upward-pointing triangles), and
mole fractions of FA (crosses).Equilibrium composition of syngas as a function of temperature obtained
by cofeeding FA to the SRM reaction at 25 bar and H2O:CH4 = 1. The Gibbs minimization method is used to obtain the
syngas equilibrium composition based on reactions and R5. Initial mole
fraction of FA relative to reactions mole fraction of water: (a) HCOOH:H2O = 0.49, (b) HCOOH:H2O = 1.00, (c) HCOOH:H2O = 1.50, and (d) HCOOH:H2O = 5.66. In all panels:
mole fractions of hydrogen (squares), mole fractions of carbon monoxide
(circles), mole fractions of methane (downward-pointing triangles),
mole fractions of water (diamonds), mole fractions of carbon dioxide
(upward-pointing triangles), and mole fractions of FA (crosses).Carrying out the SRM-FA process
at P = 25 bar changes the equilibrium composition
of the reacting system, such that higher temperatures are required
to fully reform methane and to reduce the carbon dioxide content in
the syngas. This is in agreement with the Le Chatelier’s principle[143,144] which states that an increase in the pressure leads to a change
in equilibrium composition to a new state in which fewer molecules
per mole are present. Here, the thermodynamic equilibrium is shifted
toward water, carbon dioxide, and methane (R5) at low temperatures. Therefore, higher temperatures are required
to reduce the methane and carbon dioxide concentrations in the syngas.On the basis of the results shown in Figures and 6, it is clear
that the concentrations of hydrogen and carbon monoxide can be adjusted
by changing the FA concentration in the reactant feed. To have a clear
overview of this principle in Figure , the H2:CO molar ratios in the syngas are
plotted as a function of the composition of the reactant feed in the
temperature range of T = 800–1500 K at P = 1 and 25 bar. The composition of the reactant feed was
obtained by adjusting the HCOOH:H2O molar ratios between
0.11 and 5.66. Results shown in Figures a and 7b correspond
to H2O:CH4 = 1:1 molar ratio in the reactant
feed at P = 1 and 25 bar, respectively. The results
shown in Figures c
and 7d correspond to H2O:CH4 = 2:1 molar ratio in the reactant feed at P = 1 and 25 bar, respectively.
Figure 7
Different equilibrium syngas compositions
(H2:CO ratios) are obtained by cofeeding FA to the SRM
reaction at different pressures and temperatures. (a) P = 1 bar, H2O:CH4 = 1; (b) P = 25 bar, H2O:CH4 = 1; (c) P = 1 bar, H2O:CH4 = 2; and (d) P = 25 bar, H2O:CH4 = 2. In all panels, the
initial mole fraction of FA relative to the mole fraction of water:
HCOOH:H2O = 0.11 (downward-pointing triangles), HCOOH:H2O = 0.49 (diamonds), HCOOH:H2O = 1.00 (circles),
HCOOH:H2O = 1.50 (squares), and HCOOH:H2O =
5.66 (upward-pointing triangles).
Different equilibrium syngas compositions
(H2:CO ratios) are obtained by cofeeding FA to the SRM
reaction at different pressures and temperatures. (a) P = 1 bar, H2O:CH4 = 1; (b) P = 25 bar, H2O:CH4 = 1; (c) P = 1 bar, H2O:CH4 = 2; and (d) P = 25 bar, H2O:CH4 = 2. In all panels, the
initial mole fraction of FA relative to the mole fraction of water:
HCOOH:H2O = 0.11 (downward-pointing triangles), HCOOH:H2O = 0.49 (diamonds), HCOOH:H2O = 1.00 (circles),
HCOOH:H2O = 1.50 (squares), and HCOOH:H2O =
5.66 (upward-pointing triangles).Thermochemical equilibrium calculations clearly show that reactions and R5 can be combined to produce syngas with an adjustable
H2:CO molar ratio ranging from one to three. The H2:CO molar ratio can be adjusted by changing the HCOOH:H2O:CH4 ratio in the reactant feed at different temperatures.
At high pressures, higher temperatures are required to reduce the
concentration of methane and carbon dioxide in the product syngas,
as shown in Figure . However, adjusting the H2:CO molar ratio in the syngas
can be achieved at any temperature and pressure.The results
show that by feeding FA to the SRM process, the equilibrium composition
of the product syngas can be adjusted by changing the concentration
of FA in the reactant feed. Future studies should investigate the
effect of different types of catalyst for the combined SRM-FA process
at different temperatures. The proposed method for adjusting the H2:CO ratio by using FA is not limited to the methane steam
reforming process. First, it can be used in any process where adjustment
of the H2:CO ratio is required (see Figure ). Examples of this include autothermal reforming,
partial oxidation, gas-to-liquid technologies, naphtha reforming,
biomass gasification, etc. Second, formic acid can be used for energy
storage by the use of fuel cells and formic acid reformers to generate
hydrogen, heat, and electricity.
Conclusions
To adjust the H2:CO molar ratio during syngas production,
FA decomposition can be combined with the steam reforming of methane.
The option to use FA as a syngas source is exploited by combining
the two FA decomposition reactions at high temperatures. Essentially,
FA can be considered as a combined hydrogen and carbon monoxide carrier.
Thermodynamic equilibrium calculations show that the syngas composition
can be controlled by adjusting the HCOOH:H2O:CH4 ratio in the reactant feed. It is possible to obtain different H2:CO molar ratios between 1 and 3 in the product syngas. At
higher pressures, higher temperatures are required for complete methane
conversion and reducing carbon dioxide content in the syngas. On the
basis of our results, it can be concluded that cofeeding FA to the
SRM reaction can potentially reduce the traditional SRM and WGS processes
from a two-step process to a single-step process able to produce syngas
with adjustable H2:CO ratio. The proposed SMR-FA process
based on CO2 reuse may open up a range of new applications
for formic acid.
Authors: Karaked Tedsree; Tong Li; Simon Jones; Chun Wong Aaron Chan; Kai Man Kerry Yu; Paul A J Bagot; Emmanuelle A Marquis; George D W Smith; Shik Chi Edman Tsang Journal: Nat Nanotechnol Date: 2011-04-10 Impact factor: 39.213
Authors: Antero T Laitinen; Vyomesh M Parsana; Olli Jauhiainen; Marco Huotari; Leo J P van den Broeke; Wiebren de Jong; Thijs J H Vlugt; Mahinder Ramdin Journal: Ind Eng Chem Res Date: 2021-04-07 Impact factor: 3.720