Get all zeros for the flame velocity

[Sikuang Wang]

Hi everyone, hope everything is good.

  I am trying to calculate the flame velocity Uo. No matter what condition I used, I always get all zeros for flame.Uo. I check the initial conditions and did not find any errors. Does anyone else have the same issue? Here is the code:

import cantera as ct
mech_file = "gri30.yaml"
fuel_comp = {"H2": 1.0}
To = 300
Po = ct.one_atm
width = 0.2
gas = ct.Solution(mech_file)
gas.TP = To, Po
gas.set_equivalence_ratio(1.0, fuel_comp, {"O2": 1.0, "N2": 3.76})

flame1 = ct.FreeFlame(gas, width=width)
flame1.set_refine_criteria(ratio=3, slope=0.07, curve=0.14)
flame1.solve(loglevel=1, auto=True)

flame1.Uo

[Chris N]

I believe you want to access the flame1.velocity variable as that is the velocity in the domain. In the documentation, the Uo property is only non-zero when the two-point control is enabled for counter flow diffusion flame configurations.

[ROHIT MISHRA]

Hi Sikuang,

Try the variable: flame1.velocity instead. Here is the plot for the velocity of the flame you have across the grid:

--
You received this message because you are subscribed to the Google Groups "Cantera Users' Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email to cantera-user...@googlegroups.com.
To view this discussion visit https://groups.google.com/d/msgid/cantera-users/0b14d027-a8f0-4632-b782-8a90446f7262n%40googlegroups.com.