干扰装置对热电偶套管绕流及减阻性能

摘要:针对石油裂解气管路中热电偶套管容易因旋涡脱落及受绕流阻力的影响而发生疲劳损伤的问题,使用LES对热电偶套管绕流进行数值模拟,研究套管不同插入深度对尾流场的影响以及套管安装角和干扰装置对涡脱及套管减阻性能的影响。结果表明:随着套管插入深度的增加,套管受流体冲击逐渐增强且尾流中逐渐有旋涡交替脱落;随着安装角的减小,套管阻力系数和升力系数逐渐降低,但套管所受流体冲击面积逐渐增大,综合考虑选择60°安装角,减阻效果达11.4%;增加干扰装置能大幅降低阻力系数,较基础单管所受阻力系数降低35.1%,旋涡脱落频率显著降低。
　　石油裂解气管路中的热电偶套管作为热电偶传感器的保护装置,经常发生疲劳破坏或断裂,影响石油裂解气生产的连续性。在亚音速横向流中，任何非流线型的物体尾部若有足够的拖迹边缘都会产生旋涡脱落。对悬臂式安装在管路中的热电偶套管,当旋涡在套管两侧周期性脱落时,涡脱产生的升力诱发套管振动,导致热电偶套管出现断裂。由于流体具有黏性,在套管表面会产生很大的速度梯度,当发生流动分离现象时会在套管上生成绕流阻力,损坏套管。所以,对热电偶套管绕流减阻抑振控制的研究十分重要。
　　减阻技术已受到众多学者重视。Gao等[功通过粒子图像测速技术对上游添加分流板的三维柱体绕流进行风洞实验,证明分流板的存在能够有效降低阻力系数。Islam等[°]通过布置分流板,研究了单、双分流板对方柱绕流减阻和旋涡脱落的影响,发现下分流板的存在使得圆柱绕流阻力系数和升力系数均方根臧少,旋涡脱落延时;上下分流板同时存在时边界层提前分离,抑制了旋涡脱落,改变了涡脱频率。对于流场中串列柱体,Grioni等采用尺度自适应仿真湍流模型对雷诺数Re=2X105的串列双圆柱进行研究,发现柱体中心距对柱体阻力有影响,双柱绕流之间的扰动作用会使后柱阻力减小。以热电偶套管为对象,研究不同插人深度和安装角对单热电偶套管减阻性能的影响,通过数值仿真对增加千扰装置后热电偶套管的减阻抑振性能进行计算分析。
1物理模型与边界条件
1.1几何模型
　　单套管与加干扰装置套管的几何模型如图1所示。为方便计算,对模型进行简化,将热电偶套管顶部以下区域作为计算域,套管前后管路取260mm,千扰装置与套管间距为3mm,套管插人深度为管路内壁至套管底端的距离,用h表示。
 
1.2湍流模型与网格划分
　　干扰装置为不规则体,计算所用网格为非结构网格,选用大涡模拟(largeeddysimulation,LES)方法且选择亚网格尺度应力模型(smargorinskymodel,SM),系数Cs=0.2。由于需要对套管表面气体流动进行精确计算,所以对套管壁面网格进行加密,使无量纲数y+<2,利用ICEM-CFD软件画出网格如图2所示。

1.3边界条件与重要参数确定
　　计算边界条件为:进口采用速度进口,出口为压力出口,计算域内其余壁面设为固壁无滑移边界条件,具.体参数如表1所示。通过计算可知,数值仿真计算时间步长△t=0.0001s时,既能满足计算准确性又能减少计算时长。
 
2单热电偶套管绕流减阻研究
2.1套管插入深度对流场的影响
　　通过计算得到90°安装角下套管插人深度h分别为1.5D.1D和0.5D(D为套管直径)时流场的速度云图如图3所示。当插人深度为1.5D时,由于套管两侧边界层分离,套管后会形成旋涡并在尾部两侧周期性脱落,这将使套管振动并加剧套管的损坏;当插入深度小于1.5D时，气体在管路中央汇集,无明显旋涡形成。

　　图4为套管不同插人深度对应流场的平均压力云图,可见当插人深度为1.5D时,气体高压分布区域最广、压力最大,说明此时套管承受的气体冲击力较高,并且套管前后存在的较大压力差会产生很高的绕流阻力,加剧套管损坏。当插人深度小于1.5D时,高压区域显著减少,套管前后压力差降低,阻力减小。
2.2安装角度对套管减阻性能影响
　　热电偶套管绕流及减阻性能分析采用的主要无量纲参数为阻力系数Cd、升力系数C1、斯特劳哈尔数St和压力系数Cp.阻力系数Cd是由套管在平行于流体流动方向所受的阻力无量纲化后得到的,通常以阻力系数平均值`Cd作为分析对象,具体定义为:
Cd=2Fd/ρAU<sup>2</sup>(1)
　　式中:Fd为套管平行于流体流动方向所受的阻力,N;ρ为管道内流体密度,kg/m³;A为柱体迎流面投影面积,m<sup>2</sup>;U为管道内来流速度,m/s。
　　升力系数C1是由套管在垂直于流体流动方向所受的升力无量纲化后得到的,通常以升力系数均方根值C1rms作为分析对象,具体定义为:
C1=2F1/ρAU<sup>2</sup>，(2)
　　式中F1为套管垂直于流体流动方向所受的升力,N。
　　斯特劳哈尔数St是用于描述套管尾流中旋涡脱落现象的无量纲数,具体定义为:
St=ƒsD/U,(3)
　　式中:ƒs为套管尾流涡脱频率,Hz;D为套管直径,m。.
　　压力系数Cp描述的是流场中的相对压力,具体定义为:
Cp=2(P-P0)/ρU<sup>2</sup>。(4)
　　式中:P为所求压力系数点处的静压,Pa;P0为流场内无穷远处压力,Pa。
　　单热电偶套管的安装角度(以下简称为安装角)用θ表示,先后取90°、60°和45°来研究套管安装角对减阻性能的影响。升力和阻力系数监测位置均选取套管插入深度中间高度即h/2处,不同安装角计算结果如表2所示。
 
2.2.1安装角度对阻力系数的影响
　　图5为各安装角下单热电偶套管阻力系数随计算时间的变化曲线。热电偶套管绕流可视为有限长圆柱绕流,其流场受套管自由端面下洗作用11]影响,有限长圆柱的阻力系数与升力系数均表现出无规律性且升力系数曲线的峰值较小。选取图中0.6到1s之间的阻力系数值并求其平均值`Cd,得到各安装角下单热电偶套管的`Cd分别为1.012、0.897和0.836。由此可知,降低安装角后热电偶套管受到的绕流阻力明显下降;安装角为60°时，阻力降低11.4%;安装角为45°时,阻力降低17.4%。
 
　　各安装角下单热电偶套管绕流升力系数随时间的变化曲线如图6所示。由于升力系数曲线呈正负波动，因此通常以升力系数均方根值作为分析对象。选取0.6s到1s之间的升力系数值并求其均方根值C1RMS,得到在各安装角下单热电偶套管的C1RMS分别为0.102、0.079和0.072。由此可知，降低安装角后单热电偶套管绕流升力系数均方根值有明显下降,降低程度要比阻力系数大;当安装角为60°时,降低22.5%;安装角为45°时，降低29.4%.
2.2.3安装角度对斯特劳哈尔数的影响
　　图7是升力系数功率谱密度图,表示信号功率在频域的分布情况,是对各安装角下单热电偶套管绕流升力系数进行快速傅里叶变换(fastFouriertransform,FFT)后得到的。由图7可知,当安装角分别为90°、60°和45°时,单热电偶套管绕流旋涡发放频率f。分别为200.5、215.3和220.0Hz,S,相差不大，说明改变安装角对单热电偶套管绕流旋涡发放频率的影响很小。
 
2.2.4安装角度对压力系数的影响
　　由于当前流动的雷诺数较高,摩擦阻力可忽略不计,故套管的压差阻力占绕流阻力的主要部分。不同安装角下套管的周向平均压力系数随圆周角度变化如图8所示。由图8可见,压力系数最大点均位于套管正对来流位置且最大值为1.0左右。随着安装角的减小,套管绕流剪切层分离点与套管前部压力改变不明显，但背压呈逐渐增大趋势。套管前后压力系数差在逐渐减小,说明套管前部与背部压差降低,这是热电偶套管阻力降低的一个原因。
2.2.5安装角度对压力分析的影响
　　图9为不同安装角下流场压力分布云图,不同安装角下流场压力均呈对称分布,且45°安装角下壁面高压分布区域最大,90°安装角下壁面高压分布区域最小,套管背压逐渐增加,这与上文压力系数表征一致。说明安装角的降低会增大套管的迎流面和壁面高压分布区域,降低套管的安全性。综上,通过降低套管安装角,可以降低套管绕流阻力和升力。当安装角分别为90°、60°和45°时,热电偶套管绕流阻力系数分别为1.012,0.897和0.836,升力系数均方根值分别为0.102.0.079和0.072。安装角分别为60°和45°时,虽然阻力系数和升力系数均方根值较90°安装角时均有所降低,但两种安装角下得到的参数数值差异不大,并且两种情况下斯特劳哈尔数相差不大,涡脱频率几乎没有改变,因此认为60°安装角足以达到减阻抑振的目的。同时随着安装角的降低,套管壁面局部高压区的面积逐渐增大,对套管造成的冲击也就越严重。综合考虑升力、阻力和压力,套管的最佳安装角为60°。
3干扰装置对套管减阻性能影响
　　为解决气流长时间冲击与祸致振动导致的套管破坏,为热电偶套管增加干扰装置,避免套管直接暴露在气流的冲击中。采用60°安装角对套管及千扰装置进行安装,计算模型如图1所示。除物理模型不同外,计算条件均与单热电偶套管绕流计算相同。通过与相同安装角时单热电偶套管绕流计算结果对比,分析干扰装置对套管及流场的影响。
 
3.1阻力系数
　　图10为加千扰装置前后热电偶套管绕流阻力系数随时间的变化曲线。选取0.6s到1s之间的阻力系数值并求平均值`Cd,得到60°安装角时加干扰装置前后套管表面Ca分别为0.897和0.657,即加干扰装置后的阻力系数值低于单热电偶套管。表3为各工况下套管绕流平均阻力系数值,60°安装角时加于扰装置较单套管能显著减少套管受到的绕流阻力,减阻可达26.8%.
 
 
　　表4为各安装角度下套管绕流C1RMS。由表4可知,各安装角对套管绕流C1RMS的影响与其对`Cd的影响相似,均呈逐渐减小趋势,说明60°安装角时加干扰装置能有效降低套管绕流升力,较单套管降低幅度可达26.6%。
3.3斯特劳哈尔数
　　对图11中的套管绕流升力系数历时曲线分别进行FFT,得到对应的套管升力系数功率谱密度图,如图12所示。由图12可知,60°安装角下加干扰装置热电偶前后套管绕流旋涡脱落频率ƒs分别为215.3和125.0Hz,相应的St分别为0.047和0.027,说明加干扰装置不仅能够减小升力系数和阻力系数,还能够改变套管绕流旋涡脱落频率。在主频附近有其他峰值出现,说明加干扰装置后破坏了套管尾流旋涡单一的发放频率,降低套管因共振而发生损坏的风险。
 
3.4压力系数
　　加干扰装置前后套管周向压力系数随圆周角度变化分布曲线如图13所示。由图13可见,单套管时压力系数最高点在1.0左右,加干扰装置后压力系数最高点数值远低于单套管时的数值。加干扰装置后,流体的冲击主要作用在千扰装置迎流面上,套管前缘处于低压区,套.管表面周向压力系数最高点为负值。加干扰装置前后套管表面压力系数相差0.34左右,远低于单套管计算结果,导致套管所受阻力减小。套管绕流剪切层分离点向后移动从而使流动分离延迟且压力系数整体低于单套管。
 
3.5压力分析
　　加干扰装置后的流场压力云图如图14所示。由图14可见,管路内气体的压力分布相对均匀,干扰装置.的迎流面出现了局部高压区,结构两侧出现了小范围的局部低压区,产生的局部压力梯度较小,对干扰装置产生的冲击作用有限,该冲击作用产生的压力在套管壁面上均匀分布.
 
　　表5为改变安装角度与增加干扰装置后整体计算结果。对比发现，当热电偶套管安装角分别为60°(单管).45°和60°(千扰装置)时,套管绕流`Cd、C1RMS和St较90°安装角时均有所降低,其中降低最大的是60°(干扰装置)且此时`Cd、C1RMS和St降低幅度分别达35.1%,43.1%和38.6%;套管在60°(干扰装置)安装角时要比90°、60°(单管)和45°安装角时达到的减阻抑制效果好。综上,优化安装角和引人干扰装置改变了流场的分布规律,不仅减少套管升力系数和阻力系数,还降低旋涡脱落频率,.防止发生套管结构共振现象,进一步验证了的有效性。
4结论
　　针对石油裂解气管道中热电偶套管在应用中出现的断裂损坏现象,优化了套管安装角和插人深度,并通过加装干扰装置来提高套管结构的安全性。结论如下:
1)随着套管插人深度的增加,套管受流体冲击力逐渐增强且尾流中逐渐有旋涡交替脱落。降低套管的安装角,可以减小套管绕流升力系数和阻力系数,降低套管结构的破坏程度。
2)降低套管的安装角后,会增大局部高压区的面积,引发对套管更严重的冲击。综合考虑,最佳安装角度为60°。
3)加装千扰装置后,热电偶套管表面阻力系数大幅度降低,旋涡脱落频率降低。