GreenPolyProducts/cantera
5
0
Fork 0
mirror of https://github.com/Cantera/cantera.git synced 2025-02-25 18:55:29 -06:00
Code Issues Packages Projects Releases Wiki Activity

[1D] Eliminate redundant calculation of stoich. mixture fraction

Browse Source
This commit is contained in:
Ray Speth 2024-06-06 21:08:40 -04:00 committed by Ingmar Schoegl
parent da9e29ba7d
commit c2ce9c49d9
1 changed files with 1 additions and 16 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

17
interfaces/cython/cantera/onedim.py
View File

@ -914,9 +914,6 @@ class CounterflowDiffusionFlame(FlameBase):
if data:
return
moles = lambda el: (self.gas.elemental_mass_fraction(el) /
self.gas.atomic_weight(el))
# Compute stoichiometric mixture composition
Yin_f = self.fuel_inlet.Y
self.gas.TPY = self.fuel_inlet.T, self.P, Yin_f
@ -925,12 +922,6 @@ class CounterflowDiffusionFlame(FlameBase):
u0f = mdotf / rho0f
T0f = self.fuel_inlet.T
sFuel = moles('O')
if 'C' in self.gas.element_names:
sFuel -= 2 * moles('C')
if 'H' in self.gas.element_names:
sFuel -= 0.5 * moles('H')
Yin_o = self.oxidizer_inlet.Y
self.gas.TPY = self.oxidizer_inlet.T, self.P, Yin_o
mdoto = self.oxidizer_inlet.mdot
@ -938,13 +929,7 @@ class CounterflowDiffusionFlame(FlameBase):
u0o = mdoto / rho0o
T0o = self.oxidizer_inlet.T
sOx = moles('O')
if 'C' in self.gas.element_names:
sOx -= 2 * moles('C')
if 'H' in self.gas.element_names:
sOx -= 0.5 * moles('H')
zst = 1.0 / (1 - sFuel / sOx)
zst = 1 / (1 + self.gas.stoich_air_fuel_ratio(Yin_f, Yin_o, 'mass'))
Yst = zst * Yin_f + (1.0 - zst) * Yin_o
# get adiabatic flame temperature and composition