"""
A burner-stabilized lean premixed hydrogen-oxygen flame at low pressure.
Requires: cantera >= 3.0
Keywords: combustion, 1D flow, premixed flame, saving output,
multicomponent transport
"""
from pathlib import Path
import cantera as ct
p = 0.05 * ct.one_atm
tburner = 373.0
mdot = 0.06
reactants = 'H2:1.5, O2:1, AR:7' # premixed gas composition
width = 0.5 # m
loglevel = 1 # amount of diagnostic output (0 to 5)
gas = ct.Solution('h2o2.yaml')
gas.TPX = tburner, p, reactants
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.transport_model = 'mixture-averaged'
f.solve(loglevel, auto=True)
if "native" in ct.hdf_support():
output = Path() / "burner_flame.h5"
else:
output = Path() / "burner_flame.yaml"
output.unlink(missing_ok=True)
f.save(output, name="mix", description="solution with mixture-averaged transport")
f.transport_model = 'multicomponent'
f.solve(loglevel) # don't use 'auto' on subsequent solves
f.show()
f.save(output, name="multi", description="solution with multicomponent transport")
f.save('burner_flame.csv', basis="mole", overwrite=True)