Elsevier

Combustion and Flame

Volume 234, December 2021, 111646
Combustion and Flame

Ignition delay times of n-butane and i-butane under O2/CO2 atmospheres: Shock tube experiments and kinetic model

https://doi.org/10.1016/j.combustflame.2021.111646Get rights and content

Abstract

Pressurized oxy-fuel combustion is deemed an advanced oxy-fuel combustion technique due to the lower cost and little decrease in the generating efficiency compared to conventional PC system without CO2 capture. Butane is an important composition of petroleum and natural gas. In this work, the ignition delay times (IDTs) for n-butane and i-butane under O2/CO2 atmospheres were measured in a shock tube at different equivalence ratios (Φ) and pressures. Based on our previous C1single bondC3 model Oxymech2.0 and the C4 sub-model of Aramco3.0, a chemical kinetic model Oxymech2.0 Plus was updated. Oxymech2.0 Plus was validated by the newly measured IDTs of butane under O2/CO2 atmospheres. The model was also validated by the literature experimental data for IDTs of n-butane and i-butane at O2/Ar and O2/N2 atmospheres, laminar flame speeds of n-butane and i-butane in air, and species profiles of n-butane and i-butane pyrolysis. The comparison between Oxymech2.0 Plus and Aramco 3.0 models was also conducted in detailed. The results show that the updated reactions in the study significantly improve the prediction of IDTs, laminar flame speeds and species profiles of butane in O2/CO2 atmospheres. The effects of equivalence ratios and CO2 on the IDTs of n-butane and i-butane were analyzed.

Introduction

Pressurized oxy-fuel combustion is deemed an advanced oxy-fuel combustion technique due to the lower cost and little decrease in the generating efficiency compared to conventional PC system without CO2 capture [1,2]. A 50-MWh natural gas-fired Allam Cycle plant has been built in La Porte, Texas, which generates energy by using oxy-fuel and direct-fired combustion in a high-pressure turbine [3]. Therefore, the study of ignition characteristics and chemical kinetic mechanisms for fuel combustion under pressurized oxy-fuel condition is essential for actual combustor design.

Butane is not only an important component of petroleum and natural gas, but also the smallest alkane with isomers. The ignition delay time (IDT), which can be measured in shock tubes, is a well-established feature for the validation of fuel combustion chemical kinetic models. The IDTs of n-butane and i-butane under O2/Ar or O2/N2 atmospheres measured in shock tubes have been reported by many researchers [4], [5], [6], [7], [8], [9]. The shock tube experiments and conditions of butane IDTs were collected and listed in Table 1. Ogura et al. [4] proposed a model based on Gri3.0 [10] that well reproduced their measured IDTs. Healy et al. [5,6] proposed the NUIG model based on the kinetic models developed by Curran et al. [11,12] and Bourque et al. [13] that satisfactorily predicted their experimental results. In 2019, He et al. [7] found that the NUIG [5,6], USC2.0 [14] and Sandiego [15] models slightly underestimated the IDTs and updated the NUIG model to match their experiments. Zhou et al. [16] proposed the Aramco3.0 model, which is widely adopted due to the good performance for predicting the IDTs, laminar flame speeds and species profiles of C1single bondC4 hydrocarbon and oxygenated hydrocarbon species. However, as far as we know, few studies on the IDTs of butane under O2/CO2 atmospheres exist at present.

Many studies found that widely adopted chemical kinetic models such as Aramco2.0 [17], Ranzi [18], USC 2.0 [14], FFCM-1 [19] and Gri 3.0 [10], cannot accurately predict the IDTs and laminar flame speeds of CH4, C2H6 and C3H8 under O2/CO2 atmospheres [20], [21], [22]. In our previous work, Liu et al. [20,21] and Xia et al. [22] proposed a C1single bondC3 kinetic model termed the Oxymech2.0 based on the Aramco2.0 model, which produced good prediction results of the IDTs, laminar flame speeds and species profiles of CH4, C2H6 and C3H8 under O2/N2 and O2/CO2 atmospheres. It is necessary to expand the fuels tested by the Oxymech2.0 to include butane. Therefore, the butane sub-model needs to be updated and improved in O2/CO2 atmospheres in this work.

The updated reactions in the Oxymech2.0 can be divided into three parts. The first part includes reactions containing CO2, such as CO + OH 〈=〉 CO2 + H. The second part contains third-body collision reactions. The rate constants of many third-body collision reactions, such as H + O2 (+M) 〈=〉 HO2 (+M), CH3 + CH3 (+M) 〈=〉 C2H6 (+M), and C2H4 + H (+M) 〈=〉 C2H5 (+M), were updated. The third part represents fuel H abstraction reactions. The rate constants of fuel H abstraction reactions, such as CH4 + OH 〈=〉 CH3 + H2O, C2H6 + HO2 〈=〉 C2H5 + H2O2, C2H6 + OH 〈=〉 C2H5 + H2O, C3H8 + H 〈=〉 iC3H7 + H2, C3H8 + OH 〈=〉 iC3H7 + H2O and C3H8 + H 〈=〉 nC3H7 + H2, were updated. According to the model development work in our previous work [22], the fuel decomposition and H abstraction reactions are crucial to fuel ignition. Therefore, the C4 sub-model in Aramco3.0 needs to be updated and combine with Oxymech2.0 for the prediction of butane combustion in O2/CO2 atmospheres in this work.

Besides IDTs, species profiles and laminar flame speeds are two other commonly used features for validating fuel combustion kinetic models. In 2018, Li et al. [23] measured the species profiles of C1–C4 products in n-butane and i-butane pyrolysis experiments at air conditions. They proposed a model (Li-2018) based on their previous work [24,25]. The Li-2018 model well predicted the species profiles of both n-butane and i-butane. The laminar flame speeds of butane in O2/N2 atmospheres and in air conditions at 1 atm were measured by many research groups [26], [27], [28], [29].

In this work, shock tube experiments of IDTs for n-butane and i-butane diluted in carbon dioxide were conducted at pressures of 1.0 and 10.0 atm, at three different equivalence ratios (0.5, 1.0, 2.0), and at temperatures of 1200–1600 K. An updated oxy-fuel combustion kinetic model based on the Oxymech2.0 and the C4 sub-model of Aramco3.0, was proposed and referred to as the Oxymech2.0 Plus. The Oxymech2.0 Plus was validated by the newly measured IDTs of butane under O2/CO2 atmospheres, and by the reported IDTs of butane under O2/N2 atmospheres [5,6] and O2/Ar atmospheres [4]. The model was also validated by the selected experimental data of butane laminar flame speeds and species profiles [23,[26], [27], [28], [29]]. The prediction performance of the Oxymech2.0 Plus on IDTs, laminar flame speeds and species profiles were compared with that of the Aramco3.0. The effects of equivalence ratios and CO2 on the ignition of butane were also discussed in detailed.

Section snippets

Experimental apparatus

The experiments were performed in a stainless-steel shock tube facility at Huazhong University of Science and Technology (HUST). Details of the shock tube were reported in our previous work [22]. High-purity nitrogen (99.999%) and helium (99.99%) blend was used as the driving gas of the shock tube. The experimental mixtures under different conditions are listed in Table 2. The purities of n-butane, i-butane are 99.99% and that of O2 and CO2 are 99.999%. All the gasses were purchased from Wuhan

Modeling

A detailed kinetic model was developed based on the Oxymech2.0 and the C4 sub-model in the Aramco3.0. The developed model is referred to as the Oxymech2.0 Plus. The improvements of the Oxymech2.0 Plus are summarized in three parts. The first is about the butane sub-mechanism. Many researchers such as He et al. [7] found that fuel decomposition reactions and H abstraction reactions are important to the ignition, so we focus our updating on these reactions. In the n-butane sub-mechanism, the rate

IDTs of n-butane and i-butane in O2/CO2 atmospheres

Figure 3(a) shows the measured IDTs of n-butane diluted in CO2 at the equivalence ratios of 0.5, 1.0, 2.0 and the pressures of 1 and 10 atm. The results show that the IDTs decrease with the decrease in the equivalence ratio from 2.0 to 0.5 at 1 atm. The difference of the IDTs among the three equivalence ratios significantly reduces at 10 atm, indicating a weaker effect of the equivalence ratios on the IDTs at a higher pressure.

Figure 3(b) shows the measured IDTs of i-butane diluted in CO2 at

Conclusion

The IDTs of n-butane and i-butane at O2/CO2 atmospheres were measured in a shock tube. The results show that the IDTs of n-butane and i-butane decrease with the increase in the pressure, and the difference of the IDTs among the three equivalence ratios significantly reduces at 10 atm. An updated kinetic model basing on the Oxymech2.0 and the C4 sub-model of Aramco1.3, was proposed and referred to as the Oxymech2.0 Plus. The Oxymech2.0 Plus was validated by the newly measured IDTs of butane

Declaration of Competing Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References (50)

  • U. Burke et al.

    A detailed chemical kinetic modeling, ignition delay time and jet-stirred reactor study of methanol oxidation

    Combust. Flame

    (2016)
  • E. Ranzi et al.

    Hierarchical and comparative kinetic modeling of laminar flame speeds of hydrocarbon and oxygenated fuels

    Prog. Energy Combust. Sci.

    (2012)
  • Y. Liu et al.

    Ignition delay times of ethane under O2/CO2 atmosphere at different pressures by shock tube and simulation methods

    Combust. Flame

    (2019)
  • W. Xia et al.

    Shock tube and modeling study of ignition delay times of propane under O2/CO2/Ar atmosphere

    Combust. Flame

    (2020)
  • W. Li et al.

    Experimental and kinetic modeling investigation on pyrolysis and combustion of n-butane and i-butane at various pressures

    Combust. Flame

    (2018)
  • J.H. Cai et al.

    Experimental and kinetic modeling study of 2-butanol pyrolysis and combustion

    Combust. Flame

    (2013)
  • T. Hirasawa et al.

    Determination of laminar flame speeds using digital particle image velocimetry: binary Fuel blends of ethylene, n-Butane, and toluene

    Proc. Combust. Inst.

    (2002)
  • O. Park et al.

    Chemical kinetic model uncertainty minimization through laminar flame speed measurements

    Combust. Flame

    (2016)
  • K.J. Bosschaart et al.

    The laminar burning velocity of flames propagating in mixtures of hydrocarbons and air measured with the heat flux method

    Combust. Flame

    (2004)
  • O. Pryor et al.

    Measurements and interpretation of shock tube ignition delay times in highly CO2 diluted mixtures using multiple diagnostics

    Combust. Flame

    (2017)
  • J.K. Shao et al.

    Ignition delay times of methane and hydrogen highly diluted in carbon dioxide at high pressures up to 300 atm

    Proc. Combust. Inst.

    (2019)
  • E. Ninnemann et al.

    Reflected shock-initiated ignition probed via simultaneous lateral and endwall high-speed imaging with a transparent, cylindrical test-section

    Combust. Flame

    (2021)
  • S.L. Peukert et al.

    High temperature rate constants for H/D plus n-C4H10 and i-C4H10

    Proc. Combust. Inst.

    (2015)
  • C.W. Zhou et al.

    A comprehensive experimental and modeling study of isobutene oxidation

    Combust. Flame

    (2016)
  • P.H. Stewart et al.

    Pressure and temperature dependence of reactions proceeding via a bound complex. 2. Application to 2CH3 → C2H5 + H

    Combust. Flame

    (1989)
  • Cited by (11)

    • Hydrogen/PRF skeletal mechanism study based on shock tube experiments and kinetic analysis

      2023, Fuel
      Citation Excerpt :

      See Supplementary Material). The IDT, as shown in Fig. 2, is defined by the time interval between the pressure rise caused by the reflected shock wave and the timing of the steepest slope of OH* emission caused by the mixture ignition and horizontal time axis [41–44]. The zero-dimensional, constant volume, and adiabatic reactor in Chemkin 17.0 was applied to calculate and analyze the ignition process of hydrogen/n-heptane mixture.

    • Shock tube experiments and numerical study on ignition delay times of ethane in super lean and ultra-lean combustion

      2022, Combustion and Flame
      Citation Excerpt :

      Section 5 concludes the paper. The shock tube in the present study had been reported in detail elsewhere [17,22,23], and only a brief recap is given here. The shock tube consists of two sections: a 4-m long driver section with 10 cm diameter and an 8-m long driven section.

    View all citing articles on Scopus
    View full text