PelePhysics
PelePhysics copied to clipboard
Ensure C1 thermo. properties for species
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).
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?
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.
hmmmm Cantera doesn't complain about discontinuities in LiDryer.
I was being stupid. Cantera warns us if it isn't C0. C1 is another story.