1、热物理量测试技术1 概述两相流广泛应用于热能动力工程、核能工程、低温工程以及航天领域等许多领域。所谓两相流，广义上讲是指一种物质或两种物质在不同状态下的流动，其中气体和液体一起流动称为气液两相流。对于两相流中的气液混合物，它们可以是同一种物质，即汽液(如水和水蒸气)，也可以是两种不同的物质，即气液(如水和空气混合物)。气液两相流是一个相当复杂的问题， 。在单相流中，经过一段距离之后，就会建立一个稳定的速度场。但对于两相流，例如蒸汽和水，则很难建立一个稳定的流动，因为在管道流动中有压降产生，由于此压降作用会产生液体的蒸发，所以在研究气液两相流时必须考虑两相间的传热与传质问题。两相流学科还处于半经
2、验半理论阶段，对于两相流的流动和传热规律进行研究时，除了依靠各种数学物理模型外，还要依靠实验，这就需要两者相结合从而更好地进行研究。2 两相流压降测量 1压降，即两相流通过系统时产生的压力变化，是两相流体流动过程中的一个重要参数。保持两相流体流动所需的动力以及动力系统的容量和功率就取决于压降的大小。一般说来，两相流体流动时产生的压降一般由三部分组成，即摩擦阻力压降、重位压降、加速压降，管道系统出现阀门、孔板等管件时，还需测量局部压降。目前，常用差压计或传感器来测量两相流压降。2.1 利用差压计测量压降应用差压计测量气液两相流压降的测量原理图如图 1 所示。所测压降为下部抽头的压力与上部抽头压力
3、之差。在差压计的 Z1截面上可列出压力平衡式如下：（2.1）1+(21)=2+(43)+(31)式中， 为取压管中的流体密度； 为差压计的流体密度。 由（2.1）可得：（2.2）12=(31)()+(42)由上式可知，要算出压降 的值，必须知道取压管中的流体密度 和差压计读数12 。31当管中流体不流动时：（2.3）12=(42)式中， 为两相混合物平均密度。将式（2.3）代入（2.1） 。可得两相流体静止时，差压计中读数如下：（2.4）31=(42)图 1 气液两相流系统中的压降测量从上面的方程式可知，为了从差压计得到压降，确定取压管中流体密度 是十分重要的，这意味着取压管中的流体必须为单相
4、液体或气体。因此在测量两相流压降时，需要一个装置保证取压管中永远充满液体，一般在取压管后接一个气液分离器。 图 2 带有气液分离器的测量系统1-实验段；2-气液分离器；3-取压管；4-差压计；5-温度测点；6-排气阀如图 2 所示，气液混合物进入气液分离器后分离，气相在上部，液相在下部，这样就可保证差压计取压管中全部为液体。但此时必须知道差压计中液体的温度，因为差压计中液体的密度与温度有关。测量时试验段中为气液混合物，因此必须对两侧的密度差进行修正。2.2 利用传感器测量压降利用传感器测量两相流压降有两种方式，一是利用两台在测点安装的压力传感器获得的信号相减来得到压降；二是利用差压传感器来获得
5、压降。第一种方式主要有电容式和压电式压力传感器，电子设备将两压力传感器输入的电信号相减即可换算得两测压点之间的压降。适用于需要快速时间响应的场合，但是有明显的缺点，将两个独立测出的电信号再进行相减会引起误差增大。使用此法时应对两个压力传感器进行校准，力求使输出信号能较精确地变换成所需测定的压降。第二种方式主要有磁阻式差压传感器和应变仪式差压传感器。差压传感器可以克服上一种方式的缺点。由于差压传感器传递压降时，膜片位移很小，因而导压管中流体流动量较小时，气泡不易进入导压管，测量结果更精确。当应用差压传感器测量压降时，需用导压管将测压点和差压传感器连接，此时导压管中也要全部充满液体。2.3 两相流
6、摩擦损失的计算两相流的全压损失 被定义为由下列诸要素组成：=+（2.5）式中 为位能损失， 为加速度损失， 为根据平均空隙率而由 确定的值（H为测定区间的高度）。在两相流中，直接测定的值=+(1-)是 ，而摩擦损失 是要从式（ 2.5）算出的值。所以，若 的定义不同， 的值当 然也会不一样。单向流的单位长度的摩擦损失 可表示为：（2.6）=22摩擦系数 ，对于层流：（2.7）=64= 64对于光滑管中的湍流，根据Blasius的公式，可表示为：（2.8）=0.31640.25=0.3164()0.25可是在两相流中，气液间存在相对速度。而且由于断面上有空隙率分布，以式（2.6）的形式来表示摩擦
7、损失时，速度 或比重 的定义未必只有一个，所以对应于不同定义， 的w 数值也就不同。对于 ，1948年曾用对应于 ，忽略了气液w=(0+0)， =0+0(0+0)两相间的相对速度，与单向流相同的值，例如用过0.02 2,0.024的值。=22=(0+0)220+0(0+0)=12(0+0)(0+0)（2.9）两相流摩擦系数最早作系统性的研究是在 1939 年由 Boelter3进行的。但那时只表示出的实验值，没有提出一般的整理法。以后两相流摩擦损失的整理方法离开了 的直接表示 法，构成下面叙述的两相流与单相流的摩擦损失比的表示形式，可以求出后者的值，反过来确定摩擦系数的值也是可能的。2.3.1
8、 基本解析法关于水平管中的层状流、波状流那样的分层流模型的摩擦损失 的一般整理法，最早是于1949年由Lochhart和Martinelli得出的（L-M法 4），这个方法现在还在应用。设流道断面上的压力分布是均匀的，则 等于气相部分的摩擦损失 ，液相部分的摩擦损失 ，即 =图3 分层流模型设图3所示的气液相所占的面积分别为 ， ，又把平均的绝对速度记作 ， (1-) ，则 按通常的摩擦损失公式可表示为： ， ，=122（2.10）=122式中 ， 为各相的摩擦系数， ， 为水力直径： =4() =4(1-) =(1-)=(1-)（2.11）式中 为湿周长度，而 为 中的值。现在考虑气相和液相
9、都是湍流的情况， ， 由式（2.8）（2.12）=0.3146()0.25， =0.3146()0.25把式（2.11）、（2.12）代入式（2.10），整理得：= 0.3164()0.25122=0.3164(0)0.251022()1.2511.75=0()31.75式中 中的值是只有气体单独充满管道断面流动时的假想摩擦损失 。把上述 0的关系写成为=（2.13）2=0=()31.75同样，对于液相则有：（2.14）2=0=()31.75另一方面， 与 之比为0 02=00=0.3164(0)0.251022 0.3164(0)0.251022（2.15）=(00)1.75()0.25()
10、此外，由式（2.13）、（2.14）有：（2.16）2=00=()3()-1.75由式（2.13） （2.16）可以看出，两相流的摩擦损失 、 通过 、 、 、 而与X相 关，Lochhart和Martinelli推出的 或 只表示为一个变量X的函数，提出了图4所示的由实验确定这个关系的各条曲线，在图中对应于各相气液流量，区分是湍流（t）还是层流（v）的场合，提出了气-液的（t-t），（t-v），（v-t）以及（v-v）的各条曲线。图4 两相流的摩擦损失与空隙率的整理方式2.3.2 向高压域的推广上述的L-M曲线是根据大气压附近的实验值制作的，所以它不可能对高压饱和蒸汽-饱和水两相流也原封不动
11、地适用。Martinelli和Nelson为此目的把L-M曲线推广制成图5那样的各种压力的曲线（M-N法 5），这些曲线是以下方法作出的：首先，原封不动地应用饱和蒸汽-饱和水两相流在大气压下的曲线即L-M曲线，然后由如下理论确定临界压力曲线。图5 对应于各种压力的摩擦损失比与空隙率的关系临界压力下的流体流量 （看作气体）加上流量为 的流体（看作气体，是0 0=k0在临界压力下与流量为 的流体具有同一性质的流体）时，其混合物（看作两相流）的摩0擦损失 为：（2.17）=(0+0)22=022式中 是流量为 时的摩擦系数，由式（2.8）： (0+0)=0.3164(0+0) 0.25=0.3164
12、(0)0.25(1+)0.25（2.18） =0(1+)0.25式中 是流量为 时的摩擦系数。把式（2.18）代入式（2.17），则临界压力下的 为：0 0 （2.19）2=0=(1+)1.75另一方面，对于（t-t）领域，由式（2.15）有：x=20.875=()()0.143()0.428（2.20 ）=(1-)()0.143()0.428在临界压力的情况下，因 ，故 ，由此和式（2.19）得：=， =x=00=1（2.21）=(1+1)0.875图5的临界压力下的曲线是式（2.21）的值，对于大气压下的 曲线之间各压力的曲线是参照实验值内插定出的，图中还表示了各种压力对 的关系。(1-)
13、2.3.3 一般管系压力损失的整理方法（Chisholm法 6）Chisholm导出的 与X关系的基本形式说明于下：干度为X的单位断面面积平均重量流量G的两相流的连续性方程式，对于图3.1的模型有：（2.22）x=， (1-)= (1-)（2.23）=+(1-)=式中 是平均速度， 是平均比容。而动量方程为： （2.24）x()+(1-)()=()由式（2.22） （2.24），得 为： （2.25）另=(1-)21+ 1-( + )+ 2(1-)2一方面，关于压力损失，若各相为湍流。粗糙面的场合，则摩擦系数由于 不变而与雷诺数无关，设各相单独充满管断面时的摩擦损失分别为 ， ，则00（2.2
14、6）0=022=222（2.27）0=022=(1-)222由此两式可得：（2.28）=00=1- 而两相流的摩擦损失 也可以与式（2.26）同样的形式表示：（2.29）=22若假设 的值等于单相流的值，把式（ 2.25） （2.28）代入式（2.29），整理得： （2.30）2=0= 2=1+2（2.31）2=0= (1-)2=1+12（2.32）= +在这种粗糙管的假定下， 与X或 与X的关系就由理论上求得了。在式（2.31）中，令 的曲线与图 3.2所示的L-M曲线（它是对于光滑管的）很一致。这样就可以认为=18（2.31）的形式不仅对于粗糙管，而且对于光滑管、弯头、阀门等局部损失也是使用的。此外，根据更详细的解析，C值可表示为：（2.33）=1+(1)(1)0.5 + 式中常数A由实验确定，其结果列在表1中。表13 两相流空泡率的测量 1空泡率（ ）是在任意流通截面上两相混合物中气相所占截面与总截面积之比，或称为气相的时间平均面积份额，有时也称为真实容积含气率。它是描述两相流动特性、测量两相流流量的基本参数之一：= +式中， 为液相所占截面积， 为气相所占截面积， 为两相混合物总截面积。
