############################################################## # # 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])