Note
Go to the end to download the full example code.
Burner-stabilized flame including ionized species#
A burner-stabilized premixed methane-air flame with charged species.
Requires: cantera >= 3.0, matplotlib >= 2.0
from pathlib import Path
import matplotlib.pyplot as plt
import cantera as ct
p = ct.one_atm
tburner = 600.0
reactants = 'CH4:1.0, O2:2.0, N2:7.52' # premixed gas composition
width = 0.05 # m
loglevel = 1 # amount of diagnostic output (0 to 5)
gas = ct.Solution('gri30_ion.yaml')
gas.TPX = tburner, p, reactants
mdot = 0.15 * gas.density
f = ct.BurnerFlame(gas, width=width)
f.burner.mdot = mdot
f.set_refine_criteria(ratio=3.0, slope=0.05, curve=0.1)
f.show()
f.solve(loglevel, auto=True)
f.solve(loglevel=loglevel, stage=2)
if "native" in ct.hdf_support():
output = Path() / "ion_burner_flame.h5"
else:
output = Path() / "ion_burner_flame.yaml"
output.unlink(missing_ok=True)
f.save(output, name="mix", description="solution with mixture-averaged transport")
f.save('ion_burner_flame.csv', basis="mole", overwrite=True)
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> burner <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
Mass Flux: 0.08419 kg/m^2/s
Temperature: 600 K
Mass Fractions:
O2 0.2201
CH4 0.05519
N2 0.7247
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> flame <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
Pressure: 1.013e+05 Pa
-------------------------------------------------------------------------------
z velocity spread_rate T lambda eField
-------------------------------------------------------------------------------
0 0.15 0 600 0 0
0.005 0.3748 0 1483 0 0
0.01 0.5997 0 2367 0 0
0.015 0.5997 0 2367 0 0
0.025 0.5997 0 2367 0 0
0.035 0.5997 0 2367 0 0
0.05 0.5997 0 2367 0 0
-------------------------------------------------------------------------------
z Uo H2 O2 H2O CH4
-------------------------------------------------------------------------------
0 0 0 0.2201 0 0.05519
0.005 0 0.0002251 0.1146 0.05888 0.02759
0.01 0 0.0004502 0.009096 0.1178 7.615e-17
0.015 0 0.0004502 0.009096 0.1178 7.615e-17
0.025 0 0.0004502 0.009096 0.1178 7.615e-17
0.035 0 0.0004502 0.009096 0.1178 7.615e-17
0.05 0 0.0004502 0.009096 0.1178 7.615e-17
-------------------------------------------------------------------------------
z CO CO2 N2 H O
-------------------------------------------------------------------------------
0 0 0 0.7247 0 0
0.005 0.007971 0.06317 0.7238 1.956e-05 0.000187
0.01 0.01594 0.1263 0.723 3.912e-05 0.000374
0.015 0.01594 0.1263 0.723 3.912e-05 0.000374
0.025 0.01594 0.1263 0.723 3.912e-05 0.000374
0.035 0.01594 0.1263 0.723 3.912e-05 0.000374
0.05 0.01594 0.1263 0.723 3.912e-05 0.000374
-------------------------------------------------------------------------------
z OH HO2 H2O2 C CH
-------------------------------------------------------------------------------
0 0 0 0 0 0
0.005 0.001699 6.771e-07 5.135e-08 9.678e-17 1.385e-17
0.01 0.003398 1.354e-06 1.027e-07 1.936e-16 2.77e-17
0.015 0.003398 1.354e-06 1.027e-07 1.936e-16 2.77e-17
0.025 0.003398 1.354e-06 1.027e-07 1.936e-16 2.77e-17
0.035 0.003398 1.354e-06 1.027e-07 1.936e-16 2.77e-17
0.05 0.003398 1.354e-06 1.027e-07 1.936e-16 2.77e-17
-------------------------------------------------------------------------------
z CH2 CH2(S) CH3 HCO CH2O
-------------------------------------------------------------------------------
0 0 0 0 0 0
0.005 2.861e-17 1.953e-18 1.142e-16 1.547e-09 2.113e-11
0.01 5.722e-17 3.906e-18 2.283e-16 3.094e-09 4.226e-11
0.015 5.722e-17 3.906e-18 2.283e-16 3.094e-09 4.226e-11
0.025 5.722e-17 3.906e-18 2.283e-16 3.094e-09 4.226e-11
0.035 5.722e-17 3.906e-18 2.283e-16 3.094e-09 4.226e-11
0.05 5.722e-17 3.906e-18 2.283e-16 3.094e-09 4.226e-11
-------------------------------------------------------------------------------
z CH2OH CH3O CH3OH C2H C2H2
-------------------------------------------------------------------------------
0 0 0 0 0 0
0.005 1.202e-16 2.124e-18 7.384e-18 5.169e-23 5.828e-21
0.01 2.404e-16 4.249e-18 1.477e-17 1.034e-22 1.166e-20
0.015 2.404e-16 4.249e-18 1.477e-17 1.034e-22 1.166e-20
0.025 2.404e-16 4.249e-18 1.477e-17 1.034e-22 1.166e-20
0.035 2.404e-16 4.249e-18 1.477e-17 1.034e-22 1.166e-20
0.05 2.404e-16 4.249e-18 1.477e-17 1.034e-22 1.166e-20
-------------------------------------------------------------------------------
z C2H3 C2H4 C2H5 C2H6 HCCO
-------------------------------------------------------------------------------
0 0 0 0 0 0
0.005 7.313e-26 4.62e-26 8.295e-31 4.216e-32 4.51e-19
0.01 1.463e-25 9.24e-26 1.659e-30 8.432e-32 9.021e-19
0.015 1.463e-25 9.24e-26 1.659e-30 8.432e-32 9.021e-19
0.025 1.463e-25 9.24e-26 1.659e-30 8.432e-32 9.021e-19
0.035 1.463e-25 9.24e-26 1.659e-30 8.432e-32 9.021e-19
0.05 1.463e-25 9.24e-26 1.659e-30 8.432e-32 9.021e-19
-------------------------------------------------------------------------------
z CH2CO HCCOH N NH NH2
-------------------------------------------------------------------------------
0 0 0 0 0 0
0.005 3.854e-19 5.811e-22 1.719e-08 2.663e-09 8.531e-10
0.01 7.707e-19 1.162e-21 3.439e-08 5.327e-09 1.706e-09
0.015 7.707e-19 1.162e-21 3.439e-08 5.327e-09 1.706e-09
0.025 7.707e-19 1.162e-21 3.439e-08 5.327e-09 1.706e-09
0.035 7.707e-19 1.162e-21 3.439e-08 5.327e-09 1.706e-09
0.05 7.707e-19 1.162e-21 3.439e-08 5.327e-09 1.706e-09
-------------------------------------------------------------------------------
z NH3 NNH NO NO2 N2O
-------------------------------------------------------------------------------
0 0 0 0 0 0
0.005 1.544e-09 1.158e-09 0.0018 5.49e-07 1.393e-07
0.01 3.088e-09 2.316e-09 0.0036 1.098e-06 2.787e-07
0.015 3.088e-09 2.316e-09 0.0036 1.098e-06 2.787e-07
0.025 3.088e-09 2.316e-09 0.0036 1.098e-06 2.787e-07
0.035 3.088e-09 2.316e-09 0.0036 1.098e-06 2.787e-07
0.05 3.088e-09 2.316e-09 0.0036 1.098e-06 2.787e-07
-------------------------------------------------------------------------------
z HNO CN HCN H2CN HCNN
-------------------------------------------------------------------------------
0 0 0 0 0 0
0.005 4.547e-08 2.55e-13 3.671e-11 2.022e-17 1.127e-20
0.01 9.094e-08 5.1e-13 7.342e-11 4.044e-17 2.253e-20
0.015 9.094e-08 5.1e-13 7.342e-11 4.044e-17 2.253e-20
0.025 9.094e-08 5.1e-13 7.342e-11 4.044e-17 2.253e-20
0.035 9.094e-08 5.1e-13 7.342e-11 4.044e-17 2.253e-20
0.05 9.094e-08 5.1e-13 7.342e-11 4.044e-17 2.253e-20
-------------------------------------------------------------------------------
z HCNO HOCN HNCO NCO AR
-------------------------------------------------------------------------------
0 0 0 0 0 0
0.005 4.623e-16 2.763e-12 6.82e-10 4.628e-11 0
0.01 9.246e-16 5.526e-12 1.364e-09 9.256e-11 0
0.015 9.246e-16 5.526e-12 1.364e-09 9.256e-11 0
0.025 9.246e-16 5.526e-12 1.364e-09 9.256e-11 0
0.035 9.246e-16 5.526e-12 1.364e-09 9.256e-11 0
0.05 9.246e-16 5.526e-12 1.364e-09 9.256e-11 0
-------------------------------------------------------------------------------
z C3H7 C3H8 CH2CHO CH3CHO HCO+
-------------------------------------------------------------------------------
0 0 0 0 0 0
0.005 1.74e-45 8.216e-47 2.656e-24 3.776e-25 9.147e-15
0.01 3.479e-45 1.643e-46 5.312e-24 7.553e-25 1.829e-14
0.015 3.479e-45 1.643e-46 5.312e-24 7.553e-25 1.829e-14
0.025 3.479e-45 1.643e-46 5.312e-24 7.553e-25 1.829e-14
0.035 3.479e-45 1.643e-46 5.312e-24 7.553e-25 1.829e-14
0.05 3.479e-45 1.643e-46 5.312e-24 7.553e-25 1.829e-14
-------------------------------------------------------------------------------
z H3O+ E
-------------------------------------------------------------------------------
0 0 0
0.005 5.946e-12 1.716e-16
0.01 1.189e-11 3.433e-16
0.015 1.189e-11 3.433e-16
0.025 1.189e-11 3.433e-16
0.035 1.189e-11 3.433e-16
0.05 1.189e-11 3.433e-16
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> outlet <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
-------------------------------------------------------------------------------
z
-------------------------------------------------------------------------------
0
************ Solving on 7 point grid with energy equation enabled ************
..............................................................................
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 2.848e-05 log(ss)= 4.886
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 5.068e-06 log(ss)= 4.854
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 3.653e-05 log(ss)= 5.738
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 9.143e-07 log(ss)= 7.089
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 2.604e-06 log(ss)= 6.327
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 8.341e-06 log(ss)= 6.088
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 2.505e-06 log(ss)= 5.868
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 8.026e-06 log(ss)= 5.768
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 5.142e-05 log(ss)= 5.378
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 0.0005857 log(ss)= 4.219
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 0.003753 log(ss)= 3.516
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 0.01603 log(ss)= 2.892
Attempt Newton solution of steady-state problem.
Newton steady-state solve succeeded.
Problem solved on [7] point grid(s).
..............................................................................
grid refinement disabled.
******************** Solving with grid refinement enabled ********************
..............................................................................
Attempt Newton solution of steady-state problem.
Newton steady-state solve succeeded.
Problem solved on [7] point grid(s).
..............................................................................
##############################################################################
Refining grid in flame.
New points inserted after grid points 0 1 2 3 4 5
to resolve C C2H2 C2H4 CH CH2 CH2(S) CH2CO CH2O CH2OH CH3 CH3OH CH4 CO CO2 H H2 H2O H2O2 HCN HCNO HCO HNCO HNO HO2 HOCN N N2 N2O NCO NH2 NH3 NO NO2 O O2 OH T velocity
##############################################################################
..............................................................................
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 4.005e-06 log(ss)= 5.711
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.14e-05 log(ss)= 4.597
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 3.608e-06 log(ss)= 4.695
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 3.853e-06 log(ss)= 4.7
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.646e-05 log(ss)= 4.694
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.854e-06 log(ss)= 4.683
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.188e-05 log(ss)= 4.673
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 4.008e-05 log(ss)= 4.565
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 8.455e-06 log(ss)= 4.502
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 4.514e-06 log(ss)= 4.379
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 2.712e-06 log(ss)= 4.364
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.931e-06 log(ss)= 4.344
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.031e-06 log(ss)= 4.333
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.858e-06 log(ss)= 4.33
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 2.645e-06 log(ss)= 4.309
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.059e-06 log(ss)= 4.295
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 4.524e-06 log(ss)= 4.278
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 3.221e-06 log(ss)= 4.236
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 6.878e-06 log(ss)= 4.195
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 2.448e-06 log(ss)= 4.093
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 6.972e-06 log(ss)= 4.017
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.396e-06 log(ss)= 3.904
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.342e-05 log(ss)= 3.754
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.209e-05 log(ss)= 3.641
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 5.738e-06 log(ss)= 3.91
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 4.085e-06 log(ss)= 4.068
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 2.617e-05 log(ss)= 4.391
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.863e-05 log(ss)= 4.654
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 6.632e-06 log(ss)= 4.698
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 7.082e-06 log(ss)= 4.75
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 5.672e-06 log(ss)= 4.805
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 3.634e-05 log(ss)= 5.165
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.15e-05 log(ss)= 6.88
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.637e-05 log(ss)= 5.98
Attempt Newton solution of steady-state problem.
Newton steady-state solve succeeded.
Problem solved on [13] point grid(s).
..............................................................................
##############################################################################
Refining grid in flame.
New points inserted after grid points 0 1 2 3 4
to resolve C C2H2 C2H4 CH CH2 CH2(S) CH2CO CH2O CH2OH CH3 CH3O CH3OH CH4 CO CO2 H H2 H2O H2O2 HCCO HCN HCNO HCO HNCO HNO HO2 N N2 N2O NCO NH2 NH3 NO NO2 O O2 OH T velocity
##############################################################################
..............................................................................
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 6.007e-06 log(ss)= 5.539
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 2.887e-05 log(ss)= 4.8
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.084e-06 log(ss)= 4.78
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 9.764e-07 log(ss)= 4.778
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 3.91e-07 log(ss)= 4.775
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 4.697e-07 log(ss)= 4.773
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 3.344e-07 log(ss)= 4.771
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.27e-06 log(ss)= 4.768
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.017e-06 log(ss)= 4.764
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.086e-06 log(ss)= 4.759
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.739e-06 log(ss)= 4.755
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 6.966e-07 log(ss)= 4.751
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.255e-06 log(ss)= 4.749
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.787e-06 log(ss)= 4.741
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 6.362e-07 log(ss)= 4.73
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 6.794e-07 log(ss)= 4.727
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 2.579e-06 log(ss)= 4.718
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.377e-06 log(ss)= 4.707
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.471e-06 log(ss)= 4.703
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.571e-06 log(ss)= 4.697
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 9.939e-07 log(ss)= 4.691
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.061e-06 log(ss)= 4.685
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 6.8e-06 log(ss)= 4.671
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 2.152e-06 log(ss)= 4.657
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 3.829e-07 log(ss)= 4.65
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 2.454e-06 log(ss)= 4.64
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.553e-06 log(ss)= 4.628
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.244e-06 log(ss)= 4.622
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.992e-06 log(ss)= 4.616
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 3.191e-06 log(ss)= 4.604
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 6.057e-06 log(ss)= 4.573
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.278e-06 log(ss)= 4.549
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 6.822e-07 log(ss)= 4.543
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.093e-06 log(ss)= 4.536
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.313e-06 log(ss)= 4.531
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 7.009e-07 log(ss)= 4.521
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 2.246e-06 log(ss)= 4.514
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 9.592e-06 log(ss)= 4.487
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.536e-05 log(ss)= 4.382
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 0.0001313 log(ss)= 3.402
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 0.0001869 log(ss)= 5.189
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 8.315e-06 log(ss)= 6.456
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 2.22e-06 log(ss)= 6.255
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 7.902e-07 log(ss)= 6.138
Attempt Newton solution of steady-state problem.
Newton steady-state solve succeeded.
Problem solved on [18] point grid(s).
..............................................................................
##############################################################################
Refining grid in flame.
New points inserted after grid points 0 1 2 3 4 5
to resolve C C2H2 C2H3 C2H4 C2H6 CH CH2 CH2(S) CH2CO CH2O CH2OH CH3 CH3CHO CH3O CH3OH CH4 CO CO2 H H2 H2O H2O2 HCCO HCN HCNO HCO HNCO HO2 N N2 N2O NCO NH3 NO O O2 OH T velocity
##############################################################################
..............................................................................
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 2.848e-05 log(ss)= 4.879
Attempt Newton solution of steady-state problem.
Newton steady-state solve succeeded.
Problem solved on [24] point grid(s).
..............................................................................
##############################################################################
Refining grid in flame.
New points inserted after grid points 0 1 2 3 4 5 6
to resolve C C2H C2H2 C2H3 C2H4 C2H5 C2H6 C3H7 C3H8 CH CH2 CH2(S) CH2CHO CH2CO CH2O CH2OH CH3 CH3CHO CH3O CH3OH CH4 CO CO2 H H2 H2O H2O2 HCCO HCCOH HCN HCNO HCO HNCO HO2 N N2 N2O NCO NH3 NO NO2 O O2 OH T velocity
##############################################################################
..............................................................................
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 5.339e-06 log(ss)= 5.866
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 2.281e-05 log(ss)= 5.914
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 7.306e-05 log(ss)= 5.45
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 0.000156 log(ss)= 5.206
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 0.000395 log(ss)= 4.526
Attempt Newton solution of steady-state problem.
Newton steady-state solve succeeded.
Problem solved on [31] point grid(s).
..............................................................................
##############################################################################
Refining grid in flame.
New points inserted after grid points 0 1 2 3 4 5 6 7 8
to resolve C C2H2 C2H3 C2H4 C2H5 C2H6 C3H7 C3H8 CH CH2 CH2(S) CH2CHO CH2CO CH2O CH2OH CH3 CH3CHO CH3O CH3OH CH4 CO CO2 H H2 H2O H2O2 HCCO HCCOH HCN HCNO HCO HNCO HO2 N N2 N2O NCO NH3 NO NO2 O O2 OH T velocity
##############################################################################
..............................................................................
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 4.271e-05 log(ss)= 5.197
Attempt Newton solution of steady-state problem.
Newton steady-state solve succeeded.
Problem solved on [40] point grid(s).
..............................................................................
##############################################################################
Refining grid in flame.
New points inserted after grid points 0 1 2 3 4 5 6 7 8 9 10 11 12
to resolve C C2H2 C2H3 C2H4 C2H5 C2H6 C3H7 C3H8 CH CH2 CH2(S) CH2CHO CH2CO CH2O CH2OH CH3 CH3CHO CH3O CH3OH CH4 CO CO2 H H2 H2O H2O2 HCCO HCCOH HCN HCNO HCO HNCO HO2 N2 N2O NH3 NO NO2 O O2 OH T velocity
##############################################################################
..............................................................................
Attempt Newton solution of steady-state problem.
Newton steady-state solve succeeded.
Problem solved on [53] point grid(s).
..............................................................................
##############################################################################
Refining grid in flame.
New points inserted after grid points 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17
to resolve C C2H2 C2H3 C2H4 C2H5 C2H6 C3H7 C3H8 CH CH2 CH2(S) CH2CHO CH2CO CH2O CH2OH CH3 CH3CHO CH3O CH3OH CH4 CO CO2 H H2 H2O H2O2 HCCO HCCOH HCN HCNO HCO HNCO HO2 N2 NO NO2 O O2 OH T velocity
##############################################################################
..............................................................................
Attempt Newton solution of steady-state problem.
Newton steady-state solve succeeded.
Problem solved on [71] point grid(s).
..............................................................................
##############################################################################
Refining grid in flame.
New points inserted after grid points 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
to resolve C2H2 C2H3 C2H4 C2H5 C2H6 C3H8 CH CH2 CH2(S) CH2CHO CH2CO CH2O CH2OH CH3 CH3CHO CH3O CH3OH CH4 CO H H2 H2O H2O2 HCCO HCCOH HCN HCO HNCO HO2 N2 NO NO2 O O2 OH T velocity
##############################################################################
..............................................................................
Attempt Newton solution of steady-state problem.
Newton steady-state solve succeeded.
Problem solved on [96] point grid(s).
..............................................................................
##############################################################################
Refining grid in flame.
New points inserted after grid points 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42
to resolve C2H2 C2H3 C2H4 C2H5 C2H6 C3H8 CH CH2 CH2(S) CH2CHO CH2CO CH2OH CH3 CH3CHO CH3O CH3OH H H2O2 HCCO HCO HO2 NO2 O OH
##############################################################################
..............................................................................
Attempt Newton solution of steady-state problem.
Newton steady-state solve succeeded.
Problem solved on [139] point grid(s).
..............................................................................
##############################################################################
Refining grid in flame.
New points inserted after grid points 0 1 2 3 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 68 69 70 71 72 73 74 75 76
to resolve C2H3 C2H5 C3H8 CH2 CH3O HCCO HCO
##############################################################################
..............................................................................
Attempt Newton solution of steady-state problem.
Newton steady-state solve succeeded.
Problem solved on [168] point grid(s).
..............................................................................
##############################################################################
Refining grid in flame.
New points inserted after grid points 85
to resolve CH3
##############################################################################
..............................................................................
Attempt Newton solution of steady-state problem.
Newton steady-state solve succeeded.
Problem solved on [169] point grid(s).
..............................................................................
no new points needed in flame
..............................................................................
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 0.0001139 log(ss)= 0.3763
Attempt Newton solution of steady-state problem.
Newton steady-state solve succeeded.
Problem solved on [169] point grid(s).
..............................................................................
no new points needed in flame
Temperature and Heat Release Rate#
fig, ax1 = plt.subplots()
ax1.plot(f.grid * 1000, f.heat_release_rate / 1e6, color='C4')
ax1.set_ylabel('heat release rate [MW/m³]', color='C4')
ax1.set_xlim(0, 3.0)
ax1.set(xlabel='distance from burner [mm]')
ax2 = ax1.twinx()
ax2.plot(f.grid * 1000, f.T, color='C3')
ax2.set_ylabel('temperature [K]', color='C3')
plt.show()
Major Species Profiles#
fig, ax = plt.subplots()
major = ('O2', 'CH4', 'H2O', 'CO2')
states = f.to_array()
ax.plot(states.grid * 1000, states(*major).X, label=major)
ax.set(xlabel='distance from burner [mm]', ylabel='mole fractions')
ax.set_xlim(0, 3.0)
ax.legend()
plt.show()
Ionized Species Profiles#
fig, ax = plt.subplots()
minor = ('E', 'H3O+', 'HCO+')
ax.semilogy(states.grid * 1000, states(*minor).X, label=minor, linestyle='--')
ax.set(xlabel='distance from burner [mm]', ylabel='mole fractions', )
ax.set_xlim(0, 3.0)
ax.legend()
plt.show()
Total running time of the script: (1 minutes 52.964 seconds)