Effect of Pu-Rich Agglomerate in MOX Fuel on a Lattice Calculation

The effect of Pu-rich agglomerates in U-Pu mixed oxide (MOX) fuel on a lattice calculation has been demonstrated. The Pu-rich agglomerate parameters are defined based on the measurement data of MIMAS-MOX and the focus is on the highly enriched MOX fuel in accordance with increased burnup resulting in a higher volume fraction of the Pu-rich agglomerates. The lattice calculations with a heterogeneous fuel model and a homogeneous fuel model are performed simulating the PWR 17 × 17 fuel assembly. The heterogeneous model individually treats the Pu-rich agglomerate and U-Pu matrix, whereas the homogeneous model homogenizes the compositions within the fuel pellet. A continuous-energy Monte Carlo burnup code, MVP-BURN, is used for burnup calculations up to 70 GWd/t. A statistical geometry model is applied in modeling a large number of Pu-rich agglomerates assuming that they are distributed randomly within the MOX fuel pellet. The calculated nuclear characteristics include k-inf, Pu isotopic compositions, power density and burnup of the Pu-rich agglomerates, as well as the pellet-averaged Pu compositions as a function of burnup. It is shown that the effect of Pu-rich agglomerates on the lattice calculation is negligibly small.     

MOX燃料堆芯热工特性及设计限值研究

使用MOX燃料的快堆核电站以其线功率高、燃耗高、堆芯出口温度高等特点,对堆芯热工设计提出了新的问题.本文在对MOX燃料热工性能分析的基础上,给出了主要的热工设计限值,并以电功率870 MW电站为参考,初步分析了其堆芯热工特性和设计裕量.结果表明对于MOX燃料,较高的堆芯热工参数合理可行,且具有足够的裕量.     

Statistical Uncertainty Analysis Applied to Fuel Depletion Calculations

In this paper, we present an uncertainty methodology based on a statistical approach, for assessing uncertainties in lattice code predictions of fuel composition changes with burnup due to uncertainties in the fuel geometric configuration, initial enrichment, and depletion conditions. The methodology has been applied to depletion calculations with CASMO-4 and experimental data from the ARIANE Programme to estimate the calculation uncertainties in nuclide concentration and other neutronic parameters at any time during the irradiation history. Results have shown that important information on the quality of the code's predictions can be obtained by analyzing the comparison of the code's estimates and their associated uncertainties, in the form of tolerance intervals, with experimental data and their reported errors.     

Sensitivity Study of Fuel Cost in Extended Burnup BWR Core

A sensitivity study on the fuel cost of an extended burnup BWR core has been carried out on the basis of a realistic model of discharge burnup extension. Full power operating length in months in a refueling cycle and the number of refueling batches are chosen as independent variables in the model to describe extended burnup cores of various types. The reference core for the sensitivity study adopts 9-month full power operation and 4-batch refueling scheme. The difference in the plant cost between the extended burnup core and the reference core, which is referred to as plant capacity factor (PCF) credit, is estimated and combined with the fuel cost to calculate the fuel cost with PCF credit.

The results show that the fuel cost with PCF credit decreases for the extended burnup core with stretched operating length as the burnup extends in cases of constant non-operating length in a cycle, and that it may increase for the extended burnup core with decreased batch number in cases of constant plant capacity factor. It is also suggested that the cost minimum combination of the independent variables can be found to a given discharge burnup for the extended burnup core with decreased batch number in an intermediate case between these two extreme cases. Extended burnup cores with fixed batch number tend to have a lower natural uranium requirement, but larger separative work requirement.

The economic break-even condition for the extended burnup core with decreased batch number is discussed based on the fraction of fixed part in the non-operating length, which is insensitive to the cycle length stretch-out.     

Fission Product Model for BWR Analysis with Improved Accuracy in High Burnup

A new fission product (FP) chain model has been studied to be used in a BWR lattice calculation. In attempting to establish the model, two requirements, i.e. the accuracy in predicting burnup reactivity and the easiness in practical application, are simultaneously considered. The resultant FP model consists of 81 explicit FP nuclides and two lumped pseudo nuclides having the absorption cross sections independent of burnup history and fuel composition. For the verification, extensive numerical tests covering over a wide range of operational conditions and fuel compositions have been carried out. The results indicate that the estimated errors in burnup reactivity are within 0.1Δk for exposures up to 100GWd/t. It is concluded that the present model can offer a high degree of accuracy for FP representation in BWR lattice calculation.     

混合氧化物燃料应用前景及核燃料循环战略

本文叙述了不同国家在不同核电发展条件下的核燃料循环战略初步考虑,并着重叙述了在快堆商业化进程推迟后所发展的混合氧化物(MOX)燃料的设计、制造技术及定量经济分析的基本情况,并建议我国近期应主要采用中间贮存并配以建造后处理中间试验工厂,为今后的商业后处理厂做好技术准备,同时还要跟踪 MOX 燃料技术。     

A Numerical Analysis for a BWR Fuel Assembly Drop Event

A BWR fuel assembly dropped from the crane hook during outage and clashed against the rack bottom plate of spent fuel pool. The area monitoring system indicated no radiation release, however, damage at the top of fuel channel was found in the following inspection. As fuel integrity is essential for further management, a finite element model was established to evaluate the damage condition. Several component elements including fuel rods, tie plates, and channel were set up and integrated into a full assembly. The analysis results provided the impact force on the fuel assembly and the dynamic response of each component element. The event did result in the damage of fuel channel yet fuel rods fracture was not expected. It's consistent with the inspection observations. There is therefore no concern for future operation such as interim dry storage.     

Irradiated Fuel Behavior under Power Oscillation Conditions

In order to examine high burnup fuel performance under power oscillation conditions, two tests of irradiated fuels under simulated power oscillation conditions were conducted in the Nuclear Safety Research Reactor (NSRR). Irradiated fuels at burnups of 56 and 25 GWd/tU were subjected to four to seven power oscillations, which peaked at 50 to 95 kW/m with intervals of 2 s. The power oscillations were caused by quick withdrawal and insertion of six regulating rods of the NSRR with a computerized control. Impacts of cyclic loads on the fuel performance under hypothetical unstable power oscillations arising during an anticipated transient without scram (ATWS) in boiling water reactors (BWRs) were examined in the tests. Deformation of the fuel cladding of the test rods was comparable to those observed in shorter transient tests, which simulated reactivity-initiated accidents (RIAs), at the same fuel enthalpy level up to 368 J/g. The fuel deformation was mainly caused by pellet-cladding mechanical interaction (PCMI) and was roughly proportional to the fuel enthalpy. Enhanced cladding deformation due to ratcheting by the cyclic load was not observed. Fission gas release, on the other hand, was considerably smaller than in the RIA tests, suggesting different release mechanisms in the two types of transients.     

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.     

Analysis of Core Physics Experiments on Fresh and Irradiated PWR UO2 Fuels in the REBUS Program

Critical experiments of two cores each loaded with fresh 5 × 5 test PWR-type fuel rods of ²³⁵U enrichment of 3.8 wt% or irradiated 5 × 5 test rods of rod average burnup of 55 GWd/t in the REBUS program were analyzed using diffusion, transport, and continuous-energy Monte Carlo calculation codes coupled with nuclear data libraries based on JENDL-3.2 and JENDL-3.3. Biases in effective multiplication factors k _eff's of the critical cores were about ?1:2%Δk for the diffusion calculations (JENDL-3.2), ?0:5%Δk for the transport calculations (JENDL-3.3), and ?0:5 and 0.1%Δk for the Monte Carlo calculations (JENDL-3.3 and JENDL-3.2, respectively). The measured core fission rate and Sc- or Co-activation rate distributions were generally well reproduced using the three types of calculation. The burnup reactivity determined using the measured water level reactivity coefficients was ?2:35 ± 0:07Δk/kk′. The calculated result of the Monte Carlo calculations agreed with it; however, the diffusion and transport calculations overestimated the absolute value by about 7%, which would be mainly attributed to the errors in the calculation of the reactivity caused by changing the fuel compositions from fresh fuel to irradiated fuel.     

Dependence of Fast Reactor Fuel Burnup Characteristics on Nuclear Data Libraries

In the framework of the development of burnup calculation method for commercial fast reactors, a sensitivity analysis was carried out to clarify the dependence of fuel burnup characteristics on nuclear data libraries (NDLs). The following NDLs were compared: JEF-2.2, ENDF/B-VI Release 5, JENDL-3.2 and JENDL-3.3. The NDL-dependence of material balance for main heavy metal nuclides (²³⁵U, ²³⁸U, ²³⁹Pu, ²⁴⁰Pu and ²⁴¹Pu) was small, since the number densities at the end of one-cycle burnup did not change over 1 or 2% among the above-mentioned libraries. Relatively large differences were found for minor actinide nuclides, especially for ²³⁶U, ²³⁷Np, ^242mAm, ²⁴³Am and curium isotopes. The number densities for these nuclides after burning up showed remarkable NDL-dependence over 5% through 50%. A burnup sensitivity analysis system based on the generalized perturbation theory enabled us to find out quantitatively the causative nuclides and reactions, as well as their energy regions.     

Effect of Moderator Density Distribution of Annular Flow on Fuel Assembly Neutronic Characteristics in Boiling Water Reactor Cores

The effect of the moderator density distribution of annular flow on the fuel assembly neutronic characteristics in a boiling water nuclear reactor was investigated using the SRAC95 code system. For the investigation, a model of annular flow for fuel assembly calculation was utilized. The results of the assembly calculation with the model (Method 1) and those of the fuel assembly calculation with the uniform void fraction distribution (Method 2) were compared. It was found that Method 2 underestimates the infinite multiplication factor in the fuel assembly including the gadolinia rod (type 1 assembly). This phenomenon is explained by the fact that the capture rate in the thermal energy region in gadolinia fuel is estimated to be smaller when the liquid film of annular flow at the fuel rod surface is considered. A burnup calculation was performed under the condition of a void fraction of 65% and a volumetric fraction of the liquid film in liquid phase of 1. It is found that Method 2 underestimates the infinite multiplication factor in comparison to Method 1 in the early stage of burnup, and that Method 2 becomes to overestimate the factor after a certain degree of burnup. This is because Method 2 overestimates the depletion rate of the gadolinia.     

Estimation of the Grain Boundary Gas Inventory in MIMAS/AUC MOX Fuel and Consistency with REP-Na Test Results

In the frame of its research activities on fuel safety, the French “Institut de Radioprotection et de Sûreté Nucléaire” performed the REP-Na program in the CABRI reactor devoted to the study of Reactivity Initiated Accidents. Focused on high burn-up UO2 and MOX fuel behaviour, twelve tests (8 UO2 and 4 MOX) were realized from 1993 to 2000. In all these tests, the influence of grain boundary gas was evidenced and it appeared necessary to perform some estimation of its inventory in irradiated fuel. Such evaluations are presented for the MOX MIMAS/AUC fuel, based on two different approaches: “experimental” and “theoretical.” The fission gas amount located at the grain boundaries increases with burn-up in correlation with the production, but also with the initial Pu enrichment, as soon as the agglomerates have reached the full restructuring threshold for the High Burn-up Structure. The consistency with the REP-Na test results is checked, showing that a significant cladding deformation is needed, clearly higher than for UO2 fuel in order to release all the grain boundary gas in RIA. Furthermore, to the fission gas effect, adds the helium's occluded in the irradiated fuel whose amount increases with burn-up, Pu enrichment and ²⁴¹Pu and ²⁴¹Am initial content.     

Expression of Few-Group Effective Microscopic Cross Sections in BWR Burnup Calculations

A practical method is proposed to express few-group effective microscopic cross sections for BWR burnup analysis. A set of few-group cross sections is prepared for an infinite square lattice of fuel rods as a function of the ratios of number density of nuclides such as ²³⁵ U, ²³⁸U and ²³⁹Pu, and the water quantity around a fuel rod. Spatial variation of few-group cross sections in the fuel assembly is taken into account by adjusting the water quantity around a fuel rod.

Numerical studies show that the present method can evaluate effective few-group cross sections within the accuracy of 3% in comparison with a two-dimensional integral transport calculation.     

Interpretation of the VENUS VIP-BWR program with APOLLO2.8 reference and optimized MOC schemes for 8 × 8 UO2 and MOX BWR assemblies

The interpretation of the VIP-BWR program conducted in the CEN·SCK Mol VENUS critical facility (Belgium), has been performed with the new APOLLO2.8 product and its CEA2005V4.1 library based on the JEFF3.1.1 file. Both reference SHEM-MOC (281groups without equivalence) and Optimized BWR 26G (26 groups with equivalence) schemes are used for UO₂ and MOX BWR assembly calculations. The VIP-BWR program was aimed to provide an experimental database for BWR neutronics tools in mixed Gd poisoned configurations with 8 × 8 UO₂ and MOX assemblies. The experimental conditions are relatively representative of actual industrial BWR core characteristics, at least in terms of void fraction. Measured pin-by-pin power distributions enable to exact valuable information at various interfaces. For fresh (UO₂/UO₂–Gd) and recycled UO₂ (UO₂ only) cores loadings, the information is given through the “UO₂” core. In the case of partial MOX loadings (UO₂/MOX interface), the power distributions are available through the “T-MOX” core. All critical sizes are predicted within 1 with SHEM-MOC reference calculation scheme. For UO₂ core, the (C–E) on k_eff are (95 ± 266) pcm and (203 ± 266) pcm for SHEM-MOC and Optimized scheme respectively. For MOX core, the results are (87 ± 214) pcm and (283 ± 214) pcm. The uncertainties take into account both measurement uncertainties and technological uncertainties such as enrichment, clad thicknesses, grid pitch or fuel densities.     

Prediction of Transient Boiling Transition in BWR Fuel Rod Bundles by Subchannel Analysis Program MENUETT

Transient CHF (critical heat flux) tests of a 4 X 4 rod bundle were analyzed by the subchannel analysis program MENUETT. MENUETT is based on a non-equilibrium, five equation, two-phase flow model and is available both for steady state and transient analyses. Turbulent mixing and void drift effects are taken into account to calculate cross flows in fuel rod bundles. The tendency of calculated subchannel mass fluxes and qualities agreed with experimental data. By using a critical quality correlation obtained from steady state CHF data, the position of the earliest boiling transition could be predicted regardless of non-uniform axial heat flux distributions. This transition occurrence time was predicted within a difference of 0.1～0.3 s from the experimental time. MENUETT applicability was confirmed for transient calculations predicting thermalhydraulic behavior in bundles.