1. Introduction
New solutions for the efficient generation of carbon-neutral electricity are essential to curb the climate crisis. In the framework of renewable energies, it is well known that the energy supplied by wind or solar renewable power plants presents peaks and valleys due to the weather variability [
1]. In order to face up this variable generation and following the so-called “Power to gas” (P2G) concept, hydrogen can be produced by electrolysis of water during the generation peaks [
2,
3]. Green hydrogen is then stored and distributed for a wide range of end uses (hydrogen to power). Different cost-benefit studies have concluded that green hydrogen will be a feasible and competitive energy carrier in the coming years, optimising the use of natural resources [
4,
5,
6].
One of the alternatives relies on using the generated green hydrogen for electric generation in electric power plants, as it has already been done in several projects [
5,
7,
8]. The technical report of ENT Global [
9] summarises the recent advances in different gas turbine technologies and current capabilities of gas turbines using hydrogen as fuel. However, hydrogen is a very challenging fuel due to its combustion properties and differences with respect to conventional fuels such as natural gas, propane or kerosene. The current literature shows several investigations related to the behaviour of hydrogen flames focused on the flashback phenomena, as well as the thermal NO
formation [
10,
11,
12,
13].
However, from a broader view, it would be necessary to analyse the impact of replacing hydrogen or hydrogen–natural gas mixtures in current power plants. These types of studies could be carried out through thermal engineering principles and system analysis. Recently M. Pyo et al. [
14] developed a comparative study integrating the P2G process in a combined cycle power plant (CCPP). The considered scenarios focus on hydrogen produced using a P2G process and using it directly as fuel for the CCPP. In the second case, hydrogen was converted into synthetic natural gas by capturing carbon dioxide from the exhaust gases. Projecting the cost data for 2030, their economic feasibility analysis demonstrated that the direct use of hydrogen in the CCPP would be a promising option.
The strong impact of the COVID-19 pandemic in the energy system, together with the extraordinary increase of the price of natural gas in the last few months of 2021 as shown in
Figure 1, would further support the idea of using hydrogen or natural gas/hydrogen mixtures instead of pure natural gas [
15,
16]. Furthermore, the use of carbon-free fuels would eliminate the CO
taxes, also supporting the use of hydrogen for industrial-scale combustion systems and making it a more competitive alternative.
Additionally, combustion irreversibilities should be further analysed considering detailed calculations. In that sense, K. Nishida et al. [
17] compared the entropy generation and exergy losses for premixed and diffusion flames using hydrogen and methane as fuels, considering separately the different entropy generation mechanisms by heat conduction, chemical reaction, mass diffusion and viscous dissipation. Such type of calculations could be used to feed the thermodynamic analysis, combining CFD simulations with a system analysis approach.
In this context, the present study focuses on the use of hydrogen–methane blends as fuel in a regenerative Brayton cycle, aiming to qualitatively assess the global unsteady state such as the energetic and exergetic efficiencies of the cycle and the irreversibilities of each subsystem. For this purpose, a case study was defined using real operating parameters taken from the datasheet of the industrial natural gas turbine Siemens SGT-100 [
18]. The study was implemented using a thermodynamic system analysis approach, treating the components of the system as black boxes and solving the governing equations based on energy and entropy balances. It is worth noting that although this study was carried out considering the operating conditions of a particular industrial natural gas turbine, the present results and conclusions can be extended to any other type of gas turbine.
2. Objectives and Methodology
The main purpose of the present study was to analyse the impact of using different hydrogen/methane mixtures as a thermal source in a simple regenerative Brayton cycle, represented as a set of subsystems (compressor, combustion chamber, turbine and regenerator). A thermodynamic analysis was carried out to solve the energy and entropy balances in each subsystem, obtaining the global energy and exergy efficiencies in terms of the hydrogen content in natural gas, represented by
, which defines the volumetric percentage of H
in the fuel mixture varying from 0 to 1. The net power of the cycle, the pressure ratio in the compressor and the temperature of the combustion products and exhaust gases were fixed as design constraints on the basis of the industrial natural gas turbine Siemens SGT-100 datasheet [
18].
The chemical reactions in the combustion chamber were defined following the hypothesis of complete combustion by a single-step reaction between a hydrogen–methane mixture and air. The ideal gas assumption was used to specify the properties of the working fluids, considering air passing through the compressor and the regenerator, hydrogen and methane as fuels in the combustion chamber and the resulting mixture in the combustion products (present in the turbine and regenerator).
Once the I/O conditions were fixed at each subsystem, the equations were solved through the Engineering Equation Solver (EES) code [
19] for different values of
.
3. Problem Statement
In order to approximate the present calculations to real industrial gas turbine conditions, the datasheet of the industrial natural gas turbine Siemens SGT-100 [
18] was taken as reference, characterised by a net electrical power output of 5.1 MW, exhausts temperature of 818 K and a pressure ratio of 14:1 in the compressor. Those values were fixed as well as the inlet temperature in the gas turbine at T
= 1500 K. As stated by [
20], such temperatures may not exceed the 1700 K for reliable designs. The schematic view of the considered regenerative Brayton cycle is shown in
Figure 2 including the mentioned design values.
Air is drawn into the compressor at ambient conditions, the pressure is increased up to
with a pressure ratio of 14:1 and it is heated in the regenerator using the remaining energy in the exhaust gas mixture. The preheated air is then introduced into the combustion chamber together with the gaseous fuel, which is introduced at 298 K. The hydrogen fraction (
) was varied from 0 to 1 representing pure methane and pure hydrogen fuels, respectively. The combustion products expand in the gas turbine producing the net mechanical power of the cycle
= 5.6 MW (assuming an efficiency of 90% between the turbine and the electric generator [
21]), which is the difference between the power produced in the turbine and consumed by the compressor. The isentropic efficiencies in the compressor and turbine as well as the efficiency of the regenerator were fixed at 0.85 without loss of generality. Further information about the regenerative Brayton cycle and the T-s diagram of the process can be found in the literature [
22,
23].
Once the cycle is defined, the fact of having set the net power of the cycle, the temperature of the combustion products and the exhaust temperature results in different values for the main variables at each state for different values of , as shown next.
4. Governing Equations
Defining the enthalpy and entropies at each state highlighted in red in
Figure 2 (from st1 to st7), the energy and entropy balances were solved at each subsystem. States (st1-st2-st3) were defined considering air as an ideal gas. The state (st4) accounts for an ideal gas mixture of H
and CH
and states (st5-st6-st7) were defined by an ideal gas mixture of CO
, H
O, N
and O
(combustion products).
A single-step reaction between H
0–CH
(st4) and air (st3) was considered [
24]:
The fraction
and
represent the hydrogen volume fraction in the fuel mixture and the air-excess or air–fuel equivalence ratio, respectively. Hence, the relation between the air and fuel volumes (air–fuel ratio or AFR) can be defined as [
24]:
This value was used to relate the volumetric air and fuel fluxes through Equation (
3). It can be noted that the volumetric air flow
depends on the values of
and
f. As shown below, such relation had a significant impact on the final results.
Expressing the single-step reaction depending on
and
, it is possible to define the mole fractions of the species in terms of such parameters [
24]:
Thus, enthalpies and entropies at each state were defined using the mole fractions:
The EES code was used to calculate
and
for individual species
j in terms of the temperature for the former and pressure and temperature for the latter [
19,
22]:
where
and
represent the specific heat of species
j and the ideal gas constant, respectively. Here
represents the reference pressure (101,325 Pa). It is worth to mention that (kJ/kmol) and (kJ/kmol K) units were used for enthalpies and entropies, respectively. Once the properties at each state were defined, the governing equations for the global system and each subsystem were solved.
4.1. Global System
Taking into account the global system, the net power of the cycle (5.6 MW), the energetic and exergetic efficiencies were defined as:
where
and
represent the net thermal power and the net exergy flux introduced to the combustion chamber, respectively. Temperature and pressure for the reference environment were defined as
= 101,325 Pa and
= 298 K.
4.2. Subsystems
The governing equations for each subsystem are summarised below, which complete the system of equations to be solved:
Compressor:
Regenerator:
Combustion chamber:
Low heating value of the fuel (
):
Turbine:
Those equations define the thermodynamic behaviour of the regenerative Brayton cycle considering ideal gases and the combustion process in the chamber.
It must be noted that the specific chemical exergies for hydrogen and methane (
and
) were taken from [
22]. The equations system was solved through a parametric table for 20 values of
ranging from 0 to 1.
5. Results
Results from
Figure 3a demonstrate that both energetic and exergetic efficiencies increased while increasing
. Despite the small differences of approximately 2–4%, higher efficiencies using pure hydrogen instead of pure methane were addressed, which may be important for cost reduction in large-scale power plants for long periods. The observed higher exergetic efficiencies for an increasing hydrogen content means that the higher the hydrogen content, the greater the fuel performance in terms of work production.
The main reason for the trend in
Figure 3a can be explained by the exergy destruction rates from
Figure 3b. Exergy destruction rates at the compressor, regenerator and turbine remained almost constant. Nevertheless, exergy destruction in the combustion chamber, which was the main source of irrevesibilities in the cycle (representing approximately the 25% of the introduced chemical exergy), decreased for higher values of
.
Considering Equation (
29), this trend for
can be explained through the trends of the fuel, air and product entropies, which define the entropy destruction
in the combustion chamber. This equation can also be expressed in terms of exergy destruction:
Defining the entropy destruction per unit fuel mass,
Figure 4a demonstrates that the lower exergy destruction in the chamber is related to the lower fuel mass rate for increasing
, despite the higher entropy destruction per unit of fuel mass.
It is also worth to mention that the higher the hydrogen content in the fuel, the higher the expected product temperatures due to the higher adiabatic temperature of hydrogen. Thus, since
was fixed at 1500 K, air excess
was increased for higher
to maintain constant product temperatures as shown in
Figure 4b. Accordingly, it is interesting to note that exergy destruction still remained lower for higher
despite the higher values of
.
Figure 4b also demonstrates that the air mass flux was lower for higher
despite the higher values of
. This could be explained through the governing one-step reaction from Equation (
1). Considering pure methane (
= 0) and pure hydrogen (
= 1) at stoichiometric conditions (
= 1), one can conclude that hydrogen needs less air to be completely burnt compared to methane. This fact should be considered for the selection and sizing of the compressor.
The lower mass fuel rate for higher
observed in
Figure 4a is related to several parameters. Combining Equations (
2) and (
3) with Equation (
14), an expression for
can be derived in terms of the AFR:
This property and the trend of
led to a lower thermal power in the combustion chamber as depicted in
Figure 5a. Hence, the same mechanical power of the cycle (
= 5.6 MW) was achieved with a lower thermal power output in the combustion chamber and a lower fuel mass consumption for increasing values of
. It can be noted that the fuel consumption was reduced by almost 60% at
= 1. Since the fuel costs were fixed per unit mass, this results qualitatively demonstrated that hydrogen can still be competitive even though its unitary cost is more than twice that of natural gas.
The obtained results may be complemented with the so-called Wobbe index (WI) for different values of
. The WI is commonly used to measure the interchangeability of gaseous fuels and compares the thermal output of different fuels for a given combustion chamber or burner geometry. Two fuels with the same WI would have the same thermal power output for a constant chamber geometry, pressure and temperature conditions. It is expressed in terms of the high heating values (HHV) and densities of the considered fuels and the density of air [
25]:
As shown in
Figure 5b, for
> 0.6 the WI reaches a minimum at
0.85 and returns to the value it had at 0.6 when reaching
= 1. This fact indicates that the WI value is well bounded for
> 0.6, showing that the changes it undergoes are minimal compared to the decreasing trend from
= 0 up to 0.6. Thus, a good interchangeability can be expected at
> 0.6. However, for methane rich mixtures at
< 0.6 the WI increased, concluding that for a given combustion chamber or burner geometry the value of
should be precisely defined. As stated by Y. Zhao et al. [
25], the WI cannot be used as the only parameter to measure the interchangeability between hydrogen and methane due to the different combustion characteristics mainly related to the flammability limits, flame speeds and temperatures.
Taking into account the present results, a similar thermal output (
) could be expected for
> 0.6. However, the obtained
from
Figure 5a showed a decreasing trend since the products temperature was fixed at 1500 K and
varied with
to meet such a temperature target. Therefore, it can also be said that not only the WI or the combustion characteristics should be taken into account to assess the interchangeability of hydrogen and methane in a given combustion system, but also the overall behaviour of the thermodynamic cycle.
6. Conclusions
Considering a Brayton regenerative cycle corresponding to a particular industrial natural gas turbine, a qualitative thermodynamic analysis was performed varying the hydrogen fraction in methane () up to 100% H (vol %). A system analysis was carried out so as to assess the evolution of the global energetic and exergetic efficiencies as well as the impact on the cycle behaviour. The thermodynamic analysis treated the subsystems of the gas turbine as black boxes rather than going into the details of the individual components that make one turbine different from another. This fact makes it possible to generalise the conclusions of this study to other gas turbines with similar operating conditions.
Results demonstrated that the global exergetic efficiency increased by 5% from = 0 to 1 while the thermal efficiency increased by 2%. Such a difference may be related to the higher nonlinear behaviour of the irreversibilities (second law) compared to the differences related to the energy conservation (first law).
As expected, the combustion chamber presented the highest exergy destruction rates (3500–2700 kW), which were reduced for higher values of . Such a reduction was related to the lower fuel mass consumption rate and therefore the lower thermal power in the combustion chamber when varying from 0 to 1. This fact qualitatively showed that hydrogen can still be competitive even though its cost per unit mass is more than twice that of natural gas.
Although the present study served to draw initial qualitative results with interesting conclusions, a further analysis is needed, considering more detailed equations for the turbine and compressor based on turbomachinery theory and design, and considering the different compression and expansion stages, the velocity diagrams and stator and rotor losses, among others.