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= 1.78e-06 log(ss)= 5.794
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 8.553e-06 log(ss)= 4.912
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.218e-05 log(ss)= 4.855
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.951e-05 log(ss)= 6.682
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 4.115e-06 log(ss)= 6.305
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.758e-05 log(ss)= 6.006
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 7.507e-05 log(ss)= 5.258
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 5.344e-05 log(ss)= 4.826
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 0.0001522 log(ss)= 3.497
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 0.0005779 log(ss)= 4.19
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 0.003702 log(ss)= 3.89
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 0.02109 log(ss)= 2.761
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= 3.56e-06 log(ss)= 5.379
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.425e-06 log(ss)= 4.5
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 2.03e-06 log(ss)= 4.493
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 2.89e-06 log(ss)= 4.529
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.852e-05 log(ss)= 4.686
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.483e-05 log(ss)= 4.679
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.584e-05 log(ss)= 4.603
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.127e-05 log(ss)= 4.528
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 4.756e-06 log(ss)= 4.482
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.129e-06 log(ss)= 4.471
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.808e-06 log(ss)= 4.468
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.03e-05 log(ss)= 4.434
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 2.443e-06 log(ss)= 4.396
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 2.319e-06 log(ss)= 4.371
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 9.288e-07 log(ss)= 4.355
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 3.306e-07 log(ss)= 4.348
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 2.51e-06 log(ss)= 4.341
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 3.574e-06 log(ss)= 4.313
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 4.58e-05 log(ss)= 4.186
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 3.669e-05 log(ss)= 3.932
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 0.0001045 log(ss)= 4.981
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 6.198e-06 log(ss)= 6.533
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 4.412e-06 log(ss)= 6.368
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 6.282e-06 log(ss)= 6.135
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= 3.56e-06 log(ss)= 5.584
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 2.027e-05 log(ss)= 4.806
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 3.207e-06 log(ss)= 4.822
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 6.089e-06 log(ss)= 4.814
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 4.064e-07 log(ss)= 4.786
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.447e-07 log(ss)= 4.785
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 4.634e-07 log(ss)= 4.784
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 4.399e-07 log(ss)= 4.781
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 7.046e-07 log(ss)= 4.78
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.254e-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.013e-07 log(ss)= 4.777
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 6.435e-07 log(ss)= 4.774
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 6.108e-07 log(ss)= 4.767
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 8.697e-07 log(ss)= 4.764
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 9.288e-07 log(ss)= 4.762
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 8.816e-07 log(ss)= 4.759
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.765e-07 log(ss)= 4.757
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 5.655e-07 log(ss)= 4.755
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 4.026e-07 log(ss)= 4.752
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 8.061e-08 log(ss)= 4.75
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.722e-07 log(ss)= 4.75
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 3.268e-07 log(ss)= 4.749
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 9.307e-07 log(ss)= 4.745
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 6.626e-07 log(ss)= 4.742
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 8.386e-07 log(ss)= 4.734
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.73
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.7e-06 log(ss)= 4.726
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 3.404e-07 log(ss)= 4.723
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 7.27e-07 log(ss)= 4.719
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 6.901e-07 log(ss)= 4.715
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.474e-06 log(ss)= 4.711
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 3.148e-06 log(ss)= 4.702
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.12e-06 log(ss)= 4.695
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.064e-06 log(ss)= 4.686
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 7.098e-08 log(ss)= 4.685
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 3.411e-07 log(ss)= 4.684
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 9.713e-07 log(ss)= 4.682
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.844e-06 log(ss)= 4.677
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 5.908e-06 log(ss)= 4.647
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 0.0001136 log(ss)= 4.358
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 0.0001213 log(ss)= 4.188
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 0.0001727 log(ss)= 5.289
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 1.28e-06 log(ss)= 6.555
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.836
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= 1.898e-05 log(ss)= 5.853
Attempt Newton solution of steady-state problem.
Newton steady-state solve failed.
Attempt 10 timesteps.
Final timestep info: dt= 0.0001081 log(ss)= 5.56
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 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 failed.
Attempt 10 timesteps.
Final timestep info: dt= 0.0001281 log(ss)= 4.265
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= 7.594e-05 log(ss)= 0.4135
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 47.017 seconds)