Note
Go to the end to download the full example code
Constant-pressure, adiabatic kinetics simulation#
Requires: cantera >= 2.5.0, matplotlib >= 2.0
Initializing reactor network.
Reactor 0: 12 variables.
0 sensitivity params.
Number of equations: 12
Maximum time step: 0
t [s] T [K] P [Pa] u [J/kg]
1.000e-05 1001.000 101325.000 620761.940774
2.000e-05 1001.000 101325.000 620761.940024
3.000e-05 1001.000 101325.000 620761.937443
4.000e-05 1001.000 101325.000 620761.932033
5.000e-05 1001.000 101325.000 620761.922251
6.000e-05 1001.000 101325.000 620761.905712
7.000e-05 1001.000 101325.000 620761.878725
8.000e-05 1001.000 101325.000 620761.835580
9.000e-05 1001.001 101325.000 620761.767437
1.000e-04 1001.001 101325.000 620761.660607
1.100e-04 1001.001 101325.000 620761.493864
1.200e-04 1001.002 101325.000 620761.234248
1.300e-04 1001.003 101325.000 620760.830463
1.400e-04 1001.005 101325.000 620760.202386
1.500e-04 1001.008 101325.000 620759.224189
1.600e-04 1001.012 101325.000 620757.696642
1.700e-04 1001.019 101325.000 620755.300403
1.800e-04 1001.030 101325.000 620751.514343
1.900e-04 1001.048 101325.000 620745.465936
2.000e-04 1001.076 101325.000 620735.641598
2.100e-04 1001.124 101325.000 620719.289981
2.200e-04 1001.204 101325.000 620691.109787
2.300e-04 1001.349 101325.000 620640.173038
2.400e-04 1001.623 101325.000 620542.299510
2.500e-04 1002.182 101325.000 620340.602686
2.600e-04 1003.397 101325.000 619897.667914
2.700e-04 1006.134 101325.000 618893.694074
2.800e-04 1012.324 101325.000 616618.206745
2.900e-04 1026.776 101325.000 611309.993847
Limiting global state vector component 1 (dt = 1e-05): 37.7079 > 20
2.950e-04 1040.882 101325.000 606137.126425
Limiting global state vector component 1 (dt = 1e-05): 67.2368 > 20
Limiting global state vector component 1 (dt = 5e-06): 23.8769 > 20
2.975e-04 1051.164 101325.000 602371.487295
Limiting global state vector component 1 (dt = 1e-05): 93.9215 > 20
Limiting global state vector component 1 (dt = 5e-06): 32.0613 > 20
3.000e-04 1064.775 101325.000 597393.368098
Limiting global state vector component 1 (dt = 1e-05): 136.833 > 20
Limiting global state vector component 1 (dt = 5e-06): 44.381 > 20
3.025e-04 1083.260 101325.000 590644.956051
Limiting global state vector component 1 (dt = 1e-05): 213.457 > 20
Limiting global state vector component 1 (dt = 5e-06): 64.1412 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 25.9868 > 20
3.038e-04 1095.086 101325.000 586334.638376
Limiting global state vector component 1 (dt = 1e-05): 273.708 > 20
Limiting global state vector component 1 (dt = 5e-06): 78.6964 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 31.3396 > 20
3.050e-04 1109.250 101325.000 581179.574385
Limiting global state vector component 1 (dt = 1e-05): 358.491 > 20
Limiting global state vector component 1 (dt = 5e-06): 98.1841 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 38.326 > 20
3.063e-04 1126.431 101325.000 574936.977616
Limiting global state vector component 1 (dt = 1e-05): 483.183 > 20
Limiting global state vector component 1 (dt = 5e-06): 124.957 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 47.6408 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 21.1544 > 20
3.069e-04 1136.438 101325.000 571306.548424
Limiting global state vector component 1 (dt = 1e-05): 562.674 > 20
Limiting global state vector component 1 (dt = 5e-06): 141.945 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 53.4718 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 23.6213 > 20
3.075e-04 1147.586 101325.000 567266.817260
Limiting global state vector component 1 (dt = 1e-05): 655.304 > 20
Limiting global state vector component 1 (dt = 5e-06): 161.941 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 60.2979 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 26.4955 > 20
3.081e-04 1160.060 101325.000 562752.407328
Limiting global state vector component 1 (dt = 1e-05): 759.775 > 20
Limiting global state vector component 1 (dt = 5e-06): 185.376 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 68.3113 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 29.859 > 20
3.088e-04 1174.082 101325.000 557685.008242
Limiting global state vector component 1 (dt = 1e-05): 868.917 > 20
Limiting global state vector component 1 (dt = 5e-06): 212.465 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 77.7189 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 33.8072 > 20
3.094e-04 1189.919 101325.000 551971.071480
Limiting global state vector component 1 (dt = 1e-05): 955.051 > 20
Limiting global state vector component 1 (dt = 5e-06): 242.578 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 88.7024 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 38.4451 > 20
3.100e-04 1207.888 101325.000 545499.819076
Limiting global state vector component 1 (dt = 1e-05): 976.965 > 20
Limiting global state vector component 1 (dt = 5e-06): 274.163 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 101.371 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 43.8759 > 20
Limiting global state vector component 1 (dt = 6.25e-07): 20.4741 > 20
3.103e-04 1217.787 101325.000 541940.650602
Limiting global state vector component 1 (dt = 1e-05): 948.852 > 20
Limiting global state vector component 1 (dt = 5e-06): 289.678 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 108.325 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 46.914 > 20
Limiting global state vector component 1 (dt = 6.25e-07): 21.8816 > 20
3.106e-04 1228.362 101325.000 538142.580923
Limiting global state vector component 1 (dt = 1e-05): 874.586 > 20
Limiting global state vector component 1 (dt = 5e-06): 303.948 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 115.633 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 50.1723 > 20
Limiting global state vector component 1 (dt = 6.25e-07): 23.3988 > 20
3.109e-04 1239.668 101325.000 534087.175073
Limiting global state vector component 1 (dt = 1e-05): 676.594 > 20
Limiting global state vector component 1 (dt = 5e-06): 313.817 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 123.128 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 53.6414 > 20
Limiting global state vector component 1 (dt = 6.25e-07): 25.0275 > 20
3.113e-04 1251.761 101325.000 529755.591823
Limiting global state vector component 1 (dt = 1e-05): 433.652 > 20
Limiting global state vector component 1 (dt = 5e-06): 321.25 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 130.795 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 57.309 > 20
Limiting global state vector component 1 (dt = 6.25e-07): 26.7657 > 20
3.116e-04 1264.696 101325.000 525129.174276
Limiting global state vector component 1 (dt = 1e-05): 79.251 > 20
Limiting global state vector component 1 (dt = 5e-06): 322.993 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 138.373 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 61.1381 > 20
Limiting global state vector component 1 (dt = 6.25e-07): 28.6055 > 20
3.119e-04 1278.526 101325.000 520190.330365
Limiting global state vector component 1 (dt = 1e-05): 568.175 > 20
Limiting global state vector component 1 (dt = 5e-06): 311.509 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 145.376 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 65.0625 > 20
Limiting global state vector component 1 (dt = 6.25e-07): 30.5305 > 20
3.122e-04 1293.300 101325.000 514923.786027
Limiting global state vector component 1 (dt = 1e-05): 1370.79 > 20
Limiting global state vector component 1 (dt = 5e-06): 290.015 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 151.645 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 69.0062 > 20
Limiting global state vector component 1 (dt = 6.25e-07): 32.5145 > 20
3.125e-04 1309.056 101325.000 509318.289186
Limiting global state vector component 1 (dt = 1e-05): 2132.6 > 20
Limiting global state vector component 1 (dt = 5e-06): 263.044 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 156.959 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 72.8582 > 20
Limiting global state vector component 1 (dt = 6.25e-07): 34.517 > 20
3.128e-04 1325.815 101325.000 503368.815126
Limiting global state vector component 1 (dt = 1e-05): 2977.41 > 20
Limiting global state vector component 1 (dt = 5e-06): 225.067 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 160.707 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 76.4612 > 20
Limiting global state vector component 1 (dt = 6.25e-07): 36.481 > 20
3.131e-04 1343.573 101325.000 497079.250399
Limiting global state vector component 1 (dt = 1e-05): 3494.06 > 20
Limiting global state vector component 1 (dt = 5e-06): 188.7 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 162.879 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 79.6488 > 20
Limiting global state vector component 1 (dt = 6.25e-07): 38.3331 > 20
3.134e-04 1362.296 101325.000 490465.396978
Limiting global state vector component 1 (dt = 1e-05): 3625.91 > 20
Limiting global state vector component 1 (dt = 5e-06): 156.511 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 163.198 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 82.2247 > 20
Limiting global state vector component 1 (dt = 6.25e-07): 39.9834 > 20
3.137e-04 1381.906 101325.000 483557.927461
Limiting global state vector component 1 (dt = 1e-05): 2694.11 > 20
Limiting global state vector component 1 (dt = 5e-06): 153.544 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 162.327 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 84.0278 > 20
Limiting global state vector component 1 (dt = 6.25e-07): 41.3335 > 20
Limiting global state vector component 1 (dt = 3.125e-07): 20.3737 > 20
3.139e-04 1392.008 101325.000 480008.325066
Limiting global state vector component 1 (dt = 1e-05): 1997.51 > 20
Limiting global state vector component 1 (dt = 5e-06): 159.116 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 161.33 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 84.5847 > 20
Limiting global state vector component 1 (dt = 6.25e-07): 41.8644 > 20
Limiting global state vector component 1 (dt = 3.125e-07): 20.6925 > 20
3.141e-04 1402.280 101325.000 476404.658364
Limiting global state vector component 1 (dt = 1e-05): 689.632 > 20
Limiting global state vector component 1 (dt = 5e-06): 185.461 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 160.559 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 84.9142 > 20
Limiting global state vector component 1 (dt = 6.25e-07): 42.285 > 20
Limiting global state vector component 1 (dt = 3.125e-07): 20.961 > 20
3.142e-04 1412.700 101325.000 472755.704139
Limiting global state vector component 1 (dt = 1e-05): 496.192 > 20
Limiting global state vector component 1 (dt = 5e-06): 209.53 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 159.29 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 84.9806 > 20
Limiting global state vector component 1 (dt = 6.25e-07): 42.5844 > 20
Limiting global state vector component 1 (dt = 3.125e-07): 21.1736 > 20
3.144e-04 1423.241 101325.000 469071.270992
Limiting global state vector component 1 (dt = 1e-05): 1837.57 > 20
Limiting global state vector component 1 (dt = 5e-06): 241.227 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 157.937 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 84.796 > 20
Limiting global state vector component 1 (dt = 6.25e-07): 42.7559 > 20
Limiting global state vector component 1 (dt = 3.125e-07): 21.3254 > 20
3.145e-04 1433.874 101325.000 465362.073351
Limiting global state vector component 1 (dt = 1e-05): 3674.62 > 20
Limiting global state vector component 1 (dt = 5e-06): 292.517 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 157.005 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 84.3828 > 20
Limiting global state vector component 1 (dt = 6.25e-07): 42.7955 > 20
Limiting global state vector component 1 (dt = 3.125e-07): 21.4126 > 20
3.147e-04 1444.567 101325.000 461639.553517
Limiting global state vector component 1 (dt = 1e-05): 5083.88 > 20
Limiting global state vector component 1 (dt = 5e-06): 333.712 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 155.57 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 83.7164 > 20
Limiting global state vector component 1 (dt = 6.25e-07): 42.6999 > 20
Limiting global state vector component 1 (dt = 3.125e-07): 21.4325 > 20
3.148e-04 1455.286 101325.000 457915.654531
Limiting global state vector component 1 (dt = 1e-05): 6363.31 > 20
Limiting global state vector component 1 (dt = 5e-06): 374.755 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 154.057 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 82.8278 > 20
Limiting global state vector component 1 (dt = 6.25e-07): 42.4714 > 20
Limiting global state vector component 1 (dt = 3.125e-07): 21.3842 > 20
3.150e-04 1465.999 101325.000 454202.551386
Limiting global state vector component 1 (dt = 1e-05): 7560.48 > 20
Limiting global state vector component 1 (dt = 5e-06): 417.059 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 152.587 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 81.7421 > 20
Limiting global state vector component 1 (dt = 6.25e-07): 42.115 > 20
Limiting global state vector component 1 (dt = 3.125e-07): 21.2685 > 20
3.152e-04 1476.671 101325.000 450512.352795
Limiting global state vector component 1 (dt = 1e-05): 8054.65 > 20
Limiting global state vector component 1 (dt = 5e-06): 439.927 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 150.514 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 80.4602 > 20
Limiting global state vector component 1 (dt = 6.25e-07): 41.638 > 20
Limiting global state vector component 1 (dt = 3.125e-07): 21.0879 > 20
3.153e-04 1487.268 101325.000 446856.789926
Limiting global state vector component 1 (dt = 1e-05): 7951.07 > 20
Limiting global state vector component 1 (dt = 5e-06): 445.86 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 147.947 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 79.0127 > 20
Limiting global state vector component 1 (dt = 6.25e-07): 41.0514 > 20
Limiting global state vector component 1 (dt = 3.125e-07): 20.8463 > 20
3.155e-04 1497.759 101325.000 443246.911117
Limiting global state vector component 1 (dt = 1e-05): 7206.81 > 20
Limiting global state vector component 1 (dt = 5e-06): 432.498 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 144.821 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 77.4256 > 20
Limiting global state vector component 1 (dt = 6.25e-07): 40.3684 > 20
Limiting global state vector component 1 (dt = 3.125e-07): 20.5494 > 20
3.156e-04 1508.114 101325.000 439692.802431
Limiting global state vector component 1 (dt = 1e-05): 6434.15 > 20
Limiting global state vector component 1 (dt = 5e-06): 418.039 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 141.705 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 75.7452 > 20
Limiting global state vector component 1 (dt = 6.25e-07): 39.6046 > 20
Limiting global state vector component 1 (dt = 3.125e-07): 20.204 > 20
3.158e-04 1518.308 101325.000 436203.352407
Limiting global state vector component 1 (dt = 1e-05): 4890.76 > 20
Limiting global state vector component 1 (dt = 5e-06): 377.938 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 137.793 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 73.972 > 20
Limiting global state vector component 1 (dt = 6.25e-07): 38.7748 > 20
3.161e-04 1538.126 101325.000 429447.004732
Limiting global state vector component 1 (dt = 1e-05): 1981.88 > 20
Limiting global state vector component 1 (dt = 5e-06): 296.68 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 129.936 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 70.3225 > 20
Limiting global state vector component 1 (dt = 6.25e-07): 36.9842 > 20
3.164e-04 1557.081 101325.000 423020.045368
Limiting global state vector component 1 (dt = 1e-05): 139.111 > 20
Limiting global state vector component 1 (dt = 5e-06): 229.733 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 122.485 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 66.6897 > 20
Limiting global state vector component 1 (dt = 6.25e-07): 35.1177 > 20
3.167e-04 1575.108 101325.000 416941.309538
Limiting global state vector component 1 (dt = 1e-05): 2379.3 > 20
Limiting global state vector component 1 (dt = 5e-06): 148.742 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 114.5 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 63.1631 > 20
Limiting global state vector component 1 (dt = 6.25e-07): 33.2713 > 20
3.170e-04 1592.197 101325.000 411209.708408
Limiting global state vector component 1 (dt = 1e-05): 3068.51 > 20
Limiting global state vector component 1 (dt = 5e-06): 111.69 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 108.06 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 59.8958 > 20
Limiting global state vector component 1 (dt = 6.25e-07): 31.5154 > 20
3.173e-04 1608.379 101325.000 405809.689134
Limiting global state vector component 1 (dt = 1e-05): 2876.34 > 20
Limiting global state vector component 1 (dt = 5e-06): 101.212 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 102.753 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 56.9132 > 20
Limiting global state vector component 1 (dt = 6.25e-07): 29.8881 > 20
3.177e-04 1623.713 101325.000 400717.041191
Limiting global state vector component 1 (dt = 1e-05): 2066.91 > 20
Limiting global state vector component 1 (dt = 5e-06): 112.018 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 98.5413 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 54.2243 > 20
Limiting global state vector component 1 (dt = 6.25e-07): 28.4047 > 20
3.180e-04 1638.268 101325.000 395903.831012
Limiting global state vector component 1 (dt = 1e-05): 1150.25 > 20
Limiting global state vector component 1 (dt = 5e-06): 128.955 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 94.9639 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 51.8015 > 20
Limiting global state vector component 1 (dt = 6.25e-07): 27.0634 > 20
3.183e-04 1652.118 101325.000 391341.919323
Limiting global state vector component 1 (dt = 1e-05): 584.018 > 20
Limiting global state vector component 1 (dt = 5e-06): 137.293 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 91.5116 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 49.6042 > 20
Limiting global state vector component 1 (dt = 6.25e-07): 25.8533 > 20
3.186e-04 1665.332 101325.000 387005.062694
Limiting global state vector component 1 (dt = 1e-05): 68.3256 > 20
Limiting global state vector component 1 (dt = 5e-06): 146.447 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 88.4662 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 47.6167 > 20
Limiting global state vector component 1 (dt = 6.25e-07): 24.7604 > 20
3.189e-04 1677.972 101325.000 382869.897227
Limiting global state vector component 1 (dt = 1e-05): 221.728 > 20
Limiting global state vector component 1 (dt = 5e-06): 149.973 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 85.5472 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 45.8021 > 20
Limiting global state vector component 1 (dt = 6.25e-07): 23.769 > 20
3.192e-04 1690.093 101325.000 378916.179707
Limiting global state vector component 1 (dt = 1e-05): 339.274 > 20
Limiting global state vector component 1 (dt = 5e-06): 149.019 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 82.7403 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 44.1362 > 20
Limiting global state vector component 1 (dt = 6.25e-07): 22.8651 > 20
3.195e-04 1701.741 101325.000 375126.611311
Limiting global state vector component 1 (dt = 1e-05): 390.121 > 20
Limiting global state vector component 1 (dt = 5e-06): 146.609 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 80.1037 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 42.6008 > 20
Limiting global state vector component 1 (dt = 6.25e-07): 22.0362 > 20
3.198e-04 1712.958 101325.000 371486.474041
Limiting global state vector component 1 (dt = 1e-05): 393.26 > 20
Limiting global state vector component 1 (dt = 5e-06): 143.076 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 77.6133 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 41.1782 > 20
Limiting global state vector component 1 (dt = 6.25e-07): 21.2721 > 20
3.202e-04 1723.778 101325.000 367983.219109
Limiting global state vector component 1 (dt = 1e-05): 383.099 > 20
Limiting global state vector component 1 (dt = 5e-06): 139.398 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 75.277 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 39.8548 > 20
Limiting global state vector component 1 (dt = 6.25e-07): 20.5639 > 20
3.205e-04 1734.230 101325.000 364606.079827
Limiting global state vector component 1 (dt = 1e-05): 352.141 > 20
Limiting global state vector component 1 (dt = 5e-06): 135.207 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 73.0598 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 38.618 > 20
3.211e-04 1754.135 101325.000 358194.054593
Limiting global state vector component 1 (dt = 1e-05): 279.124 > 20
Limiting global state vector component 1 (dt = 5e-06): 126.868 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 68.9815 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 36.3686 > 20
3.217e-04 1772.846 101325.000 352188.724470
Limiting global state vector component 1 (dt = 1e-05): 249.075 > 20
Limiting global state vector component 1 (dt = 5e-06): 120.393 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 65.3678 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 34.3725 > 20
3.223e-04 1790.503 101325.000 346541.357521
Limiting global state vector component 1 (dt = 1e-05): 173.832 > 20
Limiting global state vector component 1 (dt = 5e-06): 112.884 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 62.0563 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 32.5828 > 20
3.230e-04 1807.218 101325.000 341212.266130
Limiting global state vector component 1 (dt = 1e-05): 211.092 > 20
Limiting global state vector component 1 (dt = 5e-06): 109.406 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 59.1666 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 30.9729 > 20
3.236e-04 1823.087 101325.000 336168.391749
Limiting global state vector component 1 (dt = 1e-05): 192.408 > 20
Limiting global state vector component 1 (dt = 5e-06): 104.463 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 56.474 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 29.5106 > 20
3.242e-04 1838.191 101325.000 331381.744877
Limiting global state vector component 1 (dt = 1e-05): 180.84 > 20
Limiting global state vector component 1 (dt = 5e-06): 100.1 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 54.0177 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 28.1782 > 20
3.248e-04 1852.598 101325.000 326828.341096
Limiting global state vector component 1 (dt = 1e-05): 124.51 > 20
Limiting global state vector component 1 (dt = 5e-06): 94.6242 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 51.715 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 26.9571 > 20
3.255e-04 1866.369 101325.000 322487.435545
Limiting global state vector component 1 (dt = 1e-05): 156.223 > 20
Limiting global state vector component 1 (dt = 5e-06): 92.1642 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 49.673 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 25.8373 > 20
3.261e-04 1879.556 101325.000 318340.948281
Limiting global state vector component 1 (dt = 1e-05): 152.441 > 20
Limiting global state vector component 1 (dt = 5e-06): 88.8742 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 47.7557 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 24.8039 > 20
3.267e-04 1892.207 101325.000 314373.019077
Limiting global state vector component 1 (dt = 1e-05): 148.181 > 20
Limiting global state vector component 1 (dt = 5e-06): 85.7662 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 45.9742 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 23.8475 > 20
3.273e-04 1904.361 101325.000 310569.654637
Limiting global state vector component 1 (dt = 1e-05): 144.112 > 20
Limiting global state vector component 1 (dt = 5e-06): 82.861 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 44.3166 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 22.9598 > 20
3.280e-04 1916.054 101325.000 306918.444880
Limiting global state vector component 1 (dt = 1e-05): 140.37 > 20
Limiting global state vector component 1 (dt = 5e-06): 80.1497 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 42.7713 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 22.1338 > 20
3.286e-04 1927.321 101325.000 303408.331597
Limiting global state vector component 1 (dt = 1e-05): 136.721 > 20
Limiting global state vector component 1 (dt = 5e-06): 77.5996 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 41.3259 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 21.3631 > 20
3.292e-04 1938.188 101325.000 300029.418306
Limiting global state vector component 1 (dt = 1e-05): 133.007 > 20
Limiting global state vector component 1 (dt = 5e-06): 75.1856 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 39.9704 > 20
Limiting global state vector component 1 (dt = 1.25e-06): 20.6422 > 20
3.298e-04 1948.684 101325.000 296772.812500
Limiting global state vector component 1 (dt = 1e-05): 129.349 > 20
Limiting global state vector component 1 (dt = 5e-06): 72.9027 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 38.6968 > 20
3.311e-04 1968.651 101325.000 290595.204396
Limiting global state vector component 1 (dt = 1e-05): 123.068 > 20
Limiting global state vector component 1 (dt = 5e-06): 68.7472 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 36.3713 > 20
3.323e-04 1987.385 101325.000 284819.936972
Limiting global state vector component 1 (dt = 1e-05): 117.164 > 20
Limiting global state vector component 1 (dt = 5e-06): 65.0091 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 34.2969 > 20
3.336e-04 2005.025 101325.000 279400.944885
Limiting global state vector component 1 (dt = 1e-05): 111.696 > 20
Limiting global state vector component 1 (dt = 5e-06): 61.6332 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 32.4351 > 20
3.348e-04 2021.684 101325.000 274299.614076
Limiting global state vector component 1 (dt = 1e-05): 106.703 > 20
Limiting global state vector component 1 (dt = 5e-06): 58.5742 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 30.755 > 20
3.361e-04 2037.461 101325.000 269483.255683
Limiting global state vector component 1 (dt = 1e-05): 101.74 > 20
Limiting global state vector component 1 (dt = 5e-06): 55.7629 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 29.2291 > 20
3.373e-04 2052.440 101325.000 264923.950872
Limiting global state vector component 1 (dt = 1e-05): 97.7073 > 20
Limiting global state vector component 1 (dt = 5e-06): 53.229 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 27.8415 > 20
3.386e-04 2066.693 101325.000 260597.667274
Limiting global state vector component 1 (dt = 1e-05): 93.7912 > 20
Limiting global state vector component 1 (dt = 5e-06): 50.8893 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 26.5707 > 20
3.398e-04 2080.282 101325.000 256483.568908
Limiting global state vector component 1 (dt = 1e-05): 90.5624 > 20
Limiting global state vector component 1 (dt = 5e-06): 48.7605 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 25.4053 > 20
3.411e-04 2093.263 101325.000 252563.473239
Limiting global state vector component 1 (dt = 1e-05): 86.5508 > 20
Limiting global state vector component 1 (dt = 5e-06): 46.7259 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 24.3265 > 20
3.423e-04 2105.685 101325.000 248821.420321
Limiting global state vector component 1 (dt = 1e-05): 83.3533 > 20
Limiting global state vector component 1 (dt = 5e-06): 44.8751 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 23.3313 > 20
3.436e-04 2117.590 101325.000 245243.324811
Limiting global state vector component 1 (dt = 1e-05): 80.3407 > 20
Limiting global state vector component 1 (dt = 5e-06): 43.1527 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 22.4082 > 20
3.448e-04 2129.017 101325.000 241816.693903
Limiting global state vector component 1 (dt = 1e-05): 77.5088 > 20
Limiting global state vector component 1 (dt = 5e-06): 41.5463 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 21.5497 > 20
3.461e-04 2139.999 101325.000 238530.396534
Limiting global state vector component 1 (dt = 1e-05): 74.86 > 20
Limiting global state vector component 1 (dt = 5e-06): 40.0459 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 20.7493 > 20
3.473e-04 2150.567 101325.000 235374.472866
Limiting global state vector component 1 (dt = 1e-05): 72.3348 > 20
Limiting global state vector component 1 (dt = 5e-06): 38.6384 > 20
Limiting global state vector component 1 (dt = 2.5e-06): 20.001 > 20
3.486e-04 2160.748 101325.000 232339.976033
Limiting global state vector component 1 (dt = 1e-05): 69.9999 > 20
Limiting global state vector component 1 (dt = 5e-06): 37.3206 > 20
3.511e-04 2180.049 101325.000 226603.766845
Limiting global state vector component 1 (dt = 1e-05): 65.7044 > 20
Limiting global state vector component 1 (dt = 5e-06): 34.9134 > 20
3.536e-04 2198.073 101325.000 221265.915155
Limiting global state vector component 1 (dt = 1e-05): 61.8452 > 20
Limiting global state vector component 1 (dt = 5e-06): 32.7693 > 20
3.561e-04 2214.965 101325.000 216280.170135
Limiting global state vector component 1 (dt = 1e-05): 58.3649 > 20
Limiting global state vector component 1 (dt = 5e-06): 30.8479 > 20
3.586e-04 2230.844 101325.000 211607.810734
Limiting global state vector component 1 (dt = 1e-05): 55.2081 > 20
Limiting global state vector component 1 (dt = 5e-06): 29.1159 > 20
3.611e-04 2245.814 101325.000 207216.102714
Limiting global state vector component 1 (dt = 1e-05): 52.3389 > 20
Limiting global state vector component 1 (dt = 5e-06): 27.5472 > 20
3.636e-04 2259.961 101325.000 203077.131709
Limiting global state vector component 1 (dt = 1e-05): 49.718 > 20
Limiting global state vector component 1 (dt = 5e-06): 26.1196 > 20
3.661e-04 2273.362 101325.000 199166.906588
Limiting global state vector component 1 (dt = 1e-05): 47.3141 > 20
Limiting global state vector component 1 (dt = 5e-06): 24.815 > 20
3.686e-04 2286.082 101325.000 195464.661940
Limiting global state vector component 1 (dt = 1e-05): 45.106 > 20
Limiting global state vector component 1 (dt = 5e-06): 23.6183 > 20
3.711e-04 2298.178 101325.000 191952.308444
Limiting global state vector component 1 (dt = 1e-05): 43.0627 > 20
Limiting global state vector component 1 (dt = 5e-06): 22.5161 > 20
3.736e-04 2309.701 101325.000 188613.995897
Limiting global state vector component 1 (dt = 1e-05): 41.2274 > 20
Limiting global state vector component 1 (dt = 5e-06): 21.5022 > 20
3.761e-04 2320.694 101325.000 185435.760206
Limiting global state vector component 1 (dt = 1e-05): 39.4095 > 20
Limiting global state vector component 1 (dt = 5e-06): 20.5544 > 20
3.786e-04 2331.199 101325.000 182405.238946
Limiting global state vector component 1 (dt = 1e-05): 37.7702 > 20
3.836e-04 2350.877 101325.000 176744.535443
Limiting global state vector component 1 (dt = 1e-05): 34.8078 > 20
3.886e-04 2368.974 101325.000 171557.088556
Limiting global state vector component 1 (dt = 1e-05): 32.2038 > 20
3.936e-04 2385.688 101325.000 166782.363360
Limiting global state vector component 1 (dt = 1e-05): 29.8887 > 20
3.986e-04 2401.178 101325.000 162370.702984
Limiting global state vector component 1 (dt = 1e-05): 27.8273 > 20
4.036e-04 2415.579 101325.000 158280.887128
Limiting global state vector component 1 (dt = 1e-05): 25.9709 > 20
4.086e-04 2429.004 101325.000 154478.332291
Limiting global state vector component 1 (dt = 1e-05): 24.2936 > 20
4.136e-04 2441.550 101325.000 150933.746275
Limiting global state vector component 1 (dt = 1e-05): 22.7777 > 20
4.186e-04 2453.298 101325.000 147622.101947
Limiting global state vector component 1 (dt = 1e-05): 21.3793 > 20
4.236e-04 2464.320 101325.000 144521.847342
Limiting global state vector component 1 (dt = 1e-05): 20.1057 > 20
4.286e-04 2474.678 101325.000 141614.288549
4.386e-04 2493.612 101325.000 136313.940209
4.486e-04 2510.465 101325.000 131612.404198
4.586e-04 2525.524 101325.000 127424.004429
4.686e-04 2539.022 101325.000 123679.831541
4.786e-04 2551.152 101325.000 120323.487802
4.886e-04 2562.073 101325.000 117308.102855
4.986e-04 2571.921 101325.000 114594.197213
5.086e-04 2580.813 101325.000 112148.124270
5.186e-04 2588.850 101325.000 109940.914913
5.286e-04 2596.119 101325.000 107947.407740
5.386e-04 2602.698 101325.000 106145.585210
5.486e-04 2608.654 101325.000 104516.057722
5.586e-04 2614.049 101325.000 103041.659586
5.686e-04 2618.938 101325.000 101707.126205
5.786e-04 2623.367 101325.000 100498.834004
5.886e-04 2627.382 101325.000 99404.586942
5.986e-04 2631.021 101325.000 98413.440049
6.086e-04 2634.320 101325.000 97515.551241
6.186e-04 2637.311 101325.000 96702.055729
6.286e-04 2640.022 101325.000 95964.958463
6.386e-04 2642.480 101325.000 95297.041281
6.486e-04 2644.708 101325.000 94691.781381
6.586e-04 2646.728 101325.000 94143.280774
6.686e-04 2648.559 101325.000 93646.203501
6.786e-04 2650.220 101325.000 93195.720255
6.886e-04 2651.725 101325.000 92787.458286
6.986e-04 2653.089 101325.000 92417.457469
7.086e-04 2654.326 101325.000 92082.130605
7.186e-04 2655.447 101325.000 91778.227757
7.286e-04 2656.464 101325.000 91502.804247
7.386e-04 2657.385 101325.000 91253.191471
7.486e-04 2658.220 101325.000 91026.971601
7.586e-04 2658.977 101325.000 90821.953208
7.686e-04 2659.663 101325.000 90636.149790
7.786e-04 2660.285 101325.000 90467.761237
7.886e-04 2660.848 101325.000 90315.156043
7.986e-04 2661.359 101325.000 90176.855399
8.086e-04 2661.822 101325.000 90051.519048
8.186e-04 2662.242 101325.000 89937.932221
8.286e-04 2662.622 101325.000 89834.993908
8.386e-04 2662.967 101325.000 89741.706210
8.486e-04 2663.279 101325.000 89657.164694
8.586e-04 2663.562 101325.000 89580.549644
8.686e-04 2663.819 101325.000 89511.118128
8.786e-04 2664.051 101325.000 89448.196843
8.886e-04 2664.262 101325.000 89391.175542
8.986e-04 2664.453 101325.000 89339.501081
9.086e-04 2664.626 101325.000 89292.672173
9.186e-04 2664.783 101325.000 89250.234573
9.286e-04 2664.925 101325.000 89211.776602
9.386e-04 2665.054 101325.000 89176.925125
9.486e-04 2665.171 101325.000 89145.341981
9.586e-04 2665.277 101325.000 89116.720719
9.686e-04 2665.373 101325.000 89090.783616
9.786e-04 2665.459 101325.000 89067.278979
9.886e-04 2665.538 101325.000 89045.978710
9.986e-04 2665.610 101325.000 89026.676094
1.009e-03 2665.674 101325.000 89009.183800
import sys
import cantera as ct
gas = ct.Solution('h2o2.yaml')
gas.TPX = 1001.0, ct.one_atm, 'H2:2,O2:1,N2:4'
r = ct.IdealGasConstPressureReactor(gas)
sim = ct.ReactorNet([r])
sim.verbose = True
# limit advance when temperature difference is exceeded
delta_T_max = 20.
r.set_advance_limit('temperature', delta_T_max)
dt_max = 1.e-5
t_end = 100 * dt_max
states = ct.SolutionArray(gas, extra=['t'])
print('{:10s} {:10s} {:10s} {:14s}'.format(
't [s]', 'T [K]', 'P [Pa]', 'u [J/kg]'))
while sim.time < t_end:
sim.advance(sim.time + dt_max)
states.append(r.thermo.state, t=sim.time*1e3)
print('{:10.3e} {:10.3f} {:10.3f} {:14.6f}'.format(
sim.time, r.T, r.thermo.P, r.thermo.u))
# Plot the results if matplotlib is installed.
# See http://matplotlib.org/ to get it.
if '--plot' in sys.argv[1:]:
import matplotlib.pyplot as plt
plt.clf()
plt.subplot(2, 2, 1)
plt.plot(states.t, states.T)
plt.xlabel('Time (ms)')
plt.ylabel('Temperature (K)')
plt.subplot(2, 2, 2)
plt.plot(states.t, states.X[:, gas.species_index('OH')])
plt.xlabel('Time (ms)')
plt.ylabel('OH Mole Fraction')
plt.subplot(2, 2, 3)
plt.plot(states.t, states.X[:, gas.species_index('H')])
plt.xlabel('Time (ms)')
plt.ylabel('H Mole Fraction')
plt.subplot(2, 2, 4)
plt.plot(states.t, states.X[:, gas.species_index('H2')])
plt.xlabel('Time (ms)')
plt.ylabel('H2 Mole Fraction')
plt.tight_layout()
plt.show()
else:
print("To view a plot of these results, run this script with the option --plot")
Total running time of the script: (0 minutes 0.216 seconds)