Parameter studies of peak temperature rise for LMFBR blanket rod bundles

The peak temperature rise results obtained from sodium tests of a full sized, electrically heated simulated blanket rod bundle were related to operating parameters. They were found to be related to geometry, heat input gradient, coolant flow rate at constant power-to-flow ratio, and mode of circulation (forced or natural). At any fixed power-to-flow ratio, the peak temperature rise was shown to reach a maximum value at fully turbulent flow and forced convection.     

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

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

冷却剂与含格架燃料棒束多场耦合分析

燃料棒束作为压水堆燃料组件的组成部分,其热工和结构特性直接关系到反应堆的安全。本文利用ANSYS WORKBENCH软件分析了冷却剂在5×5含定位格架燃料棒束通道内流动的分布,采用冷却剂与燃料棒束多场耦合的方式研究了燃料棒束的流动传热特性和结构形变特性。结果表明：定位格架扰动冷却剂形成横向二次流并在下游棒束间形成绕流;多场耦合条件下二次流峰值速度和平均速度均小于单流场的;二次流与燃料棒的热应力使棒束发生形变,功率和流动分布的不均匀导致形变在轴向和径向的不均匀;相较于无格架情况,定位格架的存在使冷却剂的搅混流动更加明显,冷却剂对燃料棒冲击增大;在有、无定位格架两种情况下棒束形变均很小,可保持原本结构的稳定。     

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.     

基于CFD方法的铅铋冷却燃料棒束的热工水力特性分析

铅铋冷却快堆作为第4代反应堆候选之一具有安全性高等特点,研究其在正常工况下的热工水力特性具有重要意义。本文基于商用计算流体力学（CFD）软件STAR-CCM+,使用流固耦合的方法对带有绕丝结构的19棒束铅铋组件进行数值分析,探究了质量流量、功率等边界条件对组件内部流动传热特性的影响。模拟计算结果表明：CFD方法在子通道中心温度和壁面温度预测上与实验结果取得了较好的一致。同时,绕丝结构的存在使得子通道之间存在周期性的横向交混,并使得棒束表面温度呈现震荡。随质量流量的增加,子通道间横向交混增大。功率变化对通道间的横向交混速度的影响较小,冷却剂温度的横向分布无明显差异。     

Temperature distribution in mixed ThO2–UO2 fuel rods located in blanket of an inertial fusion energy breeder

Investigations of neutronic analysis and temperature distribution in fuel rods located in a blanket driven ICF (Inertial Confinement Fusion) have been performed for various mixed fuels and coolants under a first wall load of 5 MW/m². The fuel rods containing ThO₂ and UO₂ mixed by various mixing methods for achieving a flat fission power density are replaced in the blanket and cooled with different coolants; natural lithium, flibe, eutectic lithium and helium for the nuclear heat transfer. It is assumed that surface temperature of the fuel rod increases linearly from 500 °C (at top) to 700 °C (at bottom) during cooling fuel zone. Neutronic and temperature distribution calculations have been performed by MCNP4B Code and HEATING7, respectively. In the blanket fueled with pure UO₂ and cooled with helium, M (fusion energy multiplication ratio) increases to 3.9 due to uranium having higher fission cross-section than thorium. The high fission energy released in this blanket, therefore, causes proportionally increasing of temperature in the fuel rods to 823 °C. However, the M is 2.00 in the blanket fueled with pure ThO₂ and cooled with eutectic lithium because of more capture reaction than fission reaction. Maximum and minumum values of TBR (tritium breeding ratio) being one of main neutronic paremeters for a fusion reactor are 1.07 and 1.45 in the helium and the natural lithium coolant blanket, respectively. These consequences bring out that the investigated reactor can produce substantial electricity in situ during breeding fissile fuel and can be self-sufficient in the tritium required for the DT fusion driver in all cases of mixed fuels and coolant types. Quasi-constant fission power density profiles in FFB (fissile fuel breeding) zone are obtained by parabolically increasing mixture fraction of UO₂ in radial and axial directions for all coolant types. Such as, in the helium coolant blanket and the case of PMF (parabolically mixed fuel), Γ (peek-to-average fission power density ratio) of the blanket is reduced to 1.1, and the maximum temperatures of the fuel rods in radial direction of the FFB zone are also quasi-constant. At the same time, in the case of PMF, for all coolant types, the temperature profiles in the radial direction of the fuel rods rise proportionally with surface temperature from the top to the bottom of fuel rods in the axial direction. In other words, for each radial temperature profile in the axial direction, temperature differences between centerline and surface of the fuel rods are quasi-constant. According to the coolant types, these temperature diffences vary between 30 and 45 °C.     

Effects of wire-spacer shape in LMR on thermal–hydraulic performance

In this work, analyses of three-dimensional flow and convective heat transfer in wire-wrapped fuel assemblies with different shaped wire-spacers have been carried out using Reynolds-averaged Navier–Stokes equations with shear stress transport turbulence model. Three cross-sectional shapes of wire-spacer, circle, hexagon and rhombus have been tested. All the assemblies have been analyzed for single pitch of wire-spacer with periodic boundary conditions applied at inlet and outlet of the calculation domain. It is found that the assemblies exhibit the directional periodicity in radial gradients between the adjacent subchannels due to presence of wire-spacer. The overall pressure drop is highest in case of rhombus shaped wire-spacer assembly followed by hexagonal shaped. Although circular shaped wire-spacer gives lowest peak temperature as well as lowest overall temperature difference in the assembly, the rhombus shaped wire-spacer assembly gives highest Nusselt number on fuel rod surface.     

Single-phase convective heat transfer in rod bundles

The convective heat transfer for turbulent flow through rod bundles representative of nuclear fuel rods used in pressurized water reactors is examined. The rod bundles consist of a square array of parallel rods that are held on a constant pitch by support grids spaced axially along the rod bundle. Split-vane pair support grids, which create swirling flow in the rod bundle, as well as disc and standard support grids are investigated. Single-phase convective heat transfer coefficients are measured for flow downstream of support grids in a rod bundle. The rods are heated using direct resistance heating, and a bulk axial flow of air is used to cool the rods in the rod bundle. Air is used as the working fluid instead of water to reduce the power required to heat the rod bundle. Results indicate heat transfer enhancement for up to 10 hydraulic diameters downstream of the support grids. A general correlation is developed to predict the heat transfer development downstream of support grids. In addition, circumferential variations in heat transfer coefficients result in hot streaks that develop on the rods downstream of split-vane pair support grids.     

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.     

Cladding response to uniform radial compaction of clinch river breeder reactor blanket rod bundles

Neutron radiation induced growth can lead to compaction of fast reactor fuel and blanket rod bundles within their duct assemblies, resulting in local cladding stresses and coolant channel closure. The stiffness of the Clinch River Breeder Reactor (CRBRP) blanket rod bundle during uniform radial compaction was measured experimentally, and was compared to theoretical predictions which assumed no relative rod-rod motion within the bundle during compaction. The basic agreement between the measured and theoretical compaction stiffness values permitted an assessment to be made of cladding stresses and coolant channel closures in CRBRP blanket rod bundles due to rod bundle-duct differential growth.     

Heat transfer from a partially uncovered reactor core

This paper documents a model which has been developed for predicting the temperature distribution along a “flow channel” of a pressurized water reactor during simulated, uncovered core conditions. In the model, heat conduction along the fuel element, convection from the surface to the coolant, radiation exchange between the clad surface and steam, and surface exchange between adjacent fuel rods are considered. Variations of the thermophysical properties of the fuel road and of the coolant with temperature are accounted for, but oxidation of Zircaloy is not modeled. Extensive sensitivity studies on the effects of heat generation in the core, steam velocity, pressure level, uncovered core height, presence of hydrogen gas in the coolant, power skew, clad emissivity, and convective heat transfer correlations have been examined. The results show that the importance of radiation in comparison with convection increases with an increase in the fuel rod temperature, pressure, and clad emissivity.     

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.     

Dynamic Model Representing Condensate Feedwater System of BWR Plants

In order to study the effect of burst temperature on the coolant flow channel restriction, burst tests of fuel bundles were performed. Each bundle consisted of 49 rods (7×7 rods), and bursts were conducted in flowing steam. Burst temperature was changed by changing the internal gas pressure in rods. After the burst, the ballooning behavior of each rod and the degree of coolant flow area restriction in the bundle were measured.

Maximum swelling of rod occurs when the burst temperature is around α and α+β phase boundary, and this phenomenon is almost the same as that in single rod burst tests. Maximum coolant flow area restriction is also observed in this condition.     

Heat transfer in turbulent longitudinal flow through unbaffled assemblies of fuel rods

A numerical analysis of heat transfer in turbulent longitudinal flow through assemblies of unbaffled fuel rods is presented. The solution applies to triangular or rectangular arrays of fuel rods with fully developed velocity and temperature profiles, for fluids with Prandtl number 1 and « 1. In the case of liquid metals, the thermal resistance of the cladding and bond are considered, but the turbulent heat transport component is neglected. For common liquids the circumferential turbulent heat transfer is considered. Results are compared in the range of dimensionless rod spacing of 1.0–1.6. Theoretical predictions and experimental results of other authors dealing with the problem show relatively good agreement.     

Accuracy of nuclear design of fast and thermal neutron coupled core by SRAC

Supercritical-pressure light water cooled fast reactor adopts the blanket fuel assemblies with depleted uranium fuel and zirconium hydride layer in the core for negative coolant void reactivity. Thermal neutrons are generated in the core of fast reactor. It is called “fast and thermal neutron coupled core”. The purpose of the present study is to examine the accuracy of assembly and core calculations including preparation of the macroscopic cross sections with the SRAC code system for “fast and thermal neutron coupled core” in comparison with the Monte Carlo codes, MVP and MVP-BURN. Accuracy of the neutron multiplication factor and coolant void reactivity calculation has been evaluated in four types of cores of different fractions of the blanket fuel assembly with zirconium hydride rods. The conventional analysis is based on the macroscopic cross sections which are prepared with infinite lattice. The conventional SRAC calculation underestimates the neuron multiplication factor for all types of cores. Other findings are that the conventional SRAC calculation overestimates coolant void reactivity for the cores without zirconium hydride rods, and underestimates coolant void reactivity for the core of all blanket fuel assemblies having zirconium hydride rods. To overcome these problems, it has been proposed that the macroscopic cross sections of seed fuel assembly are prepared with the model that a seed fuel assembly is surrounded by blanket fuel assemblies in order to take into account the effects of the surrounding fuel assemblies. Evaluations show that accuracy of the neutron multiplication factor by the SRAC calculation can be improved by the proposed method.     

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.     

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.     

加/减速流动下棒束通道内速度分布和湍流特性研究

事故工况及海洋条件下反应堆处于非稳态工况,堆芯燃料组件内热工水力行为复杂多变,对反应堆安全提出了更高挑战,因此有必要对非稳态下燃料组件内流动换热特性开展研究。基于粒子图像测速（PIV）技术,结合远心镜头和脉冲控制器,实现对燃料组件内复杂流场的高时空分辨率、长时间的连续测量,获得了流量波动下燃料组件内时空演变的流场结构,分析了棒束通道内速度分布、湍流强度、雷诺应力等瞬时流场信息的空间演变特性。以定常流动下流场分布特性为基准,对比分析了加速度对燃料组件内空间流场分布的贡献特点。实验结果表明：加速流动提高了棒束通道内流层之间的速度梯度,抑制了横向速度和湍流强度;减速流动减弱了棒束通道内流层之间的速度梯度,提高了横向速度和湍流强度。实验结果有助于揭示燃料组件在非稳态条件下的瞬态特性,并为燃料组件的设计和优化奠定基础。     

Analytical model for calculation of the thermo hydraulic parameters in a fuel rod assembly

The paper presents the procedure of the cellular calculation of thermo hydraulic parameters of a single-phase gas flow in a fuel rod assembly. The procedure is implemented in the DARS program. The program is intended for calculation of the distribution of the gaseous coolant parameters and wall temperatures in case of arbitrary, geometrically specified, arrangement of the rods in fuel assembly and in case of arbitrary, functionally specified in space, heat release in the rods.In mathematical model the flow cross-section of the channel of intricate shape is conventionally divided to elementary cells formed by straight lines, which connect the centers of rods. Within the limits of a single cell the coolant parameters and the temperature of the corresponding part of the rod surface are assumed constant. The entire fuel assembly is viewed as a system of parallel interconnected channels.Program DARS is illustrated by calculation of a temperature mode of 85-rod assembly with spacers of wire wrapping on the rods.