堆芯瞬态核热耦合数值稳定性研究  

Research on Numerical Stability of Transient Neutronics/Thermal-Hydraulic Coupling in Reactor Core

作　　者：赵波[1,2] 曹欣荣 李权[1] 齐飞鹏[1] 黄永忠 马强[1] 陈浩 陈飞飞[1] ZHAO Bo;CAO Xinrong;LI Quan;QI Feipeng;HUANG Yongzhong;MA Qiang;CHEN Hao;CHEN Feifei(Science and Technology on Reactor System Design Technology Laboratory,Nuclear Power Institute of China,Chengdu,Sichuan Province,610213 China;Harbin Engineering University,Harbin,Heilong jiang Province,150001 China)

机构地区：[1]中国核动力研究设计院核反应堆系统设计技术重点实验室,四川成都610213 [2]哈尔滨工程大学,黑龙江哈尔滨150001

出　　处：《科技创新导报》2021年第13期90-94,共5页Science and Technology Innovation Herald

摘　　要：本文基于格林函数节块法程序TNGFM与子通道程序COBRA-EN开发了三维瞬态核热耦合计算程序,利用压水堆弹棒算例对核热耦合数值稳定性进行了分析,并提出了改进Picard迭代方法。数值模拟结果表明,减小计算时间步长可提高耦合计算数值稳定性,但计算效率较低,在物理计算外迭代中引入近似反馈可有效提高耦合收敛性能并大大减少耦合计算迭代次数。A three-dimensional transient neutronics/thermal-hydraulics coupling calculation program is developed based on Green’s function nodal method program TNGFM and subchannel thermal-hydraulic analysis program COBRA-EN.The numerical stability of neutronics/thermal-hydraulics coupling is analyzed with a PWR rod ejection example and an improved Picard method is proposed.The numerical simulation results show that reducing the time step can improve the numerical stability of coupling computation with a low efficiency.Introducing approximate feedback into external iteration of neutronics can effectively improve the coupling convergence performance and greatly reduce the number of iterations of coupling computation.

关 键 词：瞬态 核热耦合 数值稳定性 计算效率 

分 类 号：TL351.1[核科学技术—核技术及应用]