WTRLGD—A computer program for the transient analysis of waterlogged fuel rods under the RIA condition

A computer code WTRLGD has been developed to describe the transient internal pressure of a waterlogged fuel rod during power burst and also to predict the possibility of the rod failure in the mode of cladding rupture. The code predicts transient thermal behavior of the fuel rod on the basis of an assumption of axisymmetry, and thermal-hydraulic transients of the internal water on the basis of a homogeneous volume-junction model modified so as to involve the cladding deformation. Calculated transients of the rod pressure are in fairly good agreement with those measured in the NSRR experiments, simulating the fuel rod behavior under an RIA condition. The comparison between calculation and experiment verifies that the code is an effective tool for the prediction of the failure of a waterlogged fuel rod.     

基于MOOSE框架的五方程两相流分析程序开发

基于多物理场耦合框架MOOSE,采用五方程两相流模型开发了模块化程序ZEBRA,实现了高阶时间、空间离散格式两相流动传热问题的求解。采用Bartolomei开展的垂直圆管过冷沸腾实验对ZEBRA进行验证,在不同热流密度、质量流密度、压力工况下,将程序计算值与实验值进行了数值验证和计算分析。结果表明：ZEBRA中五方程模型预测值与实验值符合良好,沸腾起始点和空泡份额的预测合理,表明ZEBRA初步具备了处理两相流问题的能力。     

FARST — A computer code for the evaluation of FBR fuel rod behavior under steady-state/transient conditions

FARST, a computer code for the evaluation of fuel rod thermal and mechanical behavior under steady-state/transient conditions has been developed. The code characteristics are summarized as follows:

1. (i) FARST evaluates the fuel rod behavior under the transient conditions. The code analyzes thermal and mechanical phenomena within a fuel rod, taking into account the temperature change in coolant surrounding the fuel rod.

2. (ii) Permanent strains such as plastic, creep and swelling strains as well as thermoelastic deformations can be analyzed by using the strain increment method.

3. (iii) Axial force and contact pressure which act on the fuel stack and cladding are analyzed based on the stick/slip conditions.

4. (iv) FARST used a pellet swelling model which depends on the contact pressure between pellet and cladding, and an empirical pellet relocation model, designated as “jump relocation model”.

The code was successfully applied to analyses of the fuel rod irradiation data from pulse reactor for nuclear safety research in Cadarache (CABRI) and pulse reactor for nuclear safety research in Japan Atomic Energy Research Institute (NSRR).The code was further applied to stress analysis of a 1000 MW class large FBR plant fuel rod during transient conditions. The steady-state model which was used so far gave the conservative results for cladding stress during overpower transient, but underestimated the results for cladding stress during a rapid temperature decrease of coolant.     

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

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

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.     

Modeling of Jet Pump by RELAP4 Code in Loss-of-Coolant Accident Analysis

Waterlogged fuel rod experiments performed at the NSRR are analyzed using the computer code WTRLGD, which was devised for the analyses of thermo-dynamical behavior of a waterlogged fuel rod. The numerical results are compared with the data from the experiments in order to assess the validity of the computer code. Parameters in the analyses are volumetric fraction of water, reactor period, gap width, a pin hole and the end peaks. Thus the analyses cover almost all the waterlogged fuel rod experiments at the NSRR.

The comparison shows good agreement between the experimental results and numerical ones on the transient thermo-dynamical behaviors of fuel, such as, rod internal pressure, cladding surface temperature and cladding strain. The numerical results also quantitatively agree with the experimental data concerning the effects of the above parameters on failure threshold energy. From the above findings, the computer code is assessed to be valid enough for the analyses of the failure behavior of the waterlogged fuel rod under a reactivity initiated accident condition.     

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.     

三维耦合分析软件的落棒事故分析

作为数值反应堆中必不可少的物理和热工部分,中广核研究院有限公司开发了三维物理热工耦合分析软件,通过动态链接库技术实现了自主研发的核反应堆系统瞬态分析软件和三维核设计软件的耦合,并已与国际基准题结果对比验证。本文为耦合软件的应用,围绕华龙一号的落棒分析问题,开展不同落棒组合的耦合计算分析,并研究停堆棒组落棒和温度调节棒（R）棒组两组落棒对堆芯功率的影响。分析结果表明,非中心对称的棒组落棒事故会导致堆芯径向功率出现不对称,并使得堆芯出口回路温度不同。落棒反应性价值越大,R棒调节后的稳态功率回升相比初始稳态差异越大,DNBR公式计算值的变化趋势与功率呈现相反规律。     

压水堆燃料棒三维热-力学性能的精细模拟

压水堆燃料棒工作在复杂的辐照、热和力学环境中,对其性能进行定量评估涉及多种复杂的物理现象。目前常用的燃料性能分析程序一般对结构采用简化的轴对称假设,对辐照肿胀、辐照蠕变和高温蠕变等物理现象以及辐照-热-力等物理场之间的耦合考虑并不充分。基于ABAQUS有限元求解框架,开发了压水堆燃料棒三维热-力学性能的模拟程序,利用程序对压水堆燃料棒进行了稳态分析,以及升功率和反应性引入事故两种瞬态分析。结果表明:辐照引起燃料致密化和肿胀对燃料温度变化有重要影响;芯块应变增加主要是由裂变产物肿胀引起的;芯块几何结构导致包壳应力集中发生在芯块间的交界面处;燃料棒功率的急剧变化会加快芯块表面破裂的进程;反应性引入事故会导致芯块从内部开始破裂,并会引发芯块-包壳的接触。     

Effects of Rod Pre-Pressurization on Light Water Reactor Fuel Behavior during Reactivity Initiated Accident Conditions

Previously pressurized (pre-pressurized) fuel rod tests recently performed in the Nuclear Safety Research Reactor (NSRR) investigate the effects of initial internal pressure on fuel rod behavior during reactivity initiated accident (RIA) conditions. A single PWR type fuel rod was contained within a waterfilled, ambient temperature and ambient pressure capsule. The fuel rod was then heated by the pulsing operation of the NSRR.

Results from the tests show that the effect of pre-pressurization was significant for the fuel rods with initial internal pressure of 0.8 MPa and above, and fuel rod failure occurred from rupture of the cladding with lower threshold energy deposition for failure as the initial internal pressure was increased. The cladding rupture was governed mainly by the cladding temperature rise, not by the rod internal pressure rise during the transient. The relationships between cladding burst pressure and cladding burst temperature and between cladding strain and cladding temperature at cladding rupture obtained in the present study under an RIA condition agree with the results obtained from various in- and ex-reactor experiments under a LOCA condition, although the obtained time-averaged strain rate of the Zircaloy cladding was much greater than that in a LOCA condition.     

Thermal hydraulic and safety analysis for core conversion (HEU–LEU) of Syrian Miniature Neutron Source Reactor

The paper presents the behavior and properties analysis of the low enriched uranium fuel compared with the original high enriched uranium fuel. The MNSR reactor core was modeled with both fuel materials and the reactor behavior was studied during the steady state and abnormal conditions. The MERSAT code was used in the analysis. The steady state thermal hydraulic analysis results were compared with that obtained from the experimental results hold during commissioning the Syrian MNSR. Comparison with experimental data shows that the steady-state behavior of the HEU core was accurately predicted by the MERSAT code calculations. The validated model was then used to analyze LEU cores with two proposed UO₂ fuel pin designs. With each LEU core, the steady state and 3.77 mk rod withdrawal transient were run and the results were compared with the available published data in the literatures for the low enriched uranium fuel core. The results reveal that the low enriched uranium fuel showed a good behavior and the peak clad temperatures remain well below the clad melting temperature during reactivity insertion accident.     

Analysis on split failure of cladding of high burnup BWR rods in reactivity-initiated accident conditions by RANNS code

A computer code RANNS was developed to analyze fuel rod behaviors in the reactivity-initiated accident (RIA) conditions. RANNS performs thermal and finite-element mechanical calculation for a single rod in axis-symmetric geometry, where fuel pellet consists of 36 equal-volume ring elements and cladding metallic wall consists of eight equal-thickness ring elements and one outer oxide element. The code can calculate temperature profile inside the rod, contact pressure generated by pellet–clad mechanical interaction (PCMI), stress–strain distribution and their interactions elaborately. An experimental analysis by RANNS begins with pre-test conditions of irradiated rod which are given by the fuel performance code FEMAXI-6.In the present study, analysis was performed on the simulated RIA experiments in the “nuclear safety research reactor” (NSRR), FK-10 and FK-12, with high burnup BWR rods in a cold-start up condition, and stress–strain evolution in the PCMI process was calculated extensively. In the analysis, the pellet–clad bonding was assumed both in the heat conduction and in mechanical restraint. The calculated hoop strain increase was compared with the measured strain gauge data, and satisfactory agreement was obtained. Simulation calculations with broader power pulses anticipated in RIA of commercial BWR were carried out and the resulted cladding hoop stress was compared with the failure stress estimated by comparison of analysis with experimental data.     

Developing and modeling of the “Laguna Verde” BWR CRDA benchmark

Reactivity initiated accidents (RIA) and design basis transients are one of the most important aspects related to nuclear power reactor safety. These events are re-evaluated whenever core alterations (modifications) are made as part of the nuclear safety analysis performed to a new design. These modifications usually include, but are not limited to, power upgrades, longer cycles, new fuel assembly and control rod designs, etc. The results obtained are compared with pre-established bounding analysis values to see if the new core design fulfills the requirements of safety constraints imposed on the design. The control rod drop accident (CRDA) is the design basis transient for the reactivity events of BWR technology. The CRDA is a very localized event depending on the control rod insertion position and the fuel assemblies surrounding the control rod falling from the core. A numerical benchmark was developed based on the CRDA RIA design basis accident to further asses the performance of coupled 3D neutron kinetics/thermal-hydraulics codes. The CRDA in a BWR is a mostly neutronic driven event. This benchmark is based on a real operating nuclear power plant — unit 1 of the Laguna Verde (LV1) nuclear power plant (NPP). The definition of the benchmark is presented briefly together with the benchmark specifications. Some of the cross-sections were modified in order to make the maximum control rod worth greater than one dollar. The transient is initiated at steady-state by dropping the control rod with maximum worth at full speed. The “Laguna Verde” (LV1) BWR CRDA transient benchmark is calculated using two coupled codes: TRAC-BF1/NEM and TRAC-BF1/ENTRÉE. Neutron kinetics and thermal hydraulics models were developed for both codes. Comparison of the obtained results is presented along with some discussion of the sensitivity of results to some modeling assumptions.     

RANNS Code Analysis on the Local Mechanical Conditions of Cladding of High Burnup Fuel Rods under PCMI in RIA-Simulated Experiments in NSRR

The RANNS code analyzes behavior of a single fuel rod in reactivity-initiated accident (RIA) conditions. The code has two types of mechanical model; one-dimensional deformation model for each axial segment length of rod, and newly-developed two-dimensional local deformation model for one pellet length. Analyses were performed on the RIA-simulated experiments in the Nuclear Safety Research Reactor (NSRR), OI-10 with high burnup PWR rods, and results of cladding deformation were compared between calculations by the two models and PIE data. The pre-accident, or End-of-Life conditions of the rod were predicted by the fuel performance code FEAMXI-6. In the calculations by the two-dimensional model of RANNS, the plastic strain increases at the cladding ridges during PCMI were compared with those in between the ridges and with the PIE data, and effect of stress variance induced by local non-uniformity of strain on the crack growth was discussed.     

基于RELAP5程序的AP1000典型事故瞬态特性研究

基于最佳估算程序RELAP5/MOD3.3,对AP1000系统进行了详细的建模分析,选取冷却剂泵卡轴事故、蒸汽发生器(SG)传热管破裂事故和直接注射管线双端断裂事故作为典型事故,获得了典型事故工况下关键参数的瞬态特性和非能动系统响应特性。结果表明:对于冷却剂泵卡轴事故,一回路最大压力为16.82 MPa,燃料包壳表面温度最大值为1 299K,满足验收准则的要求;对于SG传热管破裂事故,破损SG的水体积为231.54m3,小于AP1000蒸汽发生器255.563m3的总容积;对于直接注射管线双端断裂事故,AP1000的非能动堆芯冷却系统能对一回路进行冷却和降压,并防止堆芯裸露和燃料包壳过热。     

Simulation of the chernobyl accident

An analysis of the April 26, 1986 accident at the Chernobyl-4 nuclear power plant in the Soviet Union is presented. The peak calculated core power during the accident was 550 000 MW_t. The analysis provides insights that further understanding of the plant behavior during the accident. The plant was modeled with the RELAP5/MOD2 computer code using information available in the open literature. RELAP5/MOD2 is an advanced computer code designed for best-estimate thermal-hydraulic analysis of transients in light water reactors. The Chernobyl-4 model included the reactor kinetics effects of fuel temperature, graphite temperature, core average void fraction, and automatic regulator control rod position. Preliminary calculations indicated the effects of recirculation pump coast down during performance of a test at the plant were not sufficient to initiate a reactor kinetics-driven power excursion. Another mechanism, or “trigger” is required. The accident simulation assumed the trigger was recirculation pump performance degradation caused by the onset of pump cavitation. Fuel disintegration caused by the power excursion probably led to rupture of pressure tubes. To further characterize the response of the Chernobyl-4 plant during severe accidents, simulations of an extended station blackout sequence with failure of all feedwater are also presented. For those simulations, RELAP5/MOD2 and SCDAP/MOD1 (an advanced best-estimate computer code for the prediction of reactor core behavior during a severe accident) were used. The simulations indicated that fuel rod melting was delayed significantly because the graphite acted as a heat sink.