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.     

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.     

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.     

Evaluation of high plutonia (44% PuO2) MOX as a fuel for fast breeder test reactor

Uranium plutonium mixed oxide (MOX) containing up to 30% plutonia is the conventional fuel for liquid metal cooled fast breeder reactor (LMFBR). Use of high plutonia (>30%) MOX fuel in LMFBR had been of interest but not pursued. Of late, it has regained importance for faster disposition of plutonium and also for making compact fast reactors. Some of the issues of high plutonia MOX fuels which are of concern are its chemical compatibility with liquid sodium coolant, dimensional stability and low thermal conductivity. Available literature information for MOX fuel is limited to a plutonium content of 30%. Thermodynamic assessment of mixed oxide fuels indicate that with increasing plutonia oxygen potential of the fuel increases and the fuel become more prone to chemical attack by liquid sodium coolant in case of a clad breach. In the present investigation, some of these issues of MOX fuel have been studied to evaluate this fuel for its use in fast reactor. Extensive work on the out-of-pile thermo-physical properties and fuel-coolant chemical compatibility under different simulated reactor conditions has been carried out. Results of these studies were compared with the available literature information on low plutonia MOX fuel and critically analyzed to predict in reactor behaviour of this fuel containing 44% PuO₂. The results of these out-of-pile studies have been very encouraging and helped in arriving at a suitable and achievable fuel specification for utilization of this fuel in fast breeder test reactor (FBTR). As a first step of test pin irradiation programme in FBTR, eight subassemblies of the MOX fuel are undergoing irradiation in FBTR.     

Pu recycling in a full Th-MOX PWR core. Part I: Steady state analysis

Current practice of Pu recycling in existing Light Water Reactors (LWRs) in the form of U-Pu mixed oxide fuel (MOX) is not efficient due to continuous Pu production from U-238. The use of Th-Pu mixed oxide (TOX) fuel will considerably improve Pu consumption rates because virtually no new Pu is generated from thorium. In this study, the feasibility of Pu recycling in a typical pressurized water reactor (PWR) fully loaded with TOX fuel is investigated.Detailed 3-dimensional 100% TOX and 100% MOX PWR core designs are developed. The full MOX core is considered for comparison purposes. The design stages included determination of Pu loading required to achieve 18-month fuel cycle assuming three-batch fuel management scheme, selection of poison materials, development of the core loading pattern, optimization of burnable poison loadings, evaluation of critical boron concentration requirements, estimation of reactivity coefficients, core kinetic parameters, and shutdown margin.The performance of the MOX and TOX cores under steady-state condition and during selected reactivity initiated accidents (RIAs) is compared with that of the actual uranium oxide (UOX) PWR core.Part I of this paper describes the full TOX and MOX PWR core designs and reports the results of steady state analysis. The TOX core requires a slightly higher initial Pu loading than the MOX core to achieve the target fuel cycle length. However, the TOX core exhibits superior Pu incineration capabilities.The significantly degraded worth of control materials in Pu cores is partially addressed by the use of enriched soluble boron and B₄C as a control rod absorbing material. Wet annular burnable absorber (WABA) rods are used to flatten radial power distribution. The temperature reactivity coefficients of the TOX core were found to be always negative. The TOX core has a slightly reduced, as compared to UOX core, but still sufficient shutdown margin.In the TOX core β_eff is smaller by about a factor of two in comparison to the UOX core and even lower than that of the MOX core. The combination of small β_eff and reduced control materials worth may potentially deteriorate the performance under RIA conditions and requires an additional examination. The behavior of the considered cores during the most limiting RIAs, such as rod ejection, main steam line break, and boron dilution, is further investigated and reported in Part II of the paper.     

铅基快堆氧化物燃料性能分析程序FUTURE开发与初步验证

为了对铅基快堆氧化物燃料元件稳态工况下的服役性能和行为演化进行模拟计算,本文基于串行的半隐式耦合求解方法开发了铅基快堆氧化物燃料性能分析程序FUTURE。程序采用两步分析法实现了铅基快堆氧化物燃料棒全域热力分析与局部行为模型的多物理场耦合计算。通过各计算模块与模型算例、基准公式和现有程序的对比分析,对FUTURE程序进行了各分离效应的初步验证。结果表明,FUTURE程序能准确模拟铅基快堆稳态工况条件下氧化物燃料元件内部的温度演化、结构变形、应力分布和相互作用,并实现对燃料重构、氧和钚元素的迁移、裂变气体释放和服役期内液态铅铋腐蚀等内容的计算模拟。     

Local thermal–hydraulic behaviour in tight 7-rod bundles

Advanced water-cooled reactor concepts with tight lattices have been proposed worldwide to improve the fuel utilization and the economic competitiveness. In the present work, experimental investigations were performed on thermal–hydraulic behaviour in tight hexagonal 7-rod bundles under both single-phase and two-phase conditions. Freon-12 was used as working fluid due to its convenient operating parameters. Tests were carried out under both single-phase and two-phase flow conditions. Rod surface temperatures are measured at a fixed axial elevation and in various circumferential positions. Test data with different radial power distributions are analyzed. Measured surface temperatures of unheated rods are used for the assessment of and comparison with numerical codes.In addition, numerical simulation using sub-channel analysis code MATRA and the computational fluid dynamics (CFD) code ANSYS-10 is carried out to understand the experimental data and to assess the validity of these codes in the prediction of flow and heat transfer behaviour in tight rod bundle geometries. Numerical results are compared with experimental data. A good agreement between the measured temperatures on the unheated rod surface and the CFD calculation is obtained. Both sub-channel analysis and CFD calculation indicates that the turbulent mixing in the tight rod bundle is significantly stronger than that computed with a well established correlation.     

An Apparatus for Continuously Measuring Radon-222 Exhalation from Ground

High burnup MOX and UO₂ test rods were prepared from the fuel rods irradiated in commercial BWRs. Each test rod was equipped with a fuel center thermocouple and reirradiated in the Halden boiling water reactor (HBWR) in Norway. The burnups of MOX and UO₂ test rods reached about 84GWd/tHM and 72GWd/t, respectively. Fuel temperature was measured continuously during the re-irradiation tests. Thermal conductivity change in high burnup fuel was evaluated from the results of comparison between the measured fuel temperature and the data calculated by using the fuel analysis code FEMAXI-6. The comparison results suggested that the thermal conductivity of MOX fuel pellets is comparable to that of UO₂ fuel pellets in the high burnup region around 80 GWd/t. It is probable that the impurity effect of Pu atoms gradually diminishes with increasing burnup because other factors that affect pellet thermal conductivity, such as the accumulation effect of soluble fission products and irradiation-induced defects in crystal lattice, become dominant in a high burnup region.     

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.     

CUPIDON-code decrivant le comportement thermique et mecanique d'un crayon combustible de reacteur a eau sous pression au cours d'un accident de depressurisation

In the scope of the PHEBUS experimental program to be performed in Cadarache on the behaviour of PWR's fuel assemblies under loss of coolant accidental conditions, a computer code has been developed to help designing the experimental rods and to contribute to the definition of the test runs.This code, dubbed CUPIDON, deals only with the thermal and mechanical behaviour of the rods as well as the oxidation of the cladding outside surface; it does not include any thermohydraulic subroutine. Rather, it is coupled with the RELAP code for providing necessary input data such as coolant temperatures and pressures and cladding-to-coolant heat transfer coefficients. It is restricted to a single, non irradiated, rod of short length as representing the PHEBUS experimental conditions.It is a two dimensional code, using a finite difference resolving technique. It calculates the radial thermal profile across each section of the rod, the stress and creep rate to which the cladding is submitted and the rate of formation of the oxide layer on the surface of the cladding under steady state and transient conditions. As cladding plastic strain input data, it is using the EDGAR-ZY experimental results.     

Evaluation of LWR fuel rod behavior under operational transient conditions

To evaluate the effects of fission gas flow and diffusion in the fuel-cladding gap on fuel rod thermal and mechanical behaviors in light water reactor (LWR) fuel rods under operational transient conditions, computer sub-programs which can calculate the gas flow and diffusion have been developed and integrated into the LWR fuel rod performance code BEAF. This integrated code also calculates transient temperature distribution in the fuel-pellet and cladding.The integrated code was applied to an analysis of Inter Ramp Project data, which showed that by taking into account the gas flow and diffusion effects, the calculated cladding damage indices predicted for the failed rods in the ramp test were consistent with iodine-SCC (Stress Corrosion Cracking) failure conditions which were obtained from out-of-reactor pressurized tube experiments with irradiated Zircaloy claddings. This consistency was not seen if the gas flow and diffusion effects were neglected. Evaluation were also made for the BWR 8 × 8 RJ fuel rod temperatures under power ramp conditions.     

Concept of innovative water reactor for flexible fuel cycle (FLWR)

In order to ensure sustainable energy supply in the future based on the matured light water reactor (LWR) and coming mixed oxide (MOX)-LWR technologies, a concept of innovative water reactor for flexible fuel cycle (FLWR) has been investigated in Japan Atomic Energy Research Agency (JAEA). The concept consists of two parts in the chronological sequence. The first part realizes a high conversion type core concept, which is basically intended to keep the smooth technical continuity from current LWR and coming MOX-LWR technologies without significant technical gaps. The second part represents the reduced-moderation water reactor (RMWR) core concept, which realizes a high conversion ratio over 1.0 being useful for the long-term sustainable energy supply through plutonium multiple recycling based on the well-developed LWR technologies. The key point is that the two core concepts utilize the compatible and the same size fuel assemblies, and hence, the former concept can proceed to the latter in the same reactor system, based flexibly on the future fuel cycle circumstances during the reactor operation period around 60 years.At present, reprocessed plutonium from the LWR spent fuel is to be utilized in MOX-LWR. After this stage, the first part of FLWR, i.e. the high conversion type, can be introduced as a replacement of LWR or MOX-LWR. Since the plutonium inventory of FLWR is much larger, the number of the reactor with MOX fuel will be significantly reduced compared to the MOX-LWR utilization. When the fuel cycle for plutonium multiple recycling with MOX fuel reprocessing is realized, the fuel assembly will be replaced with another type of the tight-lattice one for RMWR with different rod diameter, rod gap width and so forth even in the same reactor system, being flexibly corresponding to the fuel cycle circumstances.Investigation on the core for both the parts of the FLWR concepts has been performed, including the core conceptual design, the core characteristics under Pu multiple recycling, the thermal hydraulic investigation in the tight-lattice core, and so forth. Up to the present, promising results have been obtained.     

Conceptual Design Study of 180 MWt Small-Sized Reduced-Moderation Water Reactor Core

Conceptual design of a Small-sized Reduced-Moderation Water Reactor (S-RMWR) core, which has the thermal output of 180 MW, the conversion ratio of 1.0 and the void reactivity coefficient of negative value, has been constructed. S-RMWR is a technology demonstration reactor which also conducts material and fuel testing for commercial use of Reduced-Moderation Water Reactor (RMWR) in large-scale power plants. It has a very tight triangular fuel rod lattice and a high coolant void fraction. The RMWR core axially has two short and flat uranium plutonium mixed oxide (MOX) regions with an internal blanket region in between, in order to avoid a positive void reactivity coefficient. The MOX regions are sandwiched between upper and lower blanket regions, in order to increase a conversion ratio.

In this small reactor core, leakage of neutrons is expected to be larger than in a large core. Therefore, a core design concept different from that for a large core is necessary. Core burnup calculations and nuclear and thermal-hydraulic coupled calculations were performed in the present study with SRAC and MOSRA codes. MVP code was also used to obtain control rod worth. Because of its large neutron leakage, keeping the void reactivity coefficient negative is easier for S-RMWR than RMWR. Thus, the heights of MOX region can be taller and the plutonium enrichment can be lower than in RMWR. On the other hand, to achieve the conversion ratio of 1.0, radial blanket and stainless steel reflector assemblies are necessary, whereas they are not needed for RMWR.     

The prospect of mox fuel based PbBi cooled small nuclear power reactors

Safety performance of MOX fuel based PbBi cooled small fast power reactors has been analyzed and discussed. Though the thermal conductivity of MOX fuel is not large relative to that of nitride or metal fuel, but by proper combination of relatively small power density and relatively large natural circulation it can compensate fuel temperature decrease with coolant temperature increase smartly during unprotected loss of flow accident. Under such condition, accident analysis discussed in this paper show that under unprotected total loss of flow (ULOF) accident the reactor can survive inherently using combination of reactivity feedback. For unprotected rod run out transient over power (UTOP) accident the MOX reactor can overcome external reactivity by smaller power increase compared to that of nitride fueled reactors case. In this case doppler feedback plays much more important role compared to radial expansion component. So the MOX fueled small power reactors discussed here can survive both UTOP and ULOF accident with still enough temperature margin.     

Analysis of MOX Fuel Behavior in Halden Reactor by FEMAXI-6 Code

Capabilities of the FEMAXI-6 code to analyze the behavior of high burnup MOX fuels in LWRs have been evaluated. Coolant conditions, detailed power histories and specifications of the MIMAS-MOX fuel rods, rod 10 and rod 11, of IFA-597.4–7 irradiated in the Halden reactor were input, and calculated rod internal pressures and pellet center temperatures were compared with the measured data for the range of 0-31 MWd/kgUO₂. Some sensitivity studies were conducted mainly with respect to pellet thermal conductivity and swelling rate to investigate the changes in thermal behavior and their effects on fission gas release.

In the irradiation period up to about 23 MWd/kgUO₂, the calculated pellet center temperatures sufficiently agreed with the measured data and also the calculated rod internal pressures reproduced the tendency of an increase in the measured rod internal pressures. These results suggest that fission gas release from MOX fuels can be reasonably predicted by a diffusion process that is modeled in UO₂ pellet grains. On the other hand, the steep increase in the measured rod internal pressures observed at the power ramp around 23 MWd/kgUO₂ cannot be reproduced by FEMAXI-6 and can be regarded as the result of a relatively large amount of gas release, which possibly caused a pellet-cladding-gap closure through pellet gas-bubble swelling.     

Design of a fuel element for a lead-cooled fast reactor

The options of a lead-cooled fast reactor (LFR) of the fourth generation (GEN-IV) reactor with the electric power of 600 MW are investigated in the ELSY Project. The fuel selection, design and optimization are important steps of the project. Three types of fuel are considered as candidates: highly enriched Pu-U mixed oxide (MOX) fuel for the first core, the MOX containing between 2.5% and 5.0% of the minor actinides (MA) for next core and Pu-U-MA nitride fuel as an advanced option. Reference fuel rods with claddings made of T91 ferrite-martensitic steel and two alternative fuel assembly designs (one uses a closed hexagonal wrapper and the other is an open square variant without wrapper) have been assessed. This study focuses on the core variant with the closed hexagonal fuel assemblies. Based on the neutronic parameters provided by Monte-Carlo modeling with MCNP5 and ALEPH codes, simulations have been carried out to assess the long-term thermal-mechanical behaviour of the hottest fuel rods. A modified version of the fuel performance code FEMAXI-SCK-1, adapted for fast neutron spectrum, new fuels, cladding materials and coolant, was utilized for these calculations. The obtained results show that the fuel rods can withstand more than four effective full power years under the normal operation conditions without pellet-cladding mechanical interaction (PCMI). In a variant with solid fuel pellets, a mild PCMI can appear during the fifth year, however, it remains at an acceptable level up to the end of operation when the peak fuel pellet burnup ∼80 MW d kg⁻¹ of heavy metal (HM) and the maximum clad damage of about 82 displacements per atom (dpa) are reached. Annular pellets permit to delay PCMI for about 1 year. Based on the results of this simulation, further steps are envisioned for the optimization of the fuel rod design, aiming at achieving the fuel burnup of 100 MW d kg⁻¹ of HM.     

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

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

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.     

Artificial neural network modeling for fission gas release in LWR UO2 fuel under RIA conditions

A fission gas release (FGR) model was developed by using an artificial neural network method to predict fission gas release in UO₂ fuel under reactivity initiated accident (RIA) conditions. Based on the test data obtained in the CABRI test reactor and nuclear safety research reactor, the model takes into account the effect of the five parameters: pellet average burnup, peak fuel enthalpy, the ratio of peak fuel enthalpy to pulse width, fission gas release during base-irradiation, and grain size of a fuel pellet. The parametric study of the model, producing a physically reasonable trend of FGR for each parameter, shows that the pellet average burnup and the ratio of peak fuel enthalpy to pulse width are two of the most important parameters. Depending on the combination of input values for the five parameters, the application of the model to a fuel rod under typical RIA conditions of light water reactor produces 1.7-14.0% of FGR for the pellet average burnup ranging from 20 to 70 MW d/kg U.