覆盖亚临界到临界区的制冷工质物性快速算法Ⅰ计算方法研究.pdfVIP

中固T程热物理学会 T程热力学’j能源利用 学术会议论义 编吁：061032 覆盖亚临界到临界区的制冷工质物性快速算法1．计算方法 丁国良吴志刚 ．1二海交通大学制冷与低温工程研究所，』：海200030 Tel：+86—21-E-mail address：glding@sjtu．edu．cn 摘要为了在亚临界到临界点这个大范围内保证制冷剂物性计算的快速、稳定、可逆，提出将整个数据 区间划分成若干个子区间分别进行隐式拟合，并且保证相邻子区间连接点处的连续与光滑，以及通过 求分析解从隐式方程中得到显式物性计算式的方法。 关键词制冷工质物性快速计算临界区 1．引言 制冷系统仿真时，需要1F常多地调用工质物性的计算程序。比如对于一个划分100 个微元的三维稳态分布翅片管式换热器计算来训11，需要调用到制冷剂物性的次数就达 到约lO 7次。从制冷系统的性能仿真与优化角度讲，希望制冷剂的物性模型：1)计算速 度快；2)稳定性极好；3)关联参数之间的计算是完全可逆的；4)物性模型不仅连续， 而且需要光滑。现有的基于物性方程的计算方法【2l，在计算的速度、稳定性等方面均不 能达到要求。国际上已有的制冷剂快速物性计算方法只能用于弧临界区，不能用于非共 沸工质[3-j5J。，作者在上海交通大学已完成的简化物性方面的工作【6{l，还只能用于亚临界 循环工质，工况范围较窄。随着R410A、C02等：【质的推广应用，制冷空调装置已在采 用近临界循环、跨临界循环，因而快速物性方法应扩大到能够覆盖临界区。本文提出适 用于亚临界到临界点的快速物性方法，并分成饱和区、过热区、二相区分别介绍。应用 的例子与效果检验将在文【9】给出。 2．饱和热力性质快速计算方法 2．1隐式拟合方程 所需要计算的范嗣为图l中所示的 k．‘AICDm+lDlBIAl’。包括了饱和、二二相、 过热区。假定为了提高精度，需要划分 成m个子区域，其中饱和液体线‘AlC’ 划分成肌段依次为‘AlA2’、‘A2A3’…．、 ‘A。C’。而m段饱和气体线依次为 ‘BIB2’、‘B283’、…、‘BmC’。对丁= 其中第f段的饱和线，采用下面的拟合 方程形式 图1数据范围与分区问示意图 高等学校博{：点基金(20050248019) 基金项目：国家自然科学基金项U， 760 外的热力参数，如压力、焓等。假设有O+2)组数据国o，乃o)，乜，l’T,0，(匏，T,9，．．．，瓴。，％)’ 接拟合的精度不高，为此在拟合之前需要先做变换。 zk荆 (2) 方程(2)应很简单，如zk (zi2气ri9，…，(zi。·，正。)，(zie气死。)。 令 U=z～o· (3) ’，=r一％