制冷剂过热气体热力性质的隐式拟合方法.pdfVIP

中国工程热物理学会 编号：991070 工程热力学与能源利用学术会议 制冷剂过热气体热力性质的隐式拟合方法誊 张春路，丁国良．李灏 (上海交通大学动力与能源工程学院，上海200030) 摘要准确、快速、稳定地计算制冷剂热力性质是计算机辅助设计和工程计算的要 求。本文将隐式拟合方法与Martm-Hou方程相结合，提出制冷剂过热气体热力性质的 隐式拟合方法，该方法不仅具有形式统一、精度高、速度快、外推性和稳定性好韵优 点，而且保证了过热区与两相区之间熟力性质的连续性。 关键词制冷剂，过热气体，热力性质，拟合 1 引言 制冷剂热力性质的计算是相关领域计算机辅助设计的主要基础之一，其计算精 度、速度和稳定性将直接影响、甚至决定整个计算机辅助设计软件的计算精度、速度 和稳定性。因此，制冷剂热力性质的简化计算模型研究一直受到重视【【。】。作者在文川 中提出制冷剂饱和热力性质的隐式拟合思想，获得了很好的综合效果。然而，过热气 体的熟力性质与饱和热力性质有所不同．主要表现在两个方面。其一是过热气体热力 性质是二元函数，如果写成隐式拟合形式，将成为三元线性函数，尽管可以获得相当 多的自由项系数，但根轨迹的特性将过于复杂(从二维曲线分岔交成三维益面分岔) 而难以控制：其二是过热气体热力性质应能与已有的饱和气体热力性质具有理论上的 连续性。上述两点差异使得已有的饱和热力性质隐式拟合方法不能简单地推广到过热 气体热力性质的拟台。 为实现制冷剂过热气体热力性质在常用范围内的高效简化计算，本文借鉴了常用 的简化和适当的变形之后，获得一种适合于过热气体热力性质的三次线性拟合函数。 对多种制冷剂过热气体热力性质进行拟台，效果很好。 2制冷齐J过热气体热力性质的拟台方法 常用的过热气体热力性质计算主要涉及到两类函数．即，慨tv)和，Co,f^)。为 了从理论上保证过热气体热力性质与饱和气体热力性质的连续性，需知道相应的饱和 气体状态，故必须假定一个状态参数为已知。从制冷系统仿真角度讲，一般在涉及到 过热气体热力性质计算时，总是己知压力的(即便这～压力是迭代计算过程中的假设 ’本文研究得到国家敦委留学回国人员基金和上海交通大学科技发展基金资助 、 Ⅳ一91 487 压力)，因此奉文取压力p为已知状态参数来保证必要的连续性。 2．1 f协五砷方程 过热气体f以正力方程形式的确定可参照常见的状态方程，妇MBWR方程、 Mardn．H∞方程等，本文取用Marfia-Hou方程。 Martin-Hou方程的形式为 尺? 爿：+岛，+G ex“一哆疋) exp(一簟麓)坞+岛r+c3 ，。而+———而——+———F研—一 A、七8～T+C；exp(一^ij7T、 ‘ (V一6)。 (1) 其中，丘b为常数：疋为梅界点温度， 整理得 p一鲁=业铲+型铲 ㈦ o1％，精度相当理想。然而，在过热度为0时无法从理论上保证过热气体热力性质与 鬼和气体热力性质的连续|生．而且常数6的值也较难通过线}生最小二乘拟台获得。为 了解决上述问题，本文采用两点措施：～是将常数6取为0，尽管这样傲可能导致拟合 精度的下降；二是假设式(二)满足过热气体热刀}生质向饱和气体热力性质过渡的连续 性，从丽可以甩饱和气体热力参数取代式(2)中韵篷力项p。由此，式(2)可改写为 等一等qc专一扣姒考一争 +姒吾一吾，+“吾一吾，+《专一当， +州号专，+州号一吾，十州号一手， ㈣ 对f3)进行撮小二乘拟合，可以确定系数m(，=l，2，3，…