铅铋自然循环热工水力分析程序开发
Thermal-hydraulic Analysis Code Development for Natural Circulation of Lead-bismuth
-
摘要: 铅铋堆内冷却剂的自然循环对于反应堆的正常运行以及事故工况下的堆芯热量导出均至关重要,相关热工水力分析工作对于支持设计及安审均有重要意义。通过对铅铋堆内一回路系统内主要部件,包括堆芯、热交换器、管道等建立热工水力物理模型,开发了适用于铅铋自然循环瞬态过程模拟的热工水力分析程序,并利用铅铋自然循环回路内开展的自然循环启动实验、功率台阶影响实验等的结果进行了程序的初步验证。结果表明,程序计算得到的结果与实验结果符合较好,能够较好模拟铅铋自然循环的瞬态过程。该程序可以为铅铋堆研发过程中自然循环热工水力分析工作提供支持。Abstract: Natural circulation of coolant in the leadbismuth reactor is essentially important to remove heat from reactor core under normal or accident conditions, and the thermalhydraulic 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 leadbismuth 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 leadbismuth. 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 thermalhydraulic analysis code suitable for the simulation of transient process of natural circulation of leadbismuth is developed based on the thermalhydraulic physical models of main components in the primary loop of leadbismuth 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. Semiimplicit 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, inputoutput module, reactor component module, thermalhydraulic correlation module and solution module. Preliminary validation was evidenced by using results of startup of natural circulation experiment and power step influence experiment carried out on the natural circulation loop of leadbismuth. According to the result, there is a good agreement between the code predicting and experiment, and the transient process of natural circulation can be well simulated. Relative errors of temperature between calculation results and experimental data are within 5% during the whole transient process, and the relative errors of flowrate are almost 10%. This code can support the analysis work of natural circulation during the development of leadbismuth reactors.
-
Keywords:
- lead-bismuth eutectic ,
- natural circulation ,
- code development ,
- transient analysis
-
-
[1] MA W M, BUBELIS E, KARBOJIAN A, et al. Transient experiments from the thermal-hydraulic ADS lead bismuth loop (TALL) and comparative TRAC/AAA analysis[J]. Nuclear Engineering and Design, 2006, 236: 1422-1444. [2] MA Weimin, KARBOJIAN A, SEHGAL B R. Experimental study on natural circulation and its stability in a heavy liquid metal loop[J]. Nuclear Engineering and Design, 2007, 237: 1838-1847. [3] AGOSTINI P, BERTACCI G, GHERARDI G, et al. Natural circulation of lead-bismuth in a one-dimensional loop: Experiments and code predictions[C]∥10th International Conference on Nuclear Engineering. Arlington: American Nuclear Society, 2002. [4] BORGOHAIN A, JAISWAL B K, MAHESHWARI N K, et al. Natural circulation studies in a lead bismuth eutectic loop[J]. Progress in Nuclear Energy, 2011, 53: 308-319. [5] ZHAO P, LIU Z, YU T, et al. Code development on steadystate thermalhydraulic for small modular natural circulation leadbased fast reactor[J]. Nuclear Engineering and Design, 2020, 52: 2789-2802. [6] FIORI F, ZHOU Z W. Assessment study of RELAP5/SCDAP capability to reproduce TALL facility thermal hydraulic behavior[J]. Nuclear Engineering and Design, 2015, 295: 15-26. [7] 陆道纲,张勋,郭超. 铅铋自然循环回路热损效应分析[J]. 原子能科学技术,2015,49(9):1599-1605.
LU Daogang, ZHANG Xun, GUO Chao. Analysis on heat loss effect of leadbismuth natural circulation loop[J]. Atomic Energy Science and Technology, 2015, 49(9): 1599-1605(in Chinese).[8] IDELCHIK I E. Handbook of hydraulic resistance[M]. New York: Hemisphere Publishing Corporation, 1986: 64-65. [9] CHENG S K, TODREAS N E. Hydrodynamic models and correlations for bare and wirewrapped hexagonal rod bundles-bundle friction factors, sub-channel friction factors and mixing parameters[J]. Nuclear Engineering and Design, 1986, 92: 227-251. [10] CHENG X, TAK N. Investigation on turbulent heat transfer to leadbismuth eutectic flows in circular tubes for nuclear applications[J]. Nuclear Engineering and Design, 2006, 236: 385-393. [11] 张贵勤,肖泽军. 液态金属在环管换热研究[J]. 中国科学,1995,25(2):212-218.
ZHANG Guiqin, XIAO Zejun. Heat transfer characteristic in annulus of liquid metal sodium[J]. China Science, 1995, 25(2): 212-218(in Chinese).[12] MIKITYUK K. Heat transfer to liquid metal: Review of data and correlations for tube bundles[J]. Nuclear Engineering and Design, 2009, 239: 680-687. [13] JAEGER W. Heat transfer to liquid metals with empirical models for turbulent forced convection in various geometries[J]. Nuclear Engineering and Design, 2017, 319: 12-27. [14] 杨世铭,陶文铨. 传热学[M]. 北京:高等教育出版社,2006. [15] S750物性参数表[R]. 江苏:江苏中能化学科技股份有限公司,2019.
计量
- 文章访问数: 191
- HTML全文浏览量: 9
- PDF下载量: 233