具有肋条结构的气力输送弯管抗冲蚀优化设计

doi:10.11902/1005.4537.2022.231

具有肋条结构的气力输送弯管抗冲蚀优化设计

郭姿含, 张军, 李晖^,

集美大学海洋装备与机械工程学院福建省能源清洁利用与开发重点实验室厦门 361021

Optimal Design for Anti-erosion of Pneumatic Conveying Elbow with Rib Structure

GUO Zihan, ZHANG Jun, LI Hui^,

Fujian Province Key Laboratory of Energy Cleaning Utilization and Development, School of Marine Equipment and Mechanical Engineering, Jimei University, Xiamen 361021, China

通讯作者: 李晖，E-mail：lihui@jmu.edu.cn.，研究方向为多相流体流动

收稿日期: 2022-07-15 修回日期: 2022-08-03

基金资助:

福建省自然科学基金. 2022J01334
福建省自然科学基金. 2020J01694

Corresponding authors: LI Hui, E-mail:lihui@jmu.edu.cn

Received: 2022-07-15 Revised: 2022-08-03

Fund supported:

Natural Science Foundation of Fujian Province. 2022J01334
Natural Science Foundation of Fujian Province. 2020J01694

作者简介 About authors

郭姿含，女，1997年生，硕士生

摘要

弯管冲蚀是不可忽视的重要问题，直接关系到管路输送系统的安全运行及管道的使用寿命。针对这一问题，提出一种四边形、等腰梯形和等腰三角形3种横截面形状的肋条，分别安装在弯头外径方向不同位置，并考虑均匀安装多个肋条对冲蚀的影响。采用CFD-DPM方法对所提出的具有肋条结构的弯管抗冲蚀特性进行数值模拟。模拟结果表明：肋条安装在颗粒壁面第一次碰撞之前，一定程度上抑制冲蚀，且肋条背部形成低速逆流循环区，保护该区域壁面。3种不同横截面形状的肋条中，抗冲蚀作用最佳的为等腰三角形肋。肋深越大，保护范围越大，但影响颗粒碰撞角度，增大颗粒与其碰撞频率，并非深度越大抗冲蚀性能越佳。弯头部分均匀分布多个肋条也具有明显的抗冲蚀特性。其结论可为弯管的抗冲蚀优化设计提供新的设计方案。

关键词： 气力输送 ; 弯管 ; 肋条 ; 截面形状 ; 冲蚀 ; 数值模拟

Abstract

Elbow erosion is an important matter related directly to the safe operation of pipeline conveying system and the service life of pipeline. To relieve the harmful effect of this kind of erosion, ribs with different shapes such as quadrilateral, isosceles trapezoid and isosceles triangle sections were designed and prepared, which then were installed at different designed positions in the inner edge on the half side of the elbow with large curvature radius, and the effect of evenly installing multiple ribs on the erosion process was also considered. CFD-DPM method was used to simulate the erosion resistance of the elbow with rib structure. The results show that if the rib is installed just behind the leading edge of the particle impact on the elbow wall, the existence of rib can change the trajectory of the particles, inhibit the erosion to a certain extent, and induce the formation a low-speed counter current circulation zone on behind the rib to protect the elbow wall in this area. Among the three proposed ribs of different cross sections, the isosceles triangular ribs have the best anti-erosion effect. The greater the rib thickness is, the greater the protection range is. However, the rib thickness does affect the impact angle of particles and increase the collision frequency of particles, therefore, a proper rib thickness may be carefully selected for acquiring the better anti-erosion performance. The isosceles triangular ribs with a rib thick in 6 mm have the best anti-erosion effect at θ=25°, which is 43.63% higher than that of ordinary curved tubes. Many ribs evenly installed on the elbow also have obvious anti-erosion effect, in fact, installation of 8 ribs with isosceles triangular cross section have the best anti-erosion performance. The conclusion can provide a new reference for the anti-erosion optimization design of elbow.

Keywords： pneumatic conveying ; elbow ; rib ; section shape ; erosion ; numerical simulation

PDF (5935KB) 元数据多维度评价相关文章导出 EndNote| Ris| Bibtex 收藏本文

本文引用格式

郭姿含, 张军, 李晖. 具有肋条结构的气力输送弯管抗冲蚀优化设计. 中国腐蚀与防护学报[J], 2023, 43(3): 525-534 DOI:10.11902/1005.4537.2022.231

GUO Zihan, ZHANG Jun, LI Hui. Optimal Design for Anti-erosion of Pneumatic Conveying Elbow with Rib Structure. Journal of Chinese Society for Corrosion and Protection[J], 2023, 43(3): 525-534 DOI:10.11902/1005.4537.2022.231

ASTM国际G40-13标准认为冲蚀是由于固体表面与流体、多组分流体、撞击液体或撞击固体颗粒之间的机械相互作用而导致原始材料从固体表面逐渐流失的现象。研究表明，一个工业化国家的冲蚀磨损成本约占其国民生产总值的1%至4%^[1]。冲蚀经常发生在采矿、冶金、石油、天然气等工业生产中，特别是泵、旋风分离器、燃气轮机叶片、气力输送管道等部件中^[2~5]，其中气力输送管道中的冲蚀破坏更为严重^[6]。气力输送管道主要由水平管、垂直管和弯管组成。弯管作为改变流动方向的部件，起着连接水平管和垂直管的作用，更容易发生冲蚀。研究表明，弯管冲蚀破坏的程度约为直管段的50倍^[7]。这种冲蚀破坏不仅浪费材料，消耗能源，降低设备效率，而且加速设备故障，降低弯管的使用寿命，对管路系统的安全运行造成威胁，从而酿成更大的经济损失。目前常用计算流体力学-离散相方法 (CFD-DPM) 和计算流体力学-离散元方法 (CFD-DEM) 方法研究冲蚀。CFD-DPM仅适用于稀相流，而CFD-DEM方法在计算密相流方向有更大的优势。合理的弯头结构以提高弯管的抗冲蚀性能在工业生产中具有重要意义，研究者对弯管结构进行了大量改进以提高弯管的抗冲蚀性能。Duarte等^[8]提出一种沿水流方向扭转管壁的新颖管壁设计，这种结构会在弯头上游产生一个旋流，最大冲蚀率减小33%。San等^[9]提出一种涡流室设计，并分析流动Reynold数、颗粒Stokes数和涡流室大小对冲蚀的影响，结果表明：与普通半径的弯管相比，涡流室总能减少侵蚀。以上改变弯管结构的方法虽能一定程度上减少冲蚀，但改变了弯管原本结构，影响了管道整体布局，且冲蚀仍发生在弯管壁面上，发生冲蚀破坏时候仍需替换弯管^[10]。相比以上直接改变弯管结构的方法，内壁面增加肋条是更直接有效的方法，冲蚀破坏中只有肋条需被更换^[11]。Yao等^[12]通过实验方法分析弯管加横截面为四边形肋条的抗冲蚀性能，结果表明：肋条可提高弯管的抗冲蚀性能，矩形肋条相比方形肋条抗冲蚀性能更佳。Fan等^[13]研究了横截面为矩形、等腰三角形、等腰直角三角形3种不同截面形状肋条的抗冲蚀性能，结果表明：3种不同形状肋均提高弯管的抗冲蚀性能，等腰直角三角形肋条的抗冲蚀能力最强。以上均在横截面为正方形弯管中均匀分布多个肋条，未见圆形弯管且未考虑单个肋在弯管不同角度的抗冲蚀分布。再者，Zhu等^[14]虽将单个轴向角为5°，环形角为180°的肋条安装在90°圆形弯头外径方向不同位置上，但未考虑肋条横截面不同形状截面参数及多个肋条下的抗冲蚀特性。

为此，本文采用CFD-DPM方法，研究气力输送弯管中具有四边形、等腰梯形、等腰三角形3种横截面形状肋条的抗冲蚀特性，讨论肋条不同位置、形状、尺寸、数量下的冲蚀规律，扩展肋条结构在弯管抗冲蚀设计中的应用。其结论可直接应用于弯管的抗冲蚀优化设计。

1 理论模型

1.1 连续相模型

连续相采用Euler形式的非定常Reynolds平均Navier-Stokes方程(URANS)求解，其质量守恒方程和动量方程如下^{[9, 15,16]}：

\nabla (ρ \vec{v}) = 0

(1)

\frac{\partial (ρ \vec{v})}{\partial t} + \nabla (ρ \vec{v} \vec{v}) = - \nabla p + \nabla [μ (\nabla \vec{v} + \nabla {\vec{v}}^{T})] + ρ \vec{g} + S_{m o m}

(2)

其中，ρ表示气体连续相密度，kg/m³；v是流体速率，m/s；p为压力，Pa；g是重力加速度，m/s²；μ为气体粘度，Pa·s；S_mom是附加源项。

因URANS方程中引入了Reynolds应力项，因此需采用湍流模型来封闭控制方程^[14]。两方程k-ε模型主要包括标准k-ε模型、RNG k-ε和Realizable k-ε模型。因Realizable k-ε模型带涡旋修正，能提高剪切流和旋流的精度，故本文采用Realizable k-ε模型，其方程为^[17]：

\begin{array}{l} \frac{\partial (ρ k)}{\partial t} + \frac{\partial (ρ k u_{i})}{\partial x_{i}} = \frac{\partial}{\partial x_{j}} [(μ + \frac{μ_{t}}{σ_{k}}) \frac{\partial k}{\partial x_{j}}] + \\ G_{k} + G_{b} - ρ ε \end{array}

(3)

\begin{array}{l} \frac{\partial (ρ ε)}{\partial t} + \frac{\partial (ρ ε u_{i})}{\partial x_{j}} = \frac{\partial}{\partial x_{j}} [(μ + \frac{μ_{t}}{σ_{ε}}) \frac{\partial ε}{\partial x_{j}}] + \\ ρ C_{1} S ε - ρ C_{2} \frac{ε^{2}}{k + \sqrt[]{v ε}} + C_{1 ε} C_{3 ε} \frac{ε}{k} G_{b} \end{array}

(4)

C_{1} = m a x (0.43, \frac{η}{η + 5})

(5)

η = S \frac{k}{ε}

(6)

S = (2 S_{i j} \cdot S_{i j})^{\frac{1}{2}}

(7)

式中，G_k和G_b分别表示由平均速度梯度和升力引起的湍动能产生项；σ_k和σ_ε 分别为湍动能k和湍动耗散率ε对应的Prandtl数； $η$ 为相对应变参数； $S$ 为应变率；C_μ =0.09，C₁_ε =1.44，C₂=1.9，C₃_ε =0.9，均为标准经验值。

1.2 离散相模型

利用基于牛顿第二定律的拉格朗日方法追踪离散颗粒，颗粒在流场中受力示意图如图1所示，其控制方程如下^[16,18,19]：

图1

新窗口打开| 下载原图ZIP| 生成PPT

图1 颗粒受力示意图

Fig.1 Schematic diagram of particle stress

m_{p} \frac{d v_{p}}{d t} = \vec{F_{D}} + \vec{F_{B}} + \vec{F_{V M}} + \vec{F_{P}} + \vec{F_{G}}

(8)

\begin{matrix} \vec{F_{D}} = \frac{18 μ}{ρ_{p} d_{p}^{2}} \frac{C_{D} R e_{p}}{24} \end{matrix} (\vec{v} - \vec{v_{p}})

(9)

\vec{F_{B}} = (ρ_{p} - ρ) \vec{g} \frac{π d_{p}^{3}}{6}

(10)

\vec{F_{V M}} = \frac{π d_{p}^{3}}{12} ρ \frac{d}{d t} (\vec{v} - \vec{v_{p}})

(11)

\vec{F_{P}} = ρ \frac{π d_{p}^{3}}{6} \vec{v_{p}} \nabla \vec{u}

(12)

式中，F_D和F_B分别为颗粒所受的流体曳力和升力，F_P为连续相流场压力梯度对颗粒产生的压力梯度力，F_VM为颗粒表观质量效应产生的虚拟质量力，N；ρ_p为颗粒密度，kg/m³；d_p为颗粒的直径，m；v_p为颗粒速度，m/s；Re_p为颗粒Reynolds数，定义为：

\begin{matrix} R e_{p} = \frac{ρ d_{p} |\vec{v} - \vec{v_{p}}|}{μ} \end{matrix}

(13)

1.3 冲蚀模型与壁面碰撞方程

由于冲蚀过程复杂，受很多因素的共同影响，E/CRC冲蚀模型^[20,21]考虑了硬度和颗粒形状因素，故本文采用E/CRC冲蚀模型，其冲蚀模型如下：

\begin{matrix} R_{e r o s i o n} = \sum_{p = 1}^{N} (\frac{{\dot{m}}_{p} C {(Β Η)}^{- 0.59} F_{S} f (θ) v_{p}^{n}}{A_{f}}) \end{matrix}

(14)

f (θ) = 5.4 θ - 10.11 θ^{2} + 10.93 θ^{3} - 6.33 θ^{4} + 1.42 θ^{5}

(15)

式中，C=2.17×10^-7为经验常数；BH为靶材的布氏硬度；F_S为颗粒形状系数，球形为0.2；n为经验常数，取值2.41；A_f为受冲蚀区域面积，m²；θ为侵入角，rad；f(θ)为侵入角函数。

粒子与壁面碰撞会发生能量传递与损失，这里采用法向恢复系数e_n和切向恢复系数e_t表示^[22]：

\begin{matrix} \{\begin{matrix} e_{n} = 0.993 - 1.76 θ + 1.56 θ^{2} - 0.49 θ^{3} \\ e_{t} = 0.998 - 1.55 θ + 2.11 θ^{2} - 0.67 θ^{3} \end{matrix} \end{matrix}

(16)

2 具有肋条的弯管模型

2.1 具有肋条的弯管几何模型

利用Solidworks三维绘图软件在弯头处建立图2a所示截面为3种不同形状肋条结构。三种肋条截面的肋宽为L，肋深为H，环形角为180°，肋条安装在沿弯头曲率方向θ角度处，θ为弯头入口处到肋条前端的距离。如图2b所示为考虑弯管中安装多个肋条的情况，肋条均匀分布在90°弯头外壁面，数量为N。这些具有肋条结构的弯管仅弯头部分不同于普通弯管，垂直上升段和水平出口段与普通弯管相同。弯管材料为Al，直径D=50 mm，弯径比R/D=2，为保证管内流动充分发展，垂直入口段和水平出口段设置为10D=500 mm。

图2

新窗口打开| 下载原图ZIP| 生成PPT

图2 具有肋条的弯管几何结构

Fig.2 Elbow geometry with rib of three different cross sections (a) and multiple rib (b)

2.2 网格划分及流场条件

对普通弯管和具有肋条的弯管进行网格划分，如图3所示，采用质量较好的多面体网格进行划分，并添加边界层对流动的影响。为提高计算精确性，本文以普通弯管为例，采用最大冲蚀率为指标，进行图4所示的网格独立性验证。由图4可知，当网格数量高于25×10⁴时，最大冲蚀率值逐渐趋于稳定，故考虑计算精度和效率要求，最终网格数量为293312，网格质量较好。

图3

新窗口打开| 下载原图ZIP| 生成PPT

图3 网格划分

Fig.3 Grid division

图4

新窗口打开| 下载原图ZIP| 生成PPT

图4 网格独立性验证

Fig.4 Grid independence verification

气力输送中气固两相采用基于Euler-Lagrange方法的DPM模型。连续相为空气，密度1.18 kg/m³,粘度1.8×10^-5 Pa·s；离散相为沙粒，密度1650 kg/m³，粒径为单一粒径100 μm，颗粒质量流量为0.02 kg/s。考虑连续相和颗粒相的双向耦合作用。边界条件设置中，连续相入口为速度入口，流速30 m/s，出口为大气压压力出口，壁面采用无滑移边界条件，近壁面处采用Launder和Spalding^[23]提出的标准壁面函数；离散相为惰性颗粒垂直入口表面射入，入口出口边界条件为escape逃逸类型，壁面为reflect反射类型，假设入口处颗粒速率与连续相入口速率相同。

2.3 模型有效性验证

进行数值模拟之前，验证模拟方法及模型的有效性是必要的。采用Zhu等^[14]实验结果与本文模拟值对比分析。利用本文模拟方法模拟相同实验条件下的冲蚀率，得到图5所示管道外壁面中心线30°实验值与模拟值的对比图。可知，模拟与实验冲蚀率具有较好的一致性，验证了模拟方法及模型的可靠性。

图5

新窗口打开| 下载原图ZIP| 生成PPT

图5 实验值^[14]与模拟值对比图

Fig.5 Comparison between experimental and simulated values

3 结果与讨论

3.1 不同肋条形状对冲蚀的影响

以H=L=5 mm的3种截面肋条结构为固定条件，模拟分析3种不同横截面形状肋条安装在弯头不同角度下的冲蚀规律。图6为普通弯管在该条件下的冲蚀分布及颗粒轨迹。可知普通弯管的冲蚀形状呈“V”形分布，冲蚀最严重部位出现在沿弯头曲率方向35~60°附近，最大冲蚀率为1.1×10^-4 kg/(m²·s)。由于颗粒轨迹对冲蚀形状的形成有重要影响，可通过颗粒轨迹分析冲蚀形状。由对应的颗粒轨迹图可知颗粒在35°左右第一次与管壁碰撞形成椭圆形冲蚀严重区，之后分成两股相反对称的颗粒流，直至与壁面发生第二次碰撞，便形成图示的“V”形分布，这与Peng等^[18]和Pereira等^[24]得到的结果一致。

图6

新窗口打开| 下载原图ZIP| 生成PPT

图6 普通弯管冲蚀分布及对应颗粒轨迹

Fig.6 Erosion distribution and particle trajectory of ordinary elbow

图7分别为3种截面形状的肋条安装在弯头不同角度θ下的最大冲蚀率变化，红色点线代表普通弯管的最大冲蚀率。由图7a可知，当肋条截面形状为四边形时，肋条部分的最大冲蚀率始终低于弯管部分，随弯管趋势变化。肋条位于θ=35°之前时，均具有不同程度的抗冲蚀效果，特别是θ=15°时，弯管的最大冲蚀率为6.89×10^-5 kg/(m²·s)，相比普通弯管减小37.36%；肋条安装在θ=35°~60°时，正处于普通弯管冲蚀严重区，肋条的安装反而加大了冲蚀；θ大于60°时，颗粒逐渐流出弯管，肋条的作用逐渐削弱，不具有明显的抗冲蚀性能。由图7b可知，当肋条截面形状为等腰梯形时，肋条部分的最大冲蚀率在θ=15°~35°出现高于弯管部分的现象，因单个肋条厚度较大，便于更换，充分发挥了肋条牺牲自身保护管壁的作用，具有较好的抗冲蚀效果。肋条位于θ=40°之前时，均具有不同程度的抗冲蚀效果，特别是θ=25°时，弯管部分最大冲蚀率为6.9×10^-5 kg/(m²·s)，相比普通弯管最大冲蚀率减小37.27%；其他角度下的最大冲蚀率变化趋势与四边形肋条结构的弯管相同。由图7c可知，当肋条截面形状为等腰三角形时，肋条部分的最大冲蚀率同样在θ=15°~35°出现高于弯管部分的现象，肋条部分保护壁面免受更大冲蚀，充分发挥了其牺牲作用。肋条位于θ=35°的之前时，均具有一定程度的抗冲蚀效果，特别是θ=20°，弯管的最大冲蚀率为6.63×10^-5 kg/(m²·s)，相比普通弯管最大冲蚀率减小39.72%。之后角度最大冲蚀率变化规律与四边形肋、梯形肋变化相同。总之，3种形状肋条安装在θ=35°前均具有一定的抗冲蚀效果，综合考虑肋条作用和抗冲蚀效果，等腰三角形肋条抗冲蚀作用最佳，且肋条起到牺牲自己保护壁面的作用；四边形肋和等腰梯形肋抗冲蚀作用相近，但等腰梯形肋充分发挥本身作用，四边形肋与弯管变化趋势相同，不便于冲蚀后肋条的替换安装。

图7

新窗口打开| 下载原图ZIP| 生成PPT

图7 不同形状肋条在弯头不同角度下的最大冲蚀率

Fig.7 Maximum erosion rate of rib at different elbow angles with different shapes of quadrilateral rib (a), isosceles trapezoidal rib (b) and isosceles triangular rib (c)

因颗粒轨迹对冲蚀形状的形成有重要影响，现根据普通弯管中颗粒与壁面的碰撞角度，选取θ=0°、20°、50°、70° 4个角度不同形状肋条的冲蚀进行分析。图8为横截面不同形状肋条在不同角度θ下的冲蚀率分布。从横向看，不同横截面形状肋条在相同位置下的冲蚀形状相似。从纵向看，肋条的加入改变了原本的冲蚀分布形状。当θ为0°、20°时，具有一定的抗冲蚀性能，冲蚀率较低。当θ较大时，在肋条处和肋条前冲蚀较为严重，最大冲蚀率大于普通弯管，不具有抗冲蚀性能。随着角度的增加，颗粒逐渐流出弯头区域，肋条发挥作用减弱。

图8

新窗口打开| 下载原图ZIP| 生成PPT

图8 不同形状肋条在不同角度θ下的冲蚀率分布

Fig.8 Erosion rate distribution of rib with different shapes at different angles θ

图9为不同形状肋条在不同角度下的颗粒轨迹图和速度矢量图。结合图7可知肋条改变了颗粒的运动轨迹。当安装在颗粒与壁面第一次碰撞35°之前时，大部分颗粒首先与肋条碰撞后反弹回管道，阻止了大部分颗粒对壁面的的第一次直接碰撞。θ处于35°~60°时，大部分颗粒已与壁面发生碰撞，在肋条处再次与肋条碰撞造成肋条前和肋条处颗粒堆积，造成冲蚀率增大。θ大于60°时，颗粒已完成与壁面第一次碰撞并分成两股相反对称的颗粒流直至与壁面再次碰撞，逐渐流出弯头区域，肋条发挥作用较小，抗冲蚀作用较弱。同时肋条的加入改变了近壁面处的流场分布，均在肋条背部形成一个低速逆流循环区，保护该区域免受直接冲蚀。当安装在20°附近时，不仅抑制了颗粒与壁面的第一次直接碰撞，而且正好保护了冲蚀最严重的区域，所以是最优抗冲蚀位置。从横向看，不同形状下的颗粒轨迹略有不同，可由图10所示颗粒在3种不同横截面形状下的弹性碰撞说明其差异性。当颗粒以相同的角度和3种截面形状肋条碰撞后颗粒方向不同，与四边形肋碰撞后垂直返回计算域，等腰梯形肋的实际碰撞角度为74°，等腰三角形肋为57°。

图9

新窗口打开| 下载原图ZIP| 生成PPT

图9 不同形状肋条在不同角度θ下的颗粒轨迹和速度矢量图

Fig.9 Particle trajectory and velocity vectors of ribs with different shapes at different angles θ

图10

新窗口打开| 下载原图ZIP| 生成PPT

图10 颗粒在不同横截面形状下的碰撞示意图

Fig.10 Collision diagram of particles under different cross section shapes

3.2 不同肋条尺寸对冲蚀的影响

基于参考文献[10]中深度对循环区大小的形成有重要影响，现以肋条宽度L=5 mm为固定尺寸，分析不同肋高条件下的冲蚀变化。因四边形肋跟随弯管最大冲蚀变化，未充分发挥肋条作用，现仅讨论等腰梯形肋和等腰三角形肋。肋条安装在颗粒与壁面第一碰撞之前有更好的抗冲蚀性能，故研究10°、15°、20°、25°、30° 5个典型角度下肋条尺寸变化对冲蚀的影响。

图11为等腰梯形肋安装在10°、15°、20°、25°、30° 5种角度下尺寸变化对冲蚀的影响。随着肋条深度的增加，弯管和肋条部分的最大冲蚀率变化规律不明显，但总体呈现出肋条部分随H的增大波动增大，弯管部分在H<6 mm时，最大冲蚀率随肋条深度H的增加波动减小，H>6 mm时，弯管最大冲蚀率增加。图12为等腰三角形肋安装在10°、15°、20°、25°、30° 5种角度下尺寸变化对冲蚀的影响。肋条的变化始终跟随弯管的变化趋势。当θ为10°、15°、20°时，弯管和肋条的最大冲蚀率均随肋条深度H的增加而减小；当θ为25°、30°时，与等腰梯形肋变化趋势相同，弯管和肋条部分H<6 mm时，最大冲蚀率随肋深的增加减小，H>6 mm时，弯管随之增加。

图11

新窗口打开| 下载原图ZIP| 生成PPT

图11 等腰梯形肋在不同位置最大冲蚀率随深度的变化

Fig.11 Variation of the maximum erosion rate of isosceles trapezoidal ribs with depth at different positions of θ=10° (a), θ=15° (b), θ=20° (c), θ=25° (d) and θ=30° (e)

图12

新窗口打开| 下载原图ZIP| 生成PPT

图12 等腰三角形肋在不同位置最大冲蚀率随深度的变化

Fig.12 Variation of the maximum erosion rate of isosceles triangular ribs with depth at different positions of θ=10° (a), θ=15° (b), θ=20° (c), θ=25° (d) and θ=30° (e)

图13为肋条在不同深度下的速度矢量图。由图可知，随着肋深H的增加，肋条背部的旋涡区逐渐增大，对壁面的保护范围逐渐变大，但由于肋条深度的增加会影响颗粒与肋条碰撞角度，增加颗粒与肋条碰撞的频率，改变颗粒后续运动轨迹，从而从多方面影响冲蚀，便造成最大冲蚀率随肋条深度变化规律不明显的趋势^[10]。总之，综合考虑肋条的抗冲蚀能力与牺牲作用，肋深为6 mm的等腰三角形肋位于θ=25°时抗冲蚀作用最佳，最大冲蚀率为6.2×10^-5 kg/(m²·s)，相比普通弯管抗冲蚀性能提高43.63%。与参考文献[14]相比，该结构抗冲蚀性能更优。

图13

新窗口打开| 下载原图ZIP| 生成PPT

图13 肋条在不同深度下的速度矢量图

Fig.13 Velocity vectors of ribs at different depths

3.3 不同肋条数量对冲蚀的影响

以上一节中最佳抗冲蚀形状尺寸的肋条结构为固定条件，分析弯头部分均匀分布不同数量肋条对弯管冲蚀的影响。当肋条数为2时，肋条分布于弯管入口处和出口处，随着颗粒逐渐流出弯管，出口处肋条作用逐渐削弱，仅入口处肋条改变颗粒轨迹，与单肋条在入口处情况相似，弯管部分的冲蚀大于肋条部分的冲蚀，未充分发挥肋条自身作用，故仅从肋条数为3开始分析。图14为最大冲蚀率随肋条数量的变化。由图可知，肋条部分最大冲蚀率始终大于弯管部分，跟随弯管趋势变化，所考虑数量下随肋条数量的增加呈现“M”形变化趋势，抗冲蚀作用最佳的肋条数为7条，相比普通弯管抗冲蚀性能提高36.45%。图15为不同肋条数量下的速度矢量，肋条数量的增大减小了每个旋涡的大小，但增加了旋涡数量，对壁面起到一定的保护作用，但肋条数量增加影响了颗粒碰撞角度，增加颗粒与壁面碰撞频率，从而导致颗粒轨迹变化，冲蚀率发生变化。

图14

新窗口打开| 下载原图ZIP| 生成PPT

图14 不同肋条数量下的最大冲蚀率

Fig.14 Maximum erosion rate under different number of ribs

图15

新窗口打开| 下载原图ZIP| 生成PPT

图15 不同数量下肋条的速度矢量图

Fig.15 Velocity vector diagram of ribs in different quan-tities

4 结论

(1) 肋条位于35°之前时，改变了颗粒轨迹，阻止了颗粒与壁面的第一次直接碰撞，一定程度上抑制了冲蚀，并且肋条背部形成一个低速逆流循环区，对该区域壁面产生一定保护作用。

(2) 四边形肋最大冲蚀率始终小于弯管部分，跟随弯管趋势变化；等腰梯形肋和等腰三角形肋位于θ=15-35°时具有牺牲自身保护壁面作用。综合考虑肋条的牺牲作用及抗冲蚀性能，等腰三角形肋抗冲蚀作用最佳。

(3) 肋深越大，背部形成的循环区域越大，保护范围越大，但深度变化影响颗粒碰撞角度，增大颗粒与其碰撞频率，所以肋深并非越大越好。本文条件下肋深为6 mm的等腰三角形肋位于θ=25°抗冲蚀性能最佳，相比普通弯管抗冲蚀性能提高43.63%。

(4) 弯头部分均匀分布多个肋条也具有明显的抗冲蚀特性。肋条部分最大冲蚀率始终大于弯管，跟随弯管变化趋势，所考虑数量下随肋条数量的增加呈现“M”形趋势变化，肋条数为7的等腰三角形肋抗冲蚀性能最好。

参考文献

原文顺序

文献年度倒序

文中引用次数倒序

被引期刊影响因子

[1]

Bayer

R G

Fundamentals of wear failures

[A]. BeckerW T, ShipleyR J, eds. Failure Analysis and Prevention [M]. Hitchin: American Technical, 2002: 901

[本文引用: 1]

[2]

Song

X J

, Qi

D Z

, Xu

L J

, et al.

Numerical simulation prediction of erosion characteristics in a double-suction centrifugal pump

[J]. Processes, 2021, 9: 1483

DOI URL [本文引用: 1]

The double-suction centrifugal pumps installed along the Yellow River in China face serious sediment erosion due to the high sediment content which causes the poor operation efficiency of the pump units. The particle motion characteristics and erosion characteristics in a pump under different flow rates and different particle concentrations were numerically simulated based on the particle track model of solid-liquid two-phase flow. The results show that the flow rate has a significant effect on the particle tracks and the erosion caused by the particles in the impeller. The total erosion rate is positively correlated with the flow rate, and increases with the increase in flow rate. The vortex and secondary flow in the impeller have obvious influence on the particle trajectory, which increases the particle concentration at the trailing edge of the pressure surface and intensifies the impact erosion in this area. The particles carried by the vortex intensifies the local erosion. The particle concentration mainly affects the erosion rate, but has little effect on the erosion position. The influence of flow rate on pump erosion is greater than that of the particle properties. These results provide a reference for optimization of the design of anti-erosion blades of double-suction pumps and the regulation and operation of pumping stations.

[3]

Dehdarinejad

, Bayareh

Impact of non-uniform surface roughness on the erosion rate and performance of a cyclone separator

[J]. Chem. Eng. Sci., 2022, 249: 117351

DOI URL

[4]

Zhou

J W

, Liu

S Y

, et al.

Effects of particle shape and swirling intensity on elbow erosion in dilute-phase pneumatic conveying

[J]. Wear, 2017, 380-381: 66

DOI URL

[5]

Alkassar

, Agarwal

V K

, Pandey

R K

, et al.

Influence of particle attrition on erosive wear of bends in dilute phase pneumatic conveying

[J]. Wear, 2021, 476: 203594

DOI URL [本文引用: 1]

[6]

dos Santos

V F

, de Souza

F J

, Duarte

C A R

Reducing bend erosion with a twisted tape insert

[J]. Powder Technol., 2016, 301: 889

DOI URL [本文引用: 1]

[7]

Lin

, Lan

H Q

, Xu

Y G

, et al.

Effect of the gas-solid two-phase flow velocity on elbow erosion

[J]. J. Nat. Gas Sci. Eng., 2015, 26: 581

DOI URL [本文引用: 1]

[8]

Duarte

C A R

, de Souza

F J

Innovative pipe wall design to mitigate elbow erosion: A CFD analysis

[J]. Wear, 2017, 380-381: 176

DOI URL [本文引用: 1]

[9]

San

Y K

, Thien

, Lee Chieng Chen

Numerical study on erosion of a pipe bend with a vortex chamber

[J]. Part. Sci. Technol., 2019, 37: 200

DOI URL [本文引用: 2]

[10]

, Sun

Z Q

, Li

A J

, et al.

Design optimization of hemispherical protrusion for mitigating elbow erosion via CFD-DPM

[J]. Powder Technol., 2022, 398: 117128

DOI URL [本文引用: 3]

[11]

Fan

J R

, Yao

, Cen

K F

Antierosion in a 90° bend by particle impaction

[J]. AIChE J., 2002, 48: 1401

DOI URL [本文引用: 1]

[12]

Yao

, Zhang

B Z

, Fan

J R

An experimental investigation of a new method for protecting bends from erosion in gas-particle flows

[J]. Wear, 2000, 240: 215

DOI URL [本文引用: 1]

[13]

Fan

J R

, Luo

, Zhang

X Y

, et al.

Large eddy simulation of the anti-erosion characteristics of the ribbed-bend in gas-solid flows

[J]. J. Eng. Gas Turbines Power., 2004, 126: 672

DOI URL [本文引用: 1]

In order to find out the ribs with the best anti-erosion efficiency, the erosion effects in three 90° duct bend gas-solid flows with different ribs are investigated and compared with that in the bare bend. Three different kinds of ribs are studied, which have square cross section, isosceles right-angled triangle shape cross section and isosceles triangle shape cross section, respectively. The arrangement and the geometry dimensions of the ribs are the same. The gas phase flow field is solved by the large eddy simulation technique considering the interaction between the gas and the particles, whereas the particles are tracked by Lagrangian method. The results exhibit that the erosion rates of all the ribbed bends are greatly reduced due to the adding of the ribs compared with that of the bare bend. Among the three different kinds of ribs, the isosceles right-angled triangle ribs have the highest anti-erosion ability. The reason is that the particle impact velocity on the walls in the isosceles right-angled triangle ribbed bends is decreased to the utmost and the corresponding particle impact angle is altered adequately due to the special shape of the ribs. In addition, the rib erosion rate, the sidewall erosion rate and the bend wall erosion rate are calculated separately. It is found that the rib erosion rate is half of the wall erosion rate, and the sidewall erosion rate is so low that it can be omitted compared with the bend wall erosion rate.

[14]

Zhu

H J

, Li

Numerical analysis of mitigating elbow erosion with a rib

[J]. Powder Technol., 2018, 330: 445

DOI URL [本文引用: 5]

[15]

Zamani

, Seddighi

, Nazif

H R

Erosion of natural gas elbows due to rotating particles in turbulent gas-solid flow

[J]. J. Nat. Gas Sci. Eng., 2017, 40: 91

DOI URL [本文引用: 1]

[16]

L Y

, Wu

, Yan

, et al.

Numerical simulation of air-solid erosion in elbow with novel arc-shaped diversion erosion-inhibiting plate structure

[J]. Powder Technol., 2021, 393: 670

DOI URL [本文引用: 2]

[17]

Yang

D L

, Xing

B S

, Li

J P

, et al.

Experiment and simulation analysis of the suspension behavior of large (5-30 mm) nonspherical particles in vertical pneumatic conveying

[J]. Powder Technol., 2019, 354: 442

DOI URL [本文引用: 1]

[18]

Peng

W S

, Cao

X W

Numerical prediction of erosion distributions and solid particle trajectories in elbows for gas-solid flow

[J]. J. Nat. Gas Sci. Eng., 2016, 30: 455

DOI URL [本文引用: 2]

[19]

Zheng

, Liu

Y H

, Wang

H X

, et al.

Numerical study on improving the erosion life of ball seat for oil and gas reservoir fracturing

[J]. Eng. Fail. Anal., 2016, 60: 188

DOI URL [本文引用: 1]

[20]

Zhang

, Reuterfors

E P

, McLaury

B S

, et al.

Comparison of computed and measured particle velocities and erosion in water and air flows

[J]. Wear, 2007, 263: 330

DOI URL [本文引用: 1]

[21]

Parsi

, Najmi

, Najafifard

, et al.

A comprehensive review of solid particle erosion modeling for oil and gas wells and pipelines applications

[J]. J. Nat. Gas Sci. Eng., 2014, 21: 850

DOI URL [本文引用: 1]

[22]

Grant

, Tabakoff

Erosion prediction in turbomachinery resulting from environmental solid particles

[J]. J. Aircr., 1975, 12: 471

DOI URL [本文引用: 1]

[23]

Launder

B E

, Spalding

D B

The numerical computation of turbulent flows

[J]. Comput. Methods Appl. Mech. Eng., 1974, 3: 269

DOI URL [本文引用: 1]

[24]

Pereira

G C

, de Souza

F J

, de Moro Martins

D A

Numerical prediction of the erosion due to particles in elbows

[J]. Powder Technol., 2014, 261: 105

DOI URL [本文引用: 1]

Fundamentals of wear failures

2002

... ASTM国际G40-13标准认为冲蚀是由于固体表面与流体、多组分流体、撞击液体或撞击固体颗粒之间的机械相互作用而导致原始材料从固体表面逐渐流失的现象.研究表明，一个工业化国家的冲蚀磨损成本约占其国民生产总值的1%至4%^[1].冲蚀经常发生在采矿、冶金、石油、天然气等工业生产中，特别是泵、旋风分离器、燃气轮机叶片、气力输送管道等部件中^[2~5]，其中气力输送管道中的冲蚀破坏更为严重^[6].气力输送管道主要由水平管、垂直管和弯管组成.弯管作为改变流动方向的部件，起着连接水平管和垂直管的作用，更容易发生冲蚀.研究表明，弯管冲蚀破坏的程度约为直管段的50倍^[7].这种冲蚀破坏不仅浪费材料，消耗能源，降低设备效率，而且加速设备故障，降低弯管的使用寿命，对管路系统的安全运行造成威胁，从而酿成更大的经济损失.目前常用计算流体力学-离散相方法 (CFD-DPM) 和计算流体力学-离散元方法 (CFD-DEM) 方法研究冲蚀.CFD-DPM仅适用于稀相流，而CFD-DEM方法在计算密相流方向有更大的优势.合理的弯头结构以提高弯管的抗冲蚀性能在工业生产中具有重要意义，研究者对弯管结构进行了大量改进以提高弯管的抗冲蚀性能.Duarte等^[8]提出一种沿水流方向扭转管壁的新颖管壁设计，这种结构会在弯头上游产生一个旋流，最大冲蚀率减小33%.San等^[9]提出一种涡流室设计，并分析流动Reynold数、颗粒Stokes数和涡流室大小对冲蚀的影响，结果表明：与普通半径的弯管相比，涡流室总能减少侵蚀.以上改变弯管结构的方法虽能一定程度上减少冲蚀，但改变了弯管原本结构，影响了管道整体布局，且冲蚀仍发生在弯管壁面上，发生冲蚀破坏时候仍需替换弯管^[10].相比以上直接改变弯管结构的方法，内壁面增加肋条是更直接有效的方法，冲蚀破坏中只有肋条需被更换^[11].Yao等^[12]通过实验方法分析弯管加横截面为四边形肋条的抗冲蚀性能，结果表明：肋条可提高弯管的抗冲蚀性能，矩形肋条相比方形肋条抗冲蚀性能更佳.Fan等^[13]研究了横截面为矩形、等腰三角形、等腰直角三角形3种不同截面形状肋条的抗冲蚀性能，结果表明：3种不同形状肋均提高弯管的抗冲蚀性能，等腰直角三角形肋条的抗冲蚀能力最强.以上均在横截面为正方形弯管中均匀分布多个肋条，未见圆形弯管且未考虑单个肋在弯管不同角度的抗冲蚀分布.再者，Zhu等^[14]虽将单个轴向角为5°，环形角为180°的肋条安装在90°圆形弯头外径方向不同位置上，但未考虑肋条横截面不同形状截面参数及多个肋条下的抗冲蚀特性. ...

Numerical simulation prediction of erosion characteristics in a double-suction centrifugal pump

2021

Impact of non-uniform surface roughness on the erosion rate and performance of a cyclone separator

2022

Effects of particle shape and swirling intensity on elbow erosion in dilute-phase pneumatic conveying

2017

Influence of particle attrition on erosive wear of bends in dilute phase pneumatic conveying

2021

Reducing bend erosion with a twisted tape insert

2016

Effect of the gas-solid two-phase flow velocity on elbow erosion

2015

Innovative pipe wall design to mitigate elbow erosion: A CFD analysis

2017

Numerical study on erosion of a pipe bend with a vortex chamber

2019

... 连续相采用Euler形式的非定常Reynolds平均Navier-Stokes方程(URANS)求解，其质量守恒方程和动量方程如下^{[9, 15,16]}： ...

Design optimization of hemispherical protrusion for mitigating elbow erosion via CFD-DPM

2022

... 基于参考文献[10]中深度对循环区大小的形成有重要影响，现以肋条宽度L=5 mm为固定尺寸，分析不同肋高条件下的冲蚀变化.因四边形肋跟随弯管最大冲蚀变化，未充分发挥肋条作用，现仅讨论等腰梯形肋和等腰三角形肋.肋条安装在颗粒与壁面第一碰撞之前有更好的抗冲蚀性能，故研究10°、15°、20°、25°、30° 5个典型角度下肋条尺寸变化对冲蚀的影响. ...

... 图13为肋条在不同深度下的速度矢量图.由图可知，随着肋深H的增加，肋条背部的旋涡区逐渐增大，对壁面的保护范围逐渐变大，但由于肋条深度的增加会影响颗粒与肋条碰撞角度，增加颗粒与肋条碰撞的频率，改变颗粒后续运动轨迹，从而从多方面影响冲蚀，便造成最大冲蚀率随肋条深度变化规律不明显的趋势^[10].总之，综合考虑肋条的抗冲蚀能力与牺牲作用，肋深为6 mm的等腰三角形肋位于θ=25°时抗冲蚀作用最佳，最大冲蚀率为6.2×10^-5 kg/(m²·s)，相比普通弯管抗冲蚀性能提高43.63%.与参考文献[14]相比，该结构抗冲蚀性能更优. ...

Antierosion in a 90° bend by particle impaction

2002

An experimental investigation of a new method for protecting bends from erosion in gas-particle flows

2000

Large eddy simulation of the anti-erosion characteristics of the ribbed-bend in gas-solid flows

2004

Numerical analysis of mitigating elbow erosion with a rib

2018

... 因URANS方程中引入了Reynolds应力项，因此需采用湍流模型来封闭控制方程^[14].两方程k-ε模型主要包括标准k-ε模型、RNG k-ε和Realizable k-ε模型.因Realizable k-ε模型带涡旋修正，能提高剪切流和旋流的精度，故本文采用Realizable k-ε模型，其方程为^[17]： ...

... 进行数值模拟之前，验证模拟方法及模型的有效性是必要的.采用Zhu等^[14]实验结果与本文模拟值对比分析.利用本文模拟方法模拟相同实验条件下的冲蚀率，得到图5所示管道外壁面中心线30°实验值与模拟值的对比图.可知，模拟与实验冲蚀率具有较好的一致性，验证了模拟方法及模型的可靠性. ...

... [14]与模拟值对比图Comparison between experimental and simulated values

Fig.5

3 结果与讨论3.1 不同肋条形状对冲蚀的影响

以H=L=5 mm的3种截面肋条结构为固定条件，模拟分析3种不同横截面形状肋条安装在弯头不同角度下的冲蚀规律.图6为普通弯管在该条件下的冲蚀分布及颗粒轨迹.可知普通弯管的冲蚀形状呈“V”形分布，冲蚀最严重部位出现在沿弯头曲率方向35~60°附近，最大冲蚀率为1.1×10^-4 kg/(m²·s).由于颗粒轨迹对冲蚀形状的形成有重要影响，可通过颗粒轨迹分析冲蚀形状.由对应的颗粒轨迹图可知颗粒在35°左右第一次与管壁碰撞形成椭圆形冲蚀严重区，之后分成两股相反对称的颗粒流，直至与壁面发生第二次碰撞，便形成图示的“V”形分布，这与Peng等^[18]和Pereira等^[24]得到的结果一致. ...

Erosion of natural gas elbows due to rotating particles in turbulent gas-solid flow

2017

... 连续相采用Euler形式的非定常Reynolds平均Navier-Stokes方程(URANS)求解，其质量守恒方程和动量方程如下^{[9, 15,16]}： ...

Numerical simulation of air-solid erosion in elbow with novel arc-shaped diversion erosion-inhibiting plate structure

2021

... 连续相采用Euler形式的非定常Reynolds平均Navier-Stokes方程(URANS)求解，其质量守恒方程和动量方程如下^{[9, 15,16]}： ...

... 利用基于牛顿第二定律的拉格朗日方法追踪离散颗粒，颗粒在流场中受力示意图如图1所示，其控制方程如下^[16,18,19]： ...

Experiment and simulation analysis of the suspension behavior of large (5-30?mm) nonspherical particles in vertical pneumatic conveying

2019

Numerical prediction of erosion distributions and solid particle trajectories in elbows for gas-solid flow

2016

... 利用基于牛顿第二定律的拉格朗日方法追踪离散颗粒，颗粒在流场中受力示意图如图1所示，其控制方程如下^[16,18,19]： ...

... 以H=L=5 mm的3种截面肋条结构为固定条件，模拟分析3种不同横截面形状肋条安装在弯头不同角度下的冲蚀规律.图6为普通弯管在该条件下的冲蚀分布及颗粒轨迹.可知普通弯管的冲蚀形状呈“V”形分布，冲蚀最严重部位出现在沿弯头曲率方向35~60°附近，最大冲蚀率为1.1×10^-4 kg/(m²·s).由于颗粒轨迹对冲蚀形状的形成有重要影响，可通过颗粒轨迹分析冲蚀形状.由对应的颗粒轨迹图可知颗粒在35°左右第一次与管壁碰撞形成椭圆形冲蚀严重区，之后分成两股相反对称的颗粒流，直至与壁面发生第二次碰撞，便形成图示的“V”形分布，这与Peng等^[18]和Pereira等^[24]得到的结果一致. ...

Numerical study on improving the erosion life of ball seat for oil and gas reservoir fracturing

2016

... 利用基于牛顿第二定律的拉格朗日方法追踪离散颗粒，颗粒在流场中受力示意图如图1所示，其控制方程如下^[16,18,19]： ...

Comparison of computed and measured particle velocities and erosion in water and air flows

2007

... 由于冲蚀过程复杂，受很多因素的共同影响，E/CRC冲蚀模型^[20,21]考虑了硬度和颗粒形状因素，故本文采用E/CRC冲蚀模型，其冲蚀模型如下： ...

A comprehensive review of solid particle erosion modeling for oil and gas wells and pipelines applications

2014

... 由于冲蚀过程复杂，受很多因素的共同影响，E/CRC冲蚀模型^[20,21]考虑了硬度和颗粒形状因素，故本文采用E/CRC冲蚀模型，其冲蚀模型如下： ...

Erosion prediction in turbomachinery resulting from environmental solid particles

1975

... 粒子与壁面碰撞会发生能量传递与损失，这里采用法向恢复系数e_n和切向恢复系数e_t表示^[22]： ...

The numerical computation of turbulent flows

1974

... 气力输送中气固两相采用基于Euler-Lagrange方法的DPM模型.连续相为空气，密度1.18 kg/m³,粘度1.8×10^-5 Pa·s；离散相为沙粒，密度1650 kg/m³，粒径为单一粒径100 μm，颗粒质量流量为0.02 kg/s.考虑连续相和颗粒相的双向耦合作用.边界条件设置中，连续相入口为速度入口，流速30 m/s，出口为大气压压力出口，壁面采用无滑移边界条件，近壁面处采用Launder和Spalding^[23]提出的标准壁面函数；离散相为惰性颗粒垂直入口表面射入，入口出口边界条件为escape逃逸类型，壁面为reflect反射类型，假设入口处颗粒速率与连续相入口速率相同. ...

Numerical prediction of the erosion due to particles in elbows

2014

〈

〉