对冲式核级阀流场与压力分布的数值模拟.docxVIP

对冲式核级阀流场与压力分布的数值模拟 原子能是一种高效的能源，在世界能源的组成中发挥着不可忽视的作用。因此，核计算机化的安全问题也被提到了新的水平。核电站水力系统属复杂大型水力控制系统,系统中许多位置严禁逆流现象发生,否则可能导致严重事故,因此,系统中需大量止回元件。 止回阀是流体控制工程中的重要安全元件,它通过截止流体逆流来保证流体定向流动,从而保护其他流体机械及管线的安全。因此,在核电站水力系统中各种止回阀被大量使用,如旋启式止回阀。旋启式止回阀有一圆盘式的塞子铰接在边缘上,通过回流关闭。由于它具有设计简单、价格低、经济性好、可靠性高、日常维修方便且尺寸范围广泛等优点,该阀门在核能系统中被广泛使用,由于活塞猛然关闭撞击阀体产生的超压已被作为实验测量和模型的课题来研究。由此可见,传统止回阀的运动部件(即阀头或阀瓣)在做启闭运动时速度往往很快,这是阀门自身功能特性所必需的。但过快的启闭会在止回阀阀头附近引起瞬间的液体压力波动(水击现象),其峰值通常能达到管道正常压力的数倍以上,还会因反射水击波使系统发生严重震荡;另一方面,阀瓣的过快运动会直接引起它与阀体的撞击,在管线流量或流速很大的情况下,这种撞击严重影响止回阀寿命,甚至直接使之破坏失效。对阀门的启闭过程进行动态仿真可很好地了解阀门在开启和关闭过程中内部流场、压力分布及阀门各位置的受力情况,为阀门的设计改进与可靠性分析提供依据。Leutwyler等通过CFD软件计算扭矩和力,对碟盘阀二维模型进行了模拟;朱奇等对百万千瓦汽轮机的主调节阀在不同开度下的稳态与非稳态流场进行了分析;McElhaney从设计角度分析了止回阀的运动特性;韩旭等对对冲式止回阀全开关闭状态时各种逆流流量下阀门内部的流场、压力分布和阀头受力情况进行了稳态仿真模拟。目前,对止回阀的仿真模拟,特别是应用动网格技术对其启闭运动特性进行实时仿真的研究很少。 对冲式止回阀从阀头受力机理上对传统止回阀进行了结构改进,利用差动活塞实现阀头与流体的对冲运动。所谓对冲运动是指阀头的运动方向与流体的运动方向相反,止回阀在流体顺流时能快速开启,在流体逆流时阀头能与流体相互形成阻尼,并完成关闭运动。因此,阀头与阀座的撞击力大幅减小,流体的水击现象也得到了显著减缓。通过压差控制阀头自动运动,是对传统止回阀的技术性改进,但其导流管喉部由于孔道较小,流体流过此处时会出现射流现象,导致局部流速过大,对阀头冲击不稳定,从而影响阀门工作的稳定性。本文对该阀进行局部分析,通过动态模拟,得到阀门逆流关闭过程中导流管喉部的流场、压力等参数的分布情况,从而分析其结构性对整个阀门的影响。 1 主框架和模型的构建 1.1 流管顶部互通 对冲式止回阀的结构示于图1。该阀门由阀体、带阀座的套管、支架、导流管、阀头、活塞、轴承、后盖、整流罩和弹簧等部分组成。 当流体顺流时,与导流管喉部相通的前腔室为低压区,活塞另一侧腔室及阀门入口通道部分为高压区,因活塞直径大于阀座直径,所以作用于活塞上的合力指向入口方向,阀头向入口方向移动,阀门打开。阀头的移动方向与流体流动的方向相反。 当流体逆流时,与导流管喉部相通的前腔室为高压区,活塞另一侧腔室及阀门入口通道部分为低压区,且活塞直径大于阀座直径,所以作用于活塞上的合力指向出口方向,阀头向出口方向移动,阀门关闭。阀头的移动方向与流体的流动方向相反。 1.2 ujphxjv 本文研究的阀门内部流场为湍流,其控制方程如下。 连续性方程: ??xi(ρui)=0??xi(ρui)=0 式中:xi为坐标分量;ρ为流体密度;ui为速度分量,m·s-1。 动量方程: ??t(ρui)+??xi(ρuiuj)=-?p?xi+??xj(v??xj(ρui)-ρˉu′iu′j) 采用标准k-ε方程处理雷诺应力项-ρu′iu′j封闭问题。 湍流动能k方程: ρdkdt=??xi((μ+μtσk)?k?xi)+p-ρε 湍流耗散率ε方程: ρdεdt=??xi((μ+μtσε)?ε?xi)+C1εεkp-C2ερε2k 其中:t为时间,s;p为修正压力;μ为动力黏度;μt为湍流黏性系数;在FLUENT中,作为默认值常数,C1ε=1.44,C2ε=1.92,湍流动能k与耗散率ε的湍流普朗特数分别为σk=1.0,σε=1.3。 在本文中,方程的离散采用有限体积法的全隐格式,进出口均为压力边界条件,在模拟时给定进出口的总压差,壁面视为无滑移边界条件。 1.3 局部流场分析 本文利用Gambit对止回阀进行二维建模及网格划分,模型共包含65.79万个网格,图2示出网格数为2.68万的模型网格示意图。由于模型具有对称性,为节约计算机资源,提高运算速度,最终仅选用对称轴一侧的半个模型进行局部流场分析。在每一时间步长内,每一节点上的速度和位置均与周围的框架通