[1D] Make FreeFlame mdot guess pressure dependent

This gives a better guess for both high and low pressure flames, since velocity
is roughly constant while mass flux can vary by orders of magnitude.
This commit is contained in:
Ray Speth 2016-02-15 17:00:04 -05:00
parent 832ee88692
commit 5a08cf043f

View file

@ -407,7 +407,7 @@ class FreeFlame(FlameBase):
if not self.inlet.mdot:
# nonzero initial guess increases likelihood of convergence
self.inlet.mdot = 1.0
self.inlet.mdot = 1.0 * self.gas.density
# get adiabatic flame temperature and composition
self.gas.equilibrate('HP')