优化水和水蒸汽热力性质IAPWS—IF97计算模型的新方法.pdfVIP

PaperScience Technology 2009 Vo1．28No．6 优化水和水蒸汽热力性质 IAPWS—IF97 计算模型的新方法术 周艳明 ’ 李继庚 刘焕彬 (1．华南理工大学制浆造纸工程国家重点实验室 ，广东 广州 510640； 2．电子科技大学中山学院，广东 中山528402) 摘 要 ：为提高 IAPWS—IF97(2007．8修订版)的计算速度和精度 ，用面向对象的方法实现了水和水蒸汽热力性质 的软件并用代码优化 的方法对其计算模型进行 了优化。优化的关键方法是用合并 已知量、代数变换、删除公共子 表达式、代码外提、强度削弱对其 1、2、3区的基本方程 、导出方程和反推方程进行优化。关键在于用强度削弱的方 法对它们 的幂运算进行优化 。给出了优化规则和相应的优化算法 。优化后 ，运算次数平均减少 2．42倍 ，计算速度 平均提高 5．63倍。通过在造纸工业中的应用，表明优化后的程序具有计算速度更快、精度更高和提高了边界一致 性的优点 。 关键 词 ：IAPWS—IF97；代码优化 ；强度 削弱 ；代数 变换 ；造纸工业 中图分类号：TK1 文献标识码 ：A 文章编号 ：1671—4571(2009)06-0020-04 水和水蒸汽的热力性质广泛应用于热能、化工 、 IAPWS—IF97。首先将基本方程的系数和指数以数 冶金、建材、造纸等工业领域，其数据的准确性和可 组形式的私有成员属性存储并在类 的构造函数中初 靠性直接关系到生产设备 的设计、使用及安全性等 始化。对于基本、导出和反推方程，分别以相应的保 一 系列问题 ，而如何更快速、准确地计算它的热力性 护成员方法实现 ，成员方法命名时尽量做到顾名思 质也十分重要。 义 ，即名称与公式 的希腊名称一致并包含基本、反 虽然工业领域涌现许多采用 IAPWS—IF97 推、导出方程及 区域的意义。如 1区的基本方程可 来计算水和水蒸汽热力性质的软件 ，但鲜有对 以命名为 doublegamma—pi—tau一1(doublepi，double 其计算模型优化进行论述的。在国内相关重要期刊 tau)，导出方程可以命名为 doublegamma—d—pi—tau一 上 ，只有文 提出用 Homer法则对幂运算进行优 1(doublepi，doubletau)，其 中d表示 derivative即 化 ，本文认为这还不够，因为 Homer法则 只是 降低 导出的含义。各种属性的计算公式可以用一个公有 了幂运算的指数 ，而没有彻底消除幂运算。而提高 成员方法 实现 ，如 voidGet—parameters—pi—tau一1 IAPWS—IF97的计算速度和精度是非常重要且具有 (doubleP，doubleT，double＆v，double＆u，double 现实 的意义 。本 文根 据最 新 的 IAPWS—IF97 S，double＆ h，double＆ Cp，double Cv，double (2007．8修订版)，用面向对象的方法实现了计算水 w)。另外 ，对于该区的计算机程序验证 (Comput— 和水蒸汽热力性质的软件并用代码优化的方法对其 er—programverification)测试用例 ，可以用一个单独 计算模型进行了优化。详细地介绍了优化方法 、规 的保护成员方法实现 ，如 voidTest—Get—parameters— pitau 则和算法并简要地介绍了 IAPWS—IF97优化前后 — 一 1(void)，其它所有方程的计算机程序验证 在速度、精度和一致性方面的对比。