Search NASA⌕ Search

SEARCH · Search NASA

Results for “Spatially dependent self-shielding”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

Spatially dependent embedded self-shielding method for nonuniform temperature distribution

The embedded self-shielding method (ESSM) is simple and robust, and applicable to any geometrical configuration to account for resonance self-shielding by performing fixed source transport calculations. ESSM is very similar to the subgroup method, which has been successful in the MPACT neutronics code developed by CASL. The subgroup method introduces errors in a least-squares fitting process that is used to generate subgroup quadrature sets and cannot be applied to thermal energy groups, as ESSM introduces error for spatially dependent self-shielded cross sections inside the fuel. This study focused on integrating the conventional ESSM with an improved Stoker-Weiss method to obtain spatially dependent resonance self-shielded cross sections for a nonuniform fuel temperature distribution. Here, the proposed method would serve as a substitute for the subgroup method in high fidelity simulations using CASL’s VERA. Benchmark results show that this new, spatially dependent ESSM (SDESSM) is an appropriate candidate to replace the subgroup method.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

The Extended Embedded Self-Shielding Method in SCALE 6.3/Polaris

The SCALE transport lattice code, Polaris, has been previously developed to generate few-group homogenized cross sections for whole-core nodal diffusion simulators in which the embedded self-shielding method (ESSM) is used for resonance self-shielding calculations to process cross sections. Although the ESSM capability has been very successful in light-water reactor analysis, it may require enhancements in computational efficiency; treatment of spatially dependent resonance self-shielding effects; and handling of interrelated resonance effects among fuel, cladding, and control rod materials. Therefore, this study focuses on improving computational efficiency by using a Dancoff-based Wigner–Seitz approximation combined with a material-based resonance categorization, through which a spatially dependent ESSM capability is developed to accurately estimate self-shielded cross sections inside the fuel. Benchmark results show that the new capability significantly enhances computational efficiency and accuracy for spatially dependent local zones within the fuel and through depletion.

ESSM↗

Radiation dose rate distributions of spent fuel dry casks estimated with MAVRIC based on detailed geometry and continuous-energy models

This study presents a detailed comparison of the dose rate distributions of a dry TN–32 fuel cask with two geometry models and two cross-sectional datasets. The accuracies of radiation dose rate estimation and computational efficiencies of each geometry model with two cross-sectional data-sets are compared. The use of automated variance reduction techniques can significantly improve the computational efficiency of such a realistic, deep penetration problem that involves radiation transport from different volumetric sources, thereby eliciting only a small statistical error. Monaco with Automated Variance Reduction using Importance Calculations (MAVRIC) is a computational sequence within the SCALE 6.2 code package based on consistent adjoint driven importance sampling (CADIS), a type of automated variance reduction technique. Homogenous and full fuel assembly models are built herein, and two nuclear cross-section libraries (V7–200N47G and continuous energy) are applied in this work. Based on the detailed comparisons, we found that neutron dose rate estimation is more dependent on geometry modeling than on cross-section data. For neutron-induced gamma rays, the dose rate distribution depends on both the spatial self-shielding effect and the cross-section library. The primary gamma rays respectively contribute to the total dose rate by ~91% and ~ 99% on side and top surfaces, and the dose rate accuracy is more dependent on the cross-section library than on geometry modeling. In terms of the computation efficiency and efforts spent on geometry modeling, the homogenous fuel assembly model with the MG library can produce an acceptable dose rate distribution, but the detailed fuel assembly model with the continuous-energy library is required for more precise dose rate estimation.

61 RADIATION PROTECTION AND DOSIMETRY↗

Group Structure Optimization Using the PyGroup Code

Group constants are used as coefficients in the energy domain discretized Linear Boltzmann Equation (LBE), in which the Continuous-Energy (CE) fluxes are also integrated over energy bins and transformed into group fluxes. This procedure, based on the Multi-Group (MG) theory, is widely applied in radiation transport and diffusion simulations Careful preparation of MG cross section data, or group constants, is thus a prerequisite for success in acquiring high-fidelity LBE solutions. Generation of MG cross section data itself is a complicated multi-step task, since group constants are indispensably coupled with the solution of the LBE, i.e. the flux spectrum in the problem to be solved. Theoretically, a MG cross section library is only valid for the problem that it is built upon. Nevertheless, in practice MG cross section libraries are prepared for similar type of problems, such as thermal reactors, fast reactors, etc., to facilitate the neutronic simulations. Generally, CE point-wise cross section data set, generated by both theoretical models and experimental data, is first converted to Fine-Group (FG) libraries using a pre-defined or precalculated flux spectra, then, Broad-Group (BG) libraries are collapsed from the fine-group libraries using problem-dependent flux spectra. Proper treatments of energy and spatial self-shielding and resonance cross section are required in the CE-to-MG converting process Previously, we have developed a computer code, YGROUP, for fine-to-broad group cross section collapsing. To maintain high accuracy of the resulting BG libraries, YGROUP provides users several options on collapsing strategies and weighting schemes for group binning and evaluation of BG cross section data in users' specific problems.. An approach based on Particle Swarm Optimization was applied to automatically determine the optimal group structure. We also developed a hybrid weighting scheme to further improve BG library performance for Special Nuclear Materials (SNM) detection applications. In this work, a new computer code, PyGroup, has been developed to implement the PSO approach with an improved simulation-driven fitness function design that includes both response and forward/adjoint flux errors. Together with YGROUP, and deterministic transport code TITAN/PENMSHXP, the PyGroup code package is tested on a Dual Range Coincident Counter (DRCC) model. In this paper, we demonstrated the capability of the PyGroup code for group structure optimization using the DRCC detector model. In the future, an iteration scheme can be implemented to determine the number of broad groups based on a user-defined tolerance. we will also apply PyGroup to more complicated reactor models, and investigate physics-based fitness function designs to further reduce computational cost. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Assessment of the MPACT Resonance Data Generation Procedure

Currently, heterogeneous models are being used to generate resonance self-shielded cross-section tables as a function of background cross sections for important nuclides such as 235 U and 238 U by performing the CENTRM (Continuous Energy Transport Model) slowing down calculation with the MOC (Method of Characteristics) spatial discretization and ESSM (Embedded Self-Shielding Method) calculations to obtain background cross sections. And then the resonance self-shielded cross section tables are converted into subgroup data which are to be used in estimating problem-dependent self-shielded cross sections in MPACT (Michigan Parallel Characteristics Transport Code). Although this procedure has been developed and thus resonance data have been generated and validated by benchmark calculations, assessment has never been performed to review if the resonance data are properly generated by the procedure and utilized in MPACT. This study focuses on assessing the procedure and a proper use in MPACT.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Best Practices for Shielding Analyses of Activated Metals and Spent Resins from Reactor Operation

This report investigated best practices for performing shielding evaluations of Type B waste packages, as defined in 10 CFR Part 71 regulations on packaging and transportation of radioactive material, or packages for which the contents are not defined before loading and may include a broad range of nuclides, geometries, and non-fuel materials. The following non-fuel waste streams were analyzed: Activated metals from decommissioned commercial power reactors, including Type 304 stainless steel, reactor vessel steel, and Inconel, Control blades from boiling water reactors, Neutron-activated corrosion products on surfaces of activated metals, and Spent resins from power plant operations. Measured elemental compositions, including major constituents and impurities, for steel and Inconel samples from commercial power reactors were used in activation calculations to determine radionuclide inventories in activated metals. For a simplified cask model, 60 Co contribution to the total external package dose rate at 30 days after shutdown varied from approximately 60% to 95%, depending on the activated metal, initial cobalt impurity concentration in the metal, and the thickness of the overpack gamma shield. Its maximum contribution to the total external package dose rate of approximately 100% was reached within the time interval of 2 to 5 years after shutdown and was maintained for up to 45 to 60 years after shutdown, depending on material, initial cobalt impurity concentration, and shield thickness. Thereafter, the 60 Co contribution to external package dose rate decreased with increasing decay time. Cobalt-60 is primarily produced by neutron reactions with the cobalt impurity in steel and Inconel. Other important radionuclides in activated metals contributing to package external dose rate are radionuclides with relatively short decay times, including 51 Cr, 59 Fe, 58 Co, and 54 Mn. These radionuclides may be represented as an equivalent 60 Co activity/source because 60 Co gamma ray emissions are bounding in terms of source strength and energy to other important radionuclides identified in the analyzed activated metals. Approaches for modeling the neutron-activated corrosion products that may be attached to activated reactor components were analyzed in this report. It was demonstrated that a surface source is more conservative than a uniform volumetric source for the treatment of neutron-activated corrosion products with respect to external package dose rates. An analysis of the maximum radionuclide loadings reported on spent resins identified the radionuclides 137 Cs, 60 Co, 134 Cs, 65 Zn, and 58 Co as the primary contributors to external package dose rate. For a resin cooled for 3.08 years, the external package dose rate was entirely produced by the reported 137 Cs and 60 Co inventory. The neutron sources from actinides found on spent resins or activated metals produced negligible dose rates and may be ignored in dose rate analyses. Effects of idealized waste material, source geometry, and spatial material/source distributions on external package dose rates were determined based on dose rate results for a simplified cask model under normal conditions of transport. Type 304 stainless steel, zirconium, and aluminum with adjusted mass densities based on a maximum content weight were analyzed for material modeling effects on external package dose rate. These materials produced identical external package dose rates within the statistical uncertainties of the dose rate estimates. Among four different source geometry configurations with homogeneous material of different mass densities, uniform volumetric source distribution, and the same source strength, the geometry configuration with lowest mass density (i.e., minimum self-shielding effect among the four cases) was most conservative. Spatial source distributions that better represent localized peak 60 Co activity values were more conservative than a uniform volumetric source distribution, assuming the same weight and total source strength per package. The increase in external dose rate caused by localized source peaks can be as much as the ratio between source peak activity density to average activity density, depending on the location of the activated metal with peak activity density. Therefore, the shielding analysis may be simplified if localized peak activities and the average activity per package can be measured/determined and documented at the time of cask loading. For simplicity, external package dose rates may be determined based on average source activity and a uniform volumetric source distribution. The dose rate results from that calculation model multiplied by the ratio of peak activity density to average activity density will produce maximum dose rate values for conservative estimates.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Multimodal Nondestructive Dry Cask Basket Structure and Spent Fuel Evaluation

This NEUP-IRP assembled a team from diverse institutions and a diverse set of expertise to examine the problem of assessing the structural health of dry nuclear fuel storage casks, specifically the TN-32 design. With no physical access designed into these storage casks and the high cost and risk associated with opening casks for direct inspection, the technical difficulty of structural assessment of the basket structure and fuel assemblies is high. In response to this challenge, this IRP brought together experts in both traditional and novel NDE methodologies: active acoustics (University of Mississippi), emissions tomography (University of Florida), muon tomography (Oregon State University), and passive-active ultrasonics (University of South Carolina). Further, our industry partners in the Electric Power Research Institute (EPRI) and AREVA provided invaluable insight into technical design details for the cask, practical knowledge of operating parameters at a storage facility, and key relationships that enabled activities unforeseen at the time of the proposed work but were critical in guiding the path of the research. The team decided early in the project that we had a collective goal of laying the technical ground work for practical solutions that could be implemented in the field with reasonable cost, time, and reliability. A fortunate opportunity arose early in this IRP project to obtain access to a new TN-32 cask at Columbiana High Tech. Both acoustics teams were given access for testing concepts in the first year of the project. This ended up providing critical guidance in the design of subsequent laboratory experiments. Specifically, it became clear that the only reliable acoustic pathways to the interior were through the trunnions, base plate, and top ring of the cask. Based on these field tests and discussions about what would be of practical interest, the team as a whole decided a realistic goal was to detect catastrophic failure of a single fuel assembly, or one that was missing. It is clear from this work that the active acoustics / resonance and muon tomography methodologies would be capable of detecting interior structural integrity problems at this level. The other two areas of investigation (emissions tomography and passive/active ultrasonics) would need more testing to determine a successful implementation, but each may also have other applications such as long term monitoring and monitoring during transportation. This IRP project resulted in technical advances with applications that extend beyond the problem at hand, supported the training of 10 graduate students + 1 post doc, 21 publications, and 13 conference presentations. The remainder of this summary highlights significant activities and achievements in each of the methodologies and a brief discussion on how they can be combined to leverage strengths. Active Acoustics The active acoustics team at UM used both linear and non-linear structural resonance methods to infer information about the internal structure of a cask. Significant activities were a set of four field test measurements on a full scale TN-32 cask at the Columbiana Hi Tech facility in Greensboro, NC and the design and fabrication of a 1/6th scale model TN-32 cask for laboratory measurements. A significant set of acoustic data has been acquired on both the full scale and laboratory scale casks. It was determined early in the project that mechanical resonance techniques were not going to be sensitive to small scale (centimeter) cracks and structural defects in the basket structure or fuel assemblies. We were able to ascertain that both the energy and complexity of the resonance spectra were sensitive to “bulk state” of the fuel assemblies – an intact assembly versus a severely degraded assembly. By developing a novel combination of metrics, we are able to measure a statistically significant shift due to one out of the 32 fuel assemblies being either missing or degraded to a point of complete loss of structural coherence (i.e. “rubble”). We believe this combination of resonance metrics is new to the field and will have applications beyond the fuel storage problem. Emissions Tomography The Emissions Tomography team at the University of Florida explored the concept of using scattering of radiation emissions from the stored fuel as method of evaluating the internal structural health of the cask. Due to the need for emissive fuel assemblies and the difficulties of gaining access to storage facilities, this team primarily relied on simulations to inform feasibility. A key question was the flux estimates of neutrons with a small enough number of scattering events to retain spatial information. Some key conclusions from the transmission component of this work are that the number of scattering events falls rapidly with energy – setting a 11MeV threshold, the average number of scattering events is 8∓4. Rough estimates of scanning times using an ideal detector array would be 1 week. The team also explored the case of significant structural failure resulting in a shift of the source material (fuel). The result is that a +/- 2% mass shift results in a detectable (5.13%) shift in the correlated 2 neutron emission counts and thus may be a candidate as an initial survey. Finally, simulations for predicting the azimuthal dependence of the flux based on specific loading patterns were conducted. Due to self-shielding effects of the assemblies, a periodic pattern results and could potentially be used to identify missing or severely degraded assemblies. This work forms a basis for proper detector placement with the goal of optimizing sensitivity. Muon Tomography Muon tomography (MT) is a relatively novel method of interrogating interior structure on a large scale. One of the advantages of this technique is that the exciting energy source (cosmic muons) is ubiquitous and relatively stable in time. Challenging aspects are relatively low flux rates and appropriate detector technologies. The efforts of the muon group at Oregon State University focused on two areas: (1) simulation of various combinations of muon scattering events and detection schemes to assess sensitivities and time scales for image construction and (2) building a prototype muon detection system applicable to the dry storage cask problem as well as many other similar systems. Significant accomplishments in this area are: (1) demonstration that MT is a viable method for detecting missing or severely degraded rod assemblies, (2) a scintillation based detector with 4 panels and 32 channels of position data each were constructed, tested, and calibrated. Testing of the panels using natural muon background and no scattering objects yielded largely expected results – a straight line path for the muon through the four panels. However some anomalies (paths deviating from a line) were also detected. A redesign of the board improved the signal which greatly improved the muon tracking through the panels. Testing with lead bricks as scattering bodies showed a reduction in the coincidence rate with increasing bricks indicating increased scattering events. A refinement of the scintillator detector design resulted in a spatial image resolution of around 2 cm based on lead brick spacing tests lasting 24 hours. Passive/Active Ultrasonics This team at the University of South Carolina investigated the feasibility of combining passive ultrasonics, in which an elastic pulse is generated by a degradation process (such as a crack extending), and active ultrasonics, in which a known excitation is introduced into the structure and detected by an array of transducers. Significant accomplishments in this area include the development of elastic wave propagation simulation in relevant materials, development of diagnostic algorithms, and laboratory testing and validation using laser doppler vibrometry. This group demonstrated the applicability of various non-contact vibration detection methodologies such as air-coupled ultrasonic transducers and scanning laser doppler vibrometry. While passive acoustic emissions monitoring will likely have a role during the transportation of a fuel storage cask, the level of detailed information about the defect generating the event which can be ascertained from the detections in a complex structure remains for future work. Thermography The EPRI team contributed early in the project a series of thermography measurements in which various levels of water pooled in a mock-up canister were analyzed using commercially available infrared imaging technologies. While large levels of water (100 L, ~10 cm deep) were observable in the thermograms, more realistic levels of pooled water (1 L) were not detectable whether localized or distributed.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Evaluation of the Effect of Burnup Acceleration on UO 2 Microstructure Evolution

Accelerated fuel qualification (AFQ) is a methodology by which new nuclear fuels are developed in an accelerated time frame compared with historical fuel qualification approaches. AFQ generally relies on high-fidelity physics-based modeling and simulation tools to adequately describe fuel performance as well as on revolutionary methods to accelerate burnup accumulation and collect relevant data more quickly. This report summarizes the use of advanced fuel modeling and simulation tools to evaluate microstructures from commercially irradiated fuel and microstructures from proposed MiniFuel irradiations, in which burnup accumulation is accelerated while prototypic temperature conditions are maintained. In this milestone, we used the mesoscale fuel performance code MARMOT to model the evolution of irradiated UO 2 microstructures and their potential restructuring at high burnup. The simulation conditions were informed by BISON models of both commercially irradiated fuel and MiniFuel. A first set of simulations investigated the recrystallization behavior of fully dense microstructures and showed full recrystallization at burnups as low as 52 MWd/kgU at 950°C. However, these simulations did not account for the presence of fission gas bubbles (predicted by BISON). Therefore, a second set of simulations including fission gas bubbles was performed and indicated that at the lowest temperature considered (650°C), the porous UO 2 microstructures have the highest total Gibbs free energies and are likely to recrystallize earlier than higher temperature cases (800 and 950°C), which agrees with high-burnup fuel characterization data. The results also showed that at lower temperature (650°C), the total free energies of the PWR fuel and MiniFuel microstructures are not significantly different. However, at the highest temperature (950°C), MiniFuel microstructures have a lower free energy than that of the PWR fuel microstructure. The competing effects between the temperature-dependent grain nucleation rate and the reduction of the free energy of the microstructure at higher temperature as a result of diffusion indicated that restructuring may occur at even higher temperatures than those considered in this study. In addition to the microstructure evolution modeling efforts, the burnup gradient across a single fuel specimen was also considered. This evaluation was for the VXF-15 position of the High Flux Isotope Reactor (HFIR) using the code suite HFIRCON, which was developed to automate the workflow for evaluating targets and fuel as they are irradiated in HFIR. The burnup gradient evaluation showed a dependence on both the axial and radial locations within the specimen, with a maximum difference of 1.7 between the inner and outermost radial layers. This relationship was further supported by considering the fission product speciation with respect to location within the specimen, which showed a higher concentration of 239 Pu, 240 Pu, and 241 Pu on the outer radial locations of the specimen than the center. The findings of the burnup and speciation evaluation show that some amount of self-shielding is occurring in the specimen when irradiated in the high-flux environment of HFIR; however, this impact is more pronounced for natural uranium when compared to 6% enrichment due to the higher ratio of 238 U in the specimen. Further analyses are required to understand the sensitivity of this gradient to spatial mesh and enrichment of the specimen.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A Full-Core Resonance Self-Shielding Method Using a Continuous-Energy Quasi–One-Dimensional Slowing-Down Solution that Accounts for Temperature-Dependent Fuel Subregions and Resonance Interference

In this paper, a correction-based resonance self-shielding method is developed that allows annular subdivision of the fuel rod. The method performs the conventional iteration of the embedded self-shielding method (ESSM) without subdivision of the fuel to capture the interpin shielding effect. The resultant self-shielded cross sections are modified by correction factors incorporating the intrapin effects of radial variation of the shielded cross section, radial temperature distribution, and resonance interference. A quasi–one-dimensional slowing-down equation is developed to calculate such correction factors. The method is implemented in the DeCART code and compared with the conventional ESSM and subgroup method with benchmark MCNP results. The new method yields substantially improved results for both spatially dependent reaction rates and eigenvalues for typical pressurized water reactor pin cell cases with uniform and nonuniform fuel temperature profiles. Finally, the new method is also proved effective in treating assembly heterogeneity and complex material composition such as mixed oxide fuel, where resonance interference is much more intense.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗