IAPWS-IF97水和水蒸汽热力性质计算程序的研制.pdfVIP

东 北 电 力 大 学 学 报 第28卷第1期 Journal Of Northeast Dianli University Vo1．28，No．1 2008年2月 Natural Science Edition Feb．．2oo8 文章编号：1005—2992(2008)01—0032一o4 IAPWS．I F97水和水蒸汽热力性质计算程序的研制 赵金峰 ，吴立峰2，张荣欣3，李 勇3 (1．吉林省电力科学研究院有限公司，吉林长春，130021；3．东电一公司，辽宁铁岭112000；2．东北电力大学能源 与机械工程学院，吉林吉林，132102) 摘 要：根据国际水和水蒸汽性质协会(IAPWS)在1997年推出的一系列新的公式，简称为IAPWS— IF97公式，同时针对比焓．比熵子程序计算过程的复杂性，提出了用黄金分割法搜索压力的算法，建立了 一 个已知比焓．比熵求解压力的数学模型，进行水和水蒸气热力性质计算程序的研制。该程序的计算精 度能完全满足实际工程中的需要，且在实际应用中未发现有不收敛点。 关 键 词：水蒸汽性质；黄金分割法；计算程序 中图分类号：TK212 文献标识码：A 水和水蒸气热力性质的计算是汽轮机及动力工程相关领域的设计、试验、科学研究中的必不可少的 基础和工具。以往采用手工查取图表的方式来获得所需的参数和数据，但其速度慢，精度低。随着计算 机的应用和普及，采用编程计算水蒸汽参数成为必由之路。 IAPWS—IF97与IFC-67相比，IAPWS—IF97有效提高了热力学性质的计算精度和速度。但IAPWS— IF97公式仅仅给出了以压力，温度为自变量的基本方程和以压力，比焓或压力，比熵为自变量的反推方 程。这样要编制已知比焓，比熵求其它参数的比焓一比熵子程序就很复杂，文献[1，2]提出了采用二维迭 代算法编制，但若迭代初值选用不当，在计算某些参数点就容易发散，并且该算法程序编制复杂。因此 本文提出了基于优化算法的比焓一比熵子程序的编制方法，该方法简单适用，易于编制程序，并且精度能 完全满足工程需要，在实际使用中未发现不收敛点。 IAPWS—IF07包括了不同区域的一系列方程，其有效范围为 根据我国目前所用的参数范围，采用区域1,2和4的 图 nPWs一 的各区域和方程 公式已能满足一般工程计算和教学需要。区域3属于临界水区和蒸汽区，区域5属于超高温过热区，没 收稿日期：2007—12—15 作者简介：赵金峰(1974一)，吉林省电力科学研究院有限公司，硕士研究生，主要从事汽轮机热力试验及节能、优化运行方面的研究 工作． 第1期 赵金峰等：IAPWS—IF97水和水蒸汽热力性质计算程序的研制 33 有被包括在所编制的程序内。 对于区域1、2，根据IAPWS—IF97给定的基本方程，以温度，压力为自变量可以直接求取其他热力学 性质；根据反推方程，以压力，比焓或压力，比熵表示的区域1和2的热力学性质无需迭代即可计算。对 于区域4，可以根据基本方程由压力求饱和温度；rE可以直接利用反推方程由温度求饱和压力。这些反 推方程与所对应的基本方程在数值上非常一致。程序用语言写成子程序的形式。具体确定水和水蒸汽 热力性质的子程序及其功能见表1。 表1 确定水和水蒸汽热力性质的子程序 注：表中各子程序中压力的符号为P，单位为MPa，温度的符号为t，单位为oc，比容的符号为 ，单位m3／kg。比焓的符号为h，单位为 kJ／kg，比熵s单位为kJ／(kg·K) 子程序psk和tsk的编制比较简单，只需利用区域4的基本方程即可。子程序 按照区域1的基 本方程进行计算，子程序ptg按照区域2的基本方程进行计算。子程序pt由子程序 和子程序ptg合 成而来。子程序ph和ps的编制则稍微复杂一些，因为涉及的区域2反推方程又分为20、26和2c三个 子区域，各子区域的反推方程不同。子程序hs编制更为复杂，它所输入的已知参数为比熵和比焓，没有 直接对应的基本方程和反推方程，需迭代求解，本文