##############################################################
#
# H2 mechanism
# Developed by N. Zettervall
# FOI, Sweden, 2017-02-15
#
# Vincent-Randonnier, A., Sabelnikov, V., Ristori, A.,
# Zettervall, N., & Fureby, C. (2018).
# A combined experimental and computational study
# of the LAPCAT II supersonic combustor.
# In 22nd AIAA international space planes and
# hypersonics systems and technologies conference (p. 5208).
#
##############################################################

units(length = “cm”, time = “s”, quantity = “mol”, act_energy = “cal/mol”)

ideal_gas(name = “gas”,
elements = ” H O C N “,
species = “”” O2 N2 H2O H O OH HO2 H2
H2O2″””,
reactions = “all”,
transport = “Mix”,
initial_state = state(temperature = 300.0,
pressure = OneAtm) )

 

####################################################################################
#
# Thermodata from GRI3.0
#
####################################################################################

species(name = “H2”,
atoms = ” H:2 “,
thermo = (
NASA( [ 200.00, 1000.00], [ 2.344331120E+00, 7.980520750E-03,
-1.947815100E-05, 2.015720940E-08, -7.376117610E-12,
-9.179351730E+02, 6.830102380E-01] ),
NASA( [ 1000.00, 3500.00], [ 3.337279200E+00, -4.940247310E-05,
4.994567780E-07, -1.795663940E-10, 2.002553760E-14,
-9.501589220E+02, -3.205023310E+00] )
),
transport = gas_transport(
geom = “linear”,
diam = 2.92,
well_depth = 38.00,
polar = 0.79,
rot_relax = 280.00),
note = “TPIS78”
)

species(name = “H”,
atoms = ” H:1 “,
thermo = (
NASA( [ 200.00, 1000.00], [ 2.500000000E+00, 7.053328190E-13,
-1.995919640E-15, 2.300816320E-18, -9.277323320E-22,
2.547365990E+04, -4.466828530E-01] ),
NASA( [ 1000.00, 3500.00], [ 2.500000010E+00, -2.308429730E-11,
1.615619480E-14, -4.735152350E-18, 4.981973570E-22,
2.547365990E+04, -4.466829140E-01] )
),
transport = gas_transport(
geom = “atom”,
diam = 2.05,
well_depth = 145.00),
note = “L 7/88”
)

species(name = “O”,
atoms = ” O:1 “,
thermo = (
NASA( [ 200.00, 1000.00], [ 3.168267100E+00, -3.279318840E-03,
6.643063960E-06, -6.128066240E-09, 2.112659710E-12,
2.912225920E+04, 2.051933460E+00] ),
NASA( [ 1000.00, 3500.00], [ 2.569420780E+00, -8.597411370E-05,
4.194845890E-08, -1.001777990E-11, 1.228336910E-15,
2.921757910E+04, 4.784338640E+00] )
),
transport = gas_transport(
geom = “atom”,
diam = 2.75,
well_depth = 80.00),
note = “L 1/90”
)

species(name = “O2”,
atoms = ” O:2 “,
thermo = (
NASA( [ 200.00, 1000.00], [ 3.782456360E+00, -2.996734160E-03,
9.847302010E-06, -9.681295090E-09, 3.243728370E-12,
-1.063943560E+03, 3.657675730E+00] ),
NASA( [ 1000.00, 3500.00], [ 3.282537840E+00, 1.483087540E-03,
-7.579666690E-07, 2.094705550E-10, -2.167177940E-14,
-1.088457720E+03, 5.453231290E+00] )
),
transport = gas_transport(
geom = “linear”,
diam = 3.46,
well_depth = 107.40,
polar = 1.60,
rot_relax = 3.80),
note = “TPIS89”
)

species(name = “OH”,
atoms = ” O:1 H:1 “,
thermo = (
NASA( [ 200.00, 1000.00], [ 3.992015430E+00, -2.401317520E-03,
4.617938410E-06, -3.881133330E-09, 1.364114700E-12,
3.615080560E+03, -1.039254580E-01] ),
NASA( [ 1000.00, 3500.00], [ 3.092887670E+00, 5.484297160E-04,
1.265052280E-07, -8.794615560E-11, 1.174123760E-14,
3.858657000E+03, 4.476696100E+00] )
),
transport = gas_transport(
geom = “linear”,
diam = 2.75,
well_depth = 80.00),
note = “RUS 78”
)

species(name = “H2O”,
atoms = ” H:2 O:1 “,
thermo = (
NASA( [ 200.00, 1000.00], [ 4.198640560E+00, -2.036434100E-03,
6.520402110E-06, -5.487970620E-09, 1.771978170E-12,
-3.029372670E+04, -8.490322080E-01] ),
NASA( [ 1000.00, 3500.00], [ 3.033992490E+00, 2.176918040E-03,
-1.640725180E-07, -9.704198700E-11, 1.682009920E-14,
-3.000429710E+04, 4.966770100E+00] )
),
transport = gas_transport(
geom = “nonlinear”,
diam = 2.60,
well_depth = 572.40,
dipole = 1.84,
rot_relax = 4.00),
note = “L 8/89”
)

species(name = “HO2”,
atoms = ” H:1 O:2 “,
thermo = (
NASA( [ 200.00, 1000.00], [ 4.301798010E+00, -4.749120510E-03,
2.115828910E-05, -2.427638940E-08, 9.292251240E-12,
2.948080400E+02, 3.716662450E+00] ),
NASA( [ 1000.00, 3500.00], [ 4.017210900E+00, 2.239820130E-03,
-6.336581500E-07, 1.142463700E-10, -1.079085350E-14,
1.118567130E+02, 3.785102150E+00] )
),
transport = gas_transport(
geom = “nonlinear”,
diam = 3.46,
well_depth = 107.40,
rot_relax = 1.00),
note = “L 5/89”
)

species(name = “H2O2”,
atoms = ” H:2 O:2 “,
thermo = (
NASA( [ 200.00, 1000.00], [ 4.276112690E+00, -5.428224170E-04,
1.673357010E-05, -2.157708130E-08, 8.624543630E-12,
-1.770258210E+04, 3.435050740E+00] ),
NASA( [ 1000.00, 3500.00], [ 4.165002850E+00, 4.908316940E-03,
-1.901392250E-06, 3.711859860E-10, -2.879083050E-14,
-1.786178770E+04, 2.916156620E+00] )
),
transport = gas_transport(
geom = “nonlinear”,
diam = 3.46,
well_depth = 107.40,
rot_relax = 3.80),
note = “L 7/88”
)

species(name = “N2”,
atoms = ” N:2 “,
thermo = (
NASA( [ 300.00, 1000.00], [ 3.298677000E+00, 1.408240400E-03,
-3.963222000E-06, 5.641515000E-09, -2.444854000E-12,
-1.020899900E+03, 3.950372000E+00] ),
NASA( [ 1000.00, 5000.00], [ 2.926640000E+00, 1.487976800E-03,
-5.684760000E-07, 1.009703800E-10, -6.753351000E-15,
-9.227977000E+02, 5.980528000E+00] )
),
transport = gas_transport(
geom = “linear”,
diam = 3.62,
well_depth = 97.53,
polar = 1.76,
rot_relax = 4.00),
note = “121286”
)

#——————————————————————————-
# Reaction data
#——————————————————————————-

reaction(‘H2 + O2 => H + HO2’, [7.400000e+05, 2.43, 53500.0])
three_body_reaction( ‘ H2 + M => H + H + M’, [4.57E+19, -1.40, 105100], efficiencies=’H2:2.5 H2O:12.0 H:1.0 O:1.0 OH:1.0 O2:1.0 N2:1.0 HO2:1.0 H2O2:1.0′)
reaction( “HO2 + H2 => H2O2 + H”, [3.0e+06, 2, 21000])
reaction( “H + O2 => OH + O”, [2.45E+14, 0, 16800])
reaction( “OH + O => H + O2”, [1.2E+13, 0, 690])
reaction( “O + H2 => OH + H”, [1.80000E+10, 1, 8826])
reaction( “OH + H => O + H2”, [8.00000E+09, 1, 6760])
reaction( “H2 + OH => H2O + H”, [1.17000E+09, 1.3, 3626])
reaction( “H2O + H => H2 + OH”, [5.09000E+09, 1.3, 18588])
reaction( “OH + OH => O + H2O”, [6.00000E+08, 1.3, 0])
reaction( “O + H2O => OH + OH”, [5.90000E+09, 1.3, 17029])
three_body_reaction( “H + O2 + M => HO2 + M”, [1.8E+18, -0.8, 0],
efficiencies = ” H2:1 H2O:6.5 N2:0.4 O2:0.4 “)
reaction( “H + HO2 => OH + OH”, [1.50000E+14, 0, 1004])
reaction( “H + HO2 => H2 + O2”, [2.50000E+13, 0, 700])
reaction( “OH + HO2 => H2O + O2”, [2.00000E+13, 0, 1000])
reaction( “HO2 + HO2 => H2O2 + O2”, [8.0E+13, 0, 0])
three_body_reaction( “H2O2 + M => OH + OH + M”, [1.30000E+17, 0, 34500])
three_body_reaction( “OH + OH + M => H2O2 + M”, [9.86000E+14, 0, -5070])
reaction( “H2O2 + OH => H2O + HO2”, [1.00000E+13, 0, 1800])
reaction( “H2O + HO2 => H2O2 + OH”, [2.86000E+13, 0, 32790])
three_body_reaction( “OH + H + M => H2O + M”, [2.20000E+22, -2, 0])
three_body_reaction( “H + H + M => H2 + M”, [1.80000E+18, -1, 0])