前言
Refprop是NIST(美国国家标准与技术研究所)的一款物性数据库。其全称为Reference Fluid Thermodynamic and Transport Properties Database,即流体热力学和输运性质数据库。在做仿真时经常需要查询物性参数,当然matlab可以调用这个软件,但是调用会影响程序的运行速度,此时可以先用Refprop查询所需的物性范围,再用matlab拟合经验公式,会极大加快运行速度。
一、Refprop查询参数
点击Options,选择Units,设置物性参数的单位,一般选择SI(国际单位制即可)。
再次点击Options,选择Properties,选择你所需要的物性参数
点击Substance,选择你所需查询的物质是何种类型(纯物质,混合物等),以空气为例外,选择Predefined Mixture(预定义的混合物),选择Air
点击Calculate,饱和状态选saturation tables,我所计算的是非饱和状态下的,选择isoproperty tables,左侧是恒定值,右侧是变量,此处以压力为恒定,变温度T。
然后输入设定的P,变化T的初始值,末值,间隔,点击OK就可以计算得到物性参数随温度变化的数据库,需要注意的是得到的数据的单位和你所需要的单位可能需要转化。
二、Matlab拟合
在工作区新建变量,将Refprop计算的数据复制粘贴进去
在matlab上方的app中选择Curve Fitting
在X data中选择自变量,此处为t,Y data中选择要拟合的变量
然后可以选择要拟合的多项式是几阶的,观察右下角图形拟合得不错就行,在左下角读出各阶系数,完成拟合经验公式
此处我选择的2阶拟合已经拟合得很好了,左下角读出二阶系数,一阶系数和常数项,导热系数cond随温度t变化的拟合公式为:cond = -7.081e-9t^2+5.953e-5t+0.02685