中心型蝶阀流场的数值模拟研究

2013-09-11 黄国权 哈尔滨工程大学机电工程学院

  利用CFD软件FLUENT对中心型蝶阀流场进行数值模拟。计算模型采用不可压缩流动的雷诺时均方程组,紊流模型采用标准模型,离散方程的求解采用压力耦合方程组的半隐式方法(SIM-PLE算法)。对阀门在不同开度情况下,流场状况分别进行了数值模拟,分别得到速度场,压力场,速度矢量场。定性的给出了阀门在不同开度下重要部位的受力情况,比较直观的给出了阀门在不同工况下流道内部的速度分布,为分析阀门内部受冲击状况提供了依据。由对阀门开度的不断变化,得到了涡流的形成过程,及速度对涡流形成及扩展的影响。

1、前言

  蝶阀以其结构简单和适合用于大中口径管道的结构特点,得到了广泛的应用。应用传统的阀门设计方法已经很难满足蝶阀设计生产的需要,在这种情况下数值模拟方法应运而生,它能够从更深层次对蝶阀的内部流场状况进行模拟,将流体力学的相关理论真正应用到蝶阀设计过程中去,对蝶阀设计制造具有深远的意义。即是利用数值模拟软件FLUENT对典型蝶阀的内部流场进行数值模拟,从而为蝶阀改进设计提供依据。以期对蝶阀的设计制造起到一定的指导作用。

2、N-S方程

  黏性流体的运动方程的Navier-Stokes方程,简称N—S方程。(1)适用于可压缩黏性流体的运动方程如式(1)。(2)理想流体的运动方程—Euler方程。若不考虑流体的黏性,则由上式可得理想流体的运动方程—Euler方程如式(2)。N—S方程比较准确地描述了实际的流动,黏性流体的流动分析均可归结为对N-S方程的研究。由于其形式甚为复杂,实际上只有极少情况下可以求出精确解,故产生了通过数值求解的研究,这也是计算流体力学进行计算的最基本的方程。可以这么说,所有的流体流动问题,都是围绕对N—S方程的求解进行的。

中心型蝶阀流场的数值模拟研究(1)
中心型蝶阀流场的数值模拟研究(2)

3、中心型蝶阀

  蝶阀是用随阀杆转动的圆形板件作启闭件,往复回转90°左右来开启,关闭和调节流体通道的一种阀门。蝶阀主要用作截断阀,亦可设计成具有调节或截断兼调节的功能。目前蝶阀在低压大中口径管道上的使用越来越多。中心型蝶阀三维模型图,如图1所示。

中心型蝶阀流场的数值模拟研究

图1 中心型蝶阀三维模型

4、中心型蝶阀流场的数值模拟

  以典型中心型蝶阀DN1000为例,为保证流场的稳定性,取蝶阀及其前部管道L1=5D(D为管道直径)与其后部管道L2=10D一同作为计算域。网格划分采用了非结构混合网格技术,利用FLUNENT软件包中的前处理软件GAMBIT强大的网格划分功能,采用自适应的网格技术对流场进行调整,使其模拟出更加精细的流动。采用不可压缩流动的雷诺时均方程组,湍流模型采用标准k-ε模型;所有方程中的对流项均用二阶迎风格式离散,离散方程的求解采用压力耦合方程组的半隐式方法(SIMPLE算法),所有方程中的对流项均用二阶格式离散,所有方程的熟练残差均为0.0001,全流场计算了定常流动,得到了阀内流场的详细分布情况。利用GAMBIT建立计算模型。面的网格划分情况,如图2所示。利用Gambit建立阀门不同开度情况下的计算区域,阀转过了相应的角度即可。

阀板附近的网格划分情况

图2 阀板附近的网格划分情况

  利用FLUENT求解器求解。为了研究阀门的流场特性,按照阀门开度为100%,60%和45%的3种典型工况,将入口速度设置为1,模拟阀门全开,较大开度,半开3种典型状态,进行比较分析。

  4.1、阀门全开时流场模拟结果分析

  阀门全开,流速1m/s时压力模拟结果,如图3所示。

流速1m/s时压力图

图3 流速1m/s时压力图

  阀门全开,流速1m/s时流场流速模拟结果,如图4所示。阀门全开,流速1m/s时流场速度矢量模拟结果,如图5所示。阀门100%开度,水流速度1m/s时流场特性分析。

  (1)从压力图3可知,整个阀板和阀壁受力状况良好,受力较均匀,只是在阀板端部小区域内形成局部高压区,对阀板前部形成一定的冲击,但考虑到这种冲击发生在阀板的大尺寸方向,且上下对称,因此对蝶阀整体的影响不大。(2)从速度图4可知,流速的上下对称性较好,速度梯度比较温和,流速过度区域比较大。在阀板前端和后端形成的低速区亦在我们的意料之中,和实践的情况也较吻合,说明模拟结果可信度较高。(3)从速度矢量图5可知,整个流场基本以层流为主,整体看流场状况较为平稳。总体来看,蝶阀全开时,流速分布较均匀,整个流态相当平稳。

流速1m/s时速度图

  图4 流速1m/s时速度图

流速1m/s时速度矢量图

图5 流速1m/s时速度矢量图

  4.2、阀门60%开度时流场模拟结果分析

  阀门60%开度,流速1m/s时流场压力模拟结果,如图6所示。

流速1m/s时压力图

图6 流速1m/s时压力图

  阀门60%开度,流速1m/s时流场速度模拟结果,如图7所示。阀门60%开度,流速1m/s时流场速度矢量模拟结果,如图8所示。阀门60%开度,水流速度1m/s时流场特性分析。

流速1m/s时速度图

图7 流速1m/s时速度图

流速1m/s时速度矢量图

图8 流速1m/s时速度矢量图

  (1)从压力图6可知,阀板总体呈现前高压,后低压的状况,同时后部的低压区较为均匀,梯度不大,相比较前部的高压区又呈现出明显的上部低压,下部高压的特点,使得阀板的受力无论前后,还是上下都出现不对称的情况,不过由于此时的压力梯度并不太大,所以阀板虽然受力较复杂,但不会产生较为严重的问题。(2)从速度图7可知,阀板上下过流区域的流速过大,会对附近的管壁

  造成一定的冲击。阀板后部的层状低速区则清晰可见。(3)速度矢量图8可知,流场仍以层流为主,未见明显异常。

  4.3、阀门45%开度时流场模拟结果分析

  阀门45%开度,流速1m/s时流场压力模拟结果,如图9所示。

流速1m/s时压力图

图9 流速1m/s时压力图

  阀门45%开度,流速1m/s时流场速度模拟结果,如图10所示。

流速1m/s时速度图

图10 流速1m/s时速度图

  阀门45%开度,流速1m/s时流场速度矢量模拟结果,如图11所示。

流速1m/s时速度矢量图

图11 流速1m/s时速度矢量图

  阀门45%开度,水流速度1m/s时流场特性分析。

  (1)从压力图9上可知,以阀板为界,前面和后面分为典型的高压区和低压区,阀板的单面受力状况相对较简单,但由于前面和后面的压差较大,造成此时阀板整体的受力环境比较恶劣。(2)从速度图10可知,阀板两端过流区的流速较大,速度梯度也较大,对管壁形成一定的冲击,由于蝶阀背面存在局部低压区,从蝶阀上方越过的流体部分折向下流,从蝶阀下方流过的流体部分折向上流,在阀板背面靠下部分形成旋涡。(3)从速度矢量图11可知,涡流区依稀可见。

5、结论

  通过分析模拟结果,可以得出以下几点结论:(1)阀门全开时或者是开度较大时过流状况较好,主要表现在整体的速度梯度较小,无论对阀板还是管壁的冲击都较小,从其速度矢量场则可以看出此时整个流场以层流为主,流场状况较简单。(2)随着阀门开度的减小,首先是在阀板后部形成低压区,随着阀门开度的进一步减小,阀板边缘过流区的速度梯度明显增大,同时在阀门45%开度左右在阀板背部开始形成清晰可见的涡流区,使阀门内部的流场状况开始变的不稳定。