基于有限元方法的液态金属子通道程序开发与验证  

Development and Validation of Subchannel Analysis Program for LMR Based on FEM

作　　者：朱昕阳 陈荣华[1] 王金顺 田家豪 秋穗正[1] 苏光辉[1] 田文喜[1] ZHU Xinyang;CHEN Ronghua;WANG Jinshun;TIAN Jiahao;QIU Suizheng;SU Guanghui;TIAN Wenxi(School of Nuclear Science and Technology,Xi’an Jiaotong University,Xi’an 710049,China)

机构地区：[1]西安交通大学核科学与技术学院,陕西西安710049

出　　处：《原子能科学技术》2024年第12期2485-2496,共12页Atomic Energy Science and Technology

基　　金：国家自然科学基金(12275205)。

摘　　要：液态金属冷却快堆因其核燃料增殖以及核废料嬗变能力而备受关注,为满足我国液态金属快堆的设计建设需求,本文基于多物理场耦合平台MOOSE开发了适用于液态金属冷却快堆的子通道分析程序FLARE,并与东芝(Toshiba)37棒束实验数据、KNS 37棒束实验数据以及欧洲铅冷快堆ALFRED设计限值进行了对比验证。本程序计算结果与Toshiba实验值和KNS实验值符合良好,并在ALFRED的计算中与同类程序SACOS-LMR结果相当,初步验证了本程序具备对液态金属快堆组件内关键参数进行准确计算的能力。本程序能为液态金属快堆组件的热工水力分析提供有效的设计与分析工具。Liquid metal reactors(LMRs)involved both sodium fast reactors(SFRs)and lead fast reactors(LFRs)have gained great development over the world due to their nuclear fuel breeding and nuclear waste transmutation capability.For the purpose of meeting the design and construction requirements of LMRs in China,a subchannel program FLARE based on the MOOSE framework for LMRs was developed.The program considers four conservation equations:mass,momentum,and energy in the axial direction,along with momentum in the transverse direction.Additionally,a hybrid approach that integrates the discontinuous Galerkin(DG)and finite element method(FEM)for discretizing conservation equations was employed in FLARE PROGRAM.Advanced closure models such as pressure drop,heat transfer,and material properties of sodium,lead bismuth,and lead were incorporated.In detail,the program integrates the Upgraded Cheng and Todreas(UCTD)wall friction model,Cheng and Todreas(CT)turbulence mixing model,CT transverse heat conduction model,along with some heat transfer correlations with a wide applicability range,while considering geometric effects and flow conditions.To solve the coupled nonlinear partial differential equations,the Jacobians-Free Newton-Krylov(JFNK)algorithm was used.To validate the accuracy of the constitutive models for swirling flow effect implemented in FLARE program,Toshiba 37-pin was computationally simulated,and the results from FLARE PROGRAM were compared with experimental data as well as results from similar programs.The results demonstrate that the computational results from the FLARE program exhibit good agreement with experimental data as well as calculations from similar programs.The maximum error at high and medium flow rates is 12.6%of the inlet-outlet temperature difference,while the error remains within 10%at low flow rates.Moreover,the correctness of the transverse heat conduction correction model was validated through low-flow conditions.To further validate the computational capability of the FLARE program,the KNS 37-pin e

关 键 词：子通道程序 MOOSE 液态金属 热工水力分析 

分 类 号：TL334[核科学技术—核技术及应用]