en
×

分享给微信好友或者朋友圈

使用微信“扫一扫”功能。
作者简介:

方斌,男,硕士,工程师。主要研究方向:飞机燃油系统,管道流动。E-mail:fangbin@comac.cc

通讯作者:

方斌,E-mail:fangbin@comac.cc

中图分类号:V233.2

文献标识码:A

DOI:10.19416/j.cnki.1674-9804.2021.03.007

参考文献 1
张沁.圆管突然扩大与突然缩小的局部水头损失系数的研究[J].科学技术与工程,2014,14(28):159-62.
参考文献 2
门宝辉.管道突然扩大局部阻力系数初探[J].西北水资源与水工程,2000(1):13-6.
参考文献 3
张俊婷,杜秀丽,徐磊.基于CFD的管道流动局部水头损失系数的数值计算[J] 太原科技大学学报,2016,37(4):328-331.
参考文献 4
许自顺,咸凯,余建发,等.基于Fluent的低雷诺数突扩圆管流场的数值模拟[J].内燃机,2016(1):27-30.
参考文献 5
王战辉,张智芳,高勇,等.突扩管内部流体流动性能的数值模拟[J].化工科技,2019,27(3):6-9.
参考文献 6
张华.管道截面突变处的局部阻力系数[J].山西水利科技,1994(2):79-83.
参考文献 7
李栋浩,王文娥,葛茂生,等.突然缩小圆管局部水头损失系数试验研究[J].水利与建筑工程学报,2011,9(4):22-24.
参考文献 8
SHI Z D,LEE T S,WINOTO S H.Flow in a pipe with a ring-type obstacle[J].International Journal of Computational Fluid Dynamics,1996(6):225-237.
参考文献 9
Anupriya S,Jayanti S.Amechanistic model for expansion loss coefficient in upward vertical annular flow[J].Applied Mathematical Modelling,2018,60:552-570.
参考文献 10
王洪伟.我所理解的流体力学[M].北京:国防工业出版社,2014.
参考文献 11
高殿荣,张齐生.平面突扩管瞬态流动的有限元模拟[J].燕山大学学报,1999(4):331-334.
参考文献 12
赵海燕,贾雪松,杨士梅,等.突扩管分离流场的数值模拟[J].科学技术与工程,2009,9(17):5238-5240.
参考文献 13
姚建军,吴君强,蒋文明.突扩管内部水环流动特性的数值研究[J].炼油技术与工程,2020,50(9):46-49.
参考文献 14
贾洲侠,徐国强,邓宏武,等.亚临界压力下航空煤油RP-3动力黏度测量[J].北京航空航天大学学报,2014,40(7):934-938.
目录contents

    摘要

    管道系统中设置的阀门、限流环等十分常见,装备限流装置的管道本质就是典型的连续突变管道。突变管道流动存在十分明显的局部水头损失,因此工程设计人员在进行管道系统设计时不得不考虑因管道突变而产生的局部水头损失;目前针对单个突变管道流动局部水头损失已有较为深入的研究,并已推导出相应的计算模型,但对于连续突变管路的水头损失研究较少,工程设计中对于连续突变管路的局部水头损失计算采用基于单个突变管道水头损失的叠加计算法,由于连续突变管路流动为极其复杂的湍流,基于单个突变管道水头损失的叠加计算法而得到的连续突变管路的局部水头损失结果误差较大。鉴于此,采用理论推导与实验相结合的方法,设计四种不同管径比的管道以及23种不同进口流量工况,并以航空煤油为介质进行实验,利用实验结果进行分析并对传统的突变管道水头损失的叠加计算法进行修正,提出了连续突缩突扩管道的局部水头损失修正计算模型,其计算结果相比于传统计算方法在精度上得到了显著性提高,为工程管道系统设计中的连续突变管路局部水头损失计算提供了有力依据。

    Abstract

    Valves and limiting rings are commonly used in current pipeline system, and the nature of the pipeline equipped with limiting devices is a typical continuous mutation pipeline. There is an obvious local head loss in the flow of the abrupt change pipeline, so the engineers have to consider the local head loss caused by the abrupt change of the pipeline in the pipeline system design. At present, the local head loss of single mutation pipe flow has been widely researched, and the corresponding calculating model has been set up. However, less study was conducted for continuous mutation pipe head loss. In engineering design, the superposition calculation method based on the head loss of single mutation pipeline is adopted for the local head loss calculation of continuous mutation pipe. The error of the local head loss of the continuous mutation pipeline is large based on this superposition calculation method of head loss of a single abrupt change pipeline because the flow of the continuous mutation pipeline is very complex. In view of this, this article adopts the method of combining theoretical derivation and experiment designs four kinds of pipes with different diameter ratio and 23 kinds of inlet flow conditions, carries on the experiment with aviation kerosene as the medium, analyzes the experimental results, and modifies the traditional superposition method of calculating the local head loss. A modified calculation model for local head loss of continuous sudden contraction and expansion pipeline is proposed. Compared with the traditional calculation method and experiment data, the accuracy of of the calculation result is significantly advanced, which provides an effective method for the calculation of local head loss of continuous sudden expansion pipeline in the design of complex engineering pipeline system.

  • 0 引言

  • 管道设计在动力燃油系统、水利设施、农业以及工业管道中应用十分广泛,在工程应用中管道系统通常比较复杂,流体在复杂管道中流动存在不同程度的流动损失,包括沿程阻力损失和突变管道部位的局部阻力损失等,其中突变管道部位的局部阻力损失最为严重,而流体的水头损失一直是管道工程设计人员所关注的问题,因此对突变管道部位的局部阻力损失的研究具有重要的实际意义[1]

  • 由于突变管道流动水头损失对管道系统设计具有重要影响,目前已经有大量针对突变管道流动特性的研究。门宝辉[2]通过实验研究了突扩管道的局部阻力系数随管径比的变化情况,并用理论推导的方法表明突扩管道的局部阻力系数不仅与管径比有关,还与流体速度有关,但未指明局部阻力系数与流体速度的具体关系。张俊婷和许自顺等人[3-4]通过数值计算的方法对突扩、突缩以及渐变管道的局部阻力损失系数进行了研究,其研究结果表明突扩管道的局部水头损失系数随着雷诺数的增大而增大,随管径比值的增大而减小;但其研究结论同样未指出局部水头损失系数与雷诺数和管径比之间的具体变化关系,同时其数值计算结果并未设计实验验证,计算结果的准确性有待研究。王战辉等人[5]通过数值模拟的方法研究了突扩管道在不同管径比和进口流速下流场中压力和速度分布情况,研究结果表明,在截面突变处流动出现死区,并且突扩截面前后压差随管径比和进口速度的增大而增大。张华[6]通过理论分析得到了突变管路截面速度分布与局部阻力系数之间的关系,计算准确性较好。李栋浩等人[7]设计实验研究了雷诺数及收缩比对局部阻力系数的影响,实验数据表明在不同收缩比情况下,局部阻力系数与雷诺数的对数值成明显的线性关系。

  • 通过总结前人对突变管路流动特性的研究可以发现,目前针对于管径比以及流体进口速度对突变管道流动水头损失影响的研究较多,并且多集中在定性研究,对建立经验模型的研究较少;此外,目前对于突变管道的研究基本都是针对单个的突扩突缩管道,对于连续的突变管道的流动研究较少。

  • 鉴于此,本文将对连续的突扩突缩管道的流动特性展开研究,流体为RP-3航空煤油,分别研究四种不同管径比的连续突变管道在不同进口速度下的流动特性,并总结得到连续突变管道流动水头损失的经验模型。

  • 1 突变管道局部水头损失理论分析

  • 1.1 管道模型

  • 在实际工程应用中,管道系统上通常设置有较多限流装置,如阀门、限流孔等,这种带限流装置管道的基本特征是突变前后的管径相同,并且突变部分的长度较短;因此,本文将以这种典型的管道结构为对象展开研究,这种典型带限流装置管道的剖面结构如图1所示。

  • 图1 限流管道剖面示意图

  • 流体先以入口速度V1进入到突缩段,随后立即进入到突扩段,其中突缩前部分1-1截面和突扩后部分3-3截面的管相同,均为D1,缩小部分2-2截面的管径为D2。流体在短距离管道中连续经过突缩和突扩两级突变,流动状态为复杂的湍流,流体经过突变管道后存在不可恢复的压力损失,并且进口雷诺数越大,不可恢复的压力损失也越大[8]

  • 由于管道中的限流装置长度较短,可以忽略其中的沿程阻力损失,水头损失主要由突变截面的局部损失造成。

  • 1.2 传统计算方法

  • 目前对于单个突扩或者突缩管道流动局部水头损失的研究较多,根据传统计算方法,突变管道的水头损失计算方法为[9]:

  • h=ξv22g
    (1)
  • 其中,v为单个突变管道的进口处流体速度,ξ为突变管道的局部损失系数,对于单个突扩管道和单个突缩管道,其局部损失系数分别为[10]:

  • ξ=1-A-A2ξ=0.51-A-A
    (2)
  • 其中,ξξ分别为单个突扩和突缩管道局部损失系数,AA分别为突变截面前后的管道流通面积。

  • 对于如图1所示的连续突缩和突扩管道流动水头损失,可利用传统方法进行叠加计算,具体计算方法为:

  • ht1=ξv122g+ξv222g
    (3)
  • 其中,v2为缩小管道部分2-2截面处的平均速度。

  • 由流量连续得到:(4)

  • 并且:

  • A1A2=D12D22
    (5)
  • 将管径比表示为:

  • β=D1D2
    (6)
  • 因此,图1所示的连续突缩和突扩管道流动水头损失的传统理论计算方法可进一步表达为:

  • ht1=v124g2β4-2β2-β-2+3
    (7)
  • 1.3 理论分析

  • 对于连续突缩突扩管道整体的流动能量损失分析,首先针对图1中的1-1截面到3-3截面建立能量方程:

  • z1+p1gρ1+v122g=z3+p3gρ3+v322g+ht2
    (8)
  • 其中,z1为1-1截面相对于基准面的垂直高度;p1为流体在1-1截面处的静压;ρ1为流体在1-1截面处流体的密度;z3为3-3截面相对于基准面的垂直高度;p3为流体在3-3截面处的静压;v3为3-3截面处流体的平均速度;ρ3为流体在3-3截面处流体的密度;ht2为流体从1-1截面到3-3截面之间的流动水头损失。当管道水平放置时:

  • z1=z3
    (9)
  • 对于1-1截面到3-3截面,由流量连续得到:

  • A1v1=A3v3
    (10)
  • 对于不可压缩流动,流体密度不变,因此:

  • ρ1=ρ3
    (11)
  • 在后续的推导和计算中,流体密度统一用ρ表示。

  • 将式(9)~式(11)代入到式(8)中得到流体从截面1-1到截面3-3的流动损失为:

  • ht2=p1-p3ρg+v122g1-A12A32
    (12)
  • 根据式(12)分析可以发现,连续突变管道中的流动局部水头损失由三部分决定,分别是突变管道前后流体的静压差、流体的入口速度以及突变前后管道的流通面积之比。

  • 对于本文中分析的典型的带限流装置的突变管道,由于突变前后的管道流通面一致,流体从截面1-1到截面3-3的流动损失为:

  • ht2=p1-p3ρg
    (13)
  • 从式(13)可以发现,典型带限流装置的突变管道流动的局部水头损失由突变管道前后的静压差决定。对于单个突变管道局部水头损失,传统计算方法中使用动量方程将压差近似转化为速度及管道尺寸参数,由此叠加计算得到连续突缩突扩管道水头损失必然存在误差;此外,由于在连续突缩突扩管道中,流体经过第一个突变后会产生旋涡,流体流动状态为复杂的湍流流动[11-12],由于两个突变截面间的距离较短,在第二个突变截面入口处流体速度分布极不均匀[13],此时由传统方法计算得到的流动水头损失误差进一步扩大。

  • 对于典型带限流装置的突变管道流动水头损失分析,由于管道长度较短,管道沿程阻力损失可以忽略,水头损失主要由突变截面的局部水头损失造成。因此,为了研究连续突变管道水头损失,就必须确定突变截面前后压差的影响因素及其影响规律,突变截面前后压差的影响因素可能包括流动参数、管道结构参数等,接下来将通过实验进行研究。

  • 2 实验设计

  • 2.1 实验管道

  • 为了研究典型带限流装置的突变管道流动水头损失,设计了如图2所示的实验管道。实验管道设计有四组不同管径比;在距离管道突缩截面上游500mm处设置有压力表,在距离管道突扩截面下游500mm处和600mm处分别设置有压力表和质量流量计,并可由流量值计算得到相应位置处的流体速度;测量装置距离突变截面超过6倍管径D2,在距离突变截面6倍管径位置处流动旋涡基本消失,因此可以避免流动旋涡对测量装置精度的影响[1]

  • 图2 实验管道示意图

  • 2.2 实验工况

  • 实验在常温常压条件下进行,管径比β有四种,每种管径比又设计有23种不同流量,具体如表1所示。

  • 表1 实验工况表

  • 为了实验结果的准确性,每次实验都等到流量计数值稳定以后才开始记录数据,并且每个工况点都在稳定后持续30s,然后在30s内共记录5组数据,最后取这5组数据的平均值作为该工况下的实验值。

  • 3 实验数据分析及模型修正

  • 3.1 数据分析

  • 首先分析管道进口处航空煤油的雷诺数对连续突缩突扩管道的局部水头损失的影响。根据流量计的测量数据可计算得到管道的进口速度v1,进一步根据管道尺寸可以获得管道进口处航空煤油的雷诺数,根据雷诺数的计算方法:

  • Re=ρvLμ
    (14)
  • 其中,ρμ分别为航空煤油的密度和动力黏度,本实验在常温下进行,此时取密度为780kg/m3,动力黏度为1 038 μPa·S[14];v为管道的进口速度v1;L为管道进口部分的特征长度,即进口部分的管径D1。四种管径比的理论和实验水头损失随进口雷诺数的变化如图3所示。

  • 图3 连续突变管道水头损失随雷诺数变化

  • 由图3可以看出,连续突缩突扩管道流动局部水头损失随管道进口雷诺数的增大而增大,但并未随管径比的变化而呈现出明显的递变规律;由传统的突变管路局部损失叠加计算法得到的水头损失理论比实验值偏小。由以下方法分析四种不同管径比的连续突变管道在23种不同进口流量下的水头损失,根据式(15)计算传统方法理论值相较于实验值的平均偏差,得到表2的数据。

  • δ=123i=123 hi 实验 -hit1hi 实验
    (15)
  • 表2 传统理论局部水头损失计算误差

  • 由表2中的数据可以发现,管径比越大,理论计算值与实验值的偏差越小,即管径比越大,传统的局部损失叠加计算法误差越小。

  • 3.2 模型修正

  • 对实验数据进一步分析发现,四种不同管径比的突变管道局部水头损失实验值、理论计算值及管径比与进口雷诺数存在较为明显的变化关系,即:

  • h实验 ht1β=f(lgRe)
    (16)
  • 采用多项式拟合的方法进行分析,发现四种不同管径比的突变管道局部水头损失随进口雷诺数的变化基本呈线性函数关系,具体为:

  • f-=-0.546lgRe+9.595,R2=0.6401f==-1.001lgRe+15.786,R2=0.8524f=-1.209lgRe+18.047,R2=0.9510f=-1.245lgRe+18.424,R2=0.9524
    (17)
  • 其中R2表示的是多项式拟合的方差,方差越接近1则拟合准确度越高;反之,则偏差越大。从多项式拟合结果分析发现,一号管道的偏差相对较大,二、三及四号管道误差相对较小。因此,可以总结得到突变管道局部水头损失实验值、理论计算值及管径比随进口雷诺数变化的一般形式为:

  • f=algRe+b
    (18)
  • 其中,ab为线性函数的两个待定系数,从拟合得到的函数关系发现四种不同管径比的管道的待定系数ab变化较大,初步分析认为两个待定系数与管道的管径比相关,再次利用多项式拟合的方法分析待定系数ab与管径比的关系,得到:

  • a=1.22β2-4.5β+2.91,R2=0.9995b=-12.59β+42.5,R2=0.9732
    (19)
  • 因此,综合以上可以得到:

  • f=1.22β2-4.5β+2.91lgRe-12.59β+42.5
    (20)
  • 即:

  • h实验 =ht11.22β-4.5+2.91βlgRe-12.59+42.5β
    (21)
  • 由此可以得到修正的连续突缩突扩管道局部水头损失计算模型为:

  • h修正 =ht11.22β-4.5+2.91βlgRe-12.59+42.5β
    (22)
  • 其中修正因子为:

  • γ=1.22β-4.5+2.91β-1lgRe-12.59+42.59β-1
    (23)
  • 为了验证修正模型的准确性,利用修正模型计算连续突变管道局部水头损失,并与实验值及传统经验模型计算值进行对比验证,四种管道的对比结果分别如图4~图7所示。

  • 图4 管道一模型计算结果对比

  • 图5 管道二模型计算结果对比

  • 图6 管道三模型计算结果对比

  • 图7 管道四模型计算结果对比

  • 分析图4~图7可以发现,对于连续突缩突扩管道的局部水头损失计算,修正模型的计算结果相比于传统理论的叠加计算法精度有显著的提高,修正模型与传统理论的连续突缩突扩管道的局部水头损失计算误差对比结果如表3所示。

  • 表3 修正模型局部水头损失计算误差对比

  • 从表3可以看出,修正模型的计算结果误差均在10%以内,但修正模型计算相比于实验值均略微偏大。总之,修正模型相比于传统方法能更准确的描述连续突缩突扩管道的局部水头损失。

  • 4 结论

  • 本文对典型连续突变管道局部水头损失的研究结果表明:

  • (1)连续突缩突扩管道局部水头损失随进口雷诺数的增大而增大;

  • (2)对于连续的突变管道的局部水头损失计算,基于传统理论的突变管道局部水头损失叠加计算法误差较大,且计算结果比实际值偏小,可能的原因是传统理论的突变管道局部水头损失计算法未考虑短距离突变时湍流涡的能量损失;

  • (3)基于传统理论的突变管道局部水头损失叠加计算法的计算误差随管径比的增大而减少;

  • (4)本文结合实验数据分析及传统理论的突变管道局部水头损失叠加计算法,推导得到了连续突缩突扩管道局部水头损失计算的修正模型,该修正模型有效提高了计算精度,修正模型计算结果与实验数据对比表明修正模型误差低于10%,为典型的连续突缩突扩管道的局部水头损失计算提供了有效方法。

  • 参考文献

    • [1] 张沁.圆管突然扩大与突然缩小的局部水头损失系数的研究[J].科学技术与工程,2014,14(28):159-62.

    • [2] 门宝辉.管道突然扩大局部阻力系数初探[J].西北水资源与水工程,2000(1):13-6.

    • [3] 张俊婷,杜秀丽,徐磊.基于CFD的管道流动局部水头损失系数的数值计算[J] 太原科技大学学报,2016,37(4):328-331.

    • [4] 许自顺,咸凯,余建发,等.基于Fluent的低雷诺数突扩圆管流场的数值模拟[J].内燃机,2016(1):27-30.

    • [5] 王战辉,张智芳,高勇,等.突扩管内部流体流动性能的数值模拟[J].化工科技,2019,27(3):6-9.

    • [6] 张华.管道截面突变处的局部阻力系数[J].山西水利科技,1994(2):79-83.

    • [7] 李栋浩,王文娥,葛茂生,等.突然缩小圆管局部水头损失系数试验研究[J].水利与建筑工程学报,2011,9(4):22-24.

    • [8] SHI Z D,LEE T S,WINOTO S H.Flow in a pipe with a ring-type obstacle[J].International Journal of Computational Fluid Dynamics,1996(6):225-237.

    • [9] Anupriya S,Jayanti S.Amechanistic model for expansion loss coefficient in upward vertical annular flow[J].Applied Mathematical Modelling,2018,60:552-570.

    • [10] 王洪伟.我所理解的流体力学[M].北京:国防工业出版社,2014.

    • [11] 高殿荣,张齐生.平面突扩管瞬态流动的有限元模拟[J].燕山大学学报,1999(4):331-334.

    • [12] 赵海燕,贾雪松,杨士梅,等.突扩管分离流场的数值模拟[J].科学技术与工程,2009,9(17):5238-5240.

    • [13] 姚建军,吴君强,蒋文明.突扩管内部水环流动特性的数值研究[J].炼油技术与工程,2020,50(9):46-49.

    • [14] 贾洲侠,徐国强,邓宏武,等.亚临界压力下航空煤油RP-3动力黏度测量[J].北京航空航天大学学报,2014,40(7):934-938.

  • 参考文献

    • [1] 张沁.圆管突然扩大与突然缩小的局部水头损失系数的研究[J].科学技术与工程,2014,14(28):159-62.

    • [2] 门宝辉.管道突然扩大局部阻力系数初探[J].西北水资源与水工程,2000(1):13-6.

    • [3] 张俊婷,杜秀丽,徐磊.基于CFD的管道流动局部水头损失系数的数值计算[J] 太原科技大学学报,2016,37(4):328-331.

    • [4] 许自顺,咸凯,余建发,等.基于Fluent的低雷诺数突扩圆管流场的数值模拟[J].内燃机,2016(1):27-30.

    • [5] 王战辉,张智芳,高勇,等.突扩管内部流体流动性能的数值模拟[J].化工科技,2019,27(3):6-9.

    • [6] 张华.管道截面突变处的局部阻力系数[J].山西水利科技,1994(2):79-83.

    • [7] 李栋浩,王文娥,葛茂生,等.突然缩小圆管局部水头损失系数试验研究[J].水利与建筑工程学报,2011,9(4):22-24.

    • [8] SHI Z D,LEE T S,WINOTO S H.Flow in a pipe with a ring-type obstacle[J].International Journal of Computational Fluid Dynamics,1996(6):225-237.

    • [9] Anupriya S,Jayanti S.Amechanistic model for expansion loss coefficient in upward vertical annular flow[J].Applied Mathematical Modelling,2018,60:552-570.

    • [10] 王洪伟.我所理解的流体力学[M].北京:国防工业出版社,2014.

    • [11] 高殿荣,张齐生.平面突扩管瞬态流动的有限元模拟[J].燕山大学学报,1999(4):331-334.

    • [12] 赵海燕,贾雪松,杨士梅,等.突扩管分离流场的数值模拟[J].科学技术与工程,2009,9(17):5238-5240.

    • [13] 姚建军,吴君强,蒋文明.突扩管内部水环流动特性的数值研究[J].炼油技术与工程,2020,50(9):46-49.

    • [14] 贾洲侠,徐国强,邓宏武,等.亚临界压力下航空煤油RP-3动力黏度测量[J].北京航空航天大学学报,2014,40(7):934-938.

  • 微信公众号二维码

    手机版网站二维码

    我要投稿 投稿指南 联系我们 二维码
    TOP