Subchannel analysis of CANDU-SCWR fuel

A cooperative study has been initiated at Xi'an Jiaotong University (XJTU) with Atomic Energy of Canada Limited (AECL) to develop a subchannel code ATHAS for preliminary analyses of flow and enthalpy distributions and cladding temperatures in CANDU fuel at super-critical water conditions. The code is applicable for transient and steady-state calculations. Then the paper uses the ATHAS code to analyze CANDU-SCWR which is operating at 25.0 MPa pressure. The results show that the maximum cladding-surface temperature of CANFLEX bundle is 804.1 °C, which is below the limit of design, and it is appropriate for use in the CANDU super-critical water-cooled reactor (SCWR) based on heat-transfer analysis.     

压力管式超临界水堆极限事故下“无堆芯熔化”概念评估

在自主研发的事故分析程序SCTRAN的基础上,开发并验证了二维导热模型和辐射换热模型,并将改进后的SCTRAN应用于加拿大压力管式超临界水堆在失水事故（LOCA）叠加丧失紧急堆芯冷却系统（LOECC）事故中的堆芯安全评估,并对燃料棒到慢化剂之间的传热效率以及关键的影响因素进行了评估。计算结果表明,在LOCA叠加LOECC工况下,燃料棒到燃料通道的辐射换热和燃料棒到蒸汽的自然对流换热能够有效导出反应堆的衰变余热,最高功率的燃料组件内、外圈燃料棒的最高包壳温度分别为1278℃和1192℃,均低于不锈钢包壳的熔化温度,因此整个事故过程中不会发生堆芯熔化。      

铅铋冷却快堆含绕丝燃料组件子通道程序开发与验证

由于铅铋冷却剂流动传热现象的复杂性,准确计算铅铋冷却含绕丝燃料组件的冷却剂和包壳温度是液态金属冷却快堆燃料组件热工分析的重点。本文基于集总参数法对守恒方程进行求解,开发了适用于铅铋冷却快堆的子通道分析程序,对液态铅铋在棒束燃料组件中的摩擦阻力模型、湍流交混模型和对流换热模型进行了适用性分析,并对7棒束大涡模拟和19棒束含绕丝传热实验进行了对比验证。结果表明：包壳和冷却剂温度的最大相对误差低于5%。程序能较好完成铅铋冷却含绕丝燃料组件的热工水力计算,可为铅铋冷却快堆设计提供支持。     

Numerical Simulation of Fluid Flow and Heat Transfer of Supercritical Fluids in Fuel Bundles

A supercritical water-cooled reactor (SCWR) was proposed as a kind of generation IV reactor in order to improve the efficiency of nuclear reactors. Although investigations on the thermal-hydraulic behavior in SCWR have attracted much attention, there is still a lack of CFD study on the heat transfer of supercritical water in fuel channels. In order to understand the thermal-hydraulic behavior of supercritical fluids in nuclear reactors, the local fluid flow and heat transfer of supercritical water in a 37-element fuel bundle has been studied numerically in this work. Results show that secondary flow appears and the cladding surface temperature (CST) is very nonuniform in the fuel bundle. The maximum cladding surface temperature (MaxCST), which is an important design parameter for SCWR, can be predicted and analyzed using the CFD method. Due to a very large circumferential temperature gradient in cladding surfaces of the fuel bundle, the precise cladding temperature distributions using the CFD method is highly recommended.     

超临界水堆子通道分析

超临界水堆作为6种第4代未来堆型中唯一的水冷堆,具有一些独特的特点,受到了广泛重视。本工作以上海核工程研究设计院的常规压水堆子通道程序为基础,开发编制了适用于超临界水堆的子通道程序,并对典型带有慢化剂水棒的超临界水堆燃料组件进行了模拟计算,得到了堆芯子通道内的温度、燃料棒包壳温度、表面传热系数等参数的分布规律。此外,研究了不同超临界流体换热关系式对计算结果的影响,结果显示,各传热关系式的计算结果存在一定差异。     

超临界水冷堆的安全分析

应用RELAP5-3D程序建立了超临界水冷堆(SCWR)的稳态模型,并在此基础上,分别对SCWR的两种瞬态和两种事故工况进行了分析。汽轮机旁路系统的存在可有效维持反应堆压力,保证反应堆安全。若SCWR失去给水,在辅助给水系统启动之前,向下流的水棒可通过热传导带走堆芯热量,并向燃料通道内提供冷却剂,缓解堆芯升温。因而,向下流的水棒体现了SCWR的安全性。主泵卡轴事故由于没有惰转,最热包壳温度值最大,因而主泵惰转可有效缓解包壳温度的升高。     

Sub-channel analysis of CANDU–SCWR and review of heat-transfer correlations

The paper presents the results of sub-channel analysis of CANDU–SCWR based on a wide review of heat transfer correlations. According to comparison with experiment data at different heat flux, Bishop Correlation is selected in SUBCHAN code to analyze CANDU–SCWR fuel channel. By detailed calculation of 43 fuel rods fuel channel in CANDU–SCWR, the paper gets the conclusion that the mass flux redistribution and reduction of heat-transfer coefficient at supercritical condition caused by the steep change of coolant density will limit the power of fuel channel in CANDU–SCWR.     

Assessment of Applicability of Two-Fluid Model Code ACE-3D to Heat Transfer Test of Supercritical Water Flowing in an Annular Channel

A supercritical-water-cooled reactor (SCWR) is a high-temperature, high-pressure water cooled reactor that operates above the critical pressure of water. In order to perform efficiently the thermal design of the SCWR, it is important to assess the thermal hydraulics in rod bundles of the core. Experimental conditions of mockup tests, however, may be limited because of technical and financial reasons. Therefore, it is required to establish an analytical design technique that can extrapolate experimental data to various design conditions of the reactor. Japan Atomic Energy Agency (JAEA) has improved the three-dimensional two-fluid model analysis code ACE-3D, which was originally developed for the two-phase flow thermal hydraulics of light water reactors, to handle the thermal hydraulic properties of water in the supercritical region. In the present study, heat transfer experiments of supercritical water flowing in a vertical annular channel around a heater pin, which were performed at JAEA, were analyzed with the improved ACE-3D to assess the prediction performance of the code. As a result, it was implied that the ACE-3D code is applicable to the prediction of wall temperatures of a single rod that simulates the fuel bundle geometry of the SCWR core.     

Development of a 2D conduction code for a finned fuel element analysis

The fuel element of KMRR (Korea Multi-purpose Research Reactor) has 8 longitudinal, rectangular fins to enhance the heat transfer performance. The existence of these fins makes it difficult to analyze the heat transfer phenomena within the fuel element using the conventional one-dimensional heat conduction model. As the uncertainty in the computation of the maximum sheath temperature significantly affects the core thermal margin, a computer code, called, TEMP2D, which is based on a two-dimensional heat conduction model has been developed to deal with the finned element and validated. This computer code TEMP2D has a fully implicit numerical scheme and can solve both the steady state and transient problems such as the changes in coolant thermal-hydraulic conditions and fuel pin power. The code accuracy, which proved to be an excellent one, was verified by comparing its results with those from two widely accepted computer codes, MARC and ADINA. The result of this code calculation has been used to compute the KMRR core thermal margin and to develop a correlation for the equivalent 1D heated diameter which can reproduce the maximum cladding temperature (or heat flux) at various steady states when used in the 1D heat conduction model.     

Numerical investigation of heat transfer in upward flows of supercritical water in circular tubes and tight fuel rod bundles

Heat transfer in upward flows of supercritical water in circular tubes and in tight fuel rod bundles is numerically investigated by using the commercial CFD code STAR-CD 3.24. The objective is to have more understandings about the phenomena happening in supercritical water and for designs of supercritical water cooled reactors. Some turbulence models are selected to carry out numerical simulations and the results are compared with experimental data and other correlations to find suitable models to predict heat transfer in supercritical water. The comparisons are not only in the low bulk temperature region, but also in the high bulk temperature region. The two-layer model (Hassid and Poreh) gives a better prediction to the heat transfer than other models, and the standard k– high Re model with the standard wall function also shows an acceptable predicting capability. Three-dimensional simulations are carried out in sub-channels of tight square lattice and triangular lattice fuel rod bundles at supercritical pressure. Results show that there is a strong non-uniformity of the circumferential distribution of the cladding surface temperature, in the square lattice bundle with a small pitch-to-diameter ratio (P/D). However, it does not occur in the triangular lattice bundle with a small P/D. It is found that this phenomenon is caused by the large non-uniformity of the flow area in the cross-section of sub-channels. Some improved designs are numerically studied and proved to be effective to avoid the large circumferential temperature gradient at the cladding surface.     

超临界水冷堆堆芯子通道稳态热工分析

超临界水冷堆(SCWR)作为6种第四代未来堆型中唯一的水冷堆,冷却剂出口温度可达500℃,具有良好的经济性.本文采用改进的COBRA-IV程序对超临界水冷堆方形组件子通道进行稳态热工分析.对计算结果进行分析可知:减小慢化剂通道中给水质量流量份额和加大慢化剂通道与相邻子通道之间的热阻,可以降低热管焓升,后者还可以得到较好的慢化效果.通过热通道的传热恶化分析发现,超临界水冷堆的设计不能避免传热恶化,必须精确计算传热恶化条件下的包壳温度才能确定包壳能否保证其完整性.     

Numerical analysis of heat transfer in supercritical water cooled flow channels

Research activities are ongoing worldwide to develop nuclear power plants with a supercritical water cooled reactor (SCWR) with the purpose to achieve a high thermal efficiency and to improve their economical competitiveness. However, there is still a big deficiency in understanding and prediction of heat transfer in supercritical fluids. In this paper, heat transfer of supercritical water has been investigated in various flow channels using the computational fluid dynamics (CFD) code CFX-5.6 to provide basic knowledge of the heat transfer behaviour and to gather the first experience in the application of CFD codes to heat transfer in supercritical fluids. Three different flow channels are selected, i.e. circular tubes, the sub-channel of a square-array rod bundle and the sub-channel of a triangular-array rod bundle. The effect of mesh structures, turbulence models, as well as flow channel configurations is analysed. Based on the present results, recommendations are made on the application of turbulence models to the heat transfer of supercritical fluids in various flow channels. A new definition for the onset of heat transfer deterioration is proposed. A strong non-uniformity of heat transfer is observed in sub-channel geometries. This non-uniformity has to be taken into account in the design of fuel assemblies of SCWR.     

Numerical investigation of supercritical water-cooled nuclear reactor in horizontal rod bundles

The commercial CFD code STAR-CD v4.02 is used as a numerical simulation tool for flows in the supercritical water-cooled nuclear reactor (SCWR). The basic heat transfer element in the reactor core can be considered as round rods and rod bundles. Reactors with vertical or horizontal flow in the core can be found. In vertically oriented core, symmetric characters of flow and heat transfer can be found and two-dimensional analyses are often performed. However, in horizontally oriented core the flow and heat transfer are fully three-dimensional due to the buoyancy effect. In this paper, horizontal rods and rod bundles at SCWR conditions are studied. Special STAR-CD subroutines were developed by the authors to correctly represent the dramatic change in physical properties of the supercritical water with temperature. In the rod bundle simulations, it is found that the geometry and orientation of the rod bundle have strong effects on the wall temperature distributions and heat transfers. In one orientation the square bundle has a higher wall temperature difference than other bundles. However, when the bundles are rotated by 90° the highest wall temperature difference is found in the hexagon bundle. Similar analysis could be useful in design and safety studies to obtain optimum fuel rod arrangement in a SCWR.     

数值反应堆堆芯通道级三维热工水力程序CorTAF开发及初步验证

堆芯是核动力系统的核心部件,其完整性是反应堆安全运行的重要前提。传统核反应堆堆芯热工水力分析方法无法满足未来先进核动力系统的高精度模拟需求。本文依托开源CFD平台OpenFOAM,针对压水堆堆芯棒束结构特点建立了冷却剂流动换热模型、燃料棒导热模型和耦合换热模型,开发了一套基于有限体积法的压水堆全堆芯通道级热工水力特性分析程序CorTAF。选取GE3×3、Weiss和PNL2×6燃料组件流动换热实验开展模型验证,计算结果与实验数据基本符合,表明该程序适用于棒束燃料组件内冷却剂流动换热特性预测。本工作对压水堆堆芯安全分析工具开发具有参考和借鉴意义。     

TACR热传输系统热工水力瞬态分析

钍燃料的利用对于缓解核燃料资源短缺具有重要意义,坎杜型反应堆(Canadian Deuterium Uranium,CANDU)在堆芯布置、中子利用效率及先进燃料循环方面具有较高的灵活性,使得其在CANDU反应堆中引入钍燃料循环更具现实意义。CANDU型反应堆中钍基燃料应用关键基础技术研究是加拿大与我国正在开展的合作课题,其中开发自主的CANDU堆堆芯热工水力设计和安全分析程序是钍基燃料应用必不可少的设计工作之一。本文针对CANDU型反应堆热传输系统结构特点,采用FORTRAN程序设计语言开发了适用于CANDU型反应堆热传输系统的热工水力瞬态分析程序CANTHAC(CANDU Thermal-Hydraulic Analysis Code)。利用CANTHAC对钍基先进CANDU堆(Thorium-based Advanced CANDU Reactor,TACR)进行了瞬态分析,计算工况包括满功率稳态、无保护蒸汽发生器(Steam Generator,SG)二次侧给水温度降低事故及完全失流事故。其中,满功率稳态计算结果与清华大学设计的钍基先进CANDU堆TACR设计值吻合较好,相对误差不超过2%,在可接受范围内;无保护SG二次侧给水温度降低事故及完全失流事故在计算条件下所得的燃料温度及系统压力等关键热工水力参数均在安全限值内,满足安全准则要求。程序为模块化编程,便于移植和改进,具有一定的通用性,为进一步研究工作奠定了基础。     

Thermal and stability considerations for a supercritical water-cooled fast reactor with downward-flow channels during power-raising phase of plant startup

This paper describes study on the procedure of raising the reactor thermal power and the reactor coolant flow rate during the power-raising phase of plant startup for the supercritical water-cooled fast reactor (SWFR), which is selected as one of the Generation IV reactor concepts. Since part of the seed fuel assemblies and all the blanket fuel assemblies of the SWFR are cooled by downward flow, the feedwater from the reactor vessel inlet nozzle to the mixing plenum located below the core is distributed among these fuel assemblies and the downcomer. The flow rate distribution as the function of both the reactor thermal power and the feedwater flow rate, which are the design parameters for the power-raising phase, is obtained by the thermal hydraulic calculations. Based on the flow rate distribution, thermal analyses and thermal-hydraulic stability analyses are carried out in order to obtain the available region of the reactor thermal power and the feedwater flow rate for the power-raising phase. The criteria for the “available” region are the maximum cladding surface temperature (MCST) and the decay ratio of thermal-hydraulic stability in three “hot” channels; two seed assemblies with upward/downward flow and a blanket assembly. The effects of various heat transfer correlations and axial power distributions are also studied.     

Heat transfer to supercritical fluids flowing in channels—empirical correlations (survey)

This literature survey is devoted to the problem of heat transfer of fluids at supercritical pressures including near critical region.

The objectives are to assess the work that was done in the area of heat transfer at supercritical pressures, to understand the specifics of heat transfer at these conditions, to compare different prediction methods for supercritical heat transfer in tubes and bundles, and to choose the most reliable ones.

The comparisons showed there is a significant difference in heat transfer coefficient values calculated according to various correlations. Only some correlations show similar results, which are quite close to the experimental data for normal supercritical heat transfer in water and carbon dioxide. Also, no one correlation can accurately predict the magnitude and onset of deteriorated heat transfer.

The exhaustive literature search, which included hundreds of papers, showed that the majority of correlations were obtained in tubes and just few of them in other flow geometries including bundles.

The variations in the prediction of supercritical heat transfer are related to the significant changes in thermophysical properties near the critical and pseudocritical points. Therefore, a discussion on the general trends of various thermophysical properties at near critical and pseudocritical points is also included.

Based on several chosen correlations, the heat transfer coefficients and temperature profiles in the CANDU-X reactor cooled with supercritical water were calculated.     

CFX simulation of high temperature thermal–chemical experiment: CS28-2

A computational fluid dynamics (CFD) model of a post-blowdown fuel channel analysis for aged CANDU reactors with crept pressure tube has been developed, and validated against a high temperature thermal–chemical experiment: CS28-2. The CS28-2 experiment is one of three series of experiments to simulate the thermal–chemical behavior of a 28-element fuel channel at a high temperature and a low steam flow rate which may occur in severe accident conditions such as a LBLOCA (large break loss of coolant accident) of CANDU reactors. Pursuant to the objective of this study, the current study has focused on understanding the involved phenomena such as the thermal radiation and convection heat transfer, and the high temperature zirconium-steam reaction in a multi-ring geometry. Therefore, a zirconium-steam oxidation model based on a parabolic rate law was implemented into the CFX-10 code, which is a commercial CFD code offered from ANSYS Inc., and other heat transfer mechanisms in the 28-element fuel channel were modeled by the original CFX-10 heat transfer packages. To assess the capability of the CFX-10 code to model the thermal–chemical behavior of the 28-element fuel channel, the measured temperatures of the fuel element simulators (FES) of three fuel rings in the test bundle and the pressure tube, and the hydrogen production in the CS28-2 experiment were compared with the CFX-10 predictions.     

重水堆停堆工况下单相自然循环流动与传热分析

使用RELAP5程序建立CANDU 6型重水堆模型,对停堆工况下主热传输系统环路内的单相自然循环进行了分析研究,并推导出重水堆单相自然循环流量模型。对Vijayan模型与RELAP5程序的自然对流传热模型（Churchill-Chu和McAdams模型）进行比较计算,结果表明,Vijayan模型计算的水平壁面传热系数低于程序模型,造成包壳温度略高,而竖直壁面传热系数则无明显差别。     

带格架四棒束超临界水流动传热数值分析

棒束内超临界水流动传热是超临界水堆堆芯热工水力研究的重要内容,但对其认识还十分有限。本文针对四棒束内超临界水的流动传热现象开展数值模拟,特别分析了定位格架对棒束通道内流动和传热的影响。结果表明,采用SSG湍流模型计算所得到的棒束壁面温度和实验结果吻合良好,定位格架的存在影响下游流体的速度分布,显著提高格架下游的传热特性,交混系数有大幅上升,使得加热棒周向壁面温度分布更加平均,最高温度出现位置发生改变。