基于算符分裂、Picard和JFNK统一耦合框架COME求解不同核反应堆模型的研究  

Research on Solving Different Nuclear Reactor Models by Coupling Multiphysics Environment(COME)Based on Operator Spliting,Picard and JFNK Methods

作　　者：周夏峰 钟昌明 张杨奕 章运山 曾伟[1] 汤琪芬 强胜龙 宫兆虎[2] Zhou Xiafeng;Zhong Changming;Zhang Yangyi;Zhang Yunshan;Zeng Wei;Tang Qifen;Qiang Shenglong;Gong Zhaohu(Huazhong University of Science and Technology,Wuhan,430074,China;Nuclear Power Institute of China,Chengdu,610213,China)

机构地区：[1]华中科技大学,武汉430074 [2]中国核动力研究设计院,成都610213

出　　处：《核动力工程》2024年第5期7-18,共12页Nuclear Power Engineering

基　　金：国家自然科学基金项目(12005073)。

摘　　要：核反应堆多物理场多尺度耦合研究是核能领域研究的难点和热点,尤其是针对核反应堆具有温度、功率、密度等物理量变化剧烈、耦合工况复杂的庞大多维、强非线性耦合系统,目前的耦合计算程序时常存在收敛慢甚至不收敛等问题,这给新一代耦合计算程序开发带来诸多挑战和困难。近年来华中科技大学虚拟反应堆耦合分析实验室基于算符分裂、Picard迭代和JFNK等多种耦合方法,初步开发了统一耦合计算框架COME。本文首先详细分析了COME中的耦合方法、总体框架和通用接口等主要特点,之后基于COME分别求解了核反应堆中子输运/扩散模型、堆芯热工子通道耦合模型、系统分析程序两相流耦合模型以及复杂物理热工耦合模型等多个耦合问题,并对比不同耦合方法的收敛性和计算效率等,为提高真实复杂多物理耦合程序的计算稳定性和收敛特性提供方法指导和开发建议。The research on multi-physics and multi-scale coupling of nuclear reactor is challenging and prominent in the field of nuclear energy.Especially for the large-scale,multi-dimensional and strongly nonlinear coupling system with drastic changes in physical parameters such as temperature,power and density and complex coupling conditions,the current coupling calculation codes often have problems such as slow convergence or even non-convergence,which brings many challenges and difficulties to the development of next-generation coupling computational codes.In recent years,the Virtual Reactor Coupling Analysis Team(VRCAT)at Huazhong University of Science and Technology has developed a unified coupling computational framework called Coupling Multiphysics Environment(COME)based on various coupling methods,including operator splitting,Picard iteration,and Jacobian-free Newton-Krylov(JFNK).This paper first provides a detailed analysis of the main features of COME,including the coupling methods,overall framework,and common interface.Then,based on COME,several coupling problems such as neutron transport/diffusion model,core thermal sub-channel coupling model,two-phase flow coupling model within the system analysis code and complex physical-thermal coupling model are solved respectively,and the convergence and computational efficiency of different coupling methods are compared,so as to provide method guidance and development suggestions for improving the computational stability and convergence characteristics of real complex multi-physical coupling codes.

关 键 词：多物理场耦合 算符分裂 Picard迭代 JFNK方法 耦合框架 

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