摘要： 采用基于非结构、非交错网格的有限体积法求解用两方程模型封闭的雷诺平均N-S 方程组，对碟阀的三维分离流动进行了数值模拟。结果表明，碟阀关闭角度大于20°时，在碟阀背面产生了流动分离，形成一对大小相等方向相反的旋涡。在碟阀下游，这两个旋涡相互作用，逐渐消失。碟阀的流动阻力系数随着碟阀关闭角度的增大呈指数性增长。
　　关键词： 碟阀；数值模拟；分离流动；阻力系数
　　作者简介：诸葛伟林，男，1973年生，博士研究生，研究方向为流体机械水动力学数值计算。通讯地址：430072 湖北武汉市武汉大学动力与机械学院。
　　　　　　　刘光临  蒋劲  符向前（武汉大学，湖北武汉430072）
　　1、引言
　　碟阀的流阻特性以往主要通过实验求得［1、2］，随着计算流体力学和计算机技术的发展，用数值计算的方法得到碟阀的流阻特性已经成为可能。与实验方法相比，用计算流体力学对碟阀的流动情况进行模拟不仅简便易行，而且还可以了解碟阀内部流场的详细情况，如压力流速分布、分离流动区域等。对于指导碟阀的设计、改善其流动状况、减小流动阻力具有重要意义。
    碟阀内部的流动情况比较复杂，尤其是在小开度的情况下存在大尺度的三维分离流动。本文采用雷诺平均N-S方程组和K-ξ两方程湍流模型对碟阀内部的流动进行数值模拟。
　　2、控制方程组
　　碟阀内部流动为不可压缩三维粘性流动，可以用不可压的雷诺平均方程组求解，采用K-ξ双方程模型构成封闭的方程组。
　　（1）连续性方程：
　　（2）动量方程组：
　　（3）湍流涡粘性系数计算公式
　　（4）湍流动能K的输运方程：
　　（5）湍流耗散率ξ的输运方程：
　　（6）湍流动能和耗散率生成项Gk和Gg的计算公式：
　　本文采用了基于非交错网格的有限体积法对控制方程组进行离散，控制体积的界面位于网格各个结点的中间平面上，利用压强校正法求解各个变量。壁面附近的流动采用壁面函数进行模拟。

3、边界条件
　　选取碟阀及其前后一段管道作为计算域，如图1所示。流动方向为+Z 方向。
　　设管道进口速度uin为1.0m/s，管道进口的湍流动能k和耗散率ξ由下式给定：
　　管道进口压力由内场直接外推，出口给定参考压力，其余参数由内场直接外推得到。在管道和阀瓣壁面上给定无滑移固壁条件。
　　4、计算网格
　　碟阀的形状并不复杂，但如果采用一般的结构网格，很难保证在阀瓣和管道交接处网格的质量。因此，本文采用了非结构网格，将阀瓣的周围的网格局部加密，以确保网格质量和求解精度。图2为管道中间截面网格图。
　　5、数值模拟结果分析
　　数值模拟的结果表明，当碟阀关闭角度大于20°时，在碟阀背面开始出现旋涡流动。如图3所示。
　　碟阀关闭角度为35°时，碟阀背面的旋涡流动更加明显。图4（a）~（c）分别为XY、XZ 及YZ 截面的流速矢量图。
　　从图中可以看出，由于碟阀背面存在局部低压区，从碟阀上方越过的流体部分折向下流，从碟阀下方流过的流体部分折向上流，在XY截面上形成两个大小相等，方向相反的旋涡，在YZ截面上可观察到有回流区存在，而XZ截面上尚未形成回流。XY截面上的两个旋涡在向下游发展过程中相互作用，逐渐减弱。
　　图5（a）为碟阀后不远处管道XY截面流速矢量图，图5（b）为接近出口时管道XY截面流速矢量图。由图可见，接近管道出口时，旋涡流动已经基本消失。
　　碟阀关闭角度继续增大，在XY截面上的回流区逐渐增大，控制整个流场结构。图6为碟阀关闭角度为85°时YZ截面流速矢量图。在阀瓣的上方和下方分别形成了两个很大的回流区，回流区中心压力最低。
　　根据不同关阀角度下的管道进、出口压力降计算出的碟阀阻力系数ζ与关阀角度φ之间的关系如图7所示，图中纵坐标为对数坐标。如图，碟阀阻力系数随着关阀角度的增大呈指数性增长， 这与参考文献［2］中实测数据的变化规律一致。
　　6、结论
　　本文采用基于非结构、非交错网格的有限体积法求解用K-ξ两方程模型封闭的N-S方程组，模拟碟阀的三维流动的方法是可行的。模拟的结果表明，碟阀关闭角度大于20°时，在碟阀背面开始出现旋涡；碟阀关闭角度较小时，在与管道轴线垂直的截面上形成的一对大小相等方向相反的旋涡控制着流场结构；碟阀关闭角度较大时，碟阀背面出现的大范围回流区控制着流场结构。在碟阀下游，旋涡相互作用，逐渐减弱，接近管道出口时旋涡流动基本消失。根据不同关阀角度碟阀的数值模拟结果计算得到的流动阻力系数与碟阀关闭角度的关系为指数增长关系。
　　参考文献
　　［1］ 宋汉武，吴荣方- 碟阀的流量特性和阻力特性的试验研究［J］- 发电设备，1995，（Z2）.
　　［2］ 徐明阳，龚力群- 液控缓闭碟阀流阻特性试验研究［J］.甘肃科学学报，1997，（4）：37-40