燃料棒接触问题的新型自适应预处理方法  

Novel Adaptive Preconditioner for Contact Problem in Fuel Rod

作　　者：刘威震 刘振海[2] 辛勇[2] 叶舒愉 张世全[1] Weizhen LIU;Zhenhai LIU;Yong XIN;Shuyu YE;Shiquan ZHANG(School of Mathematics,Sichuan University,Chengdu,Sichuan 610064,China;Science and Technology on Reactor System Design Technology Laboratory,Nuclear Power Institute of China,Chengdu,Sichuan 610000,China)

机构地区：[1]四川大学数学学院,四川成都610064 [2]中国核动力研究设计院核反应堆系统设计技术重点实验室,四川成都610000

出　　处：《计算物理》2024年第5期596-606,共11页Chinese Journal of Computational Physics

基　　金：国家自然科学基金(11971335);四川省自然科学基金面上项目(2023NSFSC0075);四川大学科研基金(2020SCUNL110);四川大学-中国核动力研究设计院联合创新基金(HG2022145)资助。

摘　　要：通过分析接触前后雅可比矩阵的稀疏结构,构造基于接触压力判断准则的新型自适应预处理方法。在接触压力为零时使用完全的雅可比矩阵作为Preconditioned Jacobian-free Newton-Krylov(PJFNK)方法的预处理矩阵,在接触压力非零时使用雅可比矩阵的块对角作为预处理矩阵。数值试验表明:本文的自适应预处理方法能够有效提升模拟效率,求解总时间与Multiphysics Object-Oriented Simulation Environment(MOOSE)中现有方法相比减少了15.2%~33.1%。The contact problem between the fuel core block and the cladding is a key point in the simulation of nuclear reactor fuel rods,and the inclusion of a contact model increases the difficulty of solving the set of algebraic equations.By analyzing the sparse structure of the Jacobi matrix before and after contact,a novel adaptive preprocessing method based on the contact pressure judgment criterion is constructed.The full Jacobi matrix is used as the preprocessing matrix for the Preconditioned Jacobian-free Newton-Krylov(PJFNK)method when the contact pressure is zero,and the block diagonal of the Jacobi matrix is used as the preprocessing matrix when the contact pressure is non-zero.Numerical experiments show that the adaptive preprocessing method in this paper can effectively improve the simulation efficiency,and the total solution time is reduced by 15.2%~33.1%compared with the existing methods in Multiphysics Object-Oriented Simulation Environment(MOOSE).

关 键 词：MOOSE PJFNK 燃料棒 力学接触 自适应预处理 

分 类 号：O242.1[理学—计算数学]