Analysis of MOX fuel behavior in reduced-moderation water reactor by fuel performance code FEMAXI-RM

To assess the feasibility of the 31% Pu-MOX fuel rod design of reduced-moderation water reactor (RMWR) in terms of thermal and mechanical behaviors, a single rod assumed to be irradiated in the core of RMWR up to 106 GWd/tHM has been analyzed by a fuel performance code FEMAXI-RM which is an extended version of FEMAXI-6 code. In the analysis, design specifications of fuel rod and irradiation conditions have been input, and available models of both MOX fuel and UO₂ fuel have been used as appropriate. The results are: fission gas release is several tens of percent, rod internal pressure does not exceed the coolant pressure, and the highest fuel center temperature is 2400 K, while cladding diameter increase caused by pellet swelling is within 1% strain. These predictions suggest that the MOX fuel rod integrity will be held during irradiation in RMWR, though actual behavior of MOX pellet swelling and cladding oxidation require to be investigated in detail.     

Analysis of mechanical load on cladding induced by fuel swelling during power ramp in high burn-up rod by fuel performance code FEMAXI-6

Mechanical load on cladding induced by fuel swelling in a high burn-up BWR type rod was analyzed by a fuel performance code FEMAXI-6. The code was developed for the analysis of LWR fuel rod behaviors in normal operation and transient conditions using finite element method (FEM).During a power ramp for the high burn-up rod, instantaneous pellet swelling can significantly exceed the level that is predicted by a “steady-rate” swelling model, causing a large circumferential strain in cladding. This phenomenon was simulated by a new swelling model to take into account the fission gas bubble growth. As a result it was found that the new model can give reasonable predictions on cladding diameter expansion in comparison with PIE data. The bubble growth model assumes that the equilibrium state equation holds for a bubble under external pressure, and simultaneous solution is obtained with both bubble size determination equation and diffusion equation of fission gas atoms. In addition, a pellet-clad bonding model which has been incorporated in the code to assume solid mechanical coupling between pellet outer surface and cladding inner surface predicted the generation of bi-axial stress state in the cladding during ramp.     

Elastic and viscoplastic pellets fragmentation modeling using an axisymmetrical 1D finite element code

During first rise to power in Power Water Reactor, fuel pellets crack because of thermal expansion. The phenomena of pellet cracking and fragments relocation have a major influence on rod behaviour and especially on the cladding behaviour in the case of pellet–cladding interaction.This article presents the modeling used to take into account the fragmented state of the pellet in the EDF fuel rod thermo-mechanical code, CYRANO3^®. The aim is to simulate more realistic stress and strain fields in the pellet.The investigated method consists in adding parameters in the 1D finite elements calculations in order to integrate the multi-dimensional fragmentation effects in the axisymmetrical 1D code CYRANO3^®. These parameters modify the material behaviour by describing the fuel as an anisotropic damaged material. The modeling accounts for the opening and closing of radial pellet cracks. It has been implemented in the code for elastic and viscoplastic fuel behaviours.     

Modification of thermal hydraulic code,THEATRe, with fuel rod behavior models and its application to reactivity-initiated accident (RIA)

During reactor operation, many complex changes occur in fuel rod which affects its thermal, mechanical and material properties. These changes also affect the reactor response to the transient and accident situations. Realistic simulation of fuel rod behavior under transients such as reactivity-initiated accident (RIA) is of great significance. In this study, thermal hydraulic analysis code THEATRe (Thermal Hydraulic Engineering Analysis Tool in Real-time) has been modified by addition of fuel rod behavior models for dynamic simulation of nuclear reactor. Transient changes in gas-gap parameters were taken into account by modeling the gas-gap behavior. Thermo-mechanical behavior of fuel rod is modeled to take into account the thermal, elastic and plastic deformation. To simulate RIA, point reactor kinetics model is also incorporated in the THEATRe code. To demonstrate the transient fuel rod behavior, AP1000 reactor is modeled and three hypothetical RIA cases are simulated. The RIA is considered at three different reactor power levels, i.e. 100, 50 and 1% of nominal power. The investigated parameters are fuel temperature, cladding stress and strain, fuel and cladding thermal conductivity and heat transfer coefficient in gas-gap. Modified code calculates the fuel rod temperatures according to updated fuel, clad and gas-gap parameters at the onset of steady-state operation and during the transient. The modified code provides lower steady-state fuel temperature as compared to the original code. Stress and strain analyses indicate that the hoop and radial strain is higher at high power locations of the fuel rod; therefore, gap closure process will initially occur in the central portion of the fuel rod and it should be given more emphasis in the safety analysis of the fuel rod and nuclear reactor during accidents and transients.     

Determination of stress within ZIRCALOY cladding due to pellet-cladding interaction

A precise calculation of the stress distribution within the Zircaloy cladding of a water-cooled reactor fuel rod subjected to a power increase is a complex problem which, in general, requires a computer code to integrate the behaviour of both the fuel and cladding. This paper develops a simplified model which decouples the clad and fuel pellet analyses, by considering two extremes of fuel pellet mechanical behaviour, which lead to two widely different boundary conditions at the pellet-clad interface. An axisymmetric fuel rod code can be used to give the mean cladding hoop strain imposed by the thermal expansion of the pellet, and when the interfacial friction coefficient is 0.5, this information along with the frictional boundary condition can be used to determine the stress distribution within the cladding near a fuel pellet crack. Results from this simplified approach, which does not involve an integrated code, are used to study the growth of stress corrosion cracks within the cladding.     

兆瓦级空间气冷堆系统启停堆瞬态分析

本文针对兆瓦级高温气冷堆布雷顿循环系统,采用Fortran语言开发系统分析程序TASS,包括堆芯、透平-发电机-压气机、回热器、冷却器和热管式辐射散热器等模型。通过设计值与程序计算值对比对TASS进行验证,并利用TASS对系统启动、停堆瞬态工况进行数值模拟。结果显示,通过分两阶段、阶梯式引入正反应性和提高涡轮机械的转轴速度,堆芯流量和功率匹配良好,系统可在3.5 h内完成启动过程,达到反应堆功率3 406 kW、流量14.2 kg/s的稳态运行。系统停堆过程中,反应堆可依靠自身的非能动余热排出能力,确保芯块和包壳温度与熔点间存在较大安全裕量,实现安全停堆。     

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.     

Thermal Stress Analysis of High-Burnup LWR Fuel Pellet Pulse-Irradiated in Reactivity-Initiated Accident Conditions

For RIA-simulated experiments in the NSRR with high-burnup PWR fuel and BWR fuel, numerical analyses were performed to evaluate the temporal changes of profiles of temperature and thermal stress in pellet induced by pulse power, using the RANNS code. The pre-pulse states of rods were calculated using the fuel performance code FEMAXI-6 along the irradiation histories in commercial reactors and the results were fed to the RANNS analysis as initial conditions of the rod. One-dimensional FEM was applied to the mechanical analysis of the fuel rod, and the calculated cladding permanent strain was compared with the measured value to confirm the validity of the PCMI calculation. The calculated changes in the profiles of temperature and stress in the pellet during an early transient phase were compared with the measured data such as the internal gas pressure rise, cracks and grain structure in the post-test pellet, anddiscussed in terms of PCMI and grain separation. The analyses indicate that the pellet cracking appearances coincided with the calculated tensile stress state and that the compressive thermal stress suppresses the fission gas bubble expansion leading to grain separation.     

Current perspective on fast reactor containment

The concept of “containment” is to provide a series of physical barriers between the radioactive products of the fission process and the public. All nuclear reactors have several such barriers and LMFBRs have more than most. These barriers are, successively:

1. fuel, which retains fission products;

2. fuel cladding, which encloses the fuel;

3. sodium coolant, which absorbs fission products released through fuel caldding;

4. primary coolant boundary, which has energy absorption and leakage control capabilities;

5. containment building, hereafter referred to as containment, which provides the final engineered barrier for control of radioactive releases;

6. exclusion distance, which provides space for natural attenuation of radioactive releases before reaching the public.

These barriers, along with the design approaches and features which protect their integrity under normal and accident conditions, assure that the public is adequately protected from the potential hazards of radioactivity residing in the core. It is only in the case of hypothesized core disruptive accidents (HCDAs) that these successive barriers can be sufficiently threatened as to pose a significant threat to the public. These HCDAs involve an extremely low probability sequence of successive failures resulting in core cooling imbalances which lead to fuel overheating. Under such conditions, the fuel and cladding barriers can be lost and energy sources can be generated which threaten the primary coolant boundary and containment. This paper addresses current perspectives on containment of HCDAs with emphasis on the approach and programs in the US.     

FRED fuel behaviour code: Main models and analysis of Halden IFA-503.2 tests

The FRED fuel rod code is being developed for thermal and mechanical simulation of fast breeder reactor (FBR) and light-water reactor (LWR) fuel behaviour under base-irradiation and accident conditions. The current version of the code calculates temperature distribution in fuel rods, stress-strain condition of cladding, fuel deformation, fuel-cladding gap conductance, and fuel rod inner pressure. The code was previously evaluated in the frame of two OECD mixed plutonium-uranium oxide (MOX) fuel performance benchmarks and then integrated into PSI's FAST code system to provide the fuel rod temperatures necessary for the neutron kinetics and thermal-hydraulic modules in transient calculations. This paper briefly overviews basic models and material property database of the FRED code used to assess the fuel behaviour under steady-state conditions. In addition, the code was used to simulate the IFA-503.2 tests, performed at the Halden reactor for two PWR and twelve VVER fuel samples under base-irradiation conditions. This paper presents the results of this simulation for two cases using a code-to-data comparison of fuel centreline temperatures, internal gas pressures, and fuel elongations. This comparison has demonstrated that the code adequately describes the important physical mechanisms of the uranium oxide (UOX) fuel rod thermal performance under steady-state conditions. Future activity should be concentrated on improving the model and extending the validation range, especially to the MOX fuel steady-state and transient behaviour.     

Analysis of the degradation of a PWR control rod during a loss of coolant accident, using the ICARE2 code

Conditions leading to AIC control rod damage during a loss of coolant accident in a PWR geometry, even in absence of violation of the LOCA licensing criteria, are investigated using several versions of the ICARE2 code (IPSN). Before being applied to the reactor case, the code and the modelling procedure are validated against the out-of-pile severe fuel damage experiment CORA-5. Three particular initial configurations are considered for the subsequent control rod damage analysis: nominal control rod and guide tube geometry, zircaloy guide tube bowing with concurrent cladding thickness reduction and finally control rod cladding perforation. For each of these cases the thermal, mechanical and chemical behaviour is presented. Phenomena such as ballooning and cladding failure of fuel rods, guide tube failure, melt relocation and final fluid channel cross-section modification are described. Finally, the conclusions of numerous sensitivity studies are discussed and some suggestions are given for possible improvements of the ICARE2 code.     

Stress biaxiality in high-burnup PWR fuel cladding under reactivity-initiated accident conditions

In order to promote a better understanding of failure mechanisms of high-burnup pressurized water reactor (PWR) fuels under reactivity-initiated accident (RIA) conditions, stress biaxiality in cladding has been estimated for the pellet-cladding (PC) mechanical interaction (PCMI) phase. The estimation was based on an analysis of the transient elongations of a pellet stack and a cladding tube measured in RIA-simulating experiments in the nuclear safety research reactor (NSRR) using the RANNS code. Stress biaxiality in the high-burnup PWR fuel cladding during the PCMI phase has been estimated to be 0.7–0.8, on average, at the mid-wall of the cladding. A comparison with fresh fuel test results and a sensitivity analysis showed that the effects of burnup and pulse width on cladding stress biaxiality are less than 10% for the investigated range. The present analysis also indicated that PC friction is strong, and that the cladding constraint on pellet stack elongation is significant irrespective of burnup. Therefore, it is recommended that strong PC friction be assumed, which is similar to the mechanical bonding condition, and that fuel pellets be treated as deformable materials in models of fuel behavior during the PCMI phase.     

燃料棒破损热力学模块开发及在严重事故分析程序中的应用

包壳肿胀和破损是严重事故早期阶段的重要现象。包壳形变不仅会造成局部流动堵塞,同时,水蒸气会从破裂处进入包壳气隙,增加包壳被蒸汽氧化的表面积。广泛使用的一体化严重事故分析程序不能分析早期事故过程中燃料棒的热力学行为,判断包壳破裂也只是基于简单的参数模型。本文开发了分析燃料棒热力学行为的FRTMB模块,并集成在严重事故分析程序ISAA中。使用开发的耦合系统ISAA FRTMB分析了CAP1400反应堆直接注射（DVI）管线小破口事故过程中燃料棒的热力学行为,并预计了包壳破裂时间及相应的失效温度。计算结果整体验证了ISAA FRTMB分析瞬态事故过程中燃料棒热力学行为以及判断包壳破裂的适用性和可靠性。     

Calculational procedure for determining creep collapse of LWR fuel rods

Observed collapses in pressurized water reactor fuel rods have been attributed to the radiation enhanced creep of Zircaloy cladding into regions where separations in the fuel pellet stack have occurred. A computer code, COLAPX, has been written to determine the growth of ovality and the ultimate collapse of fuel rod cladding under reactor operating conditions. This paper describes the theoretical bases of this code, the finite element formulation used, the constitutive relations between the displacement fields and the element forces, and the radiation, temperature and stress dependent material model for creep of Zircaloy tubing. Comparisons of the creep rate predictions and of the ovality predictions with data from irradiated tubes and fuel cladding are presented.     

Transient Elongation of a Fresh Fuel Rod under Reactivity Initiated Accident Conditions

In-pile experiments of fresh fuel rods under reactivity initiated accident (RIA) conditions have been performed in the Nuclear Safety Research Reactor at the Japan Atomic Energy Research Institute in order to understand the basic pellet cladding mechanical interaction (PCMI) behavior. Rapid fuel pellet expansion due to a power excursion would cause radial and longitudinal deformation of the cladding. This PCMI could be one of the possible incipient failure modes of an embrittled cladding of a high burnup fuel under the RIA conditions.

Basic PCMI behavior was studied by measuring cladding deformation of a fresh fuel rod without complicated irradiation effects. The transient elongation measurements of the fuel with two kinds of gap width indicated not only PCMI-induced cladding elongation, but also reduction of the pellet stack displacement by the cladding constraint. In the tests under a high-pressure and high-temperature condition simulating an operation condition of BWRs, additional ridge-type cladding deformation was generated due to the axial collapse of the cladding. A preliminary analysis for interpretation of the tests was made using a computer code for the transient analysis of fuel rods, FRAP-T6.     

小型铅-铋冷却快堆提棒事故核热耦合研究

基于热工程序COBRA-YT和物理程序SKRTCH-N,利用并行虚拟机（PVM）平台开发了核热耦合工具:COBRA-YT将冷却剂密度和燃料温度等热工参数传递给物理程序,用以更新截面;SKETCH-N执行物理计算,并将功率分布反馈给热工程序;最后,应用该耦合程序分析铅-铋冷却快堆的提棒事故。计算结果显示控制棒提起后,功率迅速升高,在1.42?s后达到最大值;5?s后包壳温度达到峰值1264℃,超出了设计限值。结果表明:在提棒事故后,均一化布置堆芯的安全会在极短时间内受到严重威胁,故该堆芯应采用分区布置。      

Oxide thickness-dependent transient cladding hoop stress

A few thrice-burned Zry-4 fuel assemblies which were loaded in one of the PWRs operating in Korea were found to be failed due to PCI during a power ramp following a rector trip, while thrice-burned Zr–Nb fuel assemblies and twice-burned Zry-4 ones were intact. To investigate the effect of fuel rod oxide thickness on power ramp-induced cladding hoop stress, three intact fuel rods were selected, which include an intact twice-burned Zry-4 fuel rod, an intact thrice-burned Zr-4 fuel rod and an intact thrice-burned Zr–Nb fuel rod. With the use of a fuel performance analysis code, burnup-dependent steady-state cladding stress and ramp power-dependent cladding stresses at the power-ramped burnup were predicted for the three intact fuel rods. It was found that the cladding oxide thickness has a considerable effect on the ramp power-dependent cladding hoop stresses. In addition, the cladding maximum stress of the thrice-burned Zry-4 fuel rod with 125 μm oxide thickness exceeded an ultimate cladding tensile strength of the Zry-4 cladding when the pellet–clad friction coefficient-dependent cladding stress concentration ratio was considered. However, the thrice-burned Zr–Nb fuel rod with 50 μm oxide thickness was evaluated to have a considerable margin against the power ramp-induced PCI failure.     

Comparative effects of structural and material parameters variability on Pellet–Cladding Interaction in a PWR fuel rod

In the water pressurised nuclear reactors, the fuel rod cladding is the first barrier against radioactive isotopes release. Its integrity must be demonstrated all along the fuel rod irradiation with increasing strenuous operating conditions. This paper deals with a study made with Electricité de France (EDF) in order to improve the understanding and modelling of the thermomechanical behaviour of fuel rods under these more arduous conditions. The aim of this study is to evaluate the separate influences of structural and material parameters variability on Pellet–Cladding Interaction (PCI). The following parameters have been tested by 3D simulations for a common ramping condition: axial and radial pellet cracks numbers, pellet fragment size, relative fragments displacement, non-symmetrical configuration, pellet–pellet friction and pellet–cladding friction. The second part of the article deals with the development of a model, which gives a better assessment of cladding stress concentration near radial fuel cracks. The implementation in the 1D fuel rod EDF code CYRANO3 and the validation of the model are presented.     

基于流固热耦合的燃料元件包壳结构完整性分析

反应堆系统发生瞬态工况时,冷却剂温度的瞬间大幅度变化会对燃料元件包壳结构完整性造成冲击,危及反应堆安全。本文以某压水堆3×3燃料组件为对象,采用流固热耦合方法对冷水事故下燃料组件的流动换热特性和燃料元件包壳温度、变形及应力进行了三维精细化模拟。结果表明:定位格架能够增强燃料棒表面的对流换热强度;包壳变形时向与刚凸接触的一侧折弯,向与弹簧接触的一侧凸起;包壳与定位格架接触部位的温度和最大等效应力随事故时间不断增大,且最大等效应力超过了包壳材料的屈服强度,将发生强度失效,影响其结构完整性。本文研究可为反应堆燃料元件包壳瞬态工况下的完整性评价提供借鉴。