ROP at Stoichiometric Surface

[ahuy...@gmail.com]

Hi all,

I'm wondering if there is a way or a syntax to retrieve the rate of progress of a reaction ONLY at the stoichiometric surface of a counterflow burner simulation.

For example:

rxnA = 103

rxnB = 0

reaction_mechanism = 'gri30.cti'

gas = ct.Solution(reaction_mechanism)

width = 0.02

f = ct.CounterflowDiffusionFlame(gas, width=width)

file_name = 'extinction.xml' #solutions from counterflow simulation

f.restore(file_name, name='solution', loglevel=0)

rop_rxnA = f.net_rates_of_progress[rxnA]

Thank you in advance.

Andy