第二章R407C中单质的热物理性质计算A??1990年3月,国际能源组织(IEA)成立了>18分支机构。到1998年6月为止不到10年的时间内，纯制冷剂（R32、R125、R134a）的研究工作取得了重大进展，制作了纯制冷剂的热力学性质数据库，完成了多种热力学状态方程式，并确定为国际状态方程式。另外，所确定的国际状态方程式都符合Helmholts函数，都是建立在汽液两相范围内的高精度状态方程式。这些热力学状态方程式为了表现流体全域，虽然采用了较复杂的函数形式，但这是为了与所有热力学性质有关的实测数据不相矛盾，具有相当高的可靠性。在传质过程的研究过程中，制冷剂热物性的计算是整个传质过程计算的基础，其准确性直接影响到了整个过程的计算精度，因此寻求合适且计算不太复杂的方法来求解R32、R125和R134a的热物理性质是研究的首要问题。制冷剂的热物理性质主要包括两类：一类是热力性质，常用的有压力、比容、温度、比焓、比熵、比热容等；另一类是迁移性质，常用的有粘度和导热系数。本章将分别介绍这两类热物理性质的计算方法。§2.1制冷剂热力性质的计算§2.1.1?制冷剂基本物理量[32]1990年公布了有关原子量的值，下面的基本物理量包括3种制冷剂的摩尔质量。为了正确表现制冷剂的热力学状态曲面，临界常数（临界温度、临界压力、临界密度）对于制作状态方程式中的各种热物性相关方程式以及推算对应状态原理中的热物性是个必要的基本物理量。有关R32、R125、R134a三种HFC系纯制冷剂的的基本物理量如下表2-1所示。???表2-1制冷剂基本物理量工质摩尔质量（kg/mol）正常沸点（℃）液体密度（kg/m3）@20℃潜热（kJ/kg）@1.013bar临界温度（℃）临界压力（kPa）临界密度（kg/m3）蒸气压@20℃(bar)气体常量（J/kg.K）R3252.02-51.898039378.4578043014.74-R125120.02-48.1121916566.25361957212.169.3R134a102.03-26.51378216101.140515105.7281.49?§2.1.2?压力的计算一、饱和蒸气压INCLUDEPICTURE"mhtml:file://D:\\01hq\\专门知识\\R407C中单质的热物理性质计算_070928.mht!/Doc/data.WebNoteBooks/20060710154907/image002.gif"\*MERGEFORMATINET
[32]饱和蒸气压是最基本的平衡物性之一，也是研究其他物性参数的基础。饱和蒸气压的计算在工程中非常重要，它可以由状态方程导出，但比较复杂，故工程上通常采用直接回归的饱和蒸气压方程。饱和蒸气压方程的形式很多，本文采用Wangner方程。Wangner通过函数优化的方法提出的4参数饱和蒸气压方程是同等精度下的最优式。Wangner方程的形式为INCLUDEPICTURE"mhtml:file://D:\\01hq\\专门知识\\R407C中单质的热物理性质计算_070928.mht!/Doc/data.WebNoteBooks/20060710154907/image004.gif"\*MERGEFORMATINET
??????????????????????????（2-1）式中INCLUDEPICTURE"mhtml:file://D:\\01hq\\专门知识\\R407C中单质的热物理性质计算_070928.mht!/Doc/data.WebNoteBooks/20060710154907/image006.gif"\*MERGEFORMATINET
因为没有具体的关于单质R32、R125、R134a在Wangner方程中的系数，因此对于R32、R125、R134a的饱和蒸气压，根据每个单质选用不同的计算公式。1、???1、???R32饱和蒸气压方程[33]至今已有多位学者测量过R32的饱和蒸气压，提供了从180K至临界温度的蒸气压测量数据，根据实验中使用的样品纯度，测量的不确定度及初步拟合的偏差情况，选择了5套蒸气压数据，拟合成如下三项的蒸气压方程：INCLUDEPICTURE"mhtml:file://D:\\01hq\\专门知识\\R407C中单质的热物理性质计算_070928.mht!/Doc/data.WebNoteBooks/20060710154907/image008.gif"\*MERGEFORMATINET
???????????????????????