# and transport properties
#
gas = IdealGasMix(rxnmech, mix)
# set its state to that of the unburned gas at the burner
gas.set(T=tburner, P=p, X=comp)
f = BurnerFlame(gas=gas, grid=initial_grid)
# set the properties at the burner
f.burner.set(massflux=mdot, mole_fractions=comp, temperature=tburner)
f.set(tol=tol_ss, tol_time=tol_ts)
f.setMaxJacAge(5, 10)
f.set(energy='off')
f.init()
f.showSolution()
f.solve(loglevel, refine_grid)
f.setRefineCriteria(ratio=200.0, slope=0.0502, curve=0.1)
f.set(energy='on')
f.solve(loglevel, refine_grid)
f.save('flame1.xml')
f.showSolution()
# write the velocity, temperature, and mole fractions to a CSV file
z = f.flame.grid()
T = f.T()
u = f.u()
# and transport properties
#
gas = IdealGasMix(rxnmech, mix)
# set its state to that of the unburned gas at the burner
gas.set(T = tburner, P = p, X = comp)
f = BurnerFlame(gas = gas, grid = initial_grid)
# set the properties at the burner
f.burner.set(massflux = mdot, mole_fractions = comp, temperature = tburner)
f.set(tol = tol_ss, tol_time = tol_ts)
f.setMaxJacAge(5, 10)
f.set(energy = 'off')
f.init()
f.showSolution()
f.solve(loglevel, refine_grid)
f.setRefineCriteria(ratio = 200.0, slope = 0.0502, curve = 0.1)
f.set(energy = 'on')
f.solve(loglevel,refine_grid)
f.save('flame1.xml')
f.showSolution()
# write the velocity, temperature, and mole fractions to a CSV file
z = f.flame.grid()
T = f.T()
