1. 引言

石油是全球主要能源之一，广泛应用于发电、供暖、工业和交通。由于石油资源分布不均，需要通过各种运输方式进行输送，目前采用较多的运输方式为管道运输[1][2]。在整个石油管道运输系统中，阀门作为石油管道运输中的“咽喉”对系统的稳定运行起到关键作用。在众多阀门中，球阀因可靠性高、启闭速度快等优点在管道运输中被大规模使用[3][4]。

球阀由流道结构的不同可分为全通径式球阀和缩径式球阀，在全通径式球阀中公称直径等于通道直径，在缩径式球阀中公称直径大于通道直径。目前大多数学者针对全通径式球阀进行内部流场研究，但对缩径式球阀内部流场研究较少。Cheng[5]等基于k-ε湍流模型对全通径式V型球阀内部流场进行了数值模拟，研究了不同开度和阀芯锥角下的内部流动特性，为天然气管道高效运行提供技术支撑。王佳琪[6]等应用CFD软件对阀芯锥角分别为30˚，60˚，90˚的全通径式V型球阀进行数值模拟，对其内部流场进行了研究，为V型球阀的设计与优化提供了参考。梁潇[7]运用Fluent针对球阀管路系统，通过改变阀门开度、管道构型等因素对其启闭过程的流场变化进行研究，为提高球阀管路系统稳定性提供依据。Ma[8]等利用自定义UDF技术和k-ε湍流模型，对全通径式球阀在变转速规律和不同开度下的内部流场进行了瞬态数值模拟，为球阀选用提供理论支撑。

本文的目的在于拓展对缩径式球阀内部流场的研究，为缩径式球阀设计与优化提供参考。本文提出改变缩径式球阀的通道直径和阀门开度，利用Fluent软件模拟分析速度云图、压力云图和流通能力，对缩径式球阀内部流场进行深入研究。

2. 计算模型

2.1. 几何模型

缩径式球阀是一种常用的流体控制装置，如图1所示其主要结构由阀体、阀杆、球体等部分组成。本文针对缩径式球阀内部流场研究运用NX12.0对其进行结构简化以有利于数值模拟的收敛，为了降低进出口边界对球阀附近区域流场的影响，在阀前与阀后分别添加5D和10D的管道(D为阀门公称直径)如图2所示[9][10]。

Figure 1.Structure diagram

图1.结构图

Figure 2.Simplified 3D structure diagram

图2.三维结构简化图

2.2. 网格划分

本文将简化后的几何模型导入Design Modeler中进行流体域的抽取，利用Fluent Meshing对抽取得到的流体域进行网格划分，对阀门球体以及阀门前后范围区域流场采用加密处理。体网格类型选择poly-hexcore以确保数值模拟的准确性、计算效率和收敛性。图3为缩径式球阀在30˚开度下的网格划分情况。

Figure 3.Mesh model

图3.网格模型

2.3. 边界条件设定

本文为探究缩径式球阀内部流场分布，研究范围涵盖了阀门不同开度(30˚、45˚、60˚)和不同通道直径(66 mm、76 mm、86 mm)的变化对内部流场的影响。计算涉及的边界条件：湍流模型设置为Realizable k-ε模型,流体密度为1350 kg/m³，流体粘度为1 kg/m∙s，流体速度为5 m/s，入口设置为速度入口，出口设置为压力出口，求解方法为SIMPLE[11][12]。

2.4. 准确性验证

为确保本文选取的湍流模型及求解方法的准确性，本文将进行准确性验证，本次准确性验证主要采用对比验证的方法。主要步骤如下：首先寻找对比验证所需的文献；接着，建立必要的对比验证几何模型；然后，应用本文选取的湍流模型及求解方法，并使用文献中的实验参数进行数值模拟；最后，将模拟结果与文献中的数据进行对比，从而验证所选取的湍流模型及求解方法的准确性。

本文对姬鹏[13]的实验数据进行数值模拟，采用与其相同的流体参数和工作参数，并选用本文的湍流模型和求解方法，通过计算对比速度分布和压力分布的情况，从而验证本文湍流模型和求解方法的准确性。本次准确性验证的几何模型如图4所示。

Figure 4.Geometric model

图4.几何模型

在本次准确性验证中，速度分布和压力分布的情况如图5所示。通过对比分析可以得出，本次准确性验证中的速度分布和压力分布的情况与文献中的速度分布和压力分布的情况基本一致。这一结果表明，本文选取的湍流模型与求解方法是准确的。

Figure 5.Accuracy verification diagram

图5.准确性验证图

3. 流场仿真模拟结果及分析

3.1. 不同开度对内部流场的影响

当流体速度为5 m/s，通道直径为76 mm时，截取XY平面观察阀门不同开度下的内部流场分布情况，图6由上至下分别为阀门开度在30˚、45˚以及60˚时的压力分布情况。由图6可知，阀内最大压力在流体进口处产生，球体内的上部压力小于下部压力，由于流体的流动在球体出口处的压力逐渐减小，最小压力产生在球体出口与管道连接处的下部。当阀门开度在30˚向60˚变化时，随着阀门开度的增大，阀内最大压力呈现一个逐渐增大的过程，而阀内最小压力的变化呈逐渐减小的趋势。

Figure 6.Pressure distribution diagram

图6.压力分布图

在图7和图8中，由上至下阀门开度分别30˚、45˚以及60˚。由速度分布图(图7)可知，当流体流入管道和流出管道时流体速度较为均匀，当流体流经球体入口处和出口处时由于流道的改变，流体速度急剧增大从而形成高速使流体喷向球体内部和管道内壁，同时随着阀门开度的增大导致最大流速的增大，该现象将逐渐加剧。

Figure 7.Velocity distribution diagram

图7.速度分布图

由粒子迹线图(图8)可知，流体在流入球体时球体内部会产生湍流，随着开度增加湍流现象将逐步加剧，流体流出球体时球体出口处与管道连接处的下部会产生涡旋，随着开度增加涡旋将愈发增大。

Figure 8.Particle trace diagram

图8.粒子迹线图

3.2. 不同通道直径对内部流场的影响

在缩径式球阀设计中，通道直径的大小是一个至关重要的参数，它直接影响着阀门的流体控制性能、压力损失等。因此，进行不同通道直径下的流场研究，有助于全面了解阀门的性能特征。

图9表明阀门开度为45˚，流速在5 m/s的情况下，不同通道直径对阀门内部压力分布的影响。在图9中，由上至下通道直径分别为66 mm、76 mm以及86 mm。如图9所示，在管道入口至球体入口这一部分区域为高压区，球体部分为中压区，球体出口至管道出口这一部分为低压区。当缩径的直径由66 mm向86 mm转变的过程中，高压区和中压区的压力逐渐减小，低压区的压力逐渐增大。

在图10中，由上至下通道直径分别为66 mm、76 mm以及86 mm。如图10所示，在缩径式球阀的球体进口、球体出口以及球体内部，均产生了剧烈的速度变化，成为了出现速度梯度的主要区域。在通道直径由66 mm向86 mm变化的过程中，由于通道直径增加，阀内最大速度减小，速度梯度减小。

Figure 9.Pressure distribution diagram

图9.压力分布图

Figure 10.Velocity distribution diagram

图10.速度分布图

由图11可知，缩径式球阀由上至下通道直径分别为66 mm、76 mm以及86 mm。如图11所示，在流体流动过程中会产生湍流和涡旋，随着通道直径的增加，最大速度的减小，球体内部的湍流逐渐减缓，球体出口处与管道连接处下部产生的涡旋逐渐缩小。

Figure 11.Particle trace diagram

图11.粒子迹线图

4. 缩径式球阀流通能力

4.1. 流量系数

流量系数在阀门设计中至关重要，它反映了流体介质、阀门结构、阀门开度等因素对阀门流通能力的影响程度，因此它被视为衡量阀门流通能力的关键指标[14]。流量系数的大小直接反映了阀门的流通能力的优劣，通常来说，流量系数越大，阀门的流通能力就越好，流量系数越小则意味着阀门的流通能力越差[15][16]。阀门的流量系数理论公式[17]如下：

$\text{Δ}P=P_1-P_2$(1)

$C=Q\sqrt{\frac{\rho }{\Delta {\rm P}}}$(2)

故：

$C=Q\sqrt{\frac{\rho }{{{\rm P}}_1-{{\rm P}}_2}}$(3)

式中， ${{\rm P}}_1$为阀门进口压力，Pa； ${{\rm P}}_2$为阀门出口压力，Pa； $\text{Δ}P$为压差，Pa； $C$为流量系数，无量纲； $Q$为体积流量，m³/h；ρ为流体密度kg/m³。

根据数值模拟所得数据，对缩径式球阀的流量系数进行计算，计算结果如表1和表2所示。由表1和表2可知，在缩径式球阀的流量系数计算中，缩径式球阀的不同通道直径和开度的变化，会导致相应的流量变化。

由表1可知，当通道直径保持为76 mm不变时，阀门开度由30˚向60˚逐渐增大，流量系数将随着阀门开度的增大而减小，阀门流通能力变弱。

Table 1.Flow coefficient

表1.流量系数

由表2可知，当阀门开度保持为45˚不变时，阀门通道直径由66 mm向86 mm增加，流量系数将随着阀门通道直径的增加而增大，阀门流通能力变强。

Table 2.Flow coefficient

表2.流量系数

4.2. 流阻系数

流阻系数是衡量流体在通过阀门或管道时所产生能量损失的重要参数[18]。对流阻系数进行深入研究和分析有助于我们更好地理解阀门结构对流体流动的影响，并且可以借此优化设计以减少能量损失。优化流阻系数不仅可以改善流体控制系统的性能，还可以降低系统的能耗和运行成本。一般情况下，流阻系数越大则意味能量损失越大、流通能力越差，流阻系数越小则表明能量损失越小、流通能力越好[19][20]。阀门流阻系数的计算公式[21]如下：

$\zeta =\frac{2\text{Δ}{\rm P}}{\rho v^2}$(4)

结合式(1)则：

$\zeta =\frac{2\left(P_1-P_2\right)}{\rho v^2}$(5)

式中， ${{\rm P}}_1$为阀门进口压力，Pa； $P_2$为阀门出口压力，Pa； $\text{Δ}{\rm P}$为压差，Pa； $\zeta$为流阻系数，无量纲； $v$为流速，m/s； $\rho$为流体密度kg/m³。

结合式(5)和数值模拟所得数据，对缩径式球阀进行流阻系数计算，计算结果如表3和表4所示。根据表3中的数据可得，当通道直径保持为76 mm不变时，随着阀门开度从30˚逐渐增大至60˚，流阻系数呈现逐渐增大的趋势。这表明随着阀门开度的增加，阀门的流通能力逐渐减弱，且在开度从45˚向60˚旋转时，流阻系数急剧增大。

Table 3.Flow resistance coefficient

表3.流阻系数

由表4可知，当阀门开度保持为45˚不变时，随着通道直径从66 mm增加到86 mm，流阻系数则呈现逐渐减小的趋势。这意味着阀门通道直径的增加会降低流体通过阀门时的阻力，从而增强阀门的流通能力。

Table 4.Flow resistance coefficient

表4.流阻系数

故在设计缩径式球阀时，需要综合考虑通道直径和开度的变化，以全面评估阀门的流通能力。

5. 总结

本文通过数值模拟对缩径式球阀的内部流场进行了深入研究，并得出了以下结论。

首先，对于不同开度下的影响，研究发现随着阀门开度增大，阀内最大压力逐渐增大，最小压力逐渐减小。而阀门开度的增加导致流体速度梯度加剧，尤其在球体入口和出口处，流速急剧增大，形成高速喷流现象。此外，随着开度增加，涡旋现象和球体内部湍流现象加剧，流体流动变得更加复杂。

其次，针对不同通道直径的研究发现，通道直径的增加使阀内高压区和中压区的压力逐渐减小，低压区的压力逐渐增大。同时，通道直径的增加导致了速度梯度减小，使得球体内部湍流逐渐减缓，球体出口处产生的涡旋逐渐减小。

最后，在缩径式球阀的流量系数和流阻系数计算中，结果显示阀门的流量系数随着阀门开度增大而减小，随着通道直径增大而增大。而流阻系数则呈现相反的趋势，随着阀门开度增大而增大，随着通道直径增大而减小。

综上所述，缩径式球阀的内部流场受到阀门开度和通道直径的影响，需要综合考虑这些因素来评估阀门性能。在实际应用中，该研究结果对于优化阀门设计、提高石油管道运输系统的稳定性和效率具有重要意义。

参考文献