双管旋涡脱落的数值模拟研究.pdfVIP

核动力工程 2006年lO月 Nuclear Oct．2006 Power Engineering 双管旋涡脱落的数值模拟 卫嫒嫒，陆道纲，杨林民 (华flg电力大学核能科学与工程系) 摘要：采用任意的拉格朗日欧拉有限元方法，数值模拟双管在横向流作用下的旋涡脱落现象。模拟过程 中通过改变两管间距，揭示了两管间距与旋涡脱落频率的变化关系。数值模拟的结果和理论计算结果相一致。 验证了计算方法的有效性。 关键词：拉格朗日欧拉有限元；旋涡脱落 1前言 核电厂蒸汽发生器的传热管，尤其是其底部， 由于流体横掠管束有可能发生旋涡脱落诱发振动 算网格的排列。一共用了4249个节点和4040个 的现象。目前，对单管的旋涡脱落现象的研究已 单元。图1b描述了管子的外围情况。图2给出了 当雷诺数为1．6x104时的速度向量。 经很多，A．Yamaguchi[1】采用任意的拉格朗日欧拉 有限元方法(ALE．FEM)分析了单管的旋涡脱落， 位于下方圆管上的节点2(图lb)的压力．时间 记录如图3。由图3可以看出压力随时间的变化 并验证了ALE。FEM在解决流固耦合问题上的有 效性；但是蒸汽发生器传热管是由管排或者管群 具有一定的周期性，由图中可以近似读出这个周 构成的，有必要对双管或多管条件下的旋涡脱落 期T≈0．072s。于是旋涡脱落的频率也就是压力 现象进行研究，揭示它们之间相互作用规律。目 1 变化的周期Z=去≈13．89Hz。 前国内采用ALE．FEM分析双管的研究还很少见。 』 Takehiko Ichioka等人[21采用计算流体力学分析了另一方面，文献【3】给出了单管的旋涡脱落频 单管和管束的旋涡脱落。本文采用ALE．FEM针率理论计算公式 对双管的旋涡脱落特性进行研究。 ‘(1) 墨=ud 2ALE—FEM理论简介 式中，墨是斯特罗哈数，雷诺数为1．6x104时 在ALE．FEM用到了三个坐标系统，欧拉坐 取墨=o．2；d为管子直径；U为来流速度。取d 标系、拉格朗日坐标系、和任意拉格朗日欧拉坐 为22mm，U为1．0m／s，代入式(1)得到旋涡脱落 标系。欧拉坐标是在静止的坐标系中描述流体的 频率Z=9．10Hz。此外，由式(1)计算出圆管直径 运动；拉格朗日坐标指的是坐标跟随流体的运动 Hz。对 为2d时的理论旋涡脱落频率Z=18．20 而运动；而任意拉格朗日坐标则介于两者之间， 既不静止，也不跟随流体，而是以一种任意设定 方式运动。 在欧拉有限元中，计算的范围是固定在空间 不移动的计算网格上，它不适用于移动边界问题。 在拉格朗日有限元中，计算网格随着流体而流动， 常常会导致计算网格畸变造成继续计算的困难。 引入ALE—FEM能够有效地解决这种困难。 3双管旋涡脱落的数值模拟 a整个流场 b管子的外围 根据研究对象，本文选取如