CFD analyses in tight-lattice subchannels and seven-rod bundle geometries of a Super Fast Reactor

This paper presents CFD analyses in heat unsymmetric subchannels and heat symmetric seven-rod bundle geometries of a Super Fast Reactor (Super FR) fuel assembly using STAR-CD. The purpose of CFD analyses in heat unsymmetric subchannels is to evaluate the effect of the power differences on the heat transfer in subchannels of the Super Fast Reactor. For heat symmetric seven-rod bundles, the effects of the gap clearance between the fuel rod and the assembly wall and the displacement of the fuel rod on the circumferential temperature distributions and Maximum Cladding Surface Temperature (MCST) are analyzed. The results show that larger power difference between fuel rods gives larger circumferential temperature difference of the hottest fuel rods. Considering cross flow between edge and ordinary subchannels, 1 mm gap between the fuel rod and the assembly wall is better for small MCST although the circumferential temperature difference in edge subchannel is large. MCST increases exponentially with the displacement. The relative error of displacement should be less than 1% if the allowable increment of MCST due to displacement is less than 6 °C.     

LOCA Analysis of Super Fast Reactor

This paper describes loss of coolant accident (LOCA) analyses of the Supercritical-pressure Water-Cooled Fast Reactor (Super Fast Reactor). The features of the Super Fast Reactor are high power density and downward flow cooled fuel channels for the improvement of the economic potential of the Super Fast Reactor with high outlet steam temperature. The LOCA induces large pressure and coolant density change in the core. This change influences the flow distribution among the downward flow parallel channels. It will affect the safety of the Super Fast Reactor. LOCA analysis of Super Fast Reactor is important to understand the safety features of the Super Fast Reactor. Keeping the flow rate in the core is important for the safety of the Super Fast Reactor. In LOCA, it is difficult to maintain an adequate flow rate due to the once-through coolant cycle and the downward flow cooled fuel assemblies. Therefore, the early actuation of the Automatic Depressurization System (ADS) and reduction of the maximum linear heat generation rates of the downward flow seed fuel assemblies and Low-Pressure Core Spray (LPCS) system are necessary for the Super Fast Reactor to cool the core under LOCA. Analysis results show that the Super Fast Reactor can satisfy the safety criteria with these systems.     

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.     

Principle of rationalizing the criteria for abnormal transients of the Super LWR with fuel rod analyses

A detailed fuel rod design is carried out for the first time in the development of Supercritical-pressure Light Water Reactor (Super LWR). The fuel rod design is similar to that of LWR, consisting of UO₂ pellets, a gas plenum and a Stainless Steel Cladding. The principle of rationalizing the criteria for abnormal transients of the Super LWR is developed. The fuel rod integrities can be assured by preventing plastic strains on the cladding, preventing the cladding buckling collapse, and keeping the pellet centerline temperature below its melting point. The FEMAXI-6 fuel analysis code is used to evaluate the fuel rod integrities in abnormal transient conditions. Detailed analyses have shown that allowable limits to the maximum fuel rod power and maximum cladding temperature can be determined to assure the fuel integrities. These limits may be useful in the plant safety analyses to confirm the fuel integrities during abnormal transients.     

Feasibility study of the application of exotic pin for tight-lattice fuel assembly

A tight-lattice fuel assembly having less space for the coolant is more feasibly applied in Liquid Metal Fast Breeder Reactor (LMFBR). The thermal hydraulic constraint due to smaller coolant space can be compensated by the high heat capacity of the liquid metal coolant. A tight pin configuration provides high fuel volume fraction which eventually gives better neutronic performance for longer core lifetime. A cylindrical pin array provides less flexible arrangement for tight-lattice assembly, which results in very narrow coolant gaps connecting its neighboring subchannels. Therefore, the so-called exotic pin shape is introduced, which enable to distribute the coolant flow more uniformly, to be applied in tight-lattice bundles with sodium coolant. As Nusselt number and wall friction correlation are absent for this type of geometry, CFD calculations are performed by employing k-ε turbulent model.     

CFD analysis of thermal-hydraulic behavior in SCWR typical flow channels

Investigations on thermal-hydraulic behavior in SCWR fuel assembly have obtained a significant attention in the international SCWR community. However, there is still a lack of understanding and ability to predict the heat transfer behavior of supercritical water. In this paper, CFD analysis is carried out to study the flow and heat transfer behavior of supercritical water in sub-channels of both square and triangular rod bundles. Effect of various parameters, e.g. thermal boundary conditions and pitch-to-diameter ratio on the thermal-hydraulic behavior is investigated. Two boundary conditions, i.e., constant heat flux at the outer surface of cladding and constant heat density in the fuel pin are applied. The results show that the structure of the secondary flow mainly depends on the rod bundle configuration as well as the pitch-to-diameter ratio, whereas, the amplitude of the secondary flow is affected by the thermal boundary conditions, as well. The secondary flow is much stronger in a square lattice than that in a triangular lattice. The turbulence behavior is similar in both square and triangular lattices. The dependence of the amplitude of the turbulent velocity fluctuation across the gap on Reynolds number becomes prominent in both lattices as the pitch-to-diameter ratio increases. The effect of thermal boundary conditions on turbulent velocity fluctuation is negligibly small. For both lattices with small pitch-to-diameter ratios (P/D < 1.3), the mixing coefficient is about 0.022. Both secondary flow and turbulent mixing show unusual behavior in the vicinity of the pseudo-critical point. Further investigation is needed. A strong circumferential non-uniformity of wall temperature and heat transfer is observed in tight lattices at constant heat flux boundary conditions, especially in square lattices. In the case with constant heat density of fuel pin, the circumferential conductive heat transfer significantly reduces the non-uniformity of circumferential distribution of wall temperature and heat transfer, which is favorable for the design of SCWR fuel assemblies.     

Analysis of fuel rod behavior under normal operating conditions in Super Fast Reactor

A Super Fast Reactor is a pressure-vessel type, fast spectrum supercritical water-cooled reactor (SCWR) that is presently researched in a Japanese project. A preliminary core has been designed with 1.59E+06 W/m³ of power density [1]. In order to ensure the fuel rod integrity, the fuel rod behaviors under the normal operating conditions are analyzed using FEMAXI-6 code. Three types of the limiting fuel rods, with the maximum cladding surface temperature (MCST), maximum power peak (MPP) and maximum discharge burnup (MDB), are chosen to cover all the fuel rods in the core. The power histories of these fuel rods are taken from the neutronics calculation results in the core design. The available design range of the fuel rod design parameters, such as the initial gas plenum pressure, gas plenum length, grain size and pellet-cladding gap size, are found out in order to satisfy the following design criteria: (1) Maximum fuel centerline temperature should be less than 1900 °C. (2) Maximum cladding stress in circumstance direction should be less than 100 MPa. (3) Pressure difference on the cladding should be less than 1/3 of buckling collapse pressure. (4) Compressive stress to yield strength ratio should be less than 0.2. (5) Cumulative damage fraction (CDF) on the cladding should be less than 1.0. Finally the improved fuel rod design is proposed.     

An analytical method of evaluating the effect of rod displacement on circumferential temperature and heat flux distributions in turbulent flow through unbaffled rod bundles

An analytical method of evaluating the circumferential variations of temperature and heat flux fields inside and around a displaced fuel rod in triangular rod bundles in turbulent flow is presented with illustrative examples. The analysis consists mainly of the derivation of the simultaneous solutions of a set of heat conduction equations for fuel, cladding and coolant under the assumption of fully developed flow and heat transfer conditions. The local coolant velocity distribution, which is necessary for deriving the temperature field in coolant, is determined by solving the Navier-Stokes equation and the turbulent mixing of coolant is taken into consideration. The results show how the circumferential variations in the temperature and heat flux fields on the outer surface of the cladding increase the lower the ratio and the larger the fuel rod displacement due to thermal conduction and peripheral coolant flow velocity distribution.     

CFR600泵支承氩气空间部分三维传热数值模拟

利用流体力学软件Fluent对中国示范快堆（CFR600）氩气空间及相关结构部件进行稳态共轭传热计算,验证分析泵支承结构位于氩气空间部分的热工设计。模拟结果显示,原型设计方案泵支承结构位于氩气空间部分周向温度分布有较大温度梯度,泵支承结构靠近主容器支承径侧温度明显高于远离支承径侧,改进的结构设计减小了泵支承周向的温度梯度,为优化泵支承结构设计提供了技术支撑。     

International Solvent Extraction Conference ISEC ′90

The supercritical-water-cooled power reactor (SCPR) is expected to reduce power costs compared with those of current LWRs because of its high thermal efficiency and simple reactor system. The high thermal efficiency is obtained by supercritical pressure water cooling. The fuel cladding surface temperature increases locally due to a synergistic effect from the increased coolant temperature, the expanded flow deflection due to coolant density change and the decreased heat transfer coefficient, if the coolant flow distribution is non-uniform in the fuel assembly. Therefore, the SCPR fuel assembly is designed using a subchannel analysis code based on the SILFEED code for BWRs.

The SCPR fuel assembly has many square-shaped water rods. The fuel rods are arranged around these water rods. The fuel rod pitch and diameter are 11.2 mm and 10.2 mm, respectively. Since coolant flow distribution in the fuel assembly strongly depends on the gap width between the fuel rod and the water rod, the proper gap width is examined. Subchannel analysis shows that the coolant flow distribution becomes uniform when the gap width is 1.0 mm. The maximum fuel cladding surface temperature is lower than 600°C and the temperature margin of the fuel cladding is increased in the design.     

Heat transfer to liquid metals flowing turbulently and longitudinally through closely spaced rod bundles : Part I

This paper presents results of a theoretical study of heat transfer to liquid metals in fully developed turbulent, in-line flow through unbaffled, spacer-free rod bundles. The bundles have equilateral triangular arrangement; and the rod spacings, rod design, and ranges of independent variables covered were chosen with reference to liquid-metal-cooled nuclear reactor applications. Three different sets of thermal boundary conditions are considered: (A) uniform heat flux in the axial direction with uniform temperature in the circumferential direction, on the outer surface of the cladding; (B) uniform heat flux in both directions, on the outer surface of the cladding; and (C) uniform heat flux in both directions on the inner surface of the cladding. The results of the third set are presented in Part II.     

Effects of Coolant Flow on Light Water Reactor Fuel Behaviors during Reactivity Initiated Accident

This paper describes the in-pile experimental results to study the influences of coolant flow on fuel behaviors under reactivity initiated accident (RIA) conditions performed in the Nuclear Safety Research Reactor (NSRR). A single PWR type test fuel rod was irradiated by a large neutron pulse in the NSRR to simulate a prompt power excursion of RIA's. The effects of coolant flow were studied at a coolant flow velocity of 0.3~1.8m/s and a coolant temperature of 20~90°C under the atmospheric pressure. It was found that the cooling conditions had considerable influences on fuel thermal behaviors under prompt heat-up. The increase of coolant flow velocity and subcooling enhanced heat transfer coefficient at cladding surface during film boiling, which resulted in large decrease of maximum cladding temperature and film boiling duration, and consequently in the increase of fuel failure threshold energy. The data tendencies are summarized and the influences of coolant flow are discussed with some computer analyses.     

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.     

CARR堆芯热组件自然循环条件下特性分析

本文建立了中国先进研究堆标准燃料组件单组件的流-固耦合共轭传热CFD分析模型。通过1组稳态流量工况的分析,拟合获得燃料组件的阻力特性曲线。在堆本体CFD分析模型强迫流动工况计算结果的基础上,开展了标准燃料组件自然循环数值模拟分析。计算结果表明,在设定工况下,不仅释热能安全载出,而且可保证热组件任何位置均不会发生冷却剂泡核沸腾和流动不稳定性。计算得到了自然循环建立过程组件内冷却剂温度、燃料包壳和芯体的温度分布、热点位置以及循环流量的变化规律,为研究热组件的瞬态热工水力特性提供了理论方法和参考数据。     

A method of evaluating local turbulent diffusivity of heat in a gridded fuel subassembly

In evaluating the turbulent diffusivity of heat associated with the coolant flow past a grid spacer within an FBR fuel subassembly, a heat diffusion technique is usually employed. However, measurement of subchannel bulk coolant temperature using thermocouples usually involves difficulty due to a steep and non-linear temperature gradient in the subchannels adjacent to a heater pin.A series solution of the heat conduction equation for the coolant flow in subchannels past a grid spacer and a heated section of a dummy fuel pin was derived under a slug flow approximation where the boundary conditions on dummy fuel pins were satisfied by means of the point-matching technique. The solution may be utilized in analyzing the turbulent diffusivity of heat within subchannel coolant flow as a function of distance from a grid spacer based on the measured temperature distribution on the wall of dummy fuel pins, which may be obtained without affecting the subchannel coolant temperature.In an illustrative example, the turbulent diffusivity of heat was most exaggerated at about 50 mm beyond a grid spacer and was approximately five times larger than the corresponding diffusivity without a grid spacer.     

CFD analysis of flow field in a triangular rod bundle

The flow field was investigated in subchannels of VVER-440 pressurized water cooled reactors’ fuel assemblies (triangular lattice, P/D = 1.35). Impacts of the mesh resolution and turbulence model were studied in order to obtain guidelines for CFD calculations of VVER-440 rod bundles. Results were compared to measurement data published by Trupp and Azad in 1975. The study pointed out that RANS method with BSL Reynolds stress model using a sufficient fine grid can provide an accurate prediction for the turbulence quantities in this lattice. Applying the experiences of the sensitivity study thermal hydraulic processes were investigated in VVER-440 rod bundle sections. Based on the examinations the spacer grids have important effects on the cross flows, axial velocity and outlet temperature distribution of subchannels therefore they have to be modeled satisfactorily in CFD calculations.     

Analysis of radiant heat transfer in a BWR fuel assembly

A computer code ‘CIDER’ was developed which analyzes radiant heat transfer in a BWR fuel rod bundle under loss of coolant conditions. In the code, (1) a channel box and fuel rods are considered to be gray bodies, (2) reflection and absorption of radiation beams in the atmosphere is neglected, (3) a fuel rod is approximated by a regular polygonal rod, and (4) radiant heat flux is calculated considering circumferential temperature distribution on each fuel rod surface, which is determined from radial and circumferential heat conduction calculations in a fuel rod. It was found that the conventional model with uniform cladding temperature overestimated heat flux about 30% in a typical situation, or correspondingly underestimated the temperature rises.     

CFR600乏组件转换桶内氩气空间三维传热数值模拟

利用商业计算流体力学程序STAR?CCM+,对中国示范快堆（CFR600）乏组件转换桶内的氩气空间及相关结构部件进行数值模拟计算,验证分析正常工况和事故工况下的预埋件的热工设计。研究结果表明,正常工况下预埋件的温度峰值满足设计限值要求,但预埋件的温度分布存在一定程度的不均匀性,并且当发生事故工况时温度分布的差异性明显增大。钠液面通过辐射换热和对流换热共同对设备构件产生影响,其中辐射换热起主要作用。对于钠液面辐射换热的主要决定因素为钠液面的发射率。本工作为乏组件转换桶的正常运行及安全分析提供了技术准备。     

The contribution of non-wetting droplets to direct cooling of the fuel during PWR post-LOCA reflood

During the reflood of a Pressurised Water Reactor (PWR) following a loss of coolant accident, precursory cooling prior to the arrival of the rewetting front is of vital importance in limiting the rise in cladding temperature before rewet. This precursory cooling is achieved by a flow of superheated vapour, with entrained saturated drops, which evaporate into the vapour and act as a heat sink. In this paper we investigate a complementary mechanism; the direct cooling of the cladding by the drops themselves. Cladding temperatures are such that wetting by these droplets does not occur. On the contrary, droplets bounce off a vapour cushion formed during the ∼10 ms or so that they are in close proximity to the cladding. Using a combination of previous experimental correlations and recent CFD calculations, we estimate the rate of heat removal from the cladding surface as a result of the droplet impingement. Thus, we estimate the heat removed as a result of one impingement and estimate the total rate of heat removal by estimating the number droplets impinging on the cladding per unit surface area. The heat extracted by those droplets is found to be about 1/10 of the heat extracted by single-phase vapour under typical reflood conditions. Though there significant uncertainties in these estimates, it does seem that direct cooling by droplets, not generally incorporated in analyses of reflood, could actually be making a significant contribution to keeping cladding temperatures down to acceptable levels.     

HFETR多层环形燃料组件堵流事故研究

高通量工程试验堆（HFETR）的燃料组件采用了多层环形窄缝流道的设计来提高换热能力。然而,需要注意的是窄缝流道发生堵流的可能性较高。本文基于RELAP5程序建立了HFETR燃料组件模型,经过计算值与试验值的对比验证,结果表明该模型合理准确。基于该模型研究了堵流事故工况下热盒燃料组件的瞬态特性及其影响因素。结果表明:①当堵流比大于0.5时,随着堵流比的增加,燃料包壳与芯体峰值温度显著上升;②即使单个流道发生全部堵流,由于周围流道的冷却,燃料包壳峰值温度最大值只有218.6℃,能够保证燃料包壳的完整性;③单个流道全部堵流事故工况初期流量等参数波动较大,而在事故发生15 s后燃料组件主要参数基本稳定。