水蒸氣课程设计报告.docVIP

核工程课程设计报告 题目：水蒸气参数计算及数据处理 姓 名： 班 级： 学 号： 指导教师： 2012 年 12 月5 日 目录 一. 设计的目的和要求 1 二. 水和水蒸气热力参数计算方法及程序流程图 2 前言 2 1 IAPWS-IF97概述 2 1.1 IAPWS-IF97的结构 2 2 各区域的基本方程 1 2.1 区域1的方程 1 2.1.1 基本方程 1 2.2 区域2的方程 2 2.2.1 区域2、3的边界辅助方程 2 2.2.2 基本方程 3 2.3 区域3的方程 5 2.3.1 区域3的基本方程 5 2.4 区域4的方程 7 2.4.1饱和压力方程（基本方程） 7 2.5 区域5的方程 8 三. 水和水蒸气热力参数程序计算结果 13 四. 实例应用 15 五. 程序及计算结果 15 六. 总结 15  设计的目的和要求 设计题目 水蒸气参数计算及数据处理课程设计 课程设计的目的 通过MATLAB 把复杂的水和水蒸气的热力性质计算公式程序化，大大降低计算工作量，达到非常高的精确度。 设计内容 熟悉MATLAM软件，懂得MATLAB的一些编程方法； 使用MARLAB进行水和水蒸气热力学性质的计算。 水和水蒸气热力参数计算方法及程序流程图 前言 水和水蒸汽作为一种常规工质在动力系统中有着广泛的应用。 特别是在火力发电中准确计算出实际工况下水和蒸汽的各项属性值 对准确控制提高燃料利用率和发电系统的安全性，都具有十分重要的应用意义。准确计算水和蒸汽的热力性质是工业应用和科学研究的一个重要基础。国际水和蒸汽性质协会(IAPWS)在1997年建议开始使用新的水和蒸汽热力学性质的IAPWS工业公式(简称为IAPWS-IF97)[1-2]。该工业标准将计算公式划分为5个区域,其中第3区域也是亚临界湿蒸汽发生器[3],即直流注汽锅炉热力计算经常用到的区域。在第3区域中计算给定比熵s和给定比焓h下的压力p、温度T、比定容热容Cv、比定压热容Cp、声速w、动力粘度η和导热系数λ等物性参数的方程都是基于亥姆霍兹自由能的基本方程推导出来的。工业上使用的水和蒸汽属性间的计算一般是在温度(T，单位为 K)、压强 (p，单位 MPa)、比焓(h单位为 kJ/kg)、比熵(s，单 位 为 kJ/(kg·℃))]、比容(v单位m3/kg)之间已知其中任意两个参数， 来求解其他参数的过程 。IAPWS-IF97中的3区反推求解，提出一种牛顿迭代反推求解方法，并与线性拟合公式和迭代公式重计算精确度和计算效率两个角度进行了相应的比较。 1 IAPWS-IF97概述 1.1 IAPWS-IF97的结构 IAPWS-IF97工业标准1997包括了不同区域的一系列方程，其有效范围为 273.15K≤T≤1073.15K p≤100MPa 1073.15K≤T≤2273.15K p≤10MPa 图1.1给出了整个有效范围内IAPWS的五个区域。除区域2、31.1直接看出。区域2、3B23方程，即，描述。区域1、2g(p,T)描述，区域3用对应的亥姆赫兹自由能基本方程f(ρ,T),对应于区域4的饱和曲线使用饱和压力方程ps(T)描述，高温区5也用g(p,T)方程描述。这五个方程在图1.1中为长方形的加框文字，它们组成了所谓的基本方程。 图1.1 IAPWS的各区域和方程2 各区域的基本方程 2.1 区域1的方程 本节介绍IAPWS-97区域1的基本方程的详细情况。 2.1.1 基本方程 该区域的基本方程是吉布斯自由能g的基本方程。该方程一无量纲形式表达，即,有 (1.1)，p*=16.53MPa,T*=1386K,R=0.461526kJ/(kg·K)。表1.1给出了方程1.1的系数ni及指数Ii和Ji。 表1.1 方程（1.1）的系数和指数 i Ii Ji ni i Ii Ji ni 1 0 -2 0.1463297121367 18 2 3 -0.44141845330846×10-5 2 0 -1 -0.84548187169114 19 2 17 -0.72694996297594×10-15 3 0 0 -0.375636101 20 3 -4 -0.31679644845054×10-4 4 0 1 0.33855169168385×101 21 3 0 -0.28270797985312×10-5 5 0 2 -0.95