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

来源：核动力工程杂志2024年第5期北京时间：

作者：周夏峰, 钟昌明, 张杨奕, 章运山, 曾伟, 汤琪芬, 强胜龙, 宫兆虎

单位：1.华中科技大学，武汉，430074；2.中国核动力研究设计院，成都，610213；

摘要：核反应堆多物理场多尺度耦合研究是核能领域研究的难点和热点，尤其是针对核反应堆具有温度、功率、密度等物理量变化剧烈、耦合工况复杂的庞大多维、强非线性耦合系统，目前的耦合计算程序时常存在收敛慢甚至不收敛等问题，这给新一代耦合计算程序开发带来诸多挑战和困难。近年来华中科技大学虚拟反应堆耦合分析实验室基于算符分裂、Picard迭代和JFNK等多种耦合方法，初步开发了统一耦合计算框架COME。本文首先详细分析了COME中的耦合方法、总体框架和通用接口等主要特点，之后基于COME分别求解了核反应堆中子输运/扩散模型、堆芯热工子通道耦合模型、系统分析程序两相流耦合模型以及复杂物理热工耦合模型等多个耦合问题，并对比不同耦合方法的收敛性和计算效率等，为提高真实复杂多物理耦合程序的计算稳定性和收敛特性提供方法指导和开发建议。

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

基金资助：国家自然科学基金项目（12005073）