def test_ion_profile(self):
reactants = 'CH4:1.0, O2:2.0, N2:7.52'
p = ct.one_atm
Tburner = 400
width = 0.01
# Solution object used to compute mixture properties
self.gas = ct.Solution('ch4_ion.yaml')
self.gas.TPX = Tburner, p, reactants
self.sim = ct.IonBurnerFlame(self.gas, width=width)
self.sim.set_refine_criteria(ratio=4, slope=0.1, curve=0.1)
self.sim.burner.mdot = self.gas.density * 0.15
self.sim.transport_model = 'Ion'
self.sim.solve(loglevel=0, stage=2, auto=True)
# Regression test
self.assertNear(max(self.sim.E), 591.76, 1e-2)
self.assertNear(max(self.sim.X[self.gas.species_index('E')]), 8.024e-9, 1e-2)
def test_ion_profile(self):
reactants = 'CH4:1.0, O2:2.0, N2:7.52'
p = ct.one_atm
Tburner = 400
width = 0.01
# Solution object used to compute mixture properties
self.gas = ct.Solution('ch4_ion.cti')
self.gas.TPX = Tburner, p, reactants
self.sim = ct.IonBurnerFlame(self.gas, width=width)
self.sim.set_refine_criteria(ratio=4, slope=0.4, curve=0.6)
self.sim.burner.mdot = self.gas.density * 0.15
self.sim.transport_model = 'Ion'
# stage one
self.sim.solve(loglevel=0, auto=True)
#stage two
self.sim.solve(loglevel=0, stage=2, enable_energy=True)
# Regression test
self.assertNear(max(self.sim.E), 552.33, 1e-2)
"""
A burner-stabilized premixed methane-air flame with charged species.
Requires: cantera >= 2.5.0
"""
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.5 # 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.IonBurnerFlame(gas, width=width)
f.burner.mdot = mdot
f.set_refine_criteria(ratio=3.0, slope=0.05, curve=0.1)
f.show_solution()
f.transport_model = 'Ion'
f.solve(loglevel, auto=True)
f.solve(loglevel=loglevel, stage=2, enable_energy=True)
f.save('CH4_burner_flame.xml', 'mix', 'solution with mixture-averaged transport')
f.write_csv('CH4_burner_flame.csv', quiet=False)