基于德拜声子模型的氦-3热物性研究.docxVIP

基于德拜声子模型的氦-3热物性研究 1 热物性的需求 例如，低温下的温度具有低泡沫、低密度、高比热容量、高热导率等性质，这是用机械制冷机选择1k以下低温环境的唯一工艺选择。而这些低温制冷机的设计计算和性能优化都离不开氦-3的热力学状态参量数据。因此,提供一幅适用温度和压力区间广、精度高的氦-3的p-h图和T-s图显得十分必要。 由于过去对氦-3热物性的研究大多局限于某些狭窄温区的经验公式,没有完整的宽范围、高精度的状态方程可供使用,因此也就没有可供使用的p-h图和T-s图。这就给渴求氦-3热物性的工作带来了极大的不便。2005年,课题组基于大量氦-3的p-ρ-T、定容比热容cV和定压比热容cp数据绘制了氦-3在0.2 K～20 K温区、压力在0.000 1 MPa～20 MPa之间的p-h图和T-s图。2006年,提出了基于德拜声子模型的氦-3状态方程以及氦-3饱和曲线特征方程和熔化曲线特征方程,因而有可能基于这些热力学理论计算出氦-3宽范围的p-h图和T-s图,以填补工程应用的空缺。 2 理论基础 介绍计算中采用的基于德拜声子模型的氦-3状态方程以及氦-3饱和曲线特征方程,熔化曲线特征方程和一些热力学参量的导出关系式。 2.1 多态性能参数及其乘子法 将描述固体比热容的德拜理论用于低温流体氦-3的比热容,导出了同时适用一般可压缩液相和气相的德拜氦-3热力学状态方程。该方程是以亥姆霍兹自由能的形式给出: A=-ΘΗ0(1)+11+eC1(τ-1)?4∑i=1δi(C2i+5τ2+C2i+6τ4)1+eiC2(δ-1)(2)+3∑i=1Ci+14δi(3)+2∑i=1(4)[δi1+eC3(δ-1)(C2i+16+C2i+1711+eC4(τ-1))]+(1-ε-C5τ)21+eC6(1/δ-1)?4∑i=1[(C5τ)i-3?(C3i+19+C3i+20δ+C3i+21δ2)](5)-C34τ(6)+C35(7)(1)A=?ΘH0(1)+11+eC1(τ?1)?∑i=14δi(C2i+5τ2+C2i+6τ4)1+eiC2(δ?1)(2)+∑i=13Ci+14δi(3)+∑i=1(4)2[δi1+eC3(δ?1)(C2i+16+C2i+1711+eC4(τ?1))]+(1?ε?C5τ)21+eC6(1/δ?1)?∑i=14[(C5τ)i?3?(C3i+19+C3i+20δ+C3i+21δ2)](5)?C34τ(6)+C35(7)(1) 式中:A为亥姆霍兹自由能;Θ为德拜温度,Θ=C36δ2/3(1 +C37δ+C38δ2);对比参数δ=ρ/ρc,τ=T/Tc,其中临界温度Tc= 3.315 7 K,临界密度ρc=41.191 kg·m-3。h0为参考点的比焓,其解可由数值计算方法得到。亥姆霍兹自由能A是温度T和密度ρ的函数。所以方程(1)可以简单地表示为A=f(T,ρ)。 2.2 3.饱和曲线的特性 2.2.1 基本参数1-r的统计特性 当0.23 KT2.88 K时,采用中温区氦-3饱和蒸气压方程 ln(p/Ρa)=n∑-1akΤk+bln(Τ/Κ)(2)ln(p/Pa)=∑?1nakTk+bln(T/K)(2) 当2.88 KT3.315 7 K(临界温度)时,采用近临界点氦-3饱和蒸气压方程 pr=1-[A(1-Τr)2-α+6∑1Ci(1-Τr)i](3)pr=1?[A(1?Tr)2?α+∑16Ci(1?Tr)i](3) 式中:无量纲量pr=p/pc,Tr=T/Tc,临界指数α=0.105。 当0.002 6 KT0.23 K时,采用接近绝对零度的氦-3饱和蒸气压方程 ln(p/Ρa)=i0-r0RΤ+52ln(Τ/Κ)+ε(Τ/Κ)(4)ε(Τ/Κ)=a+b(Τ/Κ)+c(Τ/Κ)2+dln(Τ/Κ)+e/(Τ/Κ)(5)ln(p/Pa)=i0?r0RT+52ln(T/K)+ε(T/K)(4)ε(T/K)=a+b(T/K)+c(T/K)2+dln(T/K)+e/(T/K)(5) 式中:R为气体常数;r0为0 K温度下的汽化潜热,r0=20.56 J/mol;i0是从统计原理导出的化学常数,对于3He,i0=12.512 64;a=-2.838 497,b=-0.223 233,c=0.122 17,d=-0.262 89,e=-0.028 43。 2.2.2 无线性参数 在宽温区(0.003 K ～3.315 7 K)范围内,采用以下饱和密度方程: ρ/ρc=1±(c1τβ+c2τ1+β+c3τ2+β+c4τ3+β)+c5τ+c6τ2(6) 该方程引入无量纲参数τ= (Tc-T)/Tc,密度也采用对比参数形式。临界温度Tc=3.315 7 K,临界密度ρc= 41.191 kg·m-3,临界指数β=0.365 3,c