DOI： 10.11817/j.issn.1672-7207.2021.05.029

球阀固-液两相流流动特性与压力损失数值模拟

石喜，贡力，陶虎，李露，李江涛

(兰州交通大学 土木工程学院，甘肃 兰州，730070)

摘要：为了研究PVC球阀中固-液两相流流动特性与压力损失变化规律，采用Eulerian多相流模型对DN75球阀进行数值模拟，探讨球阀的流场特性、颗粒分布特征以及压力损失的变化规律，评价其影响因素。研究结果表明：球阀开度对流场影响较大，随着开度减小，阀芯内混合流流速梯度增大，在阀后形成漩涡，漩涡尺度随阀门开度减小而增大；开度不同时，球阀进口断面颗粒体积分数分布基本相似，而阀芯截面和出口截面在小开度时颗粒体积分数分布较乱；混合流压力损失随流速增大而增大，随球阀开度增大而减小，随颗粒体积分数增大而增大；而压力损失随颗粒粒径的变化趋势与阀门开度有关，当开度为70°~90°时，压力损失随着粒径增大而增大；当开度为40°~60°时，压力损失随粒径的增大减小；以显著性水平α=0.05为检验标准，多因素方差分析结果表明颗粒体积分数和颗粒粒径对压力损失的影响不显著；而进口流速和球阀开度对压力损失具有显著影响，其中球阀开度对压力损失的影响最显著。

关键词：Eulerian模型；PVC球阀；流动特性；压力损失

中图分类号：S277.9          文献标志码：A               开放科学(资源服务)标识码(OSID)

文章编号：1672-7207（2021）05-1681-11

Numerical simulation on flow characteristics and pressure loss of solid-liquid two-phase flow in ball valve

SHI　Xi, GONG　Li, TAO　Hu, LI　Lu, LI　Jiangtao

(College of Civil Engineering, Lanzhou Jiaotong University, Lanzhou 730070, China)

Abstract: To study the flow characteristics of solid-liquid two-phase flow and pressure loss of PVC ball valves in PVC ball valve, Eulerian multiphase flow model was adopted to simulate DN75 ball valve. The flow field characteristics, particle distribution and pressure loss of ball valve were discussed, and the influencing factors were evaluated. The results show that the opening of ball valve has great influence on flow field. The gradient of mixed flow velocity in the valve core increases with the decrease of valve opening, and a whirlpool forms behind the valve. The scale of the whirlpool increases with the decrease of valve opening. The distribution of particle volume fraction in the inlet sections of the valve is basically similar at different openings, while in the core and outlet sections, it is disorder at small openings. The pressure loss of mixed flow increases with the increase of flow velocity, and decreases with the increase of ball valve opening and increases with the increase of particle volume fraction. The changes of pressure loss with particle size is related to valve opening, which increases with the increase of particle size when the opening is 70°-90° and decreases when the opening is 40°-60°. Taking the significance level α=0.05 as test standard, multi-factor variance analysis shows that the influence of particle volume fraction and particle size on pressure loss is not significant. However, the flow rate and the ball valve opening have significant effects on the pressure loss, and the valve opening is the most significant factor.

Key words: Eulerian model; PVC ball valve; flow characteristics; pressure loss

含沙水流的水力输送广泛存在于水利、环保和市政等领域。在水力输送管道中，球阀作为重要的水力元件常常用于控制流量、调节压力。含沙水流在球阀中运动时，受阀芯结构的影响，其流场特性和压力损失要比单相流复杂得多，还有可能使球阀产生冲蚀磨损，影响球阀的运行安全，因此，有必要对球阀的固-液两相流特性进行探讨。人们采用试验和数值模拟方法对各类阀门两相流的特性如气-固两相流的流动特性^[1]、冲蚀特性^[2]，气-液两相流的流场特性^[3]、流量系数^[4]、压力降^[5]等进行了研究。在阀门固-液两相流方面，WHEELER等^[6]研究了含沙水流对节流阀的冲蚀特性，发现在节流阀中加入CVD金刚石涂层可提高阀门的抗冲蚀性能；HU等^[7]采用DPM模型研究了节流阀中固-液两相流的抗冲蚀磨损特性；訚耀保等^[8]采用CFD与冲蚀模型对射流管伺服阀多相流中固体颗粒物和油液的运动及冲蚀磨损进行了数值模拟研究；马艺等^[9]利用雷诺应力模型和欧拉-欧拉方法对球阀的开启压差及压力分布特性进行了研究；杨国来等^[10]采用数值模拟对某黑水调节阀在2种开度下的阀门流场进行了分析；余天赐^[11]对球阀固-液两相流的流动及磨损特性进行了研究；崔之健等^[12]对除砂器中携砂液流对针型阀的冲蚀特性进行了数值模拟。可以看出，前人对各类阀门的两相流研究以气-液两相流压力损失变化规律的研究为主，而对球阀固-液两相流压力损失变化规律的研究较少。为此，本文作者以常见的PVC球阀为研究对象，采用Eulerian双流体模型进行数值模拟，得到不同阀门开度时流场变化、颗粒分布特征，以及流速、开度、颗粒体积分数和颗粒粒径的变化对球阀压力损失的影响规律，以期为水力输送管道的设计提供参考。

1  数学模型

1.1　两相流数学模型

选用Eulerian模型进行球阀固-液两相流数值模拟。Eulerian模型属于欧拉-欧拉连续流体模型，将颗粒相作为拟流体，认为颗粒相与流体相是共同存在的连续介质，两相都在欧拉坐标系下计算^[13]。基于Eulerian模型建立的固-液两相流基本方程包括连续性方程和动量方程^[14]。

液相和固相的连续性方程分别为：

(1)

(2)

式中：和分别为液相密度和固相密度；和分别为液相体积分数和固相体积分数，；和分别为液相速度和固相速度；为拉普拉斯算子；为时间。

液相和固相的动量方程分别为：

(3)

(4)

式中：为静压；为固相压力；g为重力加速度；和分别为液相和固相的应力-应变张量；为两相间的体积作用力。

考虑相间作用，由于颗粒密度大于水流密度，因此，可将颗粒所有的虚质量力忽略^[15]。在数值计算时考虑拖曳力和升力的影响，动量传递方程如下^[16]：

(5)

式中：K_sl为相间动量交换系数；C_lift为升力系数，本文计算时取0.5。

1.2　湍流模型

采用混合标准湍流模型进行数值计算。混合标准湍流模型是单相湍流模型的一种扩展，在计算时考虑两相之间的相互作用，其方程和方程分别如下^[16]：

(6)

(7)

式中：为混合相湍流黏性系数；和分别为与湍动能和耗散率对应的Prandtl数，分别取1.0和1.3；为混合相密度；为湍流动能产生项；为混合相速度；和为经验常数，分别取1.44和1.92。

2  几何模型与计算方法

2.1　几何模型

PVC球阀的结构如图1所示，由手柄、阀杆、阀体和球体等组成，其启闭件是1个有孔的球体，与阀杆和手柄连接，通过手柄旋转绕垂直轴线的启闭件来启闭通道。

图1　PVC球阀结构示意图

Fig. 1　Structure sketch of PVC ball valve

本文数值计算采用的球阀为DN75PVC球阀，球体通道内径为57.8 mm。图2所示是球阀开度为50°的几何模型，其中，上游和下游连接管道为DN75的PVC管道，内径为67.8 mm，管道水平放置。计算模型的几何区域取上游长度=1 025 mm(>10d，d为管道的公称外径)，下游长度=1 400 mm(>15d)。为了保证计算流动的稳定性，取压点分别设置在上游10d和下游15d的位置处。计算时，液相为水(温度为20 ℃)，固相为泥沙颗粒，密度为2 650 kg/m³。相间的拖曳力选用Syamlal-O′brien模型，考虑重力的影响，重力方向取Z轴的负方向。

图2　几何模型

Fig. 2　Geometric model

2.2　网格划分及网格无关性检验

采用Gambit软件对计算区域进行网格划分。在球阀阀芯位置采用非结构网格，其余部位采用六面体结构网格。为了保证网格数量对计算结果的影响较小，需进行网格无关性检验。以进口速度为3 m/s、球阀开度为70°、颗粒体积分数为10%、粒径为0.08 mm时的球阀压力损失为指标进行网格无关性检验，结果如图3所示。从图3可以看出：当计算区域的网格数量大于45万个时，可保证计算结果的精度。在不同开度时，本文取计算区域的网格数量为47.2万~48.4万个。图4所示为网格划分图(球阀开度为50°)。

图3　网格无关性检验

Fig. 3　Mesh independence test

图4　网格划分

Fig. 4　Mesh generation

2.3　离散格式与边界条件

采用相间耦合的SIMPLE算法进行迭代求解，对近壁区流动采用标准壁面函数法进行计算。边界条件进口采用速度进口条件，连续相运动速度与颗粒相运动速度相同；出口采用自由出流边界，所有方程迭代残差收敛标准设为10^-4。

2.4　计算模型验证

由于球阀固-液两相流方面的试验数据较少，本文采用GILLIES等^[17]在水平直管段中的水沙两相流试验结果进行验证。由于直管段与球阀中水流运动规律以及产生压力损失的机理相同，而且球阀在开度较大时的流动特性与水平直管的流动特性相似，因此，经过验证后的计算模型可用于球阀固-液两相流计算。GILLIES等^[17]的试验参数如下：管道长度为10 m，内径为103 mm，泥沙颗粒密度为2 650 kg/m³，颗粒粒径为0.09 mm，颗粒的体积分数分别为19%和29%，管壁的粗糙度为2 μm。图5所示为单位长度压力损失的计算结果与试验结果。从图5可以看出：在不同颗粒体积分数和速度时，单位长度压力损失计算结果略大于试验结果，但整体上保持较好的一致性且基本接近，表明Eulerian模型能够较好地模拟固-液两相流动。

图5　压力损失数值模拟结果与试验结果对比

Fig. 5　Pressure loss comparison of numerical simulation results with experimental results

3  结果与分析

3.1　混合相流场分布

当进口流速为2 m/s，颗粒体积分数为10%，粒径为0.08 mm时，不同开度的XY截面混合流速分布和流线如图6所示。由图6可以看出：当球阀开度为90°时，混合流速分布和流线较均匀，由于阀芯通道的直径小于管道直径，导致通道内出现流速梯度，中间速度较大，阀芯出口管壁处出现局部低速区；随着球阀开度减小，阀芯通道进口混合流面积不断减小，在通道内形成较大的流速梯度，中间部位流速较大，而流速梯度随着开度减小而增大。从流线变化可以看出：当球阀开度为90°时，流线平顺；当球阀开度为70°和50°时，由于球阀呈半开状态，混合流进入阀芯后在上侧形成漩涡；而当球阀开度为40°时，阀芯上侧仍有漩涡。从阀后流线变化看，球阀开度小于70°时发生回流，回流尺度随开度减小而增大，这是阀后过流面积减小所致。可以看出，球阀开度对混合流流场分布的影响较大。

图6　混合流速度分布及流线

Fig. 6　Velocity distribution and streamlines of mixed flow

3.2　颗粒体积分数分布

当进口流速为2 m/s、颗粒体积分数为10%、粒径为0.08 mm时，不同开度XZ和YZ截面颗粒体积分数度分布以及YZ截面上Z方向的涡量如图7所示。通过涡量可分析颗粒体积分数的变化情况。从图7可以看出：在进口S1截面，由于重力的影响使颗粒在上部颗粒体积分数较低，下部颗粒体积分数较高，不同开度时S1截面的颗粒体积分数分布相似；涡量的变化与颗粒体积分数变化相反，在涡量大的区域颗粒体积分数较低，而在涡量小的区域颗粒体积分数较高，这是由于S1截面整体流动的湍动能较小，此时涡量的变化对颗粒体积分数分布影响较大，涡量大的地方湍动剧烈使颗粒分散。

在阀芯S2截面，当开度为90°和70°时，颗粒体积分数梯度较小；整体流动的湍动能较小，颗粒分布受涡量变化影响，表现为涡量大的区域颗粒体积分数较低，而涡量小的区域颗粒体积分数较高；在开度为50°和40°时，颗粒体积分数梯度较大，并且分布较乱，出现的颗粒低体积分数区域和高体积分数区域较多，此时涡量的变化也较乱。这是由于在开度为50°和40°时，阀芯内的紊动加大，加速了颗粒与水流的运动，而流通面积的减小使回流区域增多，使得颗粒在部分区域堆积形成高体积分数区，在部分区域冲散形成低体积分数区。这说明在流动紊动剧烈时，颗粒分布不受涡量的影响。

在出口S3截面，当开度为90°和70°时，颗粒体积分数分布相似；当开度为50°时，中间部位颗粒体积分数较低，出现2个低体积分数区域；而当开度为40°时，中间体积分数低，分布相对均匀。从涡量变化来看，当开度为90°时，流动湍动能较小，颗粒体积分数分布受涡量变化的影响较大；而在其分开度时，由于受阀芯扰动的影响，S3截面的颗粒体积分数分布不受涡量影响。

图7　开度不同时颗粒体积分数分布与涡量分布

Fig. 7　Distribution of particle volume fraction and vorticity at different openings

3.3　压力损失变化

3.3.1　球阀压力损失计算　

固-液两相流经过球阀时的压力损失可表示为

(8)

式中：为混合流经过球阀的压力损失；为球阀上下游2个取压点之间总压力损失；为相应直管段混合流的摩阻损失。

通过对相应直管段计算可得到相应工况下混合流的水力坡度，摩阻损失可采用下式计算：

(9)

式中：为相同工况下混合流的水力坡度；为直管段长度。

3.3.2　球阀压力损失随阀门进口流速的变化　

当颗粒粒径为0.08 mm、颗粒体积分数为10%和20%时，球阀两相流压力损失随进口流速的变化如图8所示。从图8可见：在不同开度时，球阀压力损失随流速增大而增大。这是由于随着流速增大，两相流经过球阀时加剧了水流的紊动程度，水流的剧烈紊动加速了颗粒的运动，使得液相损失和固相损失^[18]均增大，因而压力损失增大。从图8还可看出：当开度为90°和70°时，压力损失的变化幅度较小，而在开度为50°时，压力损失的变化幅度较大，这说明当开度较小时，流速对压力损失的影响更加明显。

图8　球阀两相流压力损失随进口流速的变化

Fig. 8　Variation of pressure loss of solid-liquid two-phase flow of PVC ball valves with flow velocity

3.3.3　球阀压力损失随球阀开度的变化　

DN75球阀在颗粒粒径为0.08 mm、颗粒体积分数为10%和20%时，球阀两相流压力损失随球阀开度的变化如图9所示。从图9可以看出：球阀开度对压力损失的影响明显，随着开度增大，压力损失减小；开度越小，压力损失的增加幅度越大。这是由于当开度较小时，进入球阀的流动受阻，阀芯内的流速梯度较大，湍流紊动剧烈，在阀芯和阀后部位出现较大的回流区域，使得液相损失和固相损失都增大，形成了球阀前后较大的压力损失。

图9　液相两相流压力损失随球阀开度的变化

Fig. 9　Variation of pressure loss of solid-liquid two-phase flow of PVC ball valves with ball valve opening

3.3.4　球阀压力损失随颗粒体积分数的变化

当进口流速为2 m/s、颗粒粒径为0.05 mm和0.12 mm、球阀开度为90°和50°时，液相两相流压力损失随颗粒体积分数的变化如图10所示。从图10可以看出：在其他条件不变时，压力损失随颗粒体积分数增大而增大；从变化幅度来看，颗粒体积分数对压力损失的影响不如流速和阀门开度对压力损失的影响明显。这是由于影响压力损失的固相和液相中，在颗粒体积分数较小时以液相损失为主，而固相损失较小；随着颗粒体积分数增大，固相颗粒数量增多，压力损失以固相损失为主，此时，液相损失较小。由于固相颗粒与管壁的摩擦碰撞引起的损失大于液相引起的损失，因此，当颗粒体积分数较大时，混合流的压力损失较大。

图10　液相两相流压力损失随颗粒体积分数的变化

Fig. 10　Variation of pressure loss of solid-liquid two-phase flow of PVC ball valves with particle volume fraction

3.3.5　球阀压力损失随颗粒粒径的变化　

当进口流速2 m/s、颗粒体积分数为10%和20%、开度不同时，球阀压力损失随颗粒粒径的变化如图11所示。从图11可以看出：在开度为70°~90°时，压力损失随着颗粒粒径增大略有增大；而在开度为40°~60°时，压力损失随着颗粒粒径增大而有所减小，整体上变化幅度不大；当颗粒体积分数不变时，增大颗粒粒径使颗粒的数量减小，在开度为70°~90°时，由于球阀开度较大，流动相对平顺，在阀芯和阀后产生的回流区域较少而湍动能较弱，此时，压力损失主要取决于颗粒运动(包括颗粒与颗粒、颗粒与水流的作用)；当颗粒粒径增大时，粒径大而数量少的颗粒运动所消耗的能量大于粒径小而数量多的颗粒运动所消耗的能量，因而，随着粒径增大，压力损失增大；在开度为40°~60°时，由于开度较小，进入球阀的流动紊动剧烈，在阀芯和阀后部位产生的回流区域增多，回流使颗粒与壁面之间的碰撞剧烈程度加大，此时，当颗粒粒径增大时，压力损失可能取决于颗粒与壁面之间的作用，粒径小而数量多的颗粒与壁面之间的碰撞力大于粒度大而数量少的颗粒碰撞力，可能造成小开度时压力损失随粒径的增大而减小。

图11　球阀压力损失随颗粒粒径的变化

Fig. 11　Variation of pressure loss of ball values with particle size

3.4　压力损失变化影响因素评价

为了进一步探讨不同影响因素对球阀压力损失的影响程度，采用无交互作用的单因变量多因素方差分析法进行评价，分析不同影响因素的显著程度。

单因变量多因素方差分析是研究多个因素对同一因变量的影响水平的统计方法，对某一因变量，假设有个影响因素，记为，因素有个水平()，若记该因变量为，则该多因素方差分析的数学模型为^[19]

(10)

式中：为总均值；为因素的水平效应；为随机误差，各个相互独立且服从；；；。

以DN75球阀固-液两相流压力损失为考察对象，采用F检验法对影响压力损失的4个因素(进口流速、球阀开度、颗粒体积分数和颗粒粒径)进行多因素方差分析，结果如表1所示(检验假设所有因素对压力损失没有显著性影响)。

表1　影响因素方差分析结果

Tab. 1　Results of variance analysis for influenced factors

取检验显著性水平=0.05。从表1可以看出：颗粒体积分数和颗粒粒径的相伴概率P分别为0.848和1.000，均大于0.050，表明颗粒体积分数和颗粒粒径对压力损失的影响不显著；而进口流速和球阀开度的相伴概率P均小于0.050，表明进口流速和球阀开度对压力损失具有显著影响。由于球阀开度的相伴概率远小于进口流速的相伴概率，因此，球阀开度对压力损失的影响最显著。

4  讨论

影响球阀固-液两相流压力损失的因素较多，本文主要讨论进口流速、球阀开度、颗粒体积分数和颗粒粒径对DN75PVC球阀压力损失的影响。本研究中，压力损失随流速增大而增大，这与文献[18]所得的水平管道中冰水两相流压力损失一致。在球阀与水平管道中，当速度增大时，水流与颗粒的紊动都加剧，因而损失增大。压力损失随球阀开度增大而减小，这是由于开度增大时流通能力增强，流速梯度减小，压力损失也减小。压力损失随颗粒体积分数增大而增大，这与文献[20]所得的水平管道中固-液两相流结果一致。在水平管道中，颗粒数量增多，增加了颗粒之间的碰撞概率，而且较多数量的颗粒维持悬浮需要提供的能量也多；而在球阀中，由于流体面积较小，颗粒数量较多时增加了颗粒与颗粒、颗粒与壁面的碰撞概率。颗粒粒径的变化对压力损失的影响较复杂。文献[20]得出在其他条件不变时，水平管道中的压力损失随粒径增大而增大；而文献[18]得出冰水两相流压力损失随粒径增大而减小；在球阀中，当开度为70°~90°时，压力损失随粒径增大而增大；当开度为40°~60°时，压力损失随粒径增大而减小。粒径增加可能使压力损失增加，也可能会降低，具体取决于颗粒的运动特性及流通边界。

5  结论

1) 球阀开度对混合流流场分布影响较大。随着阀门开度减小，阀芯内流速梯度增大，在阀后形成漩涡，漩涡尺度随阀门开度减小而增大；当开度不同时，进口截面颗粒体积分数分布基本相似，阀芯中心截面和出口截面在小开度时颗粒体积分数分布较乱。

2) 球阀混合流压力损失随流速增大而增大，随开度增大而减小，随颗粒体积分数增加而增大；而压力损失随颗粒粒径的变化趋势与阀门开度有关，当开度为70°~90°时，压力损失随着粒径增大而增大，当开度为40°~60°时，压力损失随粒径增大而减小。

3) 颗粒体积分数和颗粒粒径对压力损失的影响不显著，而进口流速和球阀开度对压力损失具有显著影响，其中，球阀开度对压力损失的影响最显著。

参考文献：

[1] CHEN Meihua, LU Haifeng, JIN Yong, et al. Experimental and numerical study on gas-solid two-phase flow through regulating valve of pulverized coal flow[J].Chemical Engineering Research and Design, 2020, 155: 1-11.

[2] LIN Zhe, SUN Xiwang, YU Tianci, et al. Gas-solid two-phase flow and erosion calculation of gate valve based on the CFD-DEM model[J].Powder Technology, 2020, 366: 395-407.

[3] MA Yufeng, JI Lijun, ZHANG Jiexu, et al. CFD gas-liquid simulation of oriented valve tray[J]. Chinese Journal of Chemical Engineering, 2015, 23(10): 1603-1609.

[4] BOCCARDI G, BUBBICO R, CELATA G P, et al. Two-phase flow through pressure safety valves: experimental investigation and model prediction[J]. Chemical Engineering Science, 2005, 60(19): 5284-5293.

[5] EDVARDSEN S, DORAO C A, NYDAL O J. Experimental and numerical study of two-phase pressure drop in downhole shut-in valve with unified comprehensive model formulation[J]. Journal of Natural Gas Science and Engineering, 2015, 23: 440-449.

[6] WHEELER D W, WOOD R J K, HARRISON D, et al. Application of diamond to enhance Choke valve life in erosive duties[J]. Wear,2006, 261(10): 1087-1094.

[7] HU Gang, ZHANG Peng, WANG Guorong, et al. Performance study of erosion resistance on throttle valve of managed pressure drilling[J]. Journal of Petroleum Science and Engineering, 2017, 156: 29-40.

[8] 訚耀保, 付嘉华, 金瑶兰. 射流管伺服阀前置级冲蚀磨损数值模拟[J]. 浙江大学学报(工学版), 2015, 49(12): 2252-2260.

YIN Yaobao, FU Jiahua, JIN Yaolan. Numerical simulation of erosion wear of pre-stage of jet pipe servo valve[J]. Journal of Zhejiang University(Engineering Science), 2015, 49(12): 2252-2260.

[9] 马艺, 马中强, 张生昌. 高浓度固液两相流球阀开启压差特性分析[J]. 流体机械, 2014, 42(12): 22-28.

MA Yi, MA Zhongqiang, ZHANG Shengchang, et al. Analysis of different pressure characteristics during the opening process in high-concentration solid-liquid ball valve[J]. Fluid Machinery, 2014, 42(12): 22-28.

[10] 杨国来, 张婧芳, 叶建中. 黑水调节阀固液两相流数值模拟及内流场分析[J]. 液压气动与密封, 2016 (12): 4-6.

YANG Guolai, ZHANG Jingfang, YE Jianzhong. Analysis of blackwater regulating valve in solid-liquid two-phase numerical simulation and internal flow field[J]. Hydraulic Pneumatics & Seals, 2016 (12): 4-6.

[11] 余天赐. 不同开度下球阀内部固液两相流及磨损研究[D]. 杭州: 浙江理工大学, 2019: 1-2.

YU Tianci. Research on solid-liquid two-phase flow and erosion in ball valve[D]. Hangzhou: Zhejiang Sci-Tech University, 2019: 1-2.

[12] 崔之健, 董超云, 黄凯亦, 等. 携砂液流对除砂器针型阀冲蚀作用模拟研究[J]. 西安石油大学学报(自然科学版), 2014, 29(6): 54-57, 64.

CUI Zhijian, DONG Chaoyun, HUANG Kaiyi, et al. Numerical simulation of erosion on desander needle valve for sand carrying flow[J]. Journal of Xi'an Shiyou University(Natural science Edition), 2014, 29(6): 54-57, 64.

[13] 张政, 谢灼利. 流体-固体两相流数值模拟[J]. 化工学报, 2001, 52(1): 1-12.

ZHANG Zheng, XIE Zhuoli. Numerical simulation of fluid-solid two-phase flows[J]. Journal of Chemical Industry and Engineering, 2001, 52(1): 1-12.

[14] 刘勇, 石雯. 管道横向弯曲振动下液固两相流压强损失研究[J]. 水动力学研究进展A辑, 2019, 34(1): 115-121.

LIU Yong, SHI Wen. Study on pressure loss of liquid-solid two-phase flow in the pipeline under transverse bending vibration[J]. Chinese Journal of Hydrodynamics, 2019, 34(1): 115-121.

[15] 杨佳奇, 彭泉霖, 唐颖. 水平圆管内固液两相流的压降及流型分析[J]. 水力采煤与管道运输, 2014(3): 4-7.

YANG Jiaqi, PENG Quanlin, TANG Ying. Pressure and flow pattern analysis of solid-liquid two-phase flow in horizontal pipe[J]. Hydraulic Coal Mining & Pipeline Transportation, 2014(3): 4-7.

[16] 王继红, 张腾飞, 王树刚, 等. 水平管道内固液两相流流动特性的CFD模拟[J]. 化工学报, 2011, 62(12): 3399-3404.

WANG Jihong, ZHANG Tengfei, WANG Shugang, et al. Numerial simulation of liquid-solid slurry flow in horizontal pipeline[J]. CIESC Journal, 2011, 62(12): 3399-3404.

[17] GILLIES R G, SHOOK C A, XU Jihuai. Modelling heterogenous slurry flows at high velocities[J]. Canadian Journal of Chemical Engineering, 2004, 82(5): 1060-1065.

[18] 邓义斌, 王飞显, 范世东. 管道冰水两相流动阻力特性的数值模拟[J]. 武汉理工大学学报(交通科学与工程版), 2015, 39(3): 533-536.

DENG Yibin, WANG Feixian, FAN Shidong. Numerical simulation of resistance characteristics for ice-water two-phase flow in pipes[J]. Journal of Wuhan University of Technology(Transportation Science & Engineering), 2015, 39(3): 533-536.

[19] 张启贤, 刘新平, 林坤. 海外游客特征对在华停留时间影响的多因素方差分析[J]. 佳木斯大学学报(自然科学版), 2009, 27(5): 727-731.

ZHANG Qixian, LIU Xinping, LIN Kun. The multi-factor variance analysis of the influence of overseas visitors' characteristic on the residence time in China[J]. Journal of Jiamusi University(Natural Science Edition), 2009, 27(5): 727-731.

[20] 陈小宁. 基于CFD的固液两相流流型分析及阻力损失模型的研究[D]. 昆明: 昆明理工大学, 2017: 63-67.

CHEN Xiaoning. Flow pattern analysis and resistance loss model of solid-liquid two-phase flow based on CFD[D]. Kunming: Kunming University of Science and Technology, 2017: 63-67.

(编辑  陈灿华)

收稿日期： 2020 -07 -21; 修回日期： 2020 -10 -09

基金项目(Foundation item)：国家自然科学基金资助项目(51969011, 51669010)；兰州交通大学青年科学基金资助项目(2016018) (Project(51969011, 51669010) supported by the National Natural Science Foundation of China; Project(2016018) supported by the Young Scholars Science Foundation of Lanzhou Jiaotong University)

通信作者：石喜，博士，副教授，从事工程水力学与计算水力学研究；E-mail: shixi103@163.com

引用格式： 石喜, 贡力, 陶虎, 等. 球阀固-液两相流流动特性与压力损失数值模拟[J]. 中南大学学报(自然科学版), 2021, 52(5): 1681-1691.

Citation: SHI Xi, GONG Li, TAO Hu, et al. Numerical simulation on flow characteristics and pressure loss of solid-liquid two-phase flow in ball valve[J]. Journal of Central South University(Science and Technology), 2021, 52(5): 1681-1691.