Ensure C1 thermo. properties for species

[esclapez]

The deltaT solve in LM/LMeX can fail when the species thermodynamic is not C1, with residuals bouncing +/-. This can be the case around the low temperature to high temperature transition (often at 1000K).

[marchdf]

With CEPTR, the mechanism gives a warning when it's not C1 (through cantera). So we have this check now. We could turn it into a failure I bet. Which mech are you using?

[esclapez]

The LiDryer is one of the mechanism that has the issue. I was wondering if there isn't a way to fix that within CEPTR indeed.

[marchdf]

hmmmm Cantera doesn't complain about discontinuities in LiDryer.

[marchdf]

I was being stupid. Cantera warns us if it isn't C0. C1 is another story.

[baperry2]

Here's a script that plots and checks continuity for Cp: https://github.com/baperry2/PelePhysics/blob/check-C1/Mechanisms/test_C1.py

Plot for LiDryer:

[baperry2]

test_C1.py