CoolProp在特定温压下无法计算混合气体密度的问题排查

使用Python的CoolProp库计算特定压力和温度下的混合气体密度，多数场景运行正常，但在290.15K、34.2bar至35.5bar的温压组合下会触发错误。其他压力值或温度改为289K、292K时无此问题，怀疑并非相变等物理原因导致，询问是否为CoolProp的bug。

环境：Python 3.11.8、Spyder 5.4.5、CoolProp 6.7.0

测试代码 #

from IPython import get_ipython;   
get_ipython().run_line_magic('reset',' -sf')

from CoolProp.CoolProp import PropsSI

gas_composition   =('Methane[0.9254]'
                    '&Nitrogen[0.007]'
                    '&CO2[0.008]'
                    '&Ethane[0.048]'
                    '&Propane[0.0085]'
                    '&Butane[0.0014]'
                    '&Pentane[0.0002]'
                    '&IsoButane[0.0015]')

print(PropsSI('D', 'T', 290.15, 'P', 34.2*1e5, gas_composition))
print(PropsSI('D', 'T', 290.15, 'P', 35.5*1e5, gas_composition))

# 取消注释后触发错误
# print(PropsSI('D', 'T', 290.15, 'P', 35.0*1e5, gas_composition))

正常运行输出 #

26.664741396302592
27.76256499648096

报错信息 #

取消最后一行注释后，抛出如下错误：

ValueError: No density solutions for T=290.15,p=3.5e+06,z=[ 0.422333189772, 0.00119828899686, 0.0201760036275, 0.181829678273, 0.148521663389, 0.100540045183, 0.051552331582, 0.0738487991766 ] : PropsSI("D","T",290.15,"P",3500000,"Methane[0.9254]&Nitrogen[0.007]&CO2[0.008]&Ethane[0.048]&Propane[0.0085]&Butane[0.0014]&Pentane[0.0002]&IsoButane[0.0015]")

内容的提问来源于stack exchange，提问作者Toffomat