您现在的位置: 网站首页 >  新闻资讯 > 技术文章

利用FLUENT对气体涡轮流量计内部流场进行数值模拟

利用FLUENT对气体涡轮流量计内部流场进行数值模拟,气体涡轮流量计是计量天然气、液化气、煤气等介质的速度式仪表,为了改善气体涡轮流量计的性能,为设计提供指导和方向,近年来一些学者利用CFD技术对其内部流场进行了研究。LavanteEV等[3]利用FLUENT对气体涡轮流量计内部流场进行数值模拟,并根据仿真结果解释实验过程中的现象。对前导流器引起的流量计压力损失进行数值计算和实验测量,从流动机理上解释了结构和压损之间的关系。LIZhifei等[6]利用数值模拟得到了导流器内部的速度场和压力场,并以减小压力损失为目标优化了导流器的结构。通过对气体涡轮流量计进行CFD仿真,研究不同流量下的压损值,并通过实验证明了数值模拟的有效性。对涡轮传感器内部的速度场和压力场进行了数值仿真,提出对前后导流器、叶轮叶片形状和页顶间隙的改进。
  上述研究中未涉及针对不同螺旋升角涡轮内流场的数值模拟,以及涡轮叶片螺旋升角的改变对仪表性能影响的研究。本文对安装35°和45°叶片螺旋升角涡轮的DN150型气体涡轮流量计的内流场进行数值模拟,通过模拟结果预测仪表的始动流量和压力损失,并利用实验证明预测的正确性,为涡轮叶片螺旋升角的进一步优选提供数值方法。
1数学模型及边界条件
  利用FLUENT软件对涡轮内流场进行数值模拟时,忽略天然气的密度变化,在0~1200m3/h内,介质流动速度远远小于声速(即马赫数远小于0.3),认为流体不可压缩,且假设流动中无热量交换,不考虑能量守恒方程。
1.1微分控制方程
  气体涡轮流量计内部流动为湍流黏性流动,满足连续性方程和黏性流体运动方程。
基本微分方程[9]:
连续性方程:
1.2 湍流模型选择由于雷诺应力项的加入使时均N-S方程不封闭,为了求解引入k-ε两方程湍流模型。两方程湍流模型有标准k-ε模型,Renormalization-group(RNG)k-ε模型,和可实现的k-ε模型。其中,RNGk-ε模型主要应用于旋转机械的流动问题,在大范围的湍流模拟中有较高的精度。该模型能够比较准确地模拟各种复杂流动,其中湍流黏度由下式确定:
1.3网格划分与定解条件
  根据流量计的实际工况分别在介质入口和出口处添加10倍管径的直管段,并把整个模型剖分为3个区域:入口管道,旋转区,出口管道。旋转区域又细分为涡轮转子和支架定子两个区域,定子和转子之间的耦合采用多参考MRF(MultipleReferenceFrame)模型。利用GAMBIT前处理模块对进、出口直管段采用结构化网格,而对旋转区采用非结构化网格进行划分以满足对叶轮内部复杂区域的网格描述,各块网格通过块之间的交界面拼接在一起。网格总数为30多万个四面体非结构化网格和100多万个六面体结构化网格,旋转区网格如图1所示。
  定解条件包括介质入口、出口和固壁边界的设置。入口处给定相应流量(1200m3/h)下的主流速度值;出口采用压力出口边界条件,出口压力相对大气压为0;进、出口管道内壁,支架均取无滑移固壁边界条件。叶轮部分采用旋转坐标系,给定相应流量下的叶轮转速,将叶片的吸力面和压力面以及轮毂定义为旋转壁面条件,在旋转壁面条件的定义中,按照MRF的要求,将旋转壁面的旋转速度定义为相对速度,并且相对周围流体速度为0。
2计算结果分析
2.1压力场分析
  流量计全压定义为入口全压与出口全压之差,通过全压分析能够直接反映仪表压损的大小。全压越大表明流体经过流量计后产生的压损越大,压损过大会导致流量计不能正常使用。进口全压一定时,出口全压Pout越大,则流量计的全压△P越小,压力损失越小。如图2(a)和图3(a)所示,35°涡轮出口全压要明显小于45°涡轮出口处的全压,这说明相同的工况下45°涡轮所产生的压损较小。
螺旋升角45°涡轮压力场分布
  涡轮叶片动压的分布和大小直接影响涡轮驱动力矩的大小,35°叶片所受动压明显小于45°叶片所受动压,说明在相同工况下45°螺旋升角涡轮能获得较大的驱动力矩,如图2(b)和图3(b)所示,与35°涡轮相比,较小的流量就可推动涡轮稳定旋转,从而使仪表进入线性工作区。由此可预测安装45°螺旋角涡轮的仪表能获得较小的始动流量。
2.2速度场分析
  当气体介质以充分发展的湍流经过涡轮时,35°涡轮的速度矢量方向变化较大且向壁面集中,使得与叶片直接作用产生推动力矩的速度矢量减少,如图(4a)所示,且在出口处速度衰减较大,间接说明介质流经涡轮后压损的增加,如图4(b)所示。而45°涡轮内部的速度矢量分布比较均匀,过流性较好,与叶片直接作用的速度矢量较多,产生较大的驱动力矩,如图(5a)所示,且在出口处速度衰减较小,如图(5b)所示。
螺旋升角35°涡轮压力场分布
螺旋升角45°涡轮速度场分布
气体涡轮流量计负压检定平台
3实验对比
  气体涡轮流量计的检定采用负压检测方法,如图6所示,由标准吸风装置产生负压使标准罗茨流量计和被检定的气体涡轮流量计同时测量,安装在被测仪表两端取压口处的U型管可以测量流量计进、出口处的压力,从而得到仪表的压力损失。
  利用黄金分割法选取0~1200m3/h范围8个流量点,在每一个流量点随机采集3组不同时刻的数据,包括标准罗茨流量计和被检定流量计的累积流量及其输出脉冲数,对每组数据进行算术平均得到流量点处的平均仪表系数。通过采集U型管压差装置的指示值记录每个流量点处的压力损失,检定结果如表2所示。
  利用多项式插值对表2中的数据进行密化,得到20组插值数据,通过3次B样条拟合得到仪表系数曲线和压力损失曲线。
始动流量以仪表系数进入线性区的最小流量来确定,在小流量区内安装45°螺旋升角涡轮的流量计在流量20m3/h左右即进入线性工作区;而安装35°螺旋升角涡轮的流量计则在流量150m3/h左右时才进入线性工作区,而且在线性工作区内也存在着明显的波动,如图7所示.
气体涡轮流量计系数拟合曲线图
  35°涡轮流量计在各工况点处的压损明显大于45°涡轮流量计,最大压损达到3500Pa以上,如图8所示。上述分析表明安装45°螺旋升角涡轮的流量计与安装35°螺旋升角涡轮的流量计相比具有较小的始动流量,较小的压力损失,而且仪表计量的线性度较好。
气体涡轮流量计压力损失拟合曲线图
4结论
  对螺旋升角为35°和45°的气体涡轮流量计旋转部件内流场进行数值模拟,分析描述其内部流动的压力场和速度场,预测出安装45°螺旋升角涡轮的流量计比安装35°螺旋升角涡轮的流量计具有较小的始动流量和压力损失。利用黄金分割法选取仪表流量范围内的检定点,通过仪表负压检定平台获得了仪表系数曲线和压力损失曲线,实验结果与数值仿真中的预测相吻合,表明数值模在流量计性能预测中的有效性。涡轮叶片的螺旋升角是影响仪表性能的关键参数,合理选择涡轮的叶片螺旋升角,优化结构可进一步改善仪表的性能。
点击次数:  更新时间:2018-05-30 18:20:26  【打印此页】  【关闭