Note
Go to the end to download the full example code
Opposed-flow premixed strained flame#
This script simulates a lean hydrogen-oxygen flame stabilized in a strained flowfield, with an opposed flow consisting of equilibrium products.
Requires: cantera >= 3.0
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> reactants <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
Mass Flux: 0.12 kg/m^2/s
Temperature: 373 K
Mass Fractions:
H2 0.01024
O2 0.1016
AR 0.8881
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> flame <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
Pressure: 5066 Pa
-------------------------------------------------------------------------------
z velocity spread_rate T lambda eField
-------------------------------------------------------------------------------
0 2.24 0 373 -50.19 0
0.06 -0.1336 17.8 373 -50.19 0
0.1 -1.716 29.67 1214 -50.19 0
0.14 -3.298 35.6 2055 -50.19 0
0.2 -5.671 0 2055 -50.19 0
-------------------------------------------------------------------------------
z H2 H O O2 OH
-------------------------------------------------------------------------------
0 0.01024 0 0 0.1016 0
0.06 0.01024 0 0 0.1016 0
0.1 0.005186 7.052e-06 0.0001515 0.06091 0.001018
0.14 0.0001284 1.41e-05 0.0003029 0.02019 0.002035
0.2 0.0001284 1.41e-05 0.0003029 0.02019 0.002035
-------------------------------------------------------------------------------
z H2O HO2 H2O2 AR N2
-------------------------------------------------------------------------------
0 0 0 0 0.8881 0
0.06 0 0 0 0.8881 0
0.1 0.04459 1.92e-07 6.302e-09 0.8881 0
0.14 0.08919 3.841e-07 1.26e-08 0.8881 0
0.2 0.08919 3.841e-07 1.26e-08 0.8881 0
-------------------------------------------------------------------------------
z
-------------------------------------------------------------------------------
0
0.06
0.1
0.14
0.2
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> products <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
Mass Flux: 0.06 kg/m^2/s
Temperature: 2055 K
Mass Fractions:
H2 0.0001284
H 1.41e-05
O 0.0003029
O2 0.02019
OH 0.002035
H2O 0.08919
HO2 3.841e-07
H2O2 1.26e-08
AR 0.8881
************ Solving on 5 point grid with energy equation enabled ************
..............................................................................
Attempt Newton solution of steady-state problem... failure.
Take 10 timesteps 1.602e-05 4.599
Attempt Newton solution of steady-state problem... failure.
Take 10 timesteps 0.0001825 4.948
Attempt Newton solution of steady-state problem... failure.
Take 10 timesteps 0.004676 3.846
Attempt Newton solution of steady-state problem... failure.
Take 10 timesteps 0.001873 3.949
Attempt Newton solution of steady-state problem... failure.
Take 10 timesteps 0.003999 3.885
Attempt Newton solution of steady-state problem... success.
Problem solved on [5] point grid(s).
..............................................................................
grid refinement disabled.
******************** Solving with grid refinement enabled ********************
..............................................................................
Attempt Newton solution of steady-state problem... success.
Problem solved on [5] point grid(s).
..............................................................................
##############################################################################
Refining grid in flame.
New points inserted after grid points 0 1 2 3
to resolve AR H H2 H2O H2O2 HO2 O O2 OH T spread_rate velocity
##############################################################################
..............................................................................
Attempt Newton solution of steady-state problem... success.
Problem solved on [9] point grid(s).
..............................................................................
##############################################################################
Refining grid in flame.
New points inserted after grid points 0 1 2 3 4 5 6 7
to resolve AR H H2 H2O H2O2 HO2 O O2 OH T spread_rate velocity
##############################################################################
..............................................................................
Attempt Newton solution of steady-state problem... failure.
Take 10 timesteps 0.0002563 4.894
Attempt Newton solution of steady-state problem... failure.
Take 10 timesteps 0.004379 3.636
Attempt Newton solution of steady-state problem... failure.
Take 10 timesteps 0.004676 3.879
Attempt Newton solution of steady-state problem... success.
Problem solved on [17] point grid(s).
..............................................................................
##############################################################################
Refining grid in flame.
New points inserted after grid points 0 1 2 3 8 9 10 11 12 13 14 15
to resolve AR H H2 H2O H2O2 HO2 O O2 OH T spread_rate velocity
##############################################################################
..............................................................................
Attempt Newton solution of steady-state problem... failure.
Take 10 timesteps 0.0002563 4.293
Attempt Newton solution of steady-state problem... failure.
Take 10 timesteps 0.006568 3.85
Attempt Newton solution of steady-state problem... success.
Problem solved on [29] point grid(s).
..............................................................................
##############################################################################
Refining grid in flame.
New points inserted after grid points 0 1 2 3 4 5 6 15 16 17 18 20 21 22 23 24 25 26 27
to resolve AR H H2 H2O H2O2 HO2 O O2 OH T spread_rate velocity
##############################################################################
..............................................................................
Attempt Newton solution of steady-state problem... failure.
Take 10 timesteps 0.0003844 4.273
Attempt Newton solution of steady-state problem... failure.
Take 10 timesteps 0.006568 4.135
Attempt Newton solution of steady-state problem... failure.
Take 10 timesteps 0.07482 1.681
Attempt Newton solution of steady-state problem... success.
Problem solved on [48] point grid(s).
..............................................................................
##############################################################################
Refining grid in flame.
New points inserted after grid points 5 6 7 8 9 10 11 12 13 14 15 16 17 18
to resolve AR H H2 H2O H2O2 HO2 O O2 OH T spread_rate velocity
##############################################################################
..............................................................................
Attempt Newton solution of steady-state problem... failure.
Take 10 timesteps 0.0003844 4.082
Attempt Newton solution of steady-state problem... success.
Problem solved on [62] point grid(s).
..............................................................................
##############################################################################
Refining grid in flame.
New points inserted after grid points 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 35 36
to resolve AR H H2 H2O H2O2 HO2 O O2 OH T spread_rate velocity
##############################################################################
..............................................................................
Attempt Newton solution of steady-state problem... success.
Problem solved on [79] point grid(s).
..............................................................................
##############################################################################
Refining grid in flame.
New points inserted after grid points 18 19 20 21 22 23 24 25 26 27 28 29 30 31 33 34 35 36 37 38
to resolve AR H H2 H2O H2O2 HO2 O O2 OH
##############################################################################
..............................................................................
Attempt Newton solution of steady-state problem... success.
Problem solved on [99] point grid(s).
..............................................................................
##############################################################################
Refining grid in flame.
New points inserted after grid points 25 26 27 28 35 42 43 44 45 46 50 51 52 53 54 55
to resolve AR H H2O2 HO2 O
##############################################################################
..............................................................................
Attempt Newton solution of steady-state problem... success.
Problem solved on [115] point grid(s).
..............................................................................
##############################################################################
Refining grid in flame.
New points inserted after grid points 33 41 55 78
to resolve H H2O2 HO2
##############################################################################
..............................................................................
Attempt Newton solution of steady-state problem... success.
Problem solved on [119] point grid(s).
..............................................................................
##############################################################################
Refining grid in flame.
New points inserted after grid points 44
to resolve H2O2
##############################################################################
..............................................................................
Attempt Newton solution of steady-state problem... success.
Problem solved on [120] point grid(s).
..............................................................................
no new points needed in flame
Statistics:
Grid Timesteps Functions Time Jacobians Time
7 0 1161 0.0065 44 0.0106
11 0 139 0.0014 3 0.0015
19 30 823 0.0154 26 0.0256
31 20 285 0.0090 10 0.0173
50 30 441 0.0233 13 0.0386
64 10 137 0.0095 7 0.0272
81 0 33 0.0030 2 0.0104
101 0 32 0.0036 1 0.0065
117 0 26 0.0033 1 0.0076
121 0 18 0.0024 1 0.0081
122 0 14 0.0020 1 0.0081
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> reactants <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
Mass Flux: 0.12 kg/m^2/s
Temperature: 373 K
Mass Fractions:
H2 0.01024
O2 0.1016
AR 0.8881
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> flame <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
Pressure: 5066 Pa
-------------------------------------------------------------------------------
z velocity spread_rate T lambda eField
-------------------------------------------------------------------------------
0 2.24 2.307e-19 373 -44.83 0
0.00375 2.234 1.401 373 -44.83 0
0.0075 2.219 2.802 373 -44.83 0
0.01125 2.192 4.203 373 -44.83 0
0.015 2.156 5.604 373 -44.83 0
0.01875 2.108 7.005 373 -44.83 0
0.02062 2.081 7.705 373 -44.83 0
0.0225 2.05 8.405 373 -44.83 0
0.02438 2.018 9.106 373 -44.83 0
0.02625 1.982 9.806 373 -44.83 0
0.02812 1.944 10.51 373 -44.83 0
0.03 1.903 11.21 373 -44.83 0
0.03188 1.86 11.91 373 -44.83 0
0.03281 1.837 12.26 373 -44.83 0
0.03375 1.814 12.61 373 -44.83 0
0.03469 1.79 12.96 373 -44.83 0
0.03563 1.765 13.31 373 -44.83 0
0.03656 1.74 13.66 373 -44.83 0
0.0375 1.714 14.01 373 -44.83 0
0.03797 1.701 14.18 373 -44.83 0
0.03844 1.688 14.36 373 -44.83 0
0.03891 1.674 14.53 373 -44.83 0
0.03937 1.661 14.71 373.1 -44.83 0
0.03984 1.647 14.88 373.1 -44.83 0
0.04031 1.633 15.06 373.1 -44.83 0
0.04078 1.619 15.24 373.2 -44.83 0
0.04102 1.612 15.32 373.3 -44.83 0
0.04125 1.605 15.41 373.3 -44.83 0
0.04148 1.598 15.5 373.4 -44.83 0
0.04172 1.591 15.59 373.5 -44.83 0
0.04195 1.584 15.67 373.7 -44.83 0
0.04219 1.578 15.76 373.9 -44.83 0
0.04242 1.571 15.85 374.1 -44.83 0
0.04266 1.565 15.94 374.4 -44.83 0
0.04289 1.558 16.02 374.7 -44.83 0
0.04312 1.553 16.11 375.2 -44.83 0
0.04359 1.542 16.29 376.3 -44.83 0
0.04406 1.533 16.47 378 -44.83 0
0.04453 1.526 16.64 380.4 -44.83 0
0.045 1.523 16.82 383.6 -44.83 0
0.04547 1.523 17 387.8 -44.83 0
0.0457 1.525 17.1 390.5 -44.83 0
0.04594 1.528 17.19 393.5 -44.83 0
0.04617 1.533 17.28 396.9 -44.83 0
0.04641 1.539 17.37 400.9 -44.83 0
0.04664 1.547 17.47 405.3 -44.83 0
0.04688 1.557 17.56 410.3 -44.83 0
0.04734 1.583 17.75 422 -44.83 0
0.04781 1.617 17.95 435.9 -44.83 0
0.04828 1.657 18.15 451.9 -44.83 0
0.04875 1.703 18.35 469.8 -44.83 0
0.04898 1.728 18.45 479.5 -44.83 0
0.04922 1.755 18.56 489.6 -44.83 0
0.04945 1.782 18.66 500.2 -44.83 0
0.04969 1.811 18.77 511.2 -44.83 0
0.04992 1.841 18.88 522.6 -44.83 0
0.05016 1.872 18.99 534.4 -44.83 0
0.05039 1.904 19.1 546.5 -44.83 0
0.05063 1.937 19.21 558.9 -44.83 0
0.05086 1.97 19.32 571.6 -44.83 0
0.05109 2.004 19.43 584.6 -44.83 0
0.05156 2.073 19.66 611.3 -44.83 0
0.05203 2.143 19.9 638.8 -44.83 0
0.0525 2.214 20.13 667 -44.83 0
0.05344 2.359 20.62 725.1 -44.83 0
0.05391 2.431 20.87 754.7 -44.83 0
0.05437 2.503 21.12 784.5 -44.83 0
0.05484 2.574 21.37 814.5 -44.83 0
0.05531 2.644 21.62 844.4 -44.83 0
0.05578 2.712 21.88 874 -44.83 0
0.05625 2.779 22.14 903.1 -44.83 0
0.05672 2.842 22.4 931.4 -44.83 0
0.05719 2.901 22.67 958.7 -44.83 0
0.05766 2.956 22.93 984.7 -44.83 0
0.05812 3.006 23.2 1009 -44.83 0
0.05859 3.052 23.46 1033 -44.83 0
0.05906 3.092 23.73 1054 -44.83 0
0.06 3.156 24.27 1093 -44.83 0
0.06125 3.215 24.99 1136 -44.83 0
0.0625 3.249 25.72 1172 -44.83 0
0.065 3.263 27.2 1228 -44.83 0
0.0675 3.235 28.7 1273 -44.83 0
0.07 3.18 30.21 1312 -44.83 0
0.0725 3.101 31.74 1347 -44.83 0
0.075 3.004 33.29 1378 -44.83 0
0.08 2.76 36.46 1433 -44.83 0
0.085 2.463 39.72 1483 -44.83 0
0.09 2.119 43.09 1533 -44.83 0
0.095 1.733 46.62 1587 -44.83 0
0.1 1.304 50.26 1650 -44.83 0
0.105 0.8293 53.83 1725 -44.83 0
0.11 0.3029 56.76 1811 -44.83 0
0.1125 0.02056 57.64 1853 -44.83 0
0.115 -0.2711 58.02 1893 -44.83 0
0.1162 -0.4194 58.02 1912 -44.83 0
0.1175 -0.5685 57.9 1929 -44.83 0
0.12 -0.8671 57.31 1960 -44.83 0
0.1225 -1.163 56.34 1985 -44.83 0
0.125 -1.453 55.08 2005 -44.83 0
0.1275 -1.735 53.61 2020 -44.83 0
0.13 -2.007 51.98 2031 -44.83 0
0.1325 -2.269 50.26 2039 -44.83 0
0.135 -2.52 48.47 2044 -44.83 0
0.14 -2.989 44.8 2050 -44.83 0
0.1437 -3.315 42.02 2052 -44.83 0
0.1475 -3.619 39.22 2053 -44.83 0
0.1512 -3.901 36.42 2054 -44.83 0
0.155 -4.163 33.62 2054 -44.83 0
0.1587 -4.403 30.82 2054 -44.83 0
0.1625 -4.623 28.02 2054 -44.83 0
0.1662 -4.822 25.22 2055 -44.83 0
0.17 -5 22.41 2055 -44.83 0
0.1737 -5.157 19.61 2055 -44.83 0
0.1775 -5.293 16.81 2055 -44.83 0
0.1812 -5.409 14.01 2055 -44.83 0
0.185 -5.503 11.21 2055 -44.83 0
0.1888 -5.577 8.405 2055 -44.83 0
0.1925 -5.629 5.604 2055 -44.83 0
0.1963 -5.661 2.802 2055 -44.83 0
0.2 -5.671 0 2055 -44.83 0
-------------------------------------------------------------------------------
z H2 H O O2 OH
-------------------------------------------------------------------------------
0 0.01024 -1.555e-20 3.083e-18 0.1016 1.397e-18
0.00375 0.01024 2.878e-20 2.999e-18 0.1016 1.951e-19
0.0075 0.01024 2.177e-19 3.042e-18 0.1016 3.817e-20
0.01125 0.01024 2.579e-18 4.211e-18 0.1016 2.344e-20
0.015 0.01024 3.139e-17 1.545e-17 0.1016 7.599e-20
0.01875 0.01024 3.817e-16 1.156e-16 0.1016 5.704e-19
0.02062 0.01024 2.017e-15 4.539e-16 0.1016 2.255e-18
0.0225 0.01024 9.467e-15 2.054e-15 0.1016 1.026e-17
0.02438 0.01024 4.389e-14 9.53e-15 0.1016 4.842e-17
0.02625 0.01024 2.027e-13 4.403e-14 0.1016 2.373e-16
0.02812 0.01024 9.335e-13 2.011e-13 0.1016 1.305e-15
0.03 0.01024 4.285e-12 9.056e-13 0.1016 9.218e-15
0.03188 0.01024 1.96e-11 4.015e-12 0.1016 8.642e-14
0.03281 0.01024 4.764e-11 9.461e-12 0.1016 4.018e-13
0.03375 0.01024 1.126e-10 2.381e-11 0.1016 1.922e-12
0.03469 0.01024 2.64e-10 6.126e-11 0.1016 9.32e-12
0.03563 0.01024 6.165e-10 1.581e-10 0.1016 4.46e-11
0.03656 0.01024 1.435e-09 4.057e-10 0.1016 2.062e-10
0.0375 0.01024 3.332e-09 1.03e-09 0.1016 9.011e-10
0.03797 0.01024 5.208e-09 1.685e-09 0.1016 2.007e-09
0.03844 0.01024 8.123e-09 2.815e-09 0.1016 4.466e-09
0.03891 0.01024 1.264e-08 4.744e-09 0.1016 9.849e-09
0.03937 0.01024 1.963e-08 8.005e-09 0.1016 2.138e-08
0.03984 0.01024 3.04e-08 1.344e-08 0.1016 4.54e-08
0.04031 0.01023 4.693e-08 2.237e-08 0.1016 9.381e-08
0.04078 0.01023 7.223e-08 3.676e-08 0.1016 1.874e-07
0.04102 0.01023 8.987e-08 4.705e-08 0.1016 2.635e-07
0.04125 0.01022 1.118e-07 6.028e-08 0.1016 3.692e-07
0.04148 0.01022 1.39e-07 7.719e-08 0.1016 5.142e-07
0.04172 0.01022 1.728e-07 9.865e-08 0.1016 7.106e-07
0.04195 0.01021 2.147e-07 1.257e-07 0.1016 9.73e-07
0.04219 0.01021 2.665e-07 1.595e-07 0.1016 1.319e-06
0.04242 0.0102 3.308e-07 2.016e-07 0.1016 1.771e-06
0.04266 0.0102 4.103e-07 2.535e-07 0.1016 2.351e-06
0.04289 0.01019 5.087e-07 3.171e-07 0.1016 3.088e-06
0.04312 0.01018 6.305e-07 3.945e-07 0.1016 4.015e-06
0.04359 0.01016 9.6e-07 5.975e-07 0.1016 6.541e-06
0.04406 0.01014 1.455e-06 8.749e-07 0.1015 1.016e-05
0.04453 0.01011 2.194e-06 1.244e-06 0.1015 1.518e-05
0.045 0.01007 3.289e-06 1.725e-06 0.1014 2.194e-05
0.04547 0.01003 4.897e-06 2.334e-06 0.1013 3.075e-05
0.0457 0.01 5.977e-06 2.693e-06 0.1013 3.601e-05
0.04594 0.009973 7.279e-06 3.095e-06 0.1012 4.189e-05
0.04617 0.00994 8.844e-06 3.54e-06 0.1011 4.835e-05
0.04641 0.009904 1.071e-05 4.029e-06 0.101 5.532e-05
0.04664 0.009864 1.293e-05 4.561e-06 0.1009 6.27e-05
0.04688 0.009821 1.556e-05 5.136e-06 0.1008 7.035e-05
0.04734 0.009724 2.212e-05 6.398e-06 0.1005 8.574e-05
0.04781 0.009612 3.082e-05 7.795e-06 0.1001 0.0001002
0.04828 0.009483 4.209e-05 9.333e-06 0.09962 0.0001124
0.04875 0.009338 5.628e-05 1.105e-05 0.0991 0.0001217
0.04898 0.009258 6.471e-05 1.204e-05 0.09882 0.0001252
0.04922 0.009173 7.404e-05 1.314e-05 0.09851 0.0001279
0.04945 0.009084 8.431e-05 1.44e-05 0.09819 0.0001297
0.04969 0.00899 9.556e-05 1.588e-05 0.09785 0.0001307
0.04992 0.00889 0.0001078 1.767e-05 0.09749 0.0001309
0.05016 0.008786 0.0001211 1.989e-05 0.09712 0.0001306
0.05039 0.008677 0.0001355 2.269e-05 0.09672 0.0001298
0.05063 0.008562 0.0001509 2.626e-05 0.09631 0.0001285
0.05086 0.008443 0.0001675 3.087e-05 0.09587 0.0001271
0.05109 0.008318 0.0001852 3.681e-05 0.09541 0.0001255
0.05156 0.008054 0.0002237 5.363e-05 0.09441 0.0001224
0.05203 0.007769 0.0002666 7.986e-05 0.09329 0.0001205
0.0525 0.007462 0.000314 0.0001193 0.092 0.0001213
0.05344 0.006796 0.0004199 0.0002472 0.08886 0.0001349
0.05391 0.006427 0.0004795 0.0003428 0.0869 0.0001514
0.05437 0.006039 0.0005425 0.0004628 0.08463 0.0001763
0.05484 0.005635 0.0006081 0.0006076 0.08203 0.0002106
0.05531 0.005221 0.0006752 0.0007759 0.07909 0.0002552
0.05578 0.004802 0.0007425 0.0009652 0.07584 0.0003101
0.05625 0.004384 0.0008086 0.001171 0.07234 0.0003747
0.05672 0.003976 0.0008721 0.001389 0.06867 0.0004477
0.05719 0.003584 0.0009316 0.001614 0.06492 0.000527
0.05766 0.003215 0.0009858 0.001839 0.06119 0.0006102
0.05812 0.002874 0.001034 0.002058 0.05758 0.0006947
0.05859 0.002566 0.001075 0.002268 0.05417 0.0007784
0.05906 0.002291 0.00111 0.002465 0.05102 0.0008595
0.06 0.00185 0.001157 0.002807 0.04572 0.00101
0.06125 0.001449 0.001185 0.003157 0.04057 0.001184
0.0625 0.001205 0.001184 0.003399 0.03713 0.001337
0.065 0.001017 0.001131 0.003632 0.03383 0.001589
0.0675 0.0009679 0.001061 0.003707 0.03227 0.001806
0.07 0.00096 0.000992 0.00371 0.03137 0.001996
0.0725 0.0009604 0.0009282 0.003678 0.03072 0.002165
0.075 0.0009599 0.0008696 0.003627 0.03017 0.002316
0.08 0.0009506 0.0007667 0.003493 0.02925 0.002574
0.085 0.0009313 0.000674 0.003328 0.02843 0.002797
0.09 0.0009039 0.0005875 0.003137 0.02764 0.003001
0.095 0.0008702 0.0005047 0.002923 0.02681 0.003205
0.1 0.0008317 0.0004247 0.002686 0.02585 0.003429
0.105 0.0007881 0.000348 0.002423 0.02467 0.003676
0.11 0.0007368 0.0002768 0.002129 0.02329 0.003908
0.1125 0.0007061 0.000244 0.001971 0.02258 0.003987
0.115 0.0006719 0.0002135 0.00181 0.02189 0.004031
0.1162 0.0006534 0.0001992 0.00173 0.02156 0.004036
0.1175 0.0006342 0.0001856 0.001651 0.02125 0.004032
0.12 0.0005937 0.0001605 0.001497 0.0207 0.003993
0.1225 0.0005513 0.0001381 0.001352 0.02025 0.003917
0.125 0.0005078 0.0001184 0.001217 0.01991 0.00381
0.1275 0.0004646 0.0001011 0.001094 0.01966 0.003682
0.13 0.0004227 8.619e-05 0.0009822 0.01952 0.003538
0.1325 0.000383 7.34e-05 0.0008824 0.01944 0.003388
0.135 0.0003462 6.255e-05 0.0007938 0.01943 0.003236
0.14 0.0002826 4.574e-05 0.0006465 0.01952 0.002945
0.1437 0.0002442 3.665e-05 0.0005597 0.01962 0.00275
0.1475 0.0002135 2.996e-05 0.0004911 0.01974 0.002581
0.1512 0.0001897 2.512e-05 0.0004382 0.01985 0.002441
0.155 0.0001718 2.167e-05 0.0003985 0.01994 0.00233
0.1587 0.0001587 1.926e-05 0.0003694 0.02001 0.002245
0.1625 0.0001492 1.758e-05 0.0003485 0.02007 0.002181
0.1662 0.0001425 1.644e-05 0.0003338 0.02011 0.002135
0.17 0.0001379 1.566e-05 0.0003237 0.02013 0.002103
0.1737 0.0001347 1.513e-05 0.0003167 0.02015 0.002081
0.1775 0.0001326 1.478e-05 0.0003121 0.02017 0.002065
0.1812 0.0001312 1.455e-05 0.0003089 0.02018 0.002055
0.185 0.0001302 1.44e-05 0.0003069 0.02018 0.002048
0.1888 0.0001296 1.429e-05 0.0003055 0.02019 0.002044
0.1925 0.0001292 1.423e-05 0.0003046 0.02019 0.002041
0.1963 0.0001289 1.418e-05 0.000304 0.02019 0.002039
0.2 0.0001287 1.417e-05 0.0003034 0.02019 0.002037
-------------------------------------------------------------------------------
z H2O HO2 H2O2 AR N2
-------------------------------------------------------------------------------
0 1.336e-19 1.019e-17 4.246e-18 0.8881 3.656e-71
0.00375 1.567e-18 2.318e-17 4.246e-18 0.8881 3.656e-71
0.0075 2.022e-18 1.512e-16 4.246e-18 0.8881 3.656e-71
0.01125 4.119e-18 1.512e-15 4.257e-18 0.8881 3.656e-71
0.015 2.628e-17 1.532e-14 4.403e-18 0.8881 3.657e-71
0.01875 2.709e-16 1.44e-13 6.358e-18 0.8881 3.671e-71
0.02062 1.302e-15 5.402e-13 1.616e-17 0.8881 3.727e-71
0.0225 7.421e-15 2.41e-12 8.393e-17 0.8881 4.046e-71
0.02438 4.338e-14 1.109e-11 5.463e-16 0.8881 5.72e-71
0.02625 2.521e-13 5.061e-11 3.655e-15 0.8881 1.364e-70
0.02812 1.446e-12 2.255e-10 2.422e-14 0.8881 4.486e-70
0.03 8.175e-12 9.648e-10 1.577e-13 0.8881 1.22e-69
0.03188 4.538e-11 3.865e-09 1.002e-12 0.8881 -8.271e-70
0.03281 1.254e-10 7.923e-09 3.035e-12 0.8881 -1.914e-68
0.03375 3.796e-10 1.747e-08 1.02e-11 0.8881 -1.158e-67
0.03469 1.176e-09 3.958e-08 3.501e-11 0.8881 -5.374e-67
0.03563 3.628e-09 8.97e-08 1.187e-10 0.8881 -2.206e-66
0.03656 1.103e-08 1.998e-07 3.908e-10 0.8881 -8.427e-66
0.0375 3.276e-08 4.301e-07 1.23e-09 0.8881 -3.074e-65
0.03797 5.855e-08 6.246e-07 2.234e-09 0.8881 -6.279e-65
0.03844 1.079e-07 9.184e-07 4.194e-09 0.8881 -1.343e-64
0.03891 2.012e-07 1.355e-06 7.948e-09 0.8881 -2.935e-64
0.03937 3.746e-07 1.993e-06 1.496e-08 0.8881 -6.466e-64
0.03984 6.914e-07 2.899e-06 2.764e-08 0.8881 -1.429e-63
0.04031 1.258e-06 4.146e-06 4.964e-08 0.8881 -3.156e-63
0.04078 2.247e-06 5.794e-06 8.592e-08 0.8881 -6.957e-63
0.04102 3.001e-06 6.772e-06 1.115e-07 0.8881 -1.055e-62
0.04125 4.025e-06 7.875e-06 1.442e-07 0.8881 -1.626e-62
0.04148 5.404e-06 9.096e-06 1.854e-07 0.8881 -2.532e-62
0.04172 7.254e-06 1.042e-05 2.365e-07 0.8882 -3.964e-62
0.04195 9.717e-06 1.184e-05 2.988e-07 0.8882 -6.218e-62
0.04219 1.298e-05 1.331e-05 3.733e-07 0.8882 -9.749e-62
0.04242 1.726e-05 1.479e-05 4.608e-07 0.8882 -1.525e-61
0.04266 2.285e-05 1.626e-05 5.616e-07 0.8882 -2.376e-61
0.04289 3.012e-05 1.764e-05 6.757e-07 0.8882 -3.682e-61
0.04312 3.949e-05 1.891e-05 8.025e-07 0.8882 -5.666e-61
0.04359 6.62e-05 2.092e-05 1.093e-06 0.8882 -1.267e-60
0.04406 0.0001061 2.223e-05 1.418e-06 0.8882 -2.617e-60
0.04453 0.0001641 2.287e-05 1.777e-06 0.8882 -5.08e-60
0.045 0.0002464 2.304e-05 2.175e-06 0.8882 -9.282e-60
0.04547 0.0003598 2.288e-05 2.618e-06 0.8882 -1.588e-59
0.0457 0.0004313 2.272e-05 2.86e-06 0.8882 -2.013e-59
0.04594 0.0005156 2.25e-05 3.115e-06 0.8882 -2.514e-59
0.04617 0.0006146 2.223e-05 3.381e-06 0.8882 -3.087e-59
0.04641 0.00073 2.191e-05 3.655e-06 0.8883 -3.72e-59
0.04664 0.0008638 2.153e-05 3.929e-06 0.8883 -4.401e-59
0.04688 0.001018 2.109e-05 4.191e-06 0.8883 -5.098e-59
0.04734 0.00139 2.005e-05 4.644e-06 0.8883 -6.524e-59
0.04781 0.001847 1.885e-05 4.926e-06 0.8883 -7.804e-59
0.04828 0.002391 1.754e-05 4.971e-06 0.8883 -8.821e-59
0.04875 0.003021 1.618e-05 4.753e-06 0.8883 -9.563e-59
0.04898 0.003371 1.55e-05 4.555e-06 0.8883 -9.84e-59
0.04922 0.003744 1.482e-05 4.298e-06 0.8883 -1.007e-58
0.04945 0.004142 1.416e-05 3.99e-06 0.8883 -1.026e-58
0.04969 0.004564 1.351e-05 3.645e-06 0.8883 -1.041e-58
0.04992 0.005011 1.288e-05 3.276e-06 0.8883 -1.053e-58
0.05016 0.005483 1.227e-05 2.896e-06 0.8883 -1.063e-58
0.05039 0.005982 1.17e-05 2.518e-06 0.8883 -1.07e-58
0.05063 0.006508 1.115e-05 2.153e-06 0.8883 -1.076e-58
0.05086 0.007065 1.062e-05 1.812e-06 0.8883 -1.08e-58
0.05109 0.007656 1.013e-05 1.502e-06 0.8883 -1.084e-58
0.05156 0.008932 9.223e-06 9.937e-07 0.8882 -1.089e-58
0.05203 0.01036 8.424e-06 6.235e-07 0.8881 -1.09e-58
0.0525 0.01199 7.723e-06 3.749e-07 0.888 -1.09e-58
0.05344 0.01584 6.563e-06 1.511e-07 0.8877 -1.084e-58
0.05391 0.01816 6.086e-06 1.002e-07 0.8875 -1.081e-58
0.05437 0.02079 5.662e-06 7.734e-08 0.8873 -1.076e-58
0.05484 0.02374 5.28e-06 6.998e-08 0.8872 -1.072e-58
0.05531 0.02697 4.93e-06 7.115e-08 0.887 -1.068e-58
0.05578 0.03045 4.605e-06 7.718e-08 0.8869 -1.064e-58
0.05625 0.03412 4.299e-06 8.606e-08 0.8868 -1.061e-58
0.05672 0.03789 4.009e-06 9.655e-08 0.8867 -1.058e-58
0.05719 0.04167 3.732e-06 1.077e-07 0.8867 -1.056e-58
0.05766 0.04537 3.471e-06 1.189e-07 0.8868 -1.055e-58
0.05812 0.04889 3.225e-06 1.295e-07 0.8869 -1.055e-58
0.05859 0.05217 2.997e-06 1.391e-07 0.887 -1.055e-58
0.05906 0.05516 2.788e-06 1.476e-07 0.8871 -1.055e-58
0.06 0.06013 2.435e-06 1.611e-07 0.8873 -1.057e-58
0.06125 0.06487 2.083e-06 1.736e-07 0.8876 -1.059e-58
0.0625 0.06796 1.835e-06 1.826e-07 0.8878 -1.06e-58
0.065 0.07085 1.551e-06 1.945e-07 0.8879 -1.062e-58
0.0675 0.07217 1.382e-06 1.991e-07 0.888 -1.062e-58
0.07 0.07293 1.266e-06 1.97e-07 0.888 -1.061e-58
0.0725 0.0735 1.176e-06 1.896e-07 0.8881 -1.06e-58
0.075 0.074 1.102e-06 1.789e-07 0.8881 -1.058e-58
0.08 0.0749 9.861e-07 1.535e-07 0.8881 -1.049e-58
0.085 0.07578 8.91e-07 1.27e-07 0.8881 -1.027e-58
0.09 0.07669 8.059e-07 1.02e-07 0.888 -9.807e-59
0.095 0.07767 7.222e-07 7.87e-08 0.888 -8.917e-59
0.1 0.07878 6.332e-07 5.753e-08 0.888 -7.465e-59
0.105 0.08005 5.381e-07 3.969e-08 0.888 -5.504e-59
0.11 0.08151 4.469e-07 2.663e-08 0.8881 -3.424e-59
0.1125 0.08231 4.091e-07 2.21e-08 0.8882 -2.549e-59
0.115 0.08313 3.791e-07 1.869e-08 0.8883 -1.839e-59
0.1162 0.08354 3.673e-07 1.736e-08 0.8883 -1.542e-59
0.1175 0.08395 3.577e-07 1.623e-08 0.8883 -1.284e-59
0.12 0.08474 3.442e-07 1.45e-08 0.8883 -8.71e-60
0.1225 0.08549 3.373e-07 1.332e-08 0.8883 -5.745e-60
0.125 0.08616 3.349e-07 1.254e-08 0.8883 -3.696e-60
0.1275 0.08675 3.354e-07 1.204e-08 0.8882 -2.326e-60
0.13 0.08725 3.373e-07 1.173e-08 0.8882 -1.438e-60
0.1325 0.08767 3.399e-07 1.154e-08 0.8882 -8.782e-61
0.135 0.088 3.428e-07 1.144e-08 0.8881 -5.355e-61
0.14 0.08848 3.488e-07 1.139e-08 0.8881 -1.946e-61
0.1437 0.08872 3.537e-07 1.146e-08 0.8881 -8.497e-62
0.1475 0.08888 3.587e-07 1.158e-08 0.8881 -3.583e-62
0.1512 0.08899 3.637e-07 1.174e-08 0.8881 -1.464e-62
0.155 0.08906 3.683e-07 1.191e-08 0.8881 -5.807e-63
0.1587 0.0891 3.722e-07 1.207e-08 0.8881 -2.245e-63
0.1625 0.08913 3.755e-07 1.221e-08 0.8881 -8.476e-64
0.1662 0.08915 3.78e-07 1.232e-08 0.8881 -3.135e-64
0.17 0.08917 3.799e-07 1.241e-08 0.8881 -1.139e-64
0.1737 0.08918 3.812e-07 1.247e-08 0.8881 -4.071e-65
0.1775 0.08918 3.822e-07 1.251e-08 0.8881 -1.436e-65
0.1812 0.08918 3.828e-07 1.254e-08 0.8881 -5.012e-66
0.185 0.08919 3.832e-07 1.256e-08 0.8881 -1.734e-66
0.1888 0.08919 3.835e-07 1.258e-08 0.8881 -5.957e-67
0.1925 0.08919 3.837e-07 1.259e-08 0.8881 -2.038e-67
0.1963 0.08919 3.838e-07 1.259e-08 0.8881 -6.957e-68
0.2 0.08919 3.84e-07 1.26e-08 0.8881 -2.373e-68
-------------------------------------------------------------------------------
z
-------------------------------------------------------------------------------
0
0.00375
0.0075
0.01125
0.015
0.01875
0.02062
0.0225
0.02438
0.02625
0.02812
0.03
0.03188
0.03281
0.03375
0.03469
0.03563
0.03656
0.0375
0.03797
0.03844
0.03891
0.03937
0.03984
0.04031
0.04078
0.04102
0.04125
0.04148
0.04172
0.04195
0.04219
0.04242
0.04266
0.04289
0.04312
0.04359
0.04406
0.04453
0.045
0.04547
0.0457
0.04594
0.04617
0.04641
0.04664
0.04688
0.04734
0.04781
0.04828
0.04875
0.04898
0.04922
0.04945
0.04969
0.04992
0.05016
0.05039
0.05063
0.05086
0.05109
0.05156
0.05203
0.0525
0.05344
0.05391
0.05437
0.05484
0.05531
0.05578
0.05625
0.05672
0.05719
0.05766
0.05812
0.05859
0.05906
0.06
0.06125
0.0625
0.065
0.0675
0.07
0.0725
0.075
0.08
0.085
0.09
0.095
0.1
0.105
0.11
0.1125
0.115
0.1162
0.1175
0.12
0.1225
0.125
0.1275
0.13
0.1325
0.135
0.14
0.1437
0.1475
0.1512
0.155
0.1587
0.1625
0.1662
0.17
0.1737
0.1775
0.1812
0.185
0.1888
0.1925
0.1963
0.2
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> products <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
Mass Flux: 0.06 kg/m^2/s
Temperature: 2055 K
Mass Fractions:
H2 0.0001284
H 1.41e-05
O 0.0003029
O2 0.02019
OH 0.002035
H2O 0.08919
HO2 3.841e-07
H2O2 1.26e-08
AR 0.8881
from pathlib import Path
import cantera as ct
# parameter values
p = 0.05 * ct.one_atm # pressure
T_in = 373.0 # inlet temperature
mdot_reactants = 0.12 # kg/m^2/s
mdot_products = 0.06 # kg/m^2/s
rxnmech = 'h2o2.yaml' # reaction mechanism file
comp = 'H2:1.6, O2:1, AR:7' # premixed gas composition
width = 0.2 # m
loglevel = 1 # amount of diagnostic output (0 to 5)
# Set up the problem
gas = ct.Solution(rxnmech)
# set state to that of the unburned gas at the burner
gas.TPX = T_in, p, comp
# Create the flame simulation object
sim = ct.CounterflowPremixedFlame(gas=gas, width=width)
# Set grid refinement parameters
sim.set_refine_criteria(ratio=3, slope=0.1, curve=0.2, prune=0.02)
# set the boundary flow rates
sim.reactants.mdot = mdot_reactants
sim.products.mdot = mdot_products
sim.set_initial_guess() # assume adiabatic equilibrium products
sim.show()
sim.solve(loglevel, auto=True)
if "native" in ct.hdf_support():
output = Path() / "premixed_counterflow_flame.h5"
else:
output = Path() / "premixed_counterflow_flame.yaml"
output.unlink(missing_ok=True)
sim.save(output, name="mix", description="solution with mixture-averaged transport")
# write the velocity, temperature, and mole fractions to a CSV file
sim.save("premixed_counterflow_flame.csv", basis="mole", overwrite=True)
sim.show_stats()
sim.show()
Total running time of the script: (0 minutes 0.338 seconds)