水击泄压阀流场数值模拟研究.docxVIP

水击泄压阀流场数值模拟研究 0 输油管道的水击泄压 近年来，油藏长期存在迅速发展。长输管道一个主要的潜在不安全因素是来自管道内水击的破坏。当管道内油品流速发生变化时,就会在管道内产生水击压力脉冲,流速变化越快,脉冲压力就会越高,对输油管道造成的危害越大,当产生的压力冲击大于管道及管道设备的承受压力时就会造成管道及设备的损坏,油从损坏处喷射出来,造成油品泄漏,引发环境污染甚至发生火灾爆炸事故。及时消除异常压力波,避免水击破坏最主要的安全措施是在管道上安装水击泄压阀。 由于泄压过程的复杂性和多样性,难以通过实验分析给出较为精确的结论,且难以获得内流的细节特征。本文利用Fluent软件对不同管道压力下泄压阀开启过程进行数值模拟,运用UDF技术严格控制阀门开启速度,研究了泄压阀内部的速度和声源功率级场的细节特征,以期为实验研究及工程实践提供理论参考。 1 计算模型和方法 1.1 泄压阀使用程序 本文研究的泄压阀模型来自于罗浮阀门集团生产的某型号全启式弹簧泄压阀。材料属性和操作条件均使用程序默认。泄压阀内工作介质为石油,具体参数如下: 石油密度 ρ = 0. 8 kg /m3,动力黏度 μ =0. 11 kg / ( m·s) ,操作压力为1. 01325 × 105Pa。进口边界条件为压力进口,出口为自由出流边界条件。 1.2 湍流动能方程 数学模型为二维雷诺时均Navier-Stokes方程和双方程标准k - ε 模型。标准k - ε 模型自从被Launder and Spalding提出后,就变成工程流场计算中主要工具。适用范围广、经济、合理的精度,它在工业流场和热交换模拟中有较为广泛的应用。 湍流动能方程k和扩散方程: 以上方程中Gk表示由层流速度梯度而产生的湍流动能; Gb是湍流动能; YM是由于在可压缩湍流中,过渡的扩散产生的波动; C1,C2,C3是常量; Sk和Se是用户定义的。 对于不可压黏性流体,水击泄压阀壁面定义为无滑移绝热固壁边界; 下游边界,假设压强已经恢复到来流压力,因此出口边界条件为自由出流。 1.3 网格数目的确定和网格设置 基于非稳态的计算方法,采用Fluent软件进行数值计算。速度和压力耦合采用SIMPLE算法求解,扩散项采用二阶精度的中心差分格式。为了使计算过程降低对处理器及内存的需求和具有较好的网格质量,流场采用分块耦合求解方法,对阀芯附近网格进行加密处理,并设置边界层。 对于模拟用泄压阀的模型和网格划分如图1 所示,为验证计算结果与网格数量的关联性,模拟中对网格数目为13 万、22 万和38 万三种情形进行了数值模拟和对比分析。结果表明,13 万网格与22 万网格模型对比,计算结果显示13 万网格时候结果不够精确; 将22 万网格与38 万网格模型对比,相比之下流动特征基本没有变化,表明22 万网格有足够的精度,且很大程度减小了计算量,因此采用22 万网格进行研究。 对整个流场区域采用结构网格相处理,各固体壁面为无滑移、无渗透条件。将石油考虑为常温下不可压缩流体,因而不考虑传热的影响,无需求解能量守恒方程。 选择动网格模型,并选择网格更新模型为: 弹簧光滑法和局部重构法。设置动网格区域,将阀芯所有对应的边设置为刚体,阀芯的运动规律默认为由UDF程序来进行控制。对所有区域进行初始化,迭代的时间步长设为0. 0005 s,时间步数为30,最大迭代数为2 000,设置每个时间步对数据自动保存,开始进行迭代计算。 1.4 宽带噪声源分布 采用宽带噪声源模型模拟区域的噪声源分布。因噪声没有固定的频率,在涉及宽频噪声的情形下,通过求解雷诺时均方程得到湍流参数的统计学分布,并采用声学类比的方法获得宽带噪声源的分布。声功率W和声功率级LW的表达式分别为: 式中: ρ0为流体密度,kg /m3; α0为声速,m/s; l为湍流尺度,m; u为流体速度,m/s; c为常数; W0为基准声功率,其值为0. 01 W。 2 声源功率级的分析 为研究泄压阀处于不同压力时的内流特征,对水击泄压阀内部速度分布以及声源功率级进行分析研究。假定管道压力分别为0. 1,0. 3,0. 5 MPa,经过实验研究,对应的阀芯速度分别为0. 03,0. 045,0. 060 m/s。取阀门开高5 mm和8 mm进行分析研究。 2.1 不同管道开度对泄压阀内部流场的影响 通过非定常数值模拟,可得到泄压阀开启不同阶段流场的机理,深入探讨泄压阀开启阶段内外流场的阀口机理,将会给工程实践提供具有指导意义的理论依据。 当管道压力为0. 1 MPa时,速度在泄压阀阀芯前后发生剧烈变化,工作流体经过泄压阀口时,速度快速升高,且速度梯度很大( 见图2a和图2b) 。同时由于有节流作用导致在泄压阀口处速度增加,动能增大,此时动压是全压的主要变现形式。当