Victor, You could find and
 

Victor,

You could find and try using a one-step kinetic model for methane oxidation which has positive forward reaction orders. As a quick test you could try changing the "FORD" for CH4 from "-0.3" to "0.3".

Alternatively, you could experiment with the parameters of the ode solver. ("ODECoeffs" entry of the "constant/chemistryProperties" file)

Dave

Hi David thanks for your he
 

Hi David

thanks for your help. As you suggested, I changed the FORD to a positive value and it runs. However, i'm not sure about the bias it could make to the kinetics

I will try with the ODECoeffs, and see.

BTW, the kinetic parameters are from the book of Stephen Turns

thanks again

best

V