# # This object will be used to evaluate all thermodynamic, kinetic, # 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()
# # This object will be used to evaluate all thermodynamic, kinetic, # 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
[zloc, tvalues] = getTempData('tdata.dat') # set the temperature profile to the values read in f.flame.setFixedTempProfile(zloc, tvalues) f.set(tol=tol_ss, tol_time=tol_ts) # show the initial estimate for the solution f.showSolution() # don't solve the energy equation f.set(energy='off') # first solve the flame with mixture-averaged transport properties f.setRefineCriteria(ratio=3.0, slope=0.3, curve=1) f.setMaxJacAge(50, 50) f.setTimeStep(1.0e-5, [1, 2, 5, 10, 20]) f.solve(loglevel, refine_grid) f.save('ch4_flame_fixed_T.xml', 'mixav', 'solution with mixture-averaged transport') print '\n\n switching to multicomponent transport...\n\n' gas.switchTransportModel('Multi') f.flame.setTransportModel(gas) f.setRefineCriteria(ratio=3.0, slope=0.1, curve=0.2) f.solve(loglevel, refine_grid) f.save('ch4_flame_fixed_T.xml', 'multi', 'solution with multicomponent transport')
[zloc, tvalues] = getTempData('tdata.dat') # set the temperature profile to the values read in f.flame.setFixedTempProfile(zloc, tvalues) f.set(tol = tol_ss, tol_time = tol_ts) # show the initial estimate for the solution f.showSolution() # don't solve the energy equation f.set(energy = 'off') # first solve the flame with mixture-averaged transport properties f.setRefineCriteria(ratio = 3.0, slope = 0.3, curve = 1) f.setMaxJacAge(50, 50) f.setTimeStep(1.0e-5, [1, 2, 5, 10, 20]) f.solve(loglevel, refine_grid) f.save('ch4_flame_fixed_T.xml','mixav', 'solution with mixture-averaged transport') print '\n\n switching to multicomponent transport...\n\n' gas.switchTransportModel('Multi') f.flame.setTransportModel(gas) f.setRefineCriteria(ratio = 3.0, slope = 0.1, curve = 0.2) f.solve(loglevel, refine_grid) f.save('ch4_flame_fixed_T.xml','multi', 'solution with multicomponent transport')