Abstract

There is significant interest in utilizing hydrogen as an energy carrier in a net-zero carbon energy economy. A variety of fundamental and practical issues associated with premixed hydrogen/methane combustion must be addressed, particularly flame flashback, combustion instabilities, and NOx emissions. Significantly less attention has been given to CO emissions from hydrogen blending. However, recent H2 blending engine tests have clearly noted major reductions in CO emissions, reductions which substantially exceed what would be expected based upon C atom mass balances in the fuel. Moreover, NOx, CO, and combustor operability always tradeoff with each other, and so, for example, these CO reductions can also enable operational or design modifications to reduce NOx emissions and/or wider operability windows. This paper addresses H2 blending effects on CO emissions, addressing both the direct influences on CO emissions, as well as how it influences NOx–CO tradeoffs. The paper presents computations that differentiate between the kinetic and equilibrium effects of H2 blending, as well as how H2 blending can lead to reductions in both NOx and CO by moving the Pareto front downward. Finally, this work suggests that H2 blending, in addition to its CO2 reductions, can also increase plant operational flexibility and so provide additional value propositions to an evolving electricity market.

1 Introduction

Hydrogen is being evaluated as a zero-carbon energy carrier and energy storage medium for a variety of industrial, transportation, and electric power applications [1]. A variety of demonstrations and pilots are being explored, generally with variable levels of blending with natural gas [2,3]. Natural gas, which is composed primarily of methane, is a dominant source of electrical power, as well as for heat in other commercial and industrial applications. Key emissions of concern from combustion are usually carbon monoxide (CO), oxides of nitrogen (NOx), unburned hydrocarbons, and particulate matter (PM).

A number of recent studies have addressed pollutant emissions from CH4/H2 blends. Because nonpremixed H2-blended flames are hotter than nonpremixed CH4 blends, it is critical to carefully specify what is being held constant when such a comparison is being made. In particular, thermal NO production is a linear function of residence time and O2 concentration, and an exponential function of temperature. For this reason, a number of early studies of CH4/H2 blends reported significant increases in NO emissions [2,47] with H2 blending, likely due to elevated flame temperatures leading to higher thermal NO formation rates. However, when comparing CH4/H2 blends at constant temperature, recent work from our group has shown that roughly constant NO levels are produced, with deviations due to product O atom concentrations and fuel-dependent roles of flame NO production mechanisms [8].

NOx and CO tradeoffs play one of the dominant roles in how modern combustors are designed, as well as their operational limits. For example, CO exists at superequilibrium levels in the flame, and the rate at which CO relaxes to equilibrium increases with temperature and pressure. Only if CO levels are at equilibrium does increasing temperature lead to increased CO (due to CO2 dissociating back to CO). In general, however, CO emissions are reduced with increasing pressure, temperature, and residence time. NO emissions generally exhibit the opposite tendencies. Longer residence times lead to higher NO formation and higher temperatures/pressures to elevated formation rates. As such, combustor design and operational limits are generally associated with balancing NO at high power/temperature and CO emissions at low power/temperature. Indeed, the fact that CO emissions increase rapidly with decreases in flame temperature has huge ramifications on the global gas turbine fleet operational characteristics, as these engines can only be turned down to a minimum power level (often on the order of 60% max power) without going out of compliance on their CO permit. Typical regulations (such as the U.S. Environmental Protection Agency, EPA) are defined with respect to output-based standards (e.g., the mass of NOx/CO per unit of useful shaft work), which for a given fuel can be converted to the more easily measured concentration-based standard (in units of parts per million by volume at 15% oxygen) [9]. The relationship between normalized mass output standards and concentration is itself fuel dependent, and so output-based standards must be used when comparing emissions across different fuels, which becomes particularly important for H2-blended fuels [10].

Even while both NOx and CO have such important influences on how combustors are both designed and operated, their actual health and environmental impacts are quite different. For reference, the electric power sector contributed about 1% of total U.S. CO emissions in 2023 [11]. Moreover, CO is not a driver of new regulations because the entire U.S. is currently in attainment of the CO National Ambient Air Quality Standard (NAAQS). NO2 is a regulated criteria pollutant because of its potential to cause adverse respiratory health effects and because it contributes to acid deposition. Though all areas in the U.S. are currently meeting the air quality standards for NO2, it is still a concern because it is a precursor to two other criteria pollutants, ozone, and fine particulate matter (PM2.5). As of Mar. 31, 2024, there are 46 areas in the U.S. with 115 million people which fail to meet the 2015 NAAQS for ozone [12]. Additionally, there are five areas in the U.S. with 21 million people in nonattainment of the 2012 PM2.5 NAAQS and 11 areas with 31 million people in nonattainment of the 2006 PM2.5 NAAQS. On Feb. 7, 2024, EPA issued a new PM2.5 NAAQS, and nonattainment area designations are pending for the new annual standard of 9 μg m3.

As might be expected, H2 blending should reduce CO emissions from flames due to simple mass balance considerations, as there are fewer carbon atoms in the fuel. To illustrate, Fig. 1 quantifies the reduction by showing a calculation of the equilibrium CO emissions produced from H2/CH4 blends at constant adiabatic flame temperature. The nonzero CO emissions at 100%H2 come from atmospheric CO2.

Fig. 1
Equilibrium CO production from a hydrogen–methane blend flame at “high-power” combustor inlet condition (i.e., 20 bar, 800 K, see also Table 1) for fixed adiabatic flame temperatures. Note: the secondary y-axis shows ppmvd corrected to 15%O2 based on a 0%H2 in blend.
Fig. 1
Equilibrium CO production from a hydrogen–methane blend flame at “high-power” combustor inlet condition (i.e., 20 bar, 800 K, see also Table 1) for fixed adiabatic flame temperatures. Note: the secondary y-axis shows ppmvd corrected to 15%O2 based on a 0%H2 in blend.
Close modal

However, field data from H2/CH4 blending testing show dramatically larger, even orders of magnitude larger, CO reductions than illustrated by Fig. 1 [13]. At a field test of a NYPA Brentwood GE LM6000 test, there was a 70% reduction in CO as hydrogen increased from 0% to 35% hydrogen by volume blend, which corresponds to an approximately 20% decrease in CO2 emissions. This was measured while NOx was held constant through an increase in water to fuel ratio, which tends to increase CO emissions. At the Southern Company Mitsubishi 501G testing, baseload CO emissions remained below the analyzer detection limit at the high baseload firing temperature from 0% to 20% blending volume. At 50% load, controlled to a lower firing temperature, CO emissions were reduced between 85% and 98% for an increase of hydrogen from 0% to 20% by volume, which corresponds to an approximately 7% decrease in CO2 emissions [14,15].

These observations are a key motivator behind this paper. These CO reductions could be due to two additional factors, beyond the obvious reduction in available carbon atoms with increased H2 blending. First, H2 blending can accelerate CO oxidation to CO2. For example, the important reaction CO + OH → CO2 + H would certainly be impacted by altered levels of both OH and H in flames with H2 blending. Second, H2 blending is well known to increase turbulent flame speeds, thus making the flame shorter and thereby increasing postflame residence times for CO oxidation.

The purpose of this paper is to further explore and quantify this first effect, namely, the kinetic influences of H2CH4 blending on CO emissions. In addition to their fundamental interest, this relationship will have important impacts on both the operation of existing combustors and the design of future ones. For example, a significant number of gas turbines in the U.S. are operated intermittently, shutting down and turning back on once or even twice a day, as their operational window is limited by low-load CO emissions. H2 blending thereby could expand operational envelopes and increase their utilization factors. Second, new dedicated designs for elevated H2 levels could operate with shorter residence times than current systems, thereby enabling reductions in NO emissions for a given CO emission. This paper presents the results of computations to provide further insight into these questions and to understand the role of H2 on CO relaxation rates, with attention given to fundamental influences of pressure, temperature, residence time, and H2 blending levels.

2 Numerical Methods

2.1 Test Matrix.

A test matrix was developed to individually study the effects of fuel composition, adiabatic flame temperature, and pressure conditions. The varying pressure cases were run with different preheat temperature values to simulate different applications and/or power levels of a system. The test matrix was developed by varying equivalence ratios for three power conditions to achieve constant product temperatures across fuel blends. Adiabatic flame temperatures of 1600–2000 K were analyzed across three initial power conditions: “low” power (1 bar, 300 K), “medium” power (5 bar, 650 K), and “high” power (20 bar, 800 K). The low and high-power conditions were chosen to correspond to “idle” and “full power” of an F-class industrial gas turbine, respectively. Eight hydrogen blends from 0% to 100% by volume were studied. Also, note that adiabatic flame temperatures were held constant to ±6 K of the quoted values. The full test matrix, where the value in the table indicates the equivalence ratio, is shown in Table 1. In order to capture CO emissions at the 100%H2 cases, the simulated air includes 412 ppm of atmospheric CO2 [16].

Table 1

Equivalence ratios and fuel blends utilized for the constant adiabatic flame temperature results (using HyChem product species) used in this study

%H2 fuel blend
Tad(K)0255075909599100
1 bar, 300 K16000.5610.5530.5420.5250.5070.4990.4890.487
17000.6160.6080.5960.5760.5590.5490.5380.535
18000.6720.6650.6530.6320.610.6010.590.586
19000.7320.7250.7120.690.6680.6560.6440.641
20000.7930.7880.7740.7510.7270.7150.7020.698
5 bar, 650 K16000.4130.4070.3990.3840.3680.3610.3550.352
17000.4650.4590.4490.4340.4160.4080.40.397
18000.5190.5120.5020.4840.4650.4560.4470.444
19000.5750.5680.5570.5370.5160.5060.4960.493
20000.6330.6260.6140.5930.570.5590.5470.543
20 bar, 800 K16000.3460.3420.3330.3210.3080.3010.2950.292
17000.3960.3920.3820.3680.3530.3460.3380.336
18000.450.4430.4340.4180.40.3940.3870.384
19000.5030.4960.4870.4690.4510.4440.4350.431
20000.5610.5520.5420.5230.50720.4940.4840.485
%H2 fuel blend
Tad(K)0255075909599100
1 bar, 300 K16000.5610.5530.5420.5250.5070.4990.4890.487
17000.6160.6080.5960.5760.5590.5490.5380.535
18000.6720.6650.6530.6320.610.6010.590.586
19000.7320.7250.7120.690.6680.6560.6440.641
20000.7930.7880.7740.7510.7270.7150.7020.698
5 bar, 650 K16000.4130.4070.3990.3840.3680.3610.3550.352
17000.4650.4590.4490.4340.4160.4080.40.397
18000.5190.5120.5020.4840.4650.4560.4470.444
19000.5750.5680.5570.5370.5160.5060.4960.493
20000.6330.6260.6140.5930.570.5590.5470.543
20 bar, 800 K16000.3460.3420.3330.3210.3080.3010.2950.292
17000.3960.3920.3820.3680.3530.3460.3380.336
18000.450.4430.4340.4180.40.3940.3870.384
19000.5030.4960.4870.4690.4510.4440.4350.431
20000.5610.5520.5420.5230.50720.4940.4840.485

2.2 Kinetics Calculations.

In this study, a kinetic model of a laminar, one-dimensional premixed flame was set up in ansyschemkin using the premix package [17]. The calculations were performed over sufficient axial distance to capture in-flame and postflame processes. The grid independence study performed by Breer et al. [8] was used to extend the same numerical results across the analysis of CO emissions. The mechanism implemented in this study is a HyChem kinetic model [18]. The HyChem model was chosen following the comparison conducted by Breer et al. [8] between HyChem, GRI 3.0 [19], Glarborg [20], and UCSD [21] models. The test matrix shown in Table 1 was replicated in chemkin using a parametrized sweep across the fifteen cases corresponding to five cases at each power condition. Grid independence studies were performed to confirm appropriate gridding across the cases.

2.3 Data Processing.

Results are analyzed as a function of residence time, denoted as τres, where the reference point τres=0 is defined as the location at which the value of the second derivative of OH concentration with respect to position (d2[OH]/dx2) normalized by its maximum exceeds a small threshold of 104. Further details on considerations of this reference frame are provided in Ref. [8], but the trends are insensitive to this value except for very low residence time values not considered in this paper, where one lies within the exothermic reaction zones. The residence time was calculated as shown in Eq. (1). The equation was integrated using a trapezoidal cumulative summation method [8]. In Eq. (1), “x” is the position and “V” is the axial flow velocity, and the position was found from the positional grid computed using chemkin
(1)

Here, the emissions are quantified using two different metrics—mass of emissions produced per energy released by fuel using lower heating value of the fuel, and in terms of parts-per-million by volume, dry (ppmvd) corrected to 15%O2. Because the conversion between these two quantities is H2 concentration dependent [10], the ppm values were only calculated for the 0%H2 case (if not, two lines would be needed for each hydrogen blending value). This is done in order to plot the data in both units simultaneously and provide a benchmark comparison at 0%H2. The purpose of including these concentration values on many of the secondary y-axes is that many emissions regulations are written in terms of concentrations, and many analyzers report measured values using these concentration units. However, it is emphasized that only the mass normalized metric should be used when comparing CO production values across different H2 blends as the ppm values presented are only relevant to the 0%H2 case.

As noted earlier, H2 blending should naturally lead to decreased CO emissions from pure C atom mass balance considerations, in addition to the kinetic effect of accelerating oxidation of reactants to major combustion products. To differentiate these effects and understand the kinetic effect of CO as a function of H2CH4 blend ratio, the kinetic factor COk is defined in the following equation:
(2)

This factor weighs the CO concentration by the ratio of CO2 concentrations for the CH4/H2 blended mixture and pure CH4. As such, variations in COk levels with H2 blending can be interpreted as a kinetic effect of H2 blending. If this ratio does not change with H2 blending, then variations in CO emissions with H2 blending are purely due to a variation in C atom mass in the reactants.

3 Results and Discussion

3.1 Baseline CO Trends.

Figure 2 plots typical CO trends for CH4. It shows the initial rise in CO, its peak to superequilibrium levels on the order of 104 ppm, and its relaxation toward equilibrium. While the peak superequilibrium CO levels are not strongly temperature dependent, increases in flame temperature make the flame thinner (as manifested by the width of the CO peak), increase the relaxation rate of CO toward equilibrium, and increase equilibrium concentration. While results are not shown for multiple pressure/power levels, increasing pressure decreases flame thickness, increases relaxation rates toward equilibrium, and decreases equilibrium CO levels.

Fig. 2
CO production from a methane–air flame at high-power condition (20 bar, 800 K) for three fixed adiabatic flame temperatures. Note: the secondary y-axis shows ppmvd corrected to 15%O2 based on a 0%H2 in blend.
Fig. 2
CO production from a methane–air flame at high-power condition (20 bar, 800 K) for three fixed adiabatic flame temperatures. Note: the secondary y-axis shows ppmvd corrected to 15%O2 based on a 0%H2 in blend.
Close modal

Consider next Fig. 3, which illustrates H2 blending effects. Increasing H2 levels decreases flame thickness and accelerates relaxation rates toward equilibrium, causing decreases in CO emissions at all residence time values.

Fig. 3
CO production for different fuel blends at high power (20 bar, 800 K) for fixed Tad = 1800 K. Note: the secondary y-axis shows ppmvd corrected to 15%O2 based on a 0%H2 in blend.
Fig. 3
CO production for different fuel blends at high power (20 bar, 800 K) for fixed Tad = 1800 K. Note: the secondary y-axis shows ppmvd corrected to 15%O2 based on a 0%H2 in blend.
Close modal

Having considered baseline trends, next are the more detailed characterizations. There are several useful ways to quantify and analyze this problem. These include characterizing H2 blending effects on (1) CO and COk emissions at a fixed residence time, and (2) residence time required to achieve a given CO values. Furthermore, given the intimate tradeoff between CO and NO emissions in combustor design, (3) examines NO emissions at fixed CO values. The rest of this paper summarizes these different representations in turn.

3.2 CO Emissions at Fixed Residence Time.

Start by considering CO emissions as a function of H2 blending at different residence times values. Figure 4 plots the low-power case, including the equilibrium CO levels, at residence times of 5, 15, and 25 ms. The figure shows the monotonically decreasing equilibrium CO levels with increased H2 blending noted earlier. This reduction is due to both kinetic effects and equilibrium effects. For example, the 15 ms case has essentially reached equilibrium at this condition, and so the reduction is entirely an equilibrium effect due to reduced C atoms in the fuel. In contrast, the 5 ms case is a mix of both effects.

Fig. 4
CO emissions at low power (1 bar, 300 K) as a function of fuel composition for fixed residence times and product temperature (Tad = 1800 K). Note: the secondary y-axis shows ppmvd corrected to 15%O2 based on a 0%H2 in blend.
Fig. 4
CO emissions at low power (1 bar, 300 K) as a function of fuel composition for fixed residence times and product temperature (Tad = 1800 K). Note: the secondary y-axis shows ppmvd corrected to 15%O2 based on a 0%H2 in blend.
Close modal

Similar results are observed at the other power conditions, with the key difference being that the medium power case reaches equilibrium in 3 ms and the high-power case in 1.9 ms. Figure 5 replots this data, normalizing CO emissions by their equilibrium value at that H2 blending composition. This normalization allows us to differentiate the kinetic effects of H2 blending from the fact that the overall C atom balance is also varying. If this ratio does not vary with H2 concentration, then the variations shown above in absolute CO concentrations are purely due to C atom balance considerations. If this ratio decreases with H2 concentration, then adding H2 also causes CO reductions by accelerating the oxidation of CO to CO2. The figure shows that this effect is particularly prominent at the 5 ms condition, where CO emissions are about six times their equilibrium value at 0% H2. In contrast, at 90%H2 blending, the CO emissions are 1.3 times above their equilibrium value. This CO/COeq ratio decreases with residence time, as the products approach equilibrium. Similar behavior is observed at the higher power conditions at lower residence times. Figure 6 illustrates a similar point by plotting both CO and COk emissions as a function of H2 blending. Note that both curves have the same value at 0%H2 by definition. The two curves are almost on top of each other at the lowest residence time case, 5 ms, indicating that the observed CO reductions are almost entirely due to the kinetic effect of H2 addition. In contrast, the COk curve is almost flat at 15 ms, even while CO emissions fall, indicating that the observed CO reduction is simply manifesting as a result of the drop in C atom mass in the reactants.

Fig. 5
CO/COeq at low power (1 bar, 300 K) condition as a function of fuel composition for fixed residence times and product temperature (Tad = 1800 K)
Fig. 5
CO/COeq at low power (1 bar, 300 K) condition as a function of fuel composition for fixed residence times and product temperature (Tad = 1800 K)
Close modal
Fig. 6
COk at low power (1 bar, 300 K) as a function of fuel composition for fixed residence times and product temperature (Tad = 1800 K)
Fig. 6
COk at low power (1 bar, 300 K) as a function of fuel composition for fixed residence times and product temperature (Tad = 1800 K)
Close modal

3.3 Residence Time Required to Achieve a Given CO Level.

Consider next the effect of H2 blending on the residence time required to achieve a given CO value. Figure 7 plots the pressure scaled time (P×τres) required to achieve CO levels of 1.1COeq of the pure CH4 case as a function of H2 blending. The value of 1.1 here is arbitrary and only influences the results for pure CH4, as the time to reach this equilibrium value is infinite. Picking some representative values, this required residence time is about 15 ms for the 1800 K low-power case with 100%CH4, that drops to about 9 ms with 50%H2 and to about 2 ms with 90%H2. This ratio drops by about a factor of 8 for low hydrogen content at the high-power condition, scaling almost linearly with pressure (for 90%H2 it is approximately four times lower) as seen in the pressure scaled representation of the τres. In other words, the above values become about 1.9, 1.1, and 0.5 ms with 0%, 50%, and 90%H2. These shorter residence times become particularly interesting when one begins to consider combustor design options for minimizing NOx and CO, as discussed in Sec. 3.4. We do not show results for the 95% and higher H2 cases as their highest CO peaks are lower than the equilibrium CO in the 0%H2 case.

Fig. 7
Scaled residence time value (P×τres, where P is pressure) required to achieve 1.1 ×COeq,100%CH4 as a function of fuel composition for varying product temperatures at low power (1 bar, 300 K) and high power (20 bar, 800 K)
Fig. 7
Scaled residence time value (P×τres, where P is pressure) required to achieve 1.1 ×COeq,100%CH4 as a function of fuel composition for varying product temperatures at low power (1 bar, 300 K) and high power (20 bar, 800 K)
Close modal

3.4 Interrelationship Between NO and CO Emissions.

Prior studies have examined H2 blending effects on NO emissions at a fixed residence time [8]. However, combustion systems optimized for elevated H2 levels can actually reduce NO emissions if residence times are reduced, which could potentially be done since CO levels drop with increased H2 blending. Such an option is an interesting regulatory option, whereby regulators could decrease NO emissions of the fleet by allowing designers to develop systems optimized for a certain H2 fraction. This section provides several visualizations of the relationship between NO and CO.

Figure 8 plots the relationship between NO and CO at the low-power condition, where τres is a parameter. Iso-τres values of 1.5 to 20 ms are also plotted in the figure. The different lines denote increased H2 blending. For a given H2 fraction, the figure shows the well-known tradeoff between NO and CO. In other words, increasing residence times decreases CO emissions but leads to increased NO emissions. Each of the curves shown in Fig. 8 represents a “Pareto front” or set of optimal solutions in a multi-objective problem for a fixed fuel blend. The figure also shows the monotonic shifting of the Pareto front with increased H2 blending, i.e., whereas the classical design tradeoff is one of balancing NO and CO, H2 addition allows one to decrease both. For example, at a fixed CO value, increasing H2 blending leads to reductions in NO. Similarly, at a fixed NO value, increasing H2 blending decreases CO. Similar points can be seen at the high-power condition plotted in Fig. 9, where the calculations have been run out to sufficient τres that CO has reached equilibrium by about 2 ms and does not change further, even while NO values monotonically rise with increased residence time. Figures 8 and 9 compare NO and CO under the same operating conditions. However, combustor design typically requires finding a residence time that optimizes CO at low power with NO at high-power conditions. To understand how H2 blending influences this tradeoff, Fig. 10 plots NO values at the high-power condition (800 K, 20 bar) as a function of CO emissions at the low-power condition (300 K, 1 bar) at the same residence times for Tad = 1800 K. This plot better illustrates the tradeoff in balancing NO at high temperature conditions, with CO emissions at low temperature conditions for a system designed with a given residence time. The figure again illustrates that H2 addition can provide benefits in NO, CO, or both by shifting the Pareto front down. This figure also demonstrates how H2 addition can enable increased turndown capabilities for gas turbines, by maintaining a fixed NO value at full power but being able to run to lower power conditions while maintaining the same CO levels.

Fig. 8
NO versus CO emissions in units of ng/J and ppmvd (15%O2 corrected) plotted at the same residence times for different %H2 in fuel blend at Tad = 1800 K, at low power (1 bar, 300 K). Iso-τ lines are shown for each blend. Note: the secondary x and y-axes show ppmvd corrected to 15%O2 based on a 0%H2 in blend.
Fig. 8
NO versus CO emissions in units of ng/J and ppmvd (15%O2 corrected) plotted at the same residence times for different %H2 in fuel blend at Tad = 1800 K, at low power (1 bar, 300 K). Iso-τ lines are shown for each blend. Note: the secondary x and y-axes show ppmvd corrected to 15%O2 based on a 0%H2 in blend.
Close modal
Fig. 9
NO versus CO emissions in units of ng/J and ppmvd (15%O2 corrected) plotted at the same residence times for different %H2 in fuel blend at Tad = 1800 K, at high power (20 bar, 800 K). Iso-τ lines are shown for each blend. Note: the secondary x and y-axes show ppmvd corrected to 15%O2 based on a 0%H2 in blend.
Fig. 9
NO versus CO emissions in units of ng/J and ppmvd (15%O2 corrected) plotted at the same residence times for different %H2 in fuel blend at Tad = 1800 K, at high power (20 bar, 800 K). Iso-τ lines are shown for each blend. Note: the secondary x and y-axes show ppmvd corrected to 15%O2 based on a 0%H2 in blend.
Close modal
Fig. 10
NO (ppm) at high power versus CO (ppm) at low power for Tad = 1800 K plotted at the same residence times for different %H2 in fuel blend. Iso-τ lines are shown for each blend case. Note: the secondary x and y-axes show ppmvd corrected to 15%O2 based on a 0%H2 in blend.
Fig. 10
NO (ppm) at high power versus CO (ppm) at low power for Tad = 1800 K plotted at the same residence times for different %H2 in fuel blend. Iso-τ lines are shown for each blend case. Note: the secondary x and y-axes show ppmvd corrected to 15%O2 based on a 0%H2 in blend.
Close modal

4 Conclusion

One of the most important tradeoffs in combustor design is balancing CO and NOx emissions, which generally have opposite tendencies. There has been significant attention to the question of how H2 blending will influence NOx emissions. This study has shown the important implications of H2 blending on CO emissions as well as how simultaneous consideration of NOx and CO can lead to a reconceptualization of both design and operational spaces. At a fixed residence time, H2 addition monotonically decreases CO emissions, due to both reduction in C atoms as well as kinetic effects. Similarly, H2 addition increases the rate at which superequilibrium CO formed in the flame relaxes down to postflame equilibrium levels. Finally, considering the NOx–CO tradeoff, H2 addition can lead to simultaneous reductions in both. These results have interesting implications for the utilization of H2 not only as a means for reducing the carbon intensity of gas turbine power but also expanding its operational flexibility.

Acknowledgment

This study has been partially supported by the Electric Power Research Institute (EPRI) and LCRI project titled “Investigation for Hydrogen GT DLN NOx Entitlement Curve,” and the Department of Energy (Grant No. DE-FE0032079).

Data Availability Statement

The datasets generated and supporting the findings of this article are obtainable from the corresponding author upon reasonable request.

Nomenclature

COk =

kinetic factor

MW =

molecular weight

ppmvd =

parts per million by volume, dry

Tad =

adiabatic flame temperature

Y =

mass fraction

α =

H2/CH4 fuel blend

τres =

residence time

Φ =

equivalence ratio

References

1.
O'Connor
,
J.
,
Noble
,
D.
, and
Lieuwen
,
T.
,
2022
,
Renewable Fuels
,
Cambridge University Press
,
Cambridge, UK
, Chap.
16
.
2.
[2]
Therkelsen
,
P.
,
Werts
,
T.
,
McDonell
,
V.
, and
Samuelsen
,
S.
,
2009
, “
Analysis of NOx Formation in a Hydrogen-Fueled Gas Turbine Engine
,”
ASME J. Eng. Gas Turbine Power
, 131(3), p.
031507
.10.1115/1.3028232
3.
Glanville
,
P.
,
Fridlyand
,
A.
,
Sutherland
,
B.
,
Liszka
,
M.
,
Zhao
,
Y.
,
Bingham
,
L.
, and
Jorgensen
,
K.
,
2022
, “
Impact of Hydrogen/Natural Gas Blends on Partially Premixed Combustion Equipment: NOx Emission and Operational Performance
,”
Energies
,
15
(
5
), p.
1706
.10.3390/en15051706
4.
Cellek
,
M. S.
, and
Pınarbaşı
,
A.
,
2018
, “
Investigations on Performance and Emission Characteristics of an Industrial Low Swirl Burner While Burning Natural Gas, Methane, Hydrogen-Enriched Natural Gas and Hydrogen as Fuels
,”
Int. J. Hydrogen Energy
,
43
(
2
), pp.
1194
1207
.10.1016/j.ijhydene.2017.05.107
5.
İlbaş
,
M.
, and
Yılmaz
,
İ.
,
2012
, “
Experimental Analysis of the Effects of Hydrogen Addition on Methane Combustion
,”
Int. J. Energy Res.
,
36
(
5
), pp.
643
647
.10.1002/er.1822
6.
Tuccillo
,
R.
,
Cameretti
,
M.
,
De Robbio
,
R.
,
Reale
,
F.
, and
Chiariello
,
F.
,
2019
, “
Methane-Hydrogen Blends in Micro Gas Turbines: Comparison of Different Combustor Concepts
,”
ASME
Paper No. GT2019-90229.10.1115/GT2019-90229
7.
Hoekstra
,
R.
,
Van Blarigan
,
P.
, and
Mulligan
,
N.
,
1996
, “
NOx Emissions and Efficiency of Hydrogen, Natural Gas, and Hydrogen/Natural Gas Blended Fuels
,”
SAE
Paper No. 961103.10.4271/961103
8.
Breer
,
B.
,
Rajagopalan
,
H.
,
Godbold
,
C.
,
Johnson
,
H.
,
Emerson
,
B.
,
Acharya
,
V.
,
Sun
,
W.
,
Noble
,
D.
, and
Lieuwen
,
T.
,
2023
, “
Numerical Investigation of NOx Production From Premixed Hydrogen/Methane Fuel Blends
,”
Combust. Flame
,
255
, p.
112920
.10.1016/j.combustflame.2023.112920
9.
Fellner
,
C.
,
2006
, “
Standards of Performance for Stationary Combustion Turbines
,”
Environmental Protection Agency
, Washington, DC, Report No. 06-5945.
10.
Douglas
,
C. M.
,
Shaw
,
S. L.
,
Martz
,
T. D.
,
Steele
,
R. C.
,
Noble
,
D. R.
,
Emerson
,
B. L.
, and
Lieuwen
,
T. C.
,
2022
, “
Pollutant Emissions Reporting and Performance Considerations for Hydrogen–Hydrocarbon Fuels in Gas Turbines
,”
ASME J. Eng. Gas Turbines Power
,
144
(
9
), p.
091003
.10.1115/1.4054949
11.
Environmental Protection Agency, 2024, “
National Emissions Inventory (NEI) Air Pollutant Emissions Trends Data
,”
Environmental Protection Agency
, Washington, DC, accessed Oct. 6, 2024, https://www.epa.gov/stationary-sources-air-pollution
12.
Environmental Protection Agency, 2024, “
The Green Book: Nonattainment Areas for Criteria Pollutants
,”
Environmental Protection Agency
, Washington, DC, accessed Oct. 6, 2024, https://www.epa.gov/green-book
13.
EPRI,
2023
, “
Hydrogen Blending Demonstration Synopsis: EPRI-Affiliated Testing Summary
,”
Low-Carbon Resources Initiative, Electric Power Research Institute, Palo Alto, CA
, Report No. 3002028175.
14.
Harper
,
J.
,
Cloyd
,
S.
,
Pigon
,
T.
,
Thomas
,
B.
,
Wilson
,
J.
,
Johnson
,
E.
, and
Noble
,
D.
,
2023
, “
Hydrogen Co-Firing Demonstration at Georgia Power's Plant McDonough: M501G Gas Turbine
,”
ASME
Paper No. GT2023-102660.10.1115/GT2023-102660
15.
EPRI
,
2022
, “
Southern Company Hydrogen Blending Test Report: Mitsubishi Power 501G 20% Hydrogen Test
,”
EPRI
,
Palo Alto, CA
, Report No. 3002025835.
16.
Buis
,
A.
,
2019
, “
The Atmosphere: Getting a Handle on Carbon
,”
NASA Science Editorial Team
, accessed Oct. 6, 2024, https://science.nasa.gov/earth/climate-change/greenhouse-gases/the-atmosphere-getting-a-handle-on-carbon-dioxide/
17.
Ansys,
2021
, “
Ansys Chemkin-Pro
,” Ansys, Canonsburg, PA.
18.
Park
,
J.
,
Xu
,
R.
,
Lu
,
T.
, and
Wang
,
H.
,
2022
, “
Skeletal and Reduced Model of NOx Formation in C1/A2 Jet Fuel Blend
,” accessed June 1, 2022, https://web.stanford.edu/group/haiwanglab/HyChem/pages/download.html
19.
Smith
,
G.
,
Golden
,
D.
,
Frenklach
,
M.
,
Moriarty
,
N.
,
Eit-Eneer
,
B.
,
Goldenberg
,
M.
,
Bowman
,
C.
, et al., 2024, “
Gri-Mech 3.0
,” accessed Aug. 1, 2024, http://www.me.berkeley.edu/gri_mech/
20.
Glarborg
,
P.
,
Miller
,
J. A.
,
Ruscic
,
B.
, and
Klippenstein
,
S. J.
,
2018
, “
Modeling Nitrogen Chemistry in Combustion
,”
Prog. Energy Combust. Sci.
,
67
, pp.
31
68
.10.1016/j.pecs.2018.01.002
21.
University of California at San Diego, 2024, “
University of California at San Diego, Chemical-Kinetic Mechanisms for Combustion Applications
,” University of California at San Diego, San Diego, CA, accessed Aug. 1, 2024, https://web.eng.ucsd.edu/mae/groups/combustion/mechanism.html