氦气热物性分析和总结.docxVIP

1、模拟过程的材料见文件夹“物性材料”包括铜、钢、铝

下边是氦气的物性计算过程

管内氦气认为是理想气体，满足理想气体状态方程

氦气物性需要两个参数确定，分别为温度和压力，在fluent里各物性只能是温度的函数，所以需要通过理想气体状态方程将压力P转换为温度T的函数，使得各物性最终只是温度的函数。

建模过程中的氦气区的体积通过查询为2.959*10^-3m3，即2.959L。

常压时氦气相对压力为27KPa，查询300K，该压力下氦气的密度为0.2053kg/m3折合成质量为

m=2.959*10^-3*0.2053=0.60748g

对应的物质的量为

n=m/M=0.6075/4mol=0.152mol

所以氦气的压力转换成温度的函数为P=n*R/V*T，其中R为8.314J/mol/K，所以P=427.08*T

2、氦气物性查询通过MATLAB调用refpropm软件来实现物性查询例如求氦气在300K，相对压力27Kpa下的热导率，调用格式为

通过改变T（压力已经为温度的函数），实现氦气热导率从300K到4.2K的求解对应的热导率曲线如图

进行四次多项式拟合结果如图

拟合系数为，拟合相关系数为，拟合误差分布为

拟合系数为

，拟合相关系数为

，拟合误差分布为

5到20K，拟合参数

20到200K，拟合参数

200K

200K到300K，拟合参数

拟合参数

测点温度[5.0599999

测点温度[5.0599999,5.0650001]，样

品温度[5.0500002,5.0549998]，温差约0.01K

测点温度[5.1999998,

5.2199998],样品温度[5.2199998,5.2399998]，温差约0.02

测点温度[11.5,11.525001]，

样品温度[11.35,11.375]，温差约0.15K

测点温度[16.115,16.17]，样品温

度[15.509999,15.565]，温差0.6K

速度分布，左侧能看出速度快

一些