research/cantera
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

Allow user to set flame location

Browse source 
A very small change that allows the user to set the initial location of the flame. The original hard-coded values for "locs" are retained as the default, but the user can now modify locs.
This commit is contained in:
Jeff Santner 2017-04-06 11:35:27 -05:00 • committed by Ray Speth
parent 52dbe8c007
commit 2a38b0a765
1 changed files with 7 additions and 5 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

12
interfaces/cython/cantera/onedim.py
View file

@ -408,13 +408,16 @@ class FreeFlame(FlameBase):
self.inlet.T = gas.T self.inlet.T = gas.T
self.inlet.X = gas.X self.inlet.X = gas.X
def set_initial_guess(self): def set_initial_guess(self, locs=[0.0, 0.3, 0.5, 1.0]):
""" """
Set the initial guess for the solution. The adiabatic flame Set the initial guess for the solution. The adiabatic flame
temperature and equilibrium composition are computed for the inlet gas temperature and equilibrium composition are computed for the inlet gas
composition. The temperature profile rises linearly over 20% of the composition.
domain width to Tad, then is flat. The mass fraction profiles are set
similarly. :param locs:
A list of four locations to define the temperature and mass fraction profiles.
Profiles rise linearly between the second and third location.
Locations are given as a fraction of the entire domain
""" """
super(FreeFlame, self).set_initial_guess() super(FreeFlame, self).set_initial_guess()
self.gas.TPY = self.inlet.T, self.P, self.inlet.Y self.gas.TPY = self.inlet.T, self.P, self.inlet.Y
@ -433,7 +436,6 @@ class FreeFlame(FlameBase):
Yeq = self.gas.Y Yeq = self.gas.Y
u1 = self.inlet.mdot/self.gas.density u1 = self.inlet.mdot/self.gas.density
locs = [0.0, 0.3, 0.5, 1.0]
self.set_profile('u', locs, [u0, u0, u1, u1]) self.set_profile('u', locs, [u0, u0, u1, u1])
self.set_profile('T', locs, [T0, T0, Teq, Teq]) self.set_profile('T', locs, [T0, T0, Teq, Teq])