SEARCH · Search NASA

Results for “Multigroup cross section”

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.

SourceResearch connectionFrom yearThrough year

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

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Verification of the ENDF/B-VII.1 Based MC 2 -3 Library Rev.1

The MC 2 -3 code, developed by Argonne National Laboratory under the DOE-NE NEAMS program, is a multigroup cross section generation code for fast reactor applications. Last year, the ENDF/B-VII.0 (E70) MC 2 -3 library, which has been extensively used, verified, and validated over a long period, was intensively reverified and updated to support the commercial grade dedication (CGD) requirement of the TerraPower Natrium project. This year, the ENDF/B-VII.1 (E71) MC 2 -3 library, the preliminary version of which was generated several years ago, was regenerated and rigorously verified to support the Natrium project as well as the completion of verification of the E71 library. The E71 library was verified using the process developed during the verification of the E70 library, including comparisons of cross sections with the NJOY-generated cross sections, comparisons of the resolved resonance cross sections with those using the PEDNF library, and comparison of total cross sections with the sum of partial cross sections. Additional verifications were conducted to ensure that the benchmark problem solutions with the E71 library are reasonable compared to the corresponding Monte Carlo solutions. Furthermore, the E71 gamma library was generated, which includes data for prompt gamma, delayed gamma, and delayed beta as well as neutron and gamma heating. The gamma library was verified at the level of individual isotopes. The EBR-II core solutions from MC 2 -3/ DIF3D and MCNP were compared, demonstrating that those solutions in terms of k-effective and assembly powers were in good agreement.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2019

L3:RTM.MCH.P19.06: Neutronics UQ

This report summarizes some of the work performed in the last year on the nuclear data uncertainty tasks using the multipole representation. This approach was proposed to provide a simple and integrated way of evaluating the nuclear data uncertainty in a Monte Carlo simulation and evaluating their impact on multigroup cross sections. As previously described, the windowed multipole format allows for a resonance based representation of the cross sections with direct Doppler broadening. This representation also lends itself nicely to the evaluation of sensitivities of resonance parameters and the direct perturbation of resonance parameters in the simulation environment. Both of these approaches rely on the need of conversion of the covariance data (File 32 of the ENDF library) into a pole and residues form: i.e. into a multipole covariance. Section 5 details the current pathways that have been used so far to achieve this. In order to verify the conversion process of the covariance data, the sensitivies to the resonance parameters and the sampling from these parameters, an analytical benchmark of the slowing process was developed. The benchmark provides a closed-form solution for the flux in the presence of Hydrogen and multiple resonant nuclides. Section 2 details the analytical benchmark and provides the closed-form solutions for the flux and its sensitivity to resonance parameters. A verification of the proposed methodology is performed and compared with the analytic solution to demonstrate the accuracy of the benchmark and the applicability of our embedded resonance perturbation approach in a small Monte Carlo code. Papers are currently being drafted documenting the benchmark and the generation of sensitivities. As paths to best generate multipole covariance matrices are being compared, we also explored new methods to propagate nuclear data uncertainty across the Monte Carlo neutron transport calculation. Section 4 studies the use of multipole covariance matrices to perform sensitivity analysis. Section 3.2 documents the results of our investigations into a new way of propagating nuclear data uncertainty across Monte Carlo transport simulations by sampling new epistemic cross sections for each new neutron.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2012

Development of a Statistical Sampling Method for Uncertainty Analysis with SCALE

A new statistical sampling sequence called Sampler has been developed for the SCALE code system. Random values for the input multigroup cross sections are determined by using the XSUSA program to sample uncertainty data provided in the SCALE covariance library. Using these samples, Sampler computes perturbed self-shielded cross sections and propagates the perturbed nuclear data through any specified SCALE analysis sequence, including those for criticality safety, lattice physics with depletion, and shielding calculations. Statistical analysis of the output distributions provides uncertainties and correlations in the desired responses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2012

HTTR criticality calculations with SCALE6: Studies of various geometric and unit-cell options in modeling

The fuel element of the High Temperature Engineering Test Reactor (HTTR) presents a doubly heterogeneous geometry, where tiny TRISO fuel particles dispersed in a graphite matrix form the fuel region of a cylindrical fuel rod, and a number of fuel rods together with moderator or reflector then constitute the lattice design of the core. In this study, a series of full-core HTTR criticality calculations were performed with the SCALE6 code system using various geometric and unit-cell options in order to systematically investigate their effects on neutronic analysis. Two geometric descriptions (ARRAY or HOLE) in SCALE6 can be used to construct a complicated and repeated model. The result shows that eliminating the use of HOLE in the HTTR geometric model can save the computation time by a factor of 4. Four unit-cell treatments for resonance self-shielding corrections in SCALE6 were tested to create problem-specific multigroup cross sections for the HTTR core model. Based on the same ENDF/B-VII cross-section library, their results were evaluated by comparing with continuous-energy calculations. The comparison indicates that the INFHOMMEDIUM result overestimates the system multiplication factor (k{sub eff}) by 55 mk, whereas the LATTICECELL and MULTIREGION treatments predict the k{sub eff} values with similar biases of approximately 10 mk overestimation. The DOUBLEHET result shows a more satisfactory agreement, about 4.2 mk underestimation in the k{sub eff} value. In addition, using cell-weighted cross sections instead of an explicit modeling of TRISO particles in fuel region can further reduce the computation time by a factor of 5 without sacrificing accuracy. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2012

Development of a statistical sampling method for uncertainty analysis with SCALE

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2013

An improved resonance self-shielding method for heterogeneous fast reactor assembly and core calculations

An improved resonance self-shielding method has been developed to accurately estimate the effective multigroup cross sections for heterogeneous fast reactor assembly and core calculations. In the method, the heterogeneity effect is considered by the use of isotopic escape cross sections while the resonance interference effect is accounted for through the narrow resonance approximation or slowing-down calculations for specific compositions. The isotopic escape cross sections are calculated by solving fixed-source transport equations with the method of characteristics for the whole problem domain. This method requires no pre-calculated resonance integral tables or parameters that are typically necessary in the subgroup method. Preliminary results for multi pin-cell fast reactor problems show that the escape cross sections estimated from the explicit-geometry fixed source calculations produce more accurate eigenvalue and self-shielded effective cross sections than those from conventional one-dimensional geometry models. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2013

Spatial and spectral effects in subcritical system pulsed experiments

Accurate neutronic models are needed for the interpretation of pulsed experiments in subcritical systems. In this work, the extent of spatial and spectral effects in the pulse propagation phenomena is investigated and the analysis is applied to the GUINEVERE experiment. The multigroup cross section data is generated by the Monte Carlo SERPENT code and the neutronic evolution following the source pulse is simulated by a kinetic diffusion code. The results presented show that important spatial and spectral aspects need to be properly accounted for and that a detailed energy approach may be needed to adequately capture the physical features of the system to the pulse injection. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2014

SNS Sample Activation Calculator Experimental Validation - Paper 79

The Spallation Neutron Source (SNS) at Oak Ridge National Laboratory (ORNL) uses the Sample Activation Calculator (SAC) to calculate the activation of all samples after they have been exposed to the neutron beam in any of the neutron scattering instruments at the SNS beamlines. The SAC web page takes user inputs (choice of beamline, the mass, composition and area of the sample, irradiation time, decay time, etc.) and then calculates the activation for the sample. In recent years, the SAC has been incorporated into the user proposal and sample handling process, and beamline scientists and users have noticed discrepancies in the predicted activation of their samples. As a part of the Neutronics Analysis Group at SNS, we were tasked with validating the SAC tool and diagnosing the origin of the discrepancies seen by the beamline scientists and the users. In the process of performing these measurements, we additionally found that several aspects of realistic neutron scattering operations, and thus sample exposure, cannot be captured with current configuration of the SAC tool. Typical neutron scattering samples are irradiated for multiple periods with different neutron energy ranges; the realities of accelerator operation indicate that not all of those periods will be at the same accelerator power, and scattering samples are frequently not fully illuminated by the well-collimated neutron beam. Finally, we will explore whether neutron beams, as modified by neutron guide systems with energy-dependent intensity modifications and by precise energy-band defining mechanical choppers, are suitably analyzed with multigroup cross section libraries in which there are typically very few energy groups below one eV. (authors)

43 PARTICLE ACCELERATORS↗

DOE OSTI2020

Verification of the CP2020 Library

For many years, “pencil beam” (aka “broomstick”) problems have been used in Monte Carlo neutronic code verification studies . At Los Alamos, they were used to verify the changes associated with the upgrade of NJOY / MCNP to allow continuous angular distributions (instead of discrete angular distributions) from S(α,β) scattering. Recently, “pencil beam” problems have been applied to the verification of charged particle data for MCNP. Another verification test for CP2020 was to compare the continuous energy reaction cross sections generated by ACER with multigroup reaction cross sections generated by GROUPR from the same evaluation file. This method was also used in the verification of the most recent S(α,β) neutron libraries at Los Alamos .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Analysis of SCALE Criticality and Sensitivity Calculations for Reflected HEU Cylinders [Abstract]

The SCALE code package offers several nuclear data libraries to support Monte Carlo (MC) transport, as well as MC-based derivation of $\kappa$ eff sensitivity and uncertainty (S/U) data. The CSAS sequence using the KENO MC code can utilize continuous-energy (CE) cross sections, or pre-generated multigroup (MG) cross section libraries. The use of MG libraries introduces bias into calculations in exchange for faster transport solutions. The TSUNAMI-3D sequence also utilizes KENO MC calculations. TSUNAMI-3D has two CE calculational methods: the Iterated Fission Probability (IFP) method, and the Contribution-Linked eigenvalue sensitivity/Uncertainty estimation via Tracklength importance CHaracterization (CLUTCH) method. Previous work has shown poor agreement between CLUTCH and confirmatory direct perturbation calculations in specific applications, e.g., fissionable and polyethylene reflectors

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Molten Salt Reactor Experiment Simulation using Shift/Griffin

The Department of Energy (DOE)’s NEAMS focuses its efforts on the development of advanced modeling and simulation (M&S) tools for light-water reactors (LWRs) and non–LWRs (i.e., molten salt reactors, high-temperature gas reactors, microreactors, and fast reactors). In the previous fiscal year, the Multiphysics Applications Driver Technical Area funded molten salt reactor (MSR) M&S at Oak Ridge National Laboratory (ORNL) to generate multigroup macroscopic cross sections with Shift for a MSRE 2D lattice model in Griffin. In addition, Shift’s capability to calculate gamma dose rates from activated components in the primary exchangers in a molten salt breeder reactor was also demonstrated. In fiscal year 2023, ORNL generated multigroup macroscopic cross sections using Shift for a 3D MSRE core model. MSRE depletion calculations using Griffin were also demonstrated in this fiscal year. For the depletion calculation, one-group microscopic cross sections for the 3D MSRE core were generated using Shift, and the decay transmutation library from ORIGEN was converted to an ISOXML file, which is required as input in Griffin. Several Monte Carlo codes, such as OpenMC and Serpent, were also used to benchmark and supplement multigroup cross sections generated by Shift. Multigroup libraries were generated with 8 and 20 group structures, and the study found the 8-group structure to be more accurate when comparing Griffin results to continuous energy (CE) Monte Carlo results. The average flux from CE Shift calculations is up to ~6% higher than the CE Serpent calculations because of different values applied for the energy released per fission (κ values). The average flux in the fuel salt calculated by Griffin using cross sections generated with Shift agrees well with the reference CE Shift solution; the same is valid for the corresponding Serpent results. The maximum relative error is ~6% and ~2% compared to the CE Shift and Serpent reference solutions, respectively. Meanwhile, the average flux calculated by Griffin in the graphite moderator shows a higher difference in the thermal range when compared to both reference Monte Carlo solutions; this result suggests a need for improvement in cross section generation for the graphite moderator in the thermal range in both Monte Carlo codes. Griffin depletion calculations using cross sections from Shift and Serpent were performed and compared against ORIGEN calculations, and the nuclide densities calculated by Griffin were found to be generally in agreement with those of ORIGEN. Because a different approach was taken to calculate the energy released per fission ( κ values) in Shift and Serpent, a difference in nuclide densities from differences in the average flux was observed between Griffin using Serpent and Shift cross sections. Griffin calculations with Shift cross sections produced higher average flux in the salt than with Serpent cross sections, leading to higher consumption of 235 U and higher production of 135 Xe. For time-dependent depletion calculations, cross sections were generated with Serpent, and Griffin’s results using these cross sections were compared to CE Serpent depletion results, demonstrating good agreement. The average difference in keff between Serpent and Griffin as a function of burnup is about 155 pcm. Similarly, good agreement with small differences up to ~0.5% was also noticed in the nuclide density of 235 U and 135 Xe. More details regarding the methodologies invoked to generate the cross section to make code-to-code comparisons are discussed further in this report. User feedback on Griffin and Shift capabilities that will enhance these calculations is provided in this report for future consideration. The work performed this fiscal year can be extended further for multiphysics coupling of Griffin-Pronghorn/SAM with Mole to study precursor flow and salt chemistry.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2011

The contribution of individual correlated parameters to the uncertainty of integral quantities.

Modern libraries of multigroup cross-sections often include covariances of the multigroup data. Data covariances are especially useful in calculating the variance of an important applied quantity z that can be calculated as a function of the data. Frequently, it is also of interest to know the magnitude of the contributions to this variance arising from various subsets of the data. This knowledge can be used, for example, in judging the benefits of performing specific new measurements, calculations, or evaluations on a given data subset. Here we describe a method for determining the contribution of an identified set of parameters to the variance of z that can be applied without difficulty even when the data have general correlations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2015

CEPXS

CEPXS is a multigroup-Legendre cross-section generating code. The cross sections produced by CEPXS enable coupled electron-photon transport calculations to be performed with multigroup radiation transport codes, e.g. MITS and SCEPTRE. CEPXS generates multigroup-Legendre cross sections for photons, electrons and positrons over the energy range from 100 MeV to 1.0 keV. The continuous slowing-down approximation is used for those electron interactions that result in small-energy losses. The extended transport correction is applied to the forward-peaked elastic scattering cross section for electrons. A standard multigroup-Legendre treatment is used for the other coupled electron-photon cross sections. CEPXS extracts electron cross-section information from the DATAPAC data set and photon cross-section information from Biggs-Lighthill data. The model that is used for ionization/relaxation in CEPXS is essentially the same as that employed in ITS.

Source record↗

DOE OSTI2012

Cross-section uncertainty propagation in the multigroup slowing-down equations using probability table formalism

In the frame of neutral-particle (neutron, gamma) transport, the uncertainty propagation calculation regarding the uncertainties on cross sections is often carried out without explicitly taking into account their probabilistic distribution. We investigate a new uncertainty propagation formalism where the cross-section uncertainty distributions are represented by probability tables. This technical note develops this approach for the steady-state slowing-down equation without up-scattering and in an infinite medium. This work is based on a deterministic multiband formalism that takes into account multilevel probability tables for cross sections. The first level represents the variation of cross sections versus lethargy (or energy) in each group of the multigroup lethargy mesh and thus corresponds to the classical cross-section probability tables. The higher levels represent the uncertainties on each step of the first-level cross-section probability table. This method is validated against a Monte Carlo calculation in a case of neutron slowing down in a {sup 238}U-hydrogen homogeneous mixture, showing fully consistent numerical results. The main interest of the deterministic multilevel multiband formalism is that it gives not only the mean value and the variance but also a probabilistic distribution of the fluxes. In the near future, we plan to investigate more deeply the robustness of this new approach in relation to high values of cross-section uncertainties and to introduce cross-section uncertainty correlations as well. Meanwhile, the promise of this work is its extension to the general transport steady-state equation solved by the discrete ordinates (S(N)) or Monte Carlo methods. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2025

Reduced-order modeling for efficient cross section library development in high-temperature gas reactor pebble-bed depletion analysis

Accurate modeling of running-in and equilibrium conditions in pebble-bed reactors (PBRs) requires precise microscopic multigroup neutron cross sections. In Griffin, deterministic neutronics calculations rely on multivariate interpolation over large cross section libraries, resulting in significant memory usage and performance bottlenecks. This work, together with a companion paper on Griffin integration, explores reduced-order models (ROMs) to replace interpolation with lightweight surrogates. Several ROM techniques are benchmarked, with deep neural networks (DNNs) demonstrating superior memory efficiency, scalability, and predictive accuracy. A total of 295 DNNs were trained to build a comprehensive isotope library, integrated into Griffin through a custom LibTorch interface for depletion analysis. Initial results demonstrate that DNN-based ROMs drastically reduce memory demands while preserving accuracy, enabling finer tabulations and additional state variables without overhead. In conclusion, the framework also supports online cross section generation and real-time DNN updates through transfer learning, improving fidelity by capturing self-shielding and evolving nuclide compositions during burnup.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2014

Demonstration of Full Core Reactor Depletion with MPACT

MPACT is a three-dimensional (3-D) whole core transport code capable of generating sub-pin level power distributions. This feature is accomplished by obtaining the integral transport solutions to the heterogeneous reactor problem in which the actual detailed geometrical configuration of fuel components such as the pellet and cladding is modelled explicitly.. The cross section data needed for the neutron transport calculation are obtained directly from a multigroup microscopic cross section library similar to those used in lattice physics codes. Hence MPACT involves neither a priori homogenization nor group condensation for the core spatial solution. The integral transport solution is obtained by means of the method of characteristics (MOC) that employs discrete ray tracing. Since the direct application of the 3-D MOC capability in MPACT to 3-D core configuration requires considerable amounts of memory and computing time for practical reactor applications, an alternative approximate 3-D solution method was implemented in MPACT based on a 2D-1D approach which employs planar MOC solutions in the framework of the 3-D coarse mesh finite difference (CMFD) formulation. The axial coupling is resolved by one-dimensional (1-D) diffusion solutions and the planar and axial problems are coupled through the transverse leakage. The use of a lower order 1-D solution in the axial direction is justified by the fact that most heterogeneity in the core occurs in the radial direction rather than the axial. However, it is possible to use one of the higher order methods in MPACT such as the 1-D SN kernel if improved solution accuracy is desired. The CMFD formulation which originally served as an efficient nodal formulation is used in MPACT to accelerate the whole core transport calculation. The basic mesh in the CMFD formulation is a pin cell which is much coarser than the flat source regions defined for MOC calculations. The concept of dynamic homogenization of group constants for the pin cells provides the basis of the CMFD formulation in whole core transport calculation. The intra-cell flux distribution determined from the MOC calculation is used to generate the homogenized cell constants while the MOC cell surface averaged currents are used to determine the radial nodal coupling coefficients. This dynamic implementation of the equivalence theory formulation ensures the same transport solution is obtained with CMFD as the one obtained using unaccelerated MOC calculations alone. The ability of MPACT to predict the Hot, Zero Power (HZP) condition of a reactor was reported in CASL-U-2014-0045-000. The purpose of this report is to demonstrate the capability of MPACT to simulate the operation of a reactor at Hot, Full Power (HFP) conditions during normal operation throughout a full operating cycle. The time scale for a reactor burnup cycle is typically one to two years and therefore short time phenomena such as the effect of delayed neutrons can be neglected. Instead the time dependence can be treated quasi-statically using the depletion algorithm presented in CASL-U-2013-0276-001. The full power operation requires thermal hydraulic feedback to be present during the entire simulation which will utilize the work in CASL-U-2014-0051-000. However, due to the significant computational burden of the coupled physics simulation, a simplified thermal hydraulic model is used for the analysis reported here which is principally intended to demonstrate the full core neutronics depletion. In addition to the thermal hydraulic feedback, the effect of Boron 10 depletion should be modeled in cases where major boron additions are not performed for extended periods of time. Each of these capabilities will be described in the following sections.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2016

A new multigroup method for cross-sections that vary rapidly in energy

Here, we present a numerical method for solving the time-independent thermal radiative transfer (TRT) equation or the neutron transport (NT) equation when the opacity (cross-section) varies rapidly in frequency (energy) on the microscale ε; ε corresponds to the characteristic spacing between absorption lines or resonances, and is much smaller than the macroscopic frequency (energy) variation of interest. The approach is based on a rigorous homogenization of the TRT/NT equation in the frequency (energy) variable. Discretization of the homogenized TRT/NT equation results in a multigroup-type system, and can therefore be solved by standard methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

Efficient Computation of Doppler-Broadened Elastic Scattering Kernel Moments Using Ladder-Operator Formulation

Anefficient routine for computing Legendre moments of the Doppler-broadened elastic scattering kernel, including resonance scattering effects, has been implemented in the ISOXML module of Griffin. Isotropic scattering in the center-of-mass system and the ideal gas model for target motion are assumed. A ladder-operator formulation is introduced to compute all Legendre moments from order 0 to N simultaneously, enabling near-linear scaling of computational cost with respect to the maximum Legendre order. A physics-based strategy for constructing outgoing energy grids has also been developed, in which a tailored base grid is combined with adaptive refinement to maintain accuracy while limiting the number of outgoing energy points. For energies between resonances, a constant cross-section model is employed to further reduce computational cost. In addition, a quantitative criterion is derived to determine isotope-wise cut-off incident energies based on a prescribed up-scattering probability coverage. For 238U, up to incident energies of approximately 75, 230, and 661 eV at 294, 900, and 2500 K (corresponding to a 2% up-scattering probability threshold), computation of P0 kernels requires 1–8 s and computation of P0–P5 kernels requires 0.4–4 min using a single thread, while maintaining 1–3% relative error in up-scattering probability. These results demonstrate that the proposed formulation enables accurate and computationally practical Doppler-broadened kernel generation for online multigroup cross-section production in Griffin.

Doppler-broadening↗