ndof = solution_list.size // nt
nq = ndof // nx // ny
# m = ChemicalMechanismSpec(cantera_xml='ethylene-luo.xml', group_name='ethylene-luo')
# m = ChemicalMechanismSpec(cantera_xml='coh2-hawkes.xml', group_name='coh2-hawkes')
m = ChemicalMechanismSpec(cantera_xml='heptane-liu-hewson-chen-pitsch-highT.xml', group_name='gas')
variable = 'T'
pressure = 101325.
air = m.stream(stp_air=True)
air.TP = 300., pressure
# sg = m.stream('X', 'H2:1, CO:2')
sg = m.stream('X', 'H2:1, CO:1, CH4:1, CO2:1')
fuel2 = m.copy_stream(sg)
# fuel2.TP = 300., pressure
fuel2.TP = 400., pressure
# c2h4 = m.stream('X', 'C2H4:1')
eg = m.stream('X', 'CO2:1, H2O:1, CO:0.5, H2:0.001')
c7 = m.stream('X', 'NXC7H16:1')
# fuel1 = m.copy_stream(eg)
# fuel1.TP = 300., pressure
fuel1 = m.copy_stream(c7)
fuel1.TP = 485., pressure
fuel1_name = 'C7'
fuel2_name = 'SG'
f = _Flamelet2D(m, 'equilibrium', pressure, air, fuel1, fuel2, 10., 10., grid_1=x_range, grid_2=y_range)
import numpy as np
import matplotlib.pyplot as plt
from matplotlib.colors import Normalize
m = ChemicalMechanismSpec(cantera_xml='coh2-hawkes.xml',
group_name='coh2-hawkes')
print(m.n_species, m.n_reactions)
pressure = 101325.
air = m.stream(stp_air=True)
air.TP = 300., pressure
synthesis_gas = m.stream('X', 'H2:1, CO:1')
exhaust_gas = m.stream('X', 'CO2:1, H2O:1, CO:0.5, H2:0.001')
fuel1 = m.copy_stream(synthesis_gas)
fuel1.TP = 1000., pressure
fuel2 = m.copy_stream(exhaust_gas)
fuel2.TP = 400., pressure
fuel1_name = 'SG'
fuel2_name = 'EG'
x_cp = m.stoich_mixture_fraction(fuel2, air)
y_cp = m.stoich_mixture_fraction(fuel1, air)
x_cp = x_cp if x_cp < 0.5 else 1. - x_cp
y_cp = y_cp if y_cp < 0.5 else 1. - y_cp
x_cc = 6.
y_cc = 6.