用matlab水和水蒸汽热力性质,基于Matlab的水与蒸汽热力学性质查询软件设计

水与水蒸汽由于其干净环保、无燃爆性,能量大且价格低廉等优点被广泛应用于火力发电、核电、家庭供暖等领域,是化工行业一种重要的工质。所以水与蒸汽常被用于这些领域的生产与科研、准确、快速地获取水与蒸汽的热力学性质尤为重要[1-4]。计算水与蒸汽热力学性质的工业标准公式“IAPWS-IF97”区域区分简单明确,边界一致性高,能有效提高水与蒸汽热力学性质的计算精度和速度。所以设计采用了最新的水与蒸汽的计算模型来计算其物理特性,由于水与蒸汽的物理特性繁多,计算比较繁杂,在现有诸多软件中有将所有的热力学性质算出的,内容复杂,没有针对性,软件采用Matlab图形用户界面“GUI”制作软件界面,编制出在已知温度和压力的情况下,计算出水与蒸汽的性质对应区域的比容、密度、比焓,粘度,比内能这些使用频繁的热力学参数值的界面软件[5-7]。1水与蒸汽的热力学性质计算方法在IAPWS-IF97公式中,区域1,2,5的基本方程都是比吉布斯自由能g的基础方程,其均可通过各自对应的吉布斯自由能基本方程g(P,T)表达,在这3个区域中,由其的基本方程可导出所有的热力学性质,导出的热力学关系式都与温度和压力这两个基本参数有关,对这些关系式进行编制,通过输入P和T即可求出对应区域的水与蒸汽的热力学性质值[8]。区域3的基本方程用亥姆霍兹自由能基本方程f(,T)描述,由该方程导出的热力学性质都与对比密度有关,不能通过压力和温度直接求出该区域的热力学性质,只能通过密度和温度-IF97公式也没有给出密度的计算式,因此必须迭代出密度值这一重要参数,从而求出对应的热力学性质[9]。1.1区域1,2,5的热力学关系式区域1,2,5的热力学关系式由其的基本方程比吉布斯自由能g(P,T)的基础方程导出[9],由这些关系式即可求出比容v,比内能,比焓h,其的关系式如下v=(g/p)T(1)=g-T(g/T)p-p(g/p)T(2)h=g-T(g/T)p(3)1.2热力学关系式及密度迭代算法区域3的基本方程是亥姆霍兹自由能基本方程f(,T)的基础方程[10],该方程导出该区域关于压力p,比内能u,比焓h的热力学关系式如下p=2(f/)T(4)u=f-T(f/T)(5)h=f-T(f/T)+(f/)T(6)由热力学关系式可知,不能直接通过输入压力和温度求出热力学性质,其均与对比密度有关,因此迭代出该密度值是问题的关键,首先根据压力关系式的导数结果关系式p(,)/RT=构造求根方程,再由该方程迭代出密度值这一重要参数[11],最后将密度代入热力学关系式求出热力学性质,其中比容值为密度的倒数。构造求根方程为y=p(,)RT-(n1-1+40i=2niIiIi-1Ji)(7)由该求根方程迭代求出密度的值,迭代采用割线迭代法通过Matlab编程实现,割线迭代法的Matlab程序如下function[p1,err,k,y]=secant(0,1,delta,max1)fork=1:max12=1-f(1)*(1-0)/(f(1)f(0));err=abs(2-1);0=1;1=2;y=f(1);if(err