用CFD方法研究反应堆堆芯和控制棒流致振动的力学机理.pdfVIP

用CFD方法研究反应堆堆芯和控制棒流致振动的力学机理 用CFD方法研究反应堆堆芯和控制棒流致振动的力学机理 尚智‘，赵钧2，陈硕2 1．核科学与系统工程系：2．航空航天工程系 上海交通大学机械与动力工程学院，上海，200030 摘要：将反应堆堆芯和控制棒简化为圆柱排列，用CFD方法对其中基本并列双圆柱元素进行数值研 究。利用这种方法直接计算出涡的对流和扩散。通过对流场中涡的运动物理机理的计算与描述，得 到流体运动绕流双圆柱作用力的变化，揭示出堆芯元件的流致振动力学机理。 关键词：CFD；堆芯；控制棒；流致振动；反应堆 1．引言 ．Hj涡流函数分析方法描述二维粘性不可压 核反应堆的堆芯内部的冷却剂的流动会诱 流体流动的Navier-Stokes方程，在直角坐标下 发堆芯棒束和燃料控制棒等堆芯内部部件的振 可表述为无量刚流函数涡量形式 动。这种流致振动问题在反应堆里已经引起广泛 IIJ 一十“一十V一=一V一缈 ㈣’ 的关注【l～5】，由此引发的安全性和经济性问 鲁+“警+V等=去V2缈a a)c 卸 Re 题也受到日益重视。目前对于流致振动问题做了 业+业：川， p7(2) 许多研究，有实验研究，也有理论数值研究，但 蕾+蒂2删 是这些研究的理论数值手段没有从流动涡的角 其中流函数、I，与速度分量(u。v)的关系表示为 度来研究揭示流致振动问题。 a∥ a∥ 流体流动从某种意义上可是说是涡的移动 Ⅳ=上，，=一。 (3) 砂 缸 过程，流动的涡在流动过程中，如果遇到固体壁 面，必然同固体壁面之间产生相互作用。而且流 对上述涡流函数方程的求解，许多学者提出 体流过固体构建以后，在固体构建之后都有涡街 了不同的解法，如Chorin【6】提出的随机涡法 的产生和变化，会在固体前后产生不同的流场和 求解，但此方法缺点是随流动时间的增加，所需 涡变，导致绕流拽力甚至对固体的侧向推力。如 内存和计算量会呈几何级数增大。为了更好地在 果这些力的变化剧烈，并且大小合适的话，按照 拉格朗日坐标下跟踪涡的移动和扩散， 弹性力学简化就相当于加载在固体棒状或者板 Christiansen【7】在其基础上提出了涡格子法。 状元件上的集中载荷或者分布载荷，这些载荷必 在此方法中，每一次新涡是从固体表面产生，由 然导致这些元件内部的应力变化或者微小的变 泊松方程(2)根据上一步速度场决定，每一次涡 形。这些小变化反过来又会影响固体元件外部流 量移动完成后，将涡量分配到欧拉坐标下的各个 体的涡变化和流场的改变，这种改变必然引起流 节点上，而新的速度场(或流函数)是通过在欧拉 体绕流拽力的变化，从而引起固体的振动。如果 坐标下利用有限差分方法求解泊松方程(2)来获 这种振动强度剧烈，还会引起严重的破坏事故。 得。初始速度场为势流场。这样整个计算无需对 即使是微小的振动，如果引起共振效应，那么它 非线性的涡方程直接求解，极大缩短了计算时间 的破坏也将是非常严重的。反应堆内部的燃料堆 【7】。 芯、控制棒