[1D] Use nonzero initial guess for FreeFlame mass flow rate

This helps avoids timestepping errors in the first timestep
This commit is contained in:
Ray Speth 2015-12-04 15:35:05 -05:00
parent de8348ecb2
commit a9293c3168

View file

@ -405,6 +405,10 @@ class FreeFlame(FlameBase):
u0 = self.inlet.mdot/self.gas.density
T0 = self.inlet.T
if not self.inlet.mdot:
# nonzero initial guess increases likelihood of convergence
self.inlet.mdot = 1.0
# get adiabatic flame temperature and composition
self.gas.equilibrate('HP')
Teq = self.gas.T