以计算流体动力学CFD软件为工具，通过引入Schnerr-Sauer空化模型，并结合多相流Mixture模型与RNG k-s湍流模型，对两通道非标准孔板流量计在空化发生条件下的内部流场进行数值模拟。考查等效直径比0=0.7时， 流场中蒸汽体积分数和压力分布的变化规律，入口压力对空化数的影响以及雷诺数对流出系数和压力损失的影响。 结果表明：入口压力增加到一定值时，在环隙边缘处首先发生空化现象，并且随着压力的增大，空化发生的区域增大， 空化程度加剧；空化数随入口压力的增大而减小；空化效应对流量计的流出系数的影响较大，而对压力损失的影响较 小。研究结果对孔板流量计的测量误差原因分析和提高测量精度有参考价值。

0.引言

空化是液体所特有的一种复杂的流体动力现象。当流场中某处的局部压力较低时，溶解在液体中的不凝性气体会逸出，当压力降低到对应温度下的 饱和蒸汽压时，液体开始汽化，在局部低压下液体中 瞬间形成大量空泡，这些空泡随液体流会在低压区 时生长、膨胀，而到达高压区时又会收缩、溃灭，这种 空泡爆发性生长、膨胀、收缩、溃灭的整个过程称为 水力空化现象。空化现象的发生有利有弊，目前空 化技术成功地运用在工业废水处理，饮用水消毒，选 矿等方面。在水动力学研宄领域，空化普遍出现在 核动力系统、低温热交换器、液体火箭发动机等工程领 域。当常温流体流经管路、泵、阀门、流量计等各种节 流元件时，节流压降容易导致空化的形成与发展。 空化不仅会使节流元件及下游管道被空蚀损坏、设 备效率降低，而且可能导致流量测量不准确、系统运 行不稳定。所以，对孔板流量计内流体空化流动特性 进行理论与实验研宄具有重要的工业实用价值。

赵奇等设计的两通道非标准孔板流量计与标 准孔板流量计相比，具有临界雷诺数低、永久压降低， 测量稳定性高和节能等优势。但没有考虑在入口压 力较高或流速较大的情况下，节流元件附近可能发 生的空化现象对流量计测量精度会产生影响。本文 在其设计的流量计的基础上，研宄空化数随入口压 力的变化规律以及流出系数和压力损失随雷诺数的 变化规律，并讨论空化的发生对孔板流量计测量精 度的影响，研宄工作对提高测量精度有一定的参考 价值。

1.物理模型与数学模型

1.1几何模型和网格划分

因节流元件为轴对称结构，可简化为二维模拟。 本文研宄的孔板是根据国家标准GB/T 2624—— 2006《用安装在圆形截面管道中的差压装置测量满 管流体流量》进行设计的，其几何结构如图1所 示。管道直径D=100mm，R=50mm，孔板中心孔半径 r1=17.5 mm，环孔内半径r2=38.5 mm，环孔外半径 r3=49 mm，孔板厚度E=3 mm，节流孔厚度e=1 mm，斜 角F=45°，等效直径比炉0.7。孔板上、下游的直管段 长度分别取5D和15D。

利用ICEM CFD进行网格划分，如图2所示，整 体采用四边形结构化网格，从管道两端到孔板逐渐加 密，孔板处进行局部加密，网格总体数量为176262。

1.2 数学模型

采用Schnerr-Sauer空化模型、Mixture模型与 RNG k-s湍流模型联合进行计算。Schnerr-Sauer空 化模型的蒸汽输运方程为

1.3模型参数的设置

近壁区域米用Standard wall function，压力-速度 耦合项采用SIMPLEC算法，动量和湍流动能采用一阶迎风差分格式。边界条件米用压力入口和压力出 口，进口压力的取值范围为1.013 55x105~3.5x105Pa

(绝对压力，以下均为绝对压力），出口压力取值为 0，操作压强取值为1.013 25x105Pa。湍流参数选择湍 流强度和水力直径，汽化压强取值3.166x103Pa，液相 为常温下的水，气相选择水蒸气。以上各个方程的残 差至少达到10-3，保证计算结果充分收敛。

2.数值模拟结果分析

2.1空化现象数值模拟分析

空化数是描述水力空化和空化状态的一个重要 参数，是表征空化特性的无量纲参数。其定义为

空化数的物理意义为:a=抑制空化产生的力/促 使空化出现的力。理论上讲，只要a每1就应该产生 空化，a^0.5就必然产生稳定的空化。即使在环境 压强为几十兆帕时，只要射流速度足够大，就能够 出现空化现象。但是，在实际工程应用中发现空化 数的离散度较大，用空化数来判断是否产生空化并 不准确，所以用空化数判断空化发生没有普遍应用 价值。但由于空化数相关参数容易测量、物理意义明 确，目前仍是粗略判定空化初生和空化程度的常用 方法。

图3所示为模拟得到的空化数随入口压力P 的变化趋势，空化数随入口压力的增大而减小。实际 的空化初生现象一般发生在空化数1.0~2.5之间

空化初生是空穴在极小区域内初次出现的状 态。图4所示为入口压力Pi为1.915x105Pa，入口速度为5.67m/s时，节流孔板前、后区域流体的压力云图。 可以看出在节流孔内及孔板后D/2的区域内发生压 力骤降，在0.65m处压力恢复，稳定在1.01x105Pa 附近。该压力下首次出现空化现象，由图中数据看 出，空化初生时的压力远高于蒸发压力，对应的空化 数为1.33，雷诺数为5.6x105。

图5所示为空化初生时流体中蒸汽体积分数的 等值线。从图5(a)可以看出，空化初生出现在孔板 上游端面壁面处。图5(b)为发生空化区域的局部放 大图，可以看出空化初生是在壁面上开始，在远离壁 面处蒸汽体积分数降低。

随着入口压力增加，空化范围越来越大，空化 区域内蒸汽体积分数也随着增大，当入口压力Pi为 3.5x105Pa，入口速度为8.71 m/s时，模拟所得空化数 为0.44。如图6所示，可以看出在孔板下游0.3m以 内大部分压力区域达到蒸发压力3.166x103Pa，越靠 近孔板的地方蒸汽体积分数越高。

2.2空化现象对孔板流量计测量精度的影响

赵奇等M开展的模拟研宄中没有考虑空化现象 对非标准孔板流量计测量的影响，本文通过改变流体 的不同入口压力，研宄了流出系数C和压力损失 Aw随雷诺数Re的变化情况，并对引入空化模型和 未引入空化模型的模拟结果进行对比。

雷诺数的计算公式为

本文通过改变流体的入口压力得到不同的入口 速度，计算得到不同状态下的雷诺数。流出系数是通 过孔板的实际流量值与理论流量值的比值，是一个 统计量，无法实际测出。它与管道的截面积比、取压 方式、雷诺数及管道情况等很多因素有关。在选用孔 板流量计时，首先应考虑孔板流量计的测量范围位 于流出系数为常数的范围内，以保证流量测量的稳定性。

本文通过模拟获得孔板前后的压降，根据下式进 行计算，得出流出系数。

采用D和D/2取压口取压，上游取压口的间距 为LuLi取0.9D和1.1D之间时无需对流出系数进 行校正，本文Li取1D，此处取压Pi;下游取压口的 距离为匕，因为y3=0.7，，3>0.6，所以当匕取0.49D和 0.51D之间时无需对流出系数进行校正，本文h取 0.5D，此处取压P2。其中必i、L2均为从孔板上游端面 量起。

图7所示为C随Re的变化关系。可以看出，在 两种情况下，C均随Re的增加逐渐减小，并在Re增 加到一定值后趋于常数。就工程应用而言，在选用孔 板流量计时，应确保它的流出系数落在常数区内。由 图7可知，应选择Re在1.3x105~7.2x105之间。Re = 5.6x105时，空化初生。还可以看出，在Re<7.2x105范 围内，引入空化模型的流出系数比未引入空化模型 的流出系数大；在Re >7.2 x105时，未引入空化模型 的流出系数要大。计算结果表明，在流量计测量过程 中，如果流体发生空化现象，则实际流出系数与没有 考虑空化效应的原计算值会有偏差，如果仍按原流出 系数计算流量，则会引起测量误差。当流体Re在 1.3x105~7.2x105范围内，未考虑空化现象的影响，测量值会比实际值偏小。

永久压力损失是表征装置能量消耗的技术经济 指标。压力损失按照GB/T 2624.2——2006的规定 进行计算，其公式为

图8所示为两种情况下压力损失A«与Re的 关系。模拟结果表明，在Re<7.2x105时，引入空化模 型的流量计的压力损失小于未引入空化模型的；在 Re>7.2x105时，引入空化模型的流量计的压力损失大 于未引入空化模型的。造成此种结果的原因可能如 下：在该流量计的流体流动中压力损失表现为静压 能转化为内能，该过程中，空化消耗能量为叫，汽泡 的产生使流体与管壁摩擦耗能减少量为叱。当《i>«2 时，表现为压力损失增大，对应Re>7.2x105区域；当 «1<«2时，表现为压力损失减小，对应Re<7.2x105区 域。由图可知，虽然空化的发生对流量计的压力损失 有影响，但是影响不大。

3.结束语

通过引入空化模型对两通道非标准孔板流量计 的流场进行模拟，得出以下结论：

1）随着入口压力的增加，雷诺数逐渐增大，空化 数不断减小，在低压下空化数的变化较快，在高压下 空化数的变化较慢，说明空化初生现象容易在低压 下发生。因此在进行低压、高速的流体测量时更应该 注意空化现象的发生。

2）当压力达到一定值时，空化初生发生在孔板 的上游端面靠近壁面的凸起处，如果流量计长时间 在这样的条件下使用，汽蚀作用有可能造成流量计 节流件的磨损，进而影响测量的精度。

3）空化效应对流量计的测量精度有影响，在一 定的雷诺数范围内，空化效应会引起流出系数的变 化，如果在实际测量时未考虑空化效应的影响，则会 造成流量计的测量误差。