基于MATLAB2020b的J. R.S. THOM的水强制循环沸腾压降预测模型压降倍乘系数计算器

在课堂学习中，利用Thom模型计算流体压降的系数查图会受到种种因素而导致各种令人啼笑皆非的尴尬情况，为此，我们尝试尽可能理解原作者提出的模型原理，尽可能的利用公式法，通过MATLAB R2020b，实现键入流体压力p与管道内质量干度系数χm，计算得出基于J.R.S.Thom的水强制循环沸腾压降预测模型的 r2加速压降倍乘系数，r3摩擦压降倍乘系数与r4重力压降倍乘系数，减少因人工查阅图表所带来的的误差，方便进行计算。

二、计算器模型原理

J.R.S.Thom在其论文10.1016/0017-9310(64)90002-X中给出了一个简化方案，用于计算沸水和蒸汽两相混合物循环过程中的压降的简化方案。只要通过热平衡计算算出管流出口处的流体“质量”，就可以从中估计出摩擦、加速和重力损失。

分析涵盖了15-3000 psia的压力范围和3%-100％的出口蒸汽干度χm。该方法遵循Martinelli和Nelson对水平管的分析，但已经扩展到包括更重要的垂直蒸发管的情况，

即该程序适用于0.103422Mpa-20.6844Mpa压力，3%-100%质量干度系数的环境下。

即质量干度分数大于3%的情况下。由实验认识到，质量干度和空隙率之间存在经验关系方程

这是一个折中的方法，它意味着在任何

给定的操作压力下滑移速度比是恒定的。

当整根管子的压降与绝对压力相比很小时，沸腾的混合物沿水平管流动所产生的压降由两部分组成：

(a) 混合物流经管子并汽化时因动量增加而产生的压降——加速压降。

(b) 在两相流动过程中，由于摩擦力的作用而产生的压降——摩擦压降。

在垂直管的情况下，第三部分（c）变得很重要，即由于沸腾的混合物在管内上升时重力势能的增加而产生的压降。

在实验工作中，不可能直接测量总压降的各个部分。因此，我们建立了一个理论，据此可以评估总压降中仅由加速引起的那部分。然后，这一理论被扩展到给出垂直蒸发管中的重力压降。这两部分的总和被从观察到的总压降中扣除，其余部分被认为是由摩擦引起的摩擦压降。

在250 psia以上的压力范围内，推导出的两相摩擦压降和均匀流理论预测的压降之间的差异很小，但由于惯性，特别是垂直管中的浮力效应，总压降的组成部分受到蒸汽相对于水的滑移的影响。这种影响是明显的，即使在高压下，空隙率和体积流量比之间的差异也是如此。

这样，该方法导出了一种相对简单的计算沸腾两相流压降的方法，除了在非常低的蒸汽质量情况下。

三、符号

这样，该方法导出了一种相对简单的计算沸腾两相流压降的方法，除了在非常低的蒸汽质量情况下。

各个方程中使用的符号是在文本中引入时定义的。然而，为了便于参考，它们被集中在以下列表中。其中：

下标f、g和fg表示管道中流体的一种状态，或状态的变化。因此，f表示水相，g表示蒸汽相，fg表示从水到蒸汽的变化。

下标H和S指的是计算方法。H表示该量是由均质流理论推导出来的，而s表示使用了理论中的滑移。

使用的单位是英尺ft、磅lb质量和秒s。压力的单位是psia，压力差的单位是磅力/平方英尺(pounds force/ft2)。

表1 符号表

符号	名称	单位
A	试验段的横截面积	ft²
Ag	汽相占据的截面部分	ft²
D	测试管的孔直径	ft
D’	测试管的孔直径	in
f	单相流的范宁摩擦系数	-
G	质量流量密度	lb/ft²·s
G’	质量流量密度	lb/ft²·h
g	当地重力加速度	ft/s²
gc	无量纲转换系数磅力,（lbf）到磅达（pdl）=32.2，1lbf=4.4482216152605N	-
hf	饱和温度下的水焓	Btu/lb
hi	在进口温度下的水的焓值	Btu/lb
hfg	焓值随状态变化的变化	Btu/lb
K	无量纲流量参数	-
l	测试管长度	ft
li	试验管的预热长度	ft
m	总质量流量	lb/s
mf	液相质量流量	lb/s
mg	汽相质量流量	lb/s
p	流体压强	psia
Q	沿着管子总长度向流体增加的热量	Btu/lb
R	入口与出口之间的膨胀率	-
r1	均质流加速压降的倍乘系数	ft³/lb
r2	两相流加速压降的倍乘系数	-
r3	两相流摩擦压降的倍乘系数	-
r4	两相流重力压降的倍乘系数	-
Ti	入口处水温	℉
V	流体的平均流速	ft/s
Vf与Vg	水相与汽相的平均流速	ft/s
υ	平均比体积	ft³/lb
υi	管道入口处水的比体积	ft³/lb
υf与υg	饱和水和蒸汽的比体积	ft³/lb
υfg	相变所致的比体积变化	ft³/lb
χα	混合物的面积干度或空隙率（气相所占的截面比例）	-
χm	混合物的 "质量 "干度分数	-
χv	混合物的 "体积 "流量比	-
Δp	沿管的压降	lb/ft²
Δpacl	加速压降	lb/ft²
Δpf	所有水的摩擦压降依据	lb/ft²
ΔpTPF	两相流中的摩擦压降	lb/ft²
Δpgrav	垂直向上的重力压降	lb/ft²
α	比体积比（υg/υf）	-
γ	无量纲滑移系数（见正文）	-
ρ	混合物的密度	lb/ft³

四、计算器程序算法原理

4.1质量干度系数χm分析METHODS OF ASSESSING QUALITY

由于希望了解通过出口的两相混合物的质量QUALITY，因此必须了解其表达的各种方式。

术语 "质量 "干度分数被普遍使用，其定义为

（1）

其中mg和mf分别为出口处的蒸汽和水的质量。通常情况下，χm的值是通过热平衡计算出来的。

另一种方法是按体积而不是按质量进行比较。在这种情况下，"体积 "流量比的定义为

（2）

理想情况下，这两个定义只适用于蒸汽和水相以相同速度流动的均匀混合物（滑速比S=1）。然而，“滑移”的效果是使蒸汽相比均匀的平均流速快一点，并且在管子的横截面上占据较小的部分。相反地，水相的速度比均匀的平均值慢一点，并且占据了比预期更大的横截面部分。这种状态是通过伽马射线来探测到的，伽马射线被安排在出口处扫描管子的横截面，从这些观察中第三个定义被推导出，称为 “空隙率”，其定义为