[建议]2区的一个函数没有（6.2 Supplementary Equation for the Metastable-Vapor

hrbrlz

我认为可以不要，只要已知PT求其他参数就可以了。

uesoft

这样就很简单了。把您的几个函数加入就ok了。

可是我还是很疑惑，比如说汽轮机热力循环里面，是否真的用不到？难道每次按稳态求解出pt以后，再按pt校核一下是否进入亚稳态就可以了吗？

hrbrlz

在汽轮机热力计算时用不到，汽轮机热力计算都是用稳态计算的。一般由稳态的水由于过冷才能进入亚稳态引起水滴爆裂，在汽轮机热力循环计算中不考虑过冷（即使凝结水考虑过冷度，也是计算稳态的），所以不用考虑。

hrbrlz

我所说的不是闪蒸，闪蒸的水和汽都是饱和的。

我上面说的稳态的水是湿蒸汽中的水滴，它的温度等于对应压力下的饱和温度，只有过冷才能使它的温度低于对应压力下的饱和温度（亚稳态的温度都低于饱和温度，闪蒸不可能造成这种状态），另外还需要足够的热量，才能使稳态的水滴进入亚稳态的蒸汽。

这些东西没研究过，不好下结论，这只是我个人的理解。

uesoft

过饱和蒸汽

在一定温度下，超过饱和蒸汽应有的密度而仍不液化或凝华的蒸汽（参见凝华）。说明：（1）过饱和蒸汽的密度对应于较高温度时的饱和蒸汽的密度，而前者的温度低于后者的温度，所以过饱和蒸汽也叫做“过冷蒸汽”。（2）处于过饱和状态的过饱和蒸汽并不稳定，如果受到扰动或出现蒸汽凝结时的依附物，如尘埃、带电粒子，就会部分液化或凝华而回到饱和状态。（3）自然界出现过冷蒸汽的现象是经常能见到的。高空中的水蒸气有时处于过冷状态而并不凝聚成云。喷气式飞机飞过时，喷出的微粒提供了凝聚的条件，于是在晴空中拖出长长的云带。

云室

根据蒸汽的凝结理论和过饱和现象而制成的仪器。说明：（1）威尔逊于1912年根据过冷蒸汽是一种不稳定状态（亚稳态）的性质制成云室。它的原理是在云室内制成过饱和蒸汽。当放射性物质发射的高能带电粒子经过时，蒸汽分子发生电离，过饱和蒸汽以这些离子为凝结中心而凝结成液滴。在云室中这些液滴就形成一条很窄的云带。用强光从侧面照射，就可以直接看到这种径迹。（2）云室可用来探测微观粒子的运动径迹。

[此贴子已经被作者于2007-8-15 11:27:57编辑过]

uesoft

当t >T_c (临界温度)时，对应于每一个p,可有一个v值，即相应于式（6-3）具有一个实根和两个虚根；当t = T_c(临界温度)而p=p₀ 时，v有三个相等的实根，此处为转折点。以上这些区域方程式均能与实验结果相仿佛。当t < T_c (临界温度)时，实验所得曲线与按方程式描出的曲线不能完全吻合，实验得出的水平段若按方程式描出，可得象E-P-I-Q-H曲线那样，这一段曲线的部分存在以及它的物理意义最后也由实验得到证实和阐明。气体定温被压缩至H后，如果蒸气中含杂质极少，而又不受到扰动，则可以见到气体的压力继续升高而并不开始凝结的现象，即仍保持汽态而依曲线H-Q变化其状态。因为这时蒸汽的温度已低于相应压力下的饱和温度，所以叫做过冷蒸汽。但这种状态极不稳定，如果其中有杂质，即不能出现H-Q这一层；如果受到扰动，也迅速引起凝结而回到H-I水平段所表示的状态。同样，可看到液体的延迟蒸发，即指液体沿E-P这一段变化其状态，此时叫做过热液体，也是不稳定的。中间P-I-Q这一段曲线所表示的状态则是绝对不稳定的状态，在实际上是不可能存在的。

------------------引自沈维道、郑佩芝、蒋谈安合编之<<工程热力学>>218-219

高等教育出版社、1983年第2版

[此贴子已经被作者于2007-8-15 14:04:52编辑过]

hrbrlz

ue的程序没有压力和焓求声速、粘度、普朗特数、等等的函数，在1、2、3、5区可以转成压力和温度求解，但如果是湿蒸汽就不能计算了。往往汽轮机需要计算排汽的声速，以计算排汽马赫数。

WASPCN也没有此功能。

以前可以用音速c=sqr（KRT）求解，k=1.035+0.1x,其中x是蒸汽干度。

uesoft

劳您久等，致歉。

2007年9月14日18:00，已经按照hrbrlz的要求，和WASPCN介绍的方法，发布了WaterPro7.0的一个最新版本，提供了根据压力、焓计算其他参数的功能和动态库函数，并改正了hrbrlz发现的PS()错误。增加的动态库函数有：

PH2SSP
PH2CP
PH2CV
PH2E
PH2KS
PH2ETA
PH2U
PH2RAMD
PH2PRN
PH2EPS
PH2N

计算亚稳态的函数

PT2HMV
PT2VMV
PT2SMV
PT2EMV
PT2CPMV
PT2CVMV
PT2SSPMV

下载http://www.uesoft.com/forum.php?mod=viewthread&tid=2611&extra=page%3D1

[此贴子已经被作者于2007-9-15 10:53:44编辑过]

hrbrlz

WASPCN介绍的方法计算是错误的，计算声速错误。

下面是我验证的例子。
先用压力和焓计算出对应的刚度(X)，用压力计算出对应的饱和水声速(w1)和饱和汽的声速(w2)，根据w1+(w2-w1)*X来计算湿蒸汽的声速。
1.湿蒸汽例子：
P=0.0080586MPa H=2359.73kJ/kg 用上述方法计算c=536.71m/s.
我们原来的程序计算值是382.9m/s（比较复杂的求解）。
UE上的帖子我写错了，声速简便算法是用c=sqr（KPV）求解，k=1.035+0.1x,其中x是蒸汽干度，如果用此公式计算声速是385.21m/s。

2.过热蒸汽例子：
P=0.0080586MPa H=2600kJ/kg 直接用WASPCN计算c=446.73561m/s.
用我们原来的程序算：c=446.9921m/s.
用c=sqr（KPV）求解，c=446.73561m/s，与直接用WASPCN计算一样，因原来K就是用声速计算的，只不过反算一下而已。
也就是我们原来的程序计算比较准确，在算湿蒸汽时，用w1+(w2-w1)*X计算不对。

另外当蒸汽是湿蒸汽时，用PH求其他参数也不应该先求T在通过P、T求其他参数（你们的源程序是这样的）。湿蒸汽时用相同的P不同的H求的T是一样的，都是饱和温度。如P=10MPa，H=2000和H=2500计算的T都是311℃，再用PT求其它参数时只要P一样，其它参数都是一样的，所以是错误的。

uesoft

谢谢hrbrlz指出。我们就把您和杨工当老师了，WaterPro的开发项目顾问就是你们。我现在要负责新一代三维工厂设计系统AutoPDMS和管网流体分析软件AutoPFA7.0的指导工作，只好麻烦您管理一下WaterPro，况且WaterPro可以直接用于你们的计算，负责把它搞准确了，对大家都有益。我们现在是否就按您的改？要不，我给您开发人员的电话，您直接跟她沟通好了。您直接找罗日锦，电话0731-8228804，8228803。我已给她说了此事。

[此贴子已经被作者于2007-9-21 14:30:42编辑过]

以下是引用hrbrlz在2007-9-19 14:57:00的发言：

WASPCN介绍的方法计算是错误的，计算声速错误。

下面是我验证的例子。
先用压力和焓计算出对应的刚度(X)，用压力计算出对应的饱和水声速(w1)和饱和汽的声速(w2)，根据w1+(w2-w1)*X来计算湿蒸汽的声速。
1.湿蒸汽例子：
P=0.0080586MPa H=2359.73kJ/kg 用上述方法计算c=536.71m/s.
我们原来的程序计算值是382.9m/s（比较复杂的求解）。
UE上的帖子我写错了，声速简便算法是用c=sqr（KPV）求解，k=1.035+0.1x,其中x是蒸汽干度，如果用此公式计算声速是385.21m/s。

2.过热蒸汽例子：
P=0.0080586MPa H=2600kJ/kg 直接用WASPCN计算c=446.73561m/s.
用我们原来的程序算：c=446.9921m/s.
用c=sqr（KPV）求解，c=446.73561m/s，与直接用WASPCN计算一样，因原来K就是用声速计算的，只不过反算一下而已。
也就是我们原来的程序计算比较准确，在算湿蒸汽时，用w1+(w2-w1)*X计算不对。

另外当蒸汽是湿蒸汽时，用PH求其他参数也不应该先求T在通过P、T求其他参数（你们的源程序是这样的）。湿蒸汽时用相同的P不同的H求的T是一样的，都是饱和温度。如P=10MPa，H=2000和H=2500计算的T都是311℃，再用PT求其它参数时只要P一样，其它参数都是一样的，所以是错误的。

你可能没有好好看源代码，PH求其他参数时，是先判断区域的，一旦在湿蒸汽区域（区域4），是采用

X := TH2XReg4(T,H);
S := TX2SReg4(T,X);
V := TX2VReg4(T,X);

来计算的，不会出现你所说错误。

另外，我前面提出的湿蒸汽声速计算方法，没有验证过，根据你的验证，这个算法看来是不行。我再好好研究研究，想个好的办法。

[此贴子已经被WASPCN于2008-7-16 17:56:06编辑过]

你误会了，我说的PH求“其他参数”不是指t、x、v、s，这些原来有的函数都是正确的。我说的是我上面说的声速、粘度、普朗特数等（除t、x、v、s外没有的参数）。这些参数，UE都是采用先求T在求其他参数的做法，见uewasp.cpp的766行。

uesoft

谢谢WASPCN