Added test to confirm that ck2cti.py detects missing thermo data

This commit is contained in:
Ray Speth
2012-07-26 20:49:06 +00:00
parent 0e21bfd0e0
commit db57f3a200
2 changed files with 99 additions and 0 deletions

View File

@@ -0,0 +1,84 @@
ELEMENTS
O H AR
END
SPECIES
H2 H O O2 OH H2O HO2 H2O2
AR
END
THERMO ALL
300.000 1000.000 5000.000
O L 1/90O 1 00 00 00G 200.000 3500.000 1000.000 1
2.56942078E+00-8.59741137E-05 4.19484589E-08-1.00177799E-11 1.22833691E-15 2
2.92175791E+04 4.78433864E+00 3.16826710E+00-3.27931884E-03 6.64306396E-06 3
-6.12806624E-09 2.11265971E-12 2.91222592E+04 2.05193346E+00 4
O2 TPIS89O 2 00 00 00G 200.000 3500.000 1000.000 1
3.28253784E+00 1.48308754E-03-7.57966669E-07 2.09470555E-10-2.16717794E-14 2
-1.08845772E+03 5.45323129E+00 3.78245636E+00-2.99673416E-03 9.84730201E-06 3
-9.68129509E-09 3.24372837E-12-1.06394356E+03 3.65767573E+00 4
H L 7/88H 1 00 00 00G 200.000 3500.000 1000.000 1
2.50000001E+00-2.30842973E-11 1.61561948E-14-4.73515235E-18 4.98197357E-22 2
2.54736599E+04-4.46682914E-01 2.50000000E+00 7.05332819E-13-1.99591964E-15 3
2.30081632E-18-9.27732332E-22 2.54736599E+04-4.46682853E-01 4
OH RUS 78O 1H 1 00 00G 200.000 3500.000 1000.000 1
3.09288767E+00 5.48429716E-04 1.26505228E-07-8.79461556E-11 1.17412376E-14 2
3.85865700E+03 4.47669610E+00 3.99201543E+00-2.40131752E-03 4.61793841E-06 3
-3.88113333E-09 1.36411470E-12 3.61508056E+03-1.03925458E-01 4
H2O L 8/89H 2O 1 00 00G 200.000 3500.000 1000.000 1
3.03399249E+00 2.17691804E-03-1.64072518E-07-9.70419870E-11 1.68200992E-14 2
-3.00042971E+04 4.96677010E+00 4.19864056E+00-2.03643410E-03 6.52040211E-06 3
-5.48797062E-09 1.77197817E-12-3.02937267E+04-8.49032208E-01 4
HO2 L 5/89H 1O 2 00 00G 200.000 3500.000 1000.000 1
4.01721090E+00 2.23982013E-03-6.33658150E-07 1.14246370E-10-1.07908535E-14 2
1.11856713E+02 3.78510215E+00 4.30179801E+00-4.74912051E-03 2.11582891E-05 3
-2.42763894E-08 9.29225124E-12 2.94808040E+02 3.71666245E+00 4
H2O2 L 7/88H 2O 2 00 00G 200.000 3500.000 1000.000 1
4.16500285E+00 4.90831694E-03-1.90139225E-06 3.71185986E-10-2.87908305E-14 2
-1.78617877E+04 2.91615662E+00 4.27611269E+00-5.42822417E-04 1.67335701E-05 3
-2.15770813E-08 8.62454363E-12-1.77025821E+04 3.43505074E+00 4
AR 120186AR 1 G 300.000 5000.000 1000.000 1
0.02500000E+02 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 2
-0.07453750E+04 0.04366000E+02 0.02500000E+02 0.00000000E+00 0.00000000E+00 3
0.00000000E+00 0.00000000E+00-0.07453750E+04 0.04366000E+02 4
END
REACTIONS
2O+M<=>O2+M 1.200E+17 -1.000 .00
H2/ 2.40/ H2O/15.40/ AR/ .83/
O+H+M<=>OH+M 5.000E+17 -1.000 .00
H2/2.00/ H2O/6.00/ AR/ .70/
O+H2<=>H+OH 3.870E+04 2.700 6260.00
O+HO2<=>OH+O2 2.000E+13 .000 .00
O+H2O2<=>OH+HO2 9.630E+06 2.000 4000.00
H+2O2<=>HO2+O2 2.080E+19 -1.240 .00
H+O2+H2O<=>HO2+H2O 11.26E+18 -.760 .00
H+O2+AR<=>HO2+AR 7.000E+17 -.800 .00
H+O2<=>O+OH 2.650E+16 -.6707 17041.00
2H+M<=>H2+M 1.000E+18 -1.000 .00
H2/ .00/ H2O/ .00/ AR/ .63/
2H+H2<=>2H2 9.000E+16 -.600 .00
2H+H2O<=>H2+H2O 6.000E+19 -1.250 .00
H+OH+M<=>H2O+M 2.200E+22 -2.000 .00
H2/ .73/ H2O/3.65/ AR/ .38/
H+HO2<=>O+H2O 3.970E+12 .000 671.00
H+HO2<=>O2+H2 4.480E+13 .000 1068.00
H+HO2<=>2OH 0.840E+14 .000 635.00
H+H2O2<=>HO2+H2 1.210E+07 2.000 5200.00
H+H2O2<=>OH+H2O 1.000E+13 .000 3600.00
OH+H2<=>H+H2O 2.160E+08 1.510 3430.00
2OH(+M)<=>H2O2(+M) 7.400E+13 -.370 .00
LOW / 2.300E+18 -.900 -1700.00/
TROE/ .7346 94.00 1756.00 5182.00 /
H2/2.00/ H2O/6.00/ AR/ .70/
2OH<=>O+H2O 3.570E+04 2.400 -2110.00
OH+HO2<=>O2+H2O 1.450E+13 .000 -500.00
DUPLICATE
OH+H2O2<=>HO2+H2O 2.000E+12 .000 427.00
DUPLICATE
OH+H2O2<=>HO2+H2O 1.700E+18 .000 29410.00
DUPLICATE
2HO2<=>O2+H2O2 1.300E+11 .000 -1630.00
DUPLICATE
2HO2<=>O2+H2O2 4.200E+14 .000 12000.00
DUPLICATE
OH+HO2<=>O2+H2O 0.500E+16 .000 17330.00
DUPLICATE
END

View File

@@ -20,3 +20,18 @@ class chemkinConverterTest(unittest.TestCase):
self.assertEqual(ref.speciesNames(), gas.speciesNames())
self.assertTrue((ref.reactantStoichCoeffs() == gas.reactantStoichCoeffs()).all())
def test_missingElement(self):
if os.path.exists('h2o2_missingElement.cti'):
os.remove('h2o2_missingElement.cti')
self.assertRaises(ck2cti.InputParseError,
lambda: ck2cti.convertMech('../data/h2o2_missingElement.inp',
quiet=True))
def test_missingThermo(self):
if os.path.exists('h2o2_missingThermo.cti'):
os.remove('h2o2_missingThermo.cti')
self.assertRaises(ck2cti.InputParseError,
lambda: ck2cti.convertMech('../data/h2o2_missingThermo.inp',
quiet=True))