CoolProp与REFPROP混合工质计算差异分析

CoolProp与REFPROP混合工质计算差异分析

CoolProp Thermophysical properties for the masses 项目地址: https://gitcode.com/gh_mirrors/co/CoolProp

问题背景

在使用CoolProp的REFPROP后端进行混合工质计算时，发现与REFPROP GUI及ctREFPROP计算结果存在差异。具体表现为：当使用不同的输入参数对（如温度-压力(T,p)与焓-熵(h,s)）进行计算时，CoolProp未能保持计算结果的一致性，而ctREFPROP则表现正常。

现象重现

测试案例采用O₂(21%)、N₂(74%)和H₂O(5%)的三元混合工质，在两种计算路径下进行验证：

1. CoolProp实现

import CoolProp.CoolProp as CPCP
AS = CPCP.AbstractState('REFPROP', 'O2&N2&H2O')
AS.set_mole_fractions([0.21, 0.74, 0.05])
AS.update(CPCP.PT_INPUTS, 100000.0, 287)
print(AS.hmass(), AS.smass())
AS.update(CPCP.HmassSmass_INPUTS, AS.hmass(), AS.smass())
print(f'{AS.p(), AS.T()} != 100000.0, 287')

输出结果：

305360.4799795776 6725.036383567604
(91299.94879973205, 234.65445708663927) != 100000.0, 287

2. ctREFPROP实现

from ctREFPROP.ctREFPROP import REFPROPFunctionLibrary
RP = REFPROPFunctionLibrary(os.environ['RPPREFIX'])
o = RP.REFPROPdll("Oxygen * Nitrogen * Water","TP","H;S",units,0,0,287,1e5,[0.21, 0.74, 0.05])
H, S = o.Output[0:2]
o = RP.REFPROPdll("Oxygen * Nitrogen * Water","HS","T;P",units,0,0,H,S,[0.21, 0.74, 0.05])
T, P = o.Output[0:2]
print(f'{T, P} == 100000.0, 287')

输出结果：

8924.69869464179 130.75198569170252
(100000.00000003321, 287.0000000000218) == 100000.0, 287

问题根源

经过深入分析，发现该问题源于单位制设置的不同：

1. CoolProp默认使用质量基准单位制(MASS_BASE_SI)
2. ctREFPROP示例中使用的是摩尔基准单位制(MOLAR_BASE_SI)

当在ctREFPROP中显式设置为质量基准单位制时，计算结果与CoolProp完全一致。这说明问题并非CoolProp的bug，而是单位制选择导致的差异。

技术建议

对于工程计算中的单位制选择，建议考虑以下因素：

1. 化工过程：通常更适合使用摩尔基准单位制
2. 热力循环分析：质量基准单位制更为常用
3. 混合工质特性：
   • 对于组分浓度变化大的情况，摩尔基准更直观
   • 对于流动和能量计算，质量基准更直接

性能优化建议

对于需要高性能计算的场景：

1. CoolProp优势：

   • 语法简洁易用
   • 通过AbstractState对象避免重复初始化开销
   • 内置缓存机制提升重复计算效率

2. ctREFPROP优化技巧：

   • 使用预分配内存减少动态分配开销
   • 批量处理多个状态点计算
   • 考虑使用REFPROP-manager实现并行计算

结论

本案例揭示了热力学计算中单位制选择的重要性。开发者在使用CoolProp或ctREFPROP时应当：

1. 明确理解并统一单位制设置
2. 根据应用场景选择合适的计算工具
3. 对于关键计算结果，建议通过不同工具进行交叉验证

通过正确的单位制设置和工具选择，可以确保热力学计算的准确性和可靠性。

CoolProp Thermophysical properties for the masses 项目地址: https://gitcode.com/gh_mirrors/co/CoolProp