CanteraError thrown by CVodesIntegrator::step: CVodes error encountered. Error code: -3

[shenghui qin]

dear developer

my cantera is 2.4, I use cantera to test IDT, detail is on the ipynb file. T range from 650 to 1500 K

At 933.3 K 40 bar, ignition delay time is 0.6108529137253024 ms
Reached end time 101.06 ms in 6235 steps
At 897.9 K 40 bar, ignition delay time is 0.6043271556071341 ms
Reached end time 100.78 ms in 6452 steps
At 862.5 K 40 bar, ignition delay time is 0.614169872109367 ms
Reached end time 100.10 ms in 6775 steps
At 827.1 K 40 bar, ignition delay time is 0.6704673063845971 ms

---------------------------------------------------------------------------
CanteraError                              Traceback (most recent call last)
<ipython-input-5-e57d0132bba2> in <module>
     30 for P in [40]:
     31     for i,T in enumerate(temperatures):
---> 32         ignition_delay_times[i] = get_ignition_delay(T, P)
     33     plt.semilogy(1000./temperatures, ignition_delay_times, 'o-', label='{0} bar'.format(P))
     34     for i, element in enumerate(ignition_delay_times):

<ipython-input-3-adf54ff16f0b> in get_ignition_delay(temperature, pressure, stoichiometry, plot)
     37         # the size of the step will be determined by the ODE solver
     38         # depending on how quickly things are changing.
---> 39         reactor_network.step()
     40 
     41     print("Reached end time {0:.2f} ms in {1} steps".format(times[-1]*1e3, len(times)))

interfaces/cython/cantera/reactor.pyx in cantera._cantera.ReactorNet.step()

CanteraError: 
***********************************************************************
CanteraError thrown by CVodesIntegrator::step:
CVodes error encountered. Error code: -3

At t = 0.000295852 and h = 9.29322e-16, the error test failed repeatedly or with |h| = hmin.
Components with largest weighted error estimates:
660: 141.876
657: -122.15
7: 8.71207
2: -0.00422197
664: 0.000340841
661: -0.000103909
14: 2.29653e-05
665: 1.91783e-05
663: 1.54839e-05
76: 1.78792e-06
***********************************************************************

as you can see, it fails at T = 791.7, but if I change T from 791.7 to 792.0, it works,

"Reached end time 100.19 ms in 7027 steps
At 792.0 K 40 bar, ignition delay time is 0.8989770969860696 ms"

then, next T:756.3, cvode error -3 again.

how to solve this problem?

best regards,

shenghui

[shenghui qin]

for your information:

I did some trial, 

when T is lower than 720K, it always get CVode error, but if large than 720K, it works well.

[shenghui qin]

similar to this issue,

https://groups.google.com/forum/#!searchin/cantera-users/https$3A$2F$2Fgithub.com$2FCantera$2Fncm-2019-materials$2Fblob$2Fmaster$2Fthermo_debugging.ipynb%7Csort:date/cantera-users/9kA7C5GJ-qs/BUyBPtC3BAAJ

I follow this instruction, using the debugging tool ( https://github.com/Cantera/ncm-2019-materials/blob/master/thermo_debugging.ipynb ) found too high k value from my cti file and remove them, and problem solved!