T型管道的冲蚀磨损数值模拟分析作者：彭方现闫宏伟李亚杰袁飞来源：《当代化工》2020年第03期Numerical Simulation Analysis of Erosion Wear of T-junctionPENG Fang-xian，;YAN Hong-wei，;LI Ya-jie，YUAN Fei（School of Mechanical Engineering， North University of China， Shanxi Taiyuan 030051，China）因為管道运输具有结构简单、运输量大、高效节能等一系列的优点，所以管道运输作为目前最主要的流体介质运输方式被广泛使用，成为继公路、铁路、航空、水路等运输方式后的第五大运输方式[1-4]。
但是在管道使用过程中由于流体介质中不可避免地携带的一些微小颗粒物体，会随着流体介质的流动不断地冲刷管道内壁面，对管道内壁面造成冲蚀磨损的现象，造成管道壁面破裂发生泄漏事故，产生巨大的环境污染及经济损失[5-7]。
所谓冲蚀磨损，就是指小颗粒以一定的速度和角度对靶材表面进行撞击，使靶材的表面产生损耗。
石油化工、航空航天、能源机械等领域广泛地存在冲蚀磨损现象，这也是设备失效的主要原因之一，对企业的安全生产造成重大影响，越来越引起国内外学者的普遍关注。
近年来，国内外一些有关专家学者采用科学实验或数值模拟分析的研究方法对不同类型下的油气管道的冲蚀问题进行了大量的深入研究。
Vigolo， D[8]等将实验结果与理论粒子追踪模型预测的轨迹进行比较，确定了黏性边界层在延迟粒子和降低与基板的碰撞速率方面的作用。
随着计算机的发展，计算流体动力学软件也逐步应用于流体仿真计算中。
Thiana A等[9，10]在弯头外半径和几种不同的流动条件模拟了最大侵蚀速率、侵蚀位置。
闫宏伟等[11]研究了颗粒直径、入口速度及质量流量分别对弯头冲蚀速率的影响。
许留云等[12]得出了弯管的弯曲外侧70°～90°之间的位置为冲蚀磨损最严重的部位。
三通管件在化工行业中的应用也比较广泛，三通的主要功能是为了引导管内流体介质分流或者合流，改变了其内部流体介质的流动方向，三通部件通过使管道中的流体介质在分流或者合流过程中产生了较大扰动，进一步促使管内流动介质发生摩擦、碰撞，并加剧冲蚀现象的形成，因此受到冲蚀磨损的程度比直管严重50倍。
其更容易发生泄漏危害。
例如，我国迪那气田曾发生泄漏事故的原因，即为，三通发生穿孔而导致泄漏此外。
陈宇等[13，14]研究了异面三通管在不同工况下的冲蚀磨损规律。
梁金川等[15]分析了T型管道内部的流体介质的速度场和压力场。
何兴建等[16]分析了T型弯头的冲蚀情况，得出了在T型弯头盲管区域产生了漩涡，充当保护垫作用，能减缓壁面冲蚀。
;除了三通管道本身结构、环境、材料因素的影响外，三通管道也会受到流体介质的速度、颗粒质量流量及颗粒形状等一系列的流体自身特征的影响。
基于此，本文分析了流体入口速度、颗粒质量流量及颗粒的形状三种因素对管道冲蚀磨损的影响作用。
根据流体动力学理论，流动状态划分的主要依据是雷诺数，当雷诺数超过2;300时，流动状态为湍流。
其公式为：经过计算管道内的雷诺数R e>6 000，因此管内的流体运动状态为湍流。
在流动状态中湍流是一种特别复杂的状态，无法全面准确地描述湍流状态。
Fluent仿真模拟软件中有Laminar、k-epsilon和k-omega等模型。
为了在保证计算精度的前提下，本文主要选用标准k-epsilon模型作为湍流的物理模型，因为Laminar模型是单方程模型，虽然其计算速度快，但不能保证其计算精度;k-omega模型的计算精度虽然和k-epsilon模型一样，但不易得到它的有关参数。
k-epsilon湍流物理公式[17，18]见式（2）、（3）。
1.2 ;冲蚀磨损模型因为流体中固体小颗粒所占体积分数很小，所以选择DPM（Discrete Phase Model）模型来计算管道壁面的冲蚀磨损速率。
在DPM模型中不考虑砂粒之间的相互影响，忽略砂粒本身的变化的影响[19，20]，在该工程模型中，冲蚀磨损速率可用公式定义如下：式中：R erosion;—磨损碰撞速率，kg/（m2·s）;2 ;管道模型2.1 ;管道模型及物性参数T型三通管的几何模型如图1所示，其中D=50 mm;L1=1 000 mm;L2=500 mm。
模型中连续相设置为天然气;离散相设置为砂粒，其密度为2 650 kg/m3。
端口1、2为入口，两处天然气初始速度相同且在两管相贯处充分混合，经出口端流出。
2.2 ;网格划分和边界条件设置利用Ansys meshing对T型管进行四面体网格划分，并在主管和支管相连接处进行局部网格加密，由于在管道近壁面处需要考虑流体层流影响，同时在管道壁面设置膨胀层网格。
如图2所示。
经过网格无关性检验，当网格数量数为393 801时最大冲蚀速率不再随着网格数量的增加而发生较大的变化。
在入口边界条件设置中选择速度入口边界条件，操作压力选择默认值，连续相出口边界条件设置为选择outflow;离散相壁面边界条件设置为“reflect”，出口处设置为“escape”;选取SIMPLE算法，采用二阶迎风格式。
3 ;影响因素分析在实际工程应用中，最大冲蚀率是考核管道安全的一个重要参数，T型三通管的冲蚀速率云图如图3所示，可知，在此流体流向的工况下T型管交汇处的迎流壁面处为T型管最容易发生冲蚀的部位，在迎流面偏向出口处区域为冲蚀速率最大的位置，而非入口1正对着的区域，因为由于入口2与入口1的流体在交汇处复合而使流体速度方向向出口方向处发生了偏移，故而使得冲蚀区域偏向出口区域，冲蚀速率也随着入口速度的增大而逐渐增大。
因为入口1的速度和入口2的速度大小一致，所以模拟得到的冲蚀云图基本不变，即冲蚀易发生的具体位置保持不变。
3.1 ;流体速度大小与冲蚀速率的关系管道内的流体的速度是影响冲蚀速率的重要因素之一。
颗粒的质量流量设置为0.005kg/s，颗粒的圆球度设置为1，其他参数不变，只改变天然气的速率大小进行数值模拟仿真，表1为不同的天然气入口速度下对应的最大冲蚀速率数值。
颗粒的质量流量设置为0.005 kg/s，颗粒的圆球度设置为1，其他参数不变，只改变天然气的速率大小进行数值模拟仿真，表1为不同的天然气入口速度下对应的最大冲蚀速率数值。
图4为最大冲蚀速率和天然气速度的关系曲线图，由其可知，最大冲蚀速率随着天然气入口速度的持续增大而不断增大，造成这样的主要原因之一是当管道内的流体速度较低时，由于固体小颗粒自身的速度较小，所对管道内壁的撞击冲量较小，故而冲蚀结果不明显;随着管道流体流动速度的不断增大，颗粒的动能由于气体的作用而增大，从而導致砂粒对管道管壁撞击时的冲量增大，对的管道冲蚀作用也就更加严重。
3.2 ;颗粒的质量流量与冲蚀速率的关系颗粒质量流量体现了单位时间内撞击管道壁面的砂粒数量的多少。
流体初始速度设置为2 m/s，圆球度设置为1，只改变颗粒的质量流量的大小进行数值模拟仿真。
表2为不同的质量流量对应的最大冲蚀速率值。
由图5可知，冲蚀速率和颗粒的质量流量呈正比例关系，即当颗粒的质量流量从0.001 kg/s增大到0.007 kg/s时，冲蚀速率呈线性增大;当质量流量从0.007 kg/s增大到0.009 kg/s时，最大冲蚀速率几乎不变，但当颗粒的质量流量大于0.009 kg/s时，最大冲蚀速率快速增大，呈指数式增长。
造成该结果的原因是，随着质量流量的增大，单位时间内对壁面冲击的固体小颗粒的数量增多，对壁面的冲击力度增大，所以冲蚀速率增大。
颗粒的形状也是影响管道冲蚀速率的一个重要因素，当固体小颗粒的形状越不规则时，造成的冲蚀现象越严重。
体积相同的圆球的外表面积与颗粒的外表面积之比，称为该颗粒的圆球度。
圆球度为无因次量，由于体积相同而形状不同的物体中，圆球的外表面积为最小，所以颗粒的圆球度恒小于1。
圆球度的不同，说明颗粒形状与圆球的差异程度。
当圆球度越接近0时，其形状越不规则，若其值越接近于1，则其形状也就越接近于圆球，当颗粒的形状系数为1时，其形状为标准的球形。
流体初始速度设置为2 m/s，颗粒的质量流量为0.005 kg/s，逐渐增大颗粒的圆球度的数值进行数值模拟仿真。
表3为不同的颗粒的圆球度数值对应的冲蚀速率值。
由颗粒圆球度数和冲蚀速率变化关系曲线图6所示，可知，颗粒圆球度在逐渐变大的过程中，管道壁面冲蚀速率却逐渐变小。
当颗粒圆球度由0.1增大到0.4时，冲蚀速率的数值由2.86×10-9/（kg·m2·s-1）减小到9.53×10-10/（kg·m2·s-1），减小幅度较大;当颗粒圆球度在0.4之后，随着颗粒的圆球度数值的增大，冲蚀速率呈缓慢减小的趋势。
因为当砂粒圆球度数值由0到1变化时，颗粒形状逐渐接近于球状，对壁面的冲击的方式发生变化，由切削为主的磨损变为以研磨为主的磨损，造成的冲蚀破坏也相应地减小了。
应用控制变量法针对T型三通管的冲蚀影响因素进行分析，得出了T型管道易磨损破坏的位置，可为石化管道的防护和检测定位泄漏提供合理的建议，通过设定不同的影响参数，得到的冲蚀规律为：（1）在该工况下T型三通管的迎流壁面处最易发生冲蚀磨损，最大冲蚀速率发生在偏向出口的位置处，而不是正对迎流壁面处。
（2）冲蚀速率随流体速率增大而增大;（3）蚀速率和砂粒的质量流量呈正比例关系，即当颗粒的质量流量从0.001 kg/s增大0.007 kg/s时，冲蚀速率呈线性增大;当质量流量从0.007 kg/s增大到0.009 kg/s时，冲蚀速率几乎不变，但当颗粒的质量流量大于0.009 kg/s时，冲蚀速率快速增大，呈指数式增长;（4）当砂粒圆球度系数由0.1增加到0.4时，冲蚀速率的数值减小幅度较大，呈指数式减小;当颗粒的圆球度为0.4之后时随着砂粒的圆球度的增大，冲蚀速率呈缓慢减小的趋势。
参考文献：。