铅铋自然循环热工水力分析程序开发  

Thermal-hydraulic Analysis Code Development for Natural Circulation of Lead-bismuth

作　　者：程道喜[1] 杜开文[1] 于松 齐晓光[1] 翟伟明[1] 周平[1] 马晓瑶 张学锋[1] CHENG Daoxi;DU Kaiwen;YU Song;QI Xiaoguang;ZHAI Weiming;ZHOU Ping;MA Xiaoyao;ZHANG Xuefeng(China Institute of Atomic Energy,Beijing 102413,China)

机构地区：[1]中国原子能科学研究院,北京102413

出　　处：《原子能科学技术》2022年第9期1898-1905,共8页Atomic Energy Science and Technology

摘　　要：铅铋堆内冷却剂的自然循环对于反应堆的正常运行以及事故工况下的堆芯热量导出均至关重要,相关热工水力分析工作对于支持设计及安审均有重要意义。通过对铅铋堆内一回路系统内主要部件,包括堆芯、热交换器、管道等建立热工水力物理模型,开发了适用于铅铋自然循环瞬态过程模拟的热工水力分析程序,并利用铅铋自然循环回路内开展的自然循环启动实验、功率台阶影响实验等的结果进行了程序的初步验证。结果表明,程序计算得到的结果与实验结果符合较好,能够较好模拟铅铋自然循环的瞬态过程。该程序可以为铅铋堆研发过程中自然循环热工水力分析工作提供支持。Natural circulation of coolant in the lead-bismuth reactor is essentially important to remove heat from reactor core under normal or accident conditions,and the thermal-hydraulic analysis of natural circulation is a significant work for the designing and licensing.Recently,some code development and validation work for natural circulation analysis of lead-bismuth had been carried out around the world.Most of these codes are first developed based on conventional fluids,in which some models are not adequate or necessary for lead-bismuth.Owing to the unique features of this coolant and lack of operation experiences,analysis code for natural circulation analysis is still under developing,and a lot of verification work still needs to be done.A thermal-hydraulic analysis code suitable for the simulation of transient process of natural circulation of lead-bismuth is developed based on the thermal-hydraulic physical models of main components in the primary loop of lead-bismuth reactors,such as core,heat exchanger,pipe and so on.The single-phase flow model was considered in the code and the momentum equation of the control volume is integrated along the closed natural circulation path to obtain the system momentum equation.Semi-implicit method was used to discretize the momentum and energy equations.In natural circulation analysis,system momentum equation was coupled with energy equation,and these nonlinear and coupled equations are solved iteratively.The temperature distribution of previous time step was used to solve the system momentum equation and the mass flowrate in the flow path was obtained which could be used to solve the energy equation of reactor components along the closed natural circulation path from core inlet.Repeat this process until error of temperature and flowrate between last iteration was less than 10^(-5).This code was programmed with FORTRAN 95 and adoptd modular structure which includes physical properties module,input-output module,reactor component module,thermal-hydraulic correlation module and

关 键 词：铅铋合金 自然循环 程序开发 瞬态分析 

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