4ad372b0b2
At the critical point dP/drho -> 0, so solver_rho_Tp is ill-conditioned and
a tight pressure residual does not imply a tight density. PropsSI('Dmass',
'T', Tcrit, 'P', Pcrit, 'CO2') returned 466.207 kg/m^3 instead of the true
467.600 kg/m^3. Add a short-circuit at the top of PT_flash (mirroring the
existing guard in PQ_flash) that detects exact-critical inputs within 1e-10
relative on both T and p and returns the tabulated critical density.
Co-authored-by: Claude Opus 4.7 (1M context) <noreply@anthropic.com>
Welcome to CoolProp
CoolProp is a thermophysical property database and wrappers for a selection of programming environments. It offers similar functionality to REFPROP, but CoolProp is open-source and free. It was originally developed by Ian Bell, at the time a post-doc at the University of Liege, in Liege, Belgium.
-
CoolProp has flexible licensing terms: Commercial - ok! Academic? - ok!
-
For Python, get the latest release via
pip install coolprop
-
... other binaries are available from SourceForge
-
There is also a bleeding edge nightly build of the development version.
-
The documentation is available for the latest release and the development version
-
For any kind of question regarding CoolProp and its usage, you can ask the CoolProp Discussions
-
... you might also find answers in our FAQ
-
If you found a bug or have an issue that requires the developers to become active, please file an issue in our issue tracker
-
Contributions to this project are welcomed and encouraged! If you wish to contribute bug fixes, patches, or new features, wrappers, or material properties, please submit a Pull Request with your code.
-
If you are new to Git and Github, please see the CoolProp Wiki for guidance on becoming a contributor to the project.
-
Have a look at the coverity stats
