DOE OSTI2016
Oak Ridge National Laboratory (ORNL) and the US Nuclear Regulatory Commission (NRC) have initiated a multiyear project to investigate application of burnup credit (BUC) to boiling-water reactor (BWR) fuel in storage and transportation casks and spent fuel pools (SFPs). Ideally, a single combination of all operating parameters including moderator density and axial burnup profiles, control blade exposure, and others can be identified to allow a simple, bounding approach to BWR BUC. One aspect of this investigation involves examining the effect of control blade insertion during operation on discharged assembly reactivity. In the course of these calculations, an apparent source convergence problem was discovered for fuel that has experienced partial blade insertion for the entire lifetime of the assembly. The source convergence problem appears not to impact k{sub eff}, but it manifests in predictions of the axial fission density profile. This phenomenon is of interest in identifying areas of concern with BWR BUC since it has not been studied as thoroughly as BUC for pressurized-water (PWR) reactor fuel. There may also be interest in this configuration for further study of a realistic source Monte Carlo convergence problem as opposed to some earlier manufactured source convergence problems through careful selection of initial conditions. This summary presents the depletion conditions used to generate the problem, a description of the models, a background based on well-behaved calculations, the identification of the problem, and attempts made to resolve the problem. Only a brief description of the depletion conditions and models is possible within the limited scope of this summary. Finally, conclusions are presented. This summary identifies a source convergence problem in the process of investigating BUC for BWR spent fuel. The depletion scenario which created the depleted fuel and burnable absorber compositions included a control blade inserted in the bottom 14 nodes throughout the entire depletion. This insertion is slightly over half the length of the assembly, and it appears to create two regions of nearly equal reactivity. As expected, one region is near the top of the assembly, and the other region is near the middle of the assembly. This node includes significant {sup 239}Pu generation as a result of rodded operation, but it is also the top node of the full lattice zone. While the calculated k{sub eff} values appear to be statistically equivalent, the axial neutron production profile varies significantly among independent simulations of the same model. A range of standard approaches has been used to attempt to resolve the source convergence problem and determine the true axial neutron production profile. This profile is of interest in cask criticality safety analysis because it influences the interaction among assemblies loaded in the cask. The forthcoming SCALE 6.2 release will include source convergence testing based on Shannon entropy, which may assist in identifying these types of problems in the future. It should be noted that this source convergence problem results from the depletion conditions and the final distribution of fissile and absorber material within the discharged fuel assembly. The results, while entirely generated within the SCALE code system, are representative of any Monte Carlo program attempting to calculate k{sub eff} for such a system. The most important conclusion of this work is that the BUC criticality safety community must be aware of the potential for this type of source convergence problem. Analysts must approach such partially rodded scenarios with caution and critically examine the results. Fortunately, it is extremely unlikely that BWR assemblies would experience the bladed operation necessary to create this particular source convergence problem. These results may also have some bearing on the analysis of PWR BUC for new plant designs such as AP1000 that incorporate significant rodded operations. In a broader context, the source convergence problem identified here could have been magnified if it had occurred in a Monte Carlo depletion simulation. It may be possible to identify this type of problem in a depletion calculation by comparing the apparently converged fluxes between steps to ensure that no oscillations or other serious shifts are occurring. The depletion calculations performed in this work were performed with a two-dimensional deterministic code and are therefore not subject to this problem. These results may also raise questions about the validity of 2D depletion models with extruded geometry in the 3D k{sub eff} calculation. (authors)
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗