Search NASA⌕ Search

SEARCH · Search NASA

Results for “nuclear transport codes”

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.

At least 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

VARI3D & PERSENT: Perturbation and Sensitivity Analysis

The nodal diffusion method is one of the most widely used approaches in modern reactor analysis. In the nodal diffusion method, a coarse multi-group set of “homogenized” parameters is constructed such that the complex geometry of a reactor core along with the energy dependence of neutron and gamma ray cross sections in a nuclear reactor are conserved in the simpler geometry. The homogenization is typically done on a fuel assembly level as is the case in the DIF3D code developed at Argonne National Laboratory. The nodal methodology is used primarily to predict fuel cycle behavior of nuclear systems of which there is a substantial amount of validation in the literature. Another use of the nodal method is to obtain reactivity coefficients and kinetics parameters for use in a safety analysis of a given nuclear reactor. While there are many ways to obtain reactivity worth and kinetics parameters, the work presented in this manuscript is unique as it provides the user with the ability to compute reactivity worths, kinetics parameters, and cross section sensitivities with a Cartesian and hexagonal geometry based transport code. This manuscript serves as a single manual for two separate codes: VARI3D and PERSENT. The VARI3D code (VARIational 3D) is based upon the classic finite difference diffusion theory solver available in DIF3D. The PERSENT code (PERturbation and SENitivity for Transport) is based upon the variational nodal method employed in DIF3D termed VARIANT. The VARIANT solver was added to DIF3D in 1995 and has seen continued development and use for the last 18 years. Because VARI3D primarily uses deprecated coding practices, rather than incorporating the perturbation and sensitivity treatments for transport within VARI3D, a new coding development was built using modern Fortran coding. The primary purpose of this manual is to describe the theory behind PERSENT (and by convenience, that of VARI3D) and discuss the input and output of PERSENT along with giving potential users an idea of how to use it. While this manuscript does describe the input and output of VARI3D, the PERSENT code is intended to be the replacement capability of VARI3D as PERSENT can generate nearly identical (if not superior) diffusion theory results. In this manuscript, the relevant aspects of generalized perturbation theory and exact perturbation theory that apply to both VARI3D and PERSENT are covered. The input and output of VARI3D is displayed by excerpting several of the example problems. Similarly, the input and output of PERSENT is displayed along with tips on how best to use the code. Note that the input and output of the inhomogeneous solver wrapped around DIF3D (DIF3D_IFS) is also discussed as it is needed to carry out some of the sensitivities in PERSENT such as reaction rate ratios. This manuscript describes several perturbation and sensitivity problems, and the results computed using PERSENT. From these sections, potential users should find that PERSENT provides not only the typical tables of numbers desired in perturbation and sensitivity analysis work, but also can visually plot the result for a more thorough understanding of the space and energy distribution (Section 5). Overall, PERSENT is observed to produce accurate reactivity worths and sensitivities for the displayed set of test problems and clearly demonstrates the need to have a transport-based sensitivity capability as evident from the thousands of percent errors observed in the 21-group hexagonal fast reactor problem (covered in Section 7). The uncertainty calculation capability is described in Section 3 and demonstrated in Section 7.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Updates to Relevance Vector Machine: Multiclass Classification, Variable Selection, and Proof-of-Concept Application to Safeguards Fresh Fuel Verification using List-Mode Neutron Collar Data

To expand the capabilities of safeguards authorities to verify the integrity of fresh fuel assemblies, Oak Ridge National Laboratory has retrofit the existing electronics of the JCC-71 uranium neutron coincidence collar, which contains 18 3 He neutron detectors and an external 241 AmLi(α, n) neutron interrogation source arranged to surround a fresh nuclear fuel assembly. The new electronics system allows analysts to record list-mode neutron multiplicity data in addition to the singles and doubles rates that are currently measured. Based on previous proof-of-concept research, analysis of these new data will identify off-normal fuel configurations in an assembly and characterize or localize the specific partial fuel defects. The purpose of this report it to document the analysis algorithm development and then to demonstrate its capability for the safeguards verification of fresh fuel assemblies using list mode neutron collar data. To analyze the complex list-mode data collected with the upgraded uranium neutron collar, multivariate classification algorithms are being developed using a novel classification method, the relevance vector machine. This approach may be applied to multiclass problems to estimate the probability that test data belongs to one of many possible classes of data. In addition, our method identifies the most useful variables/channels for making predictions, which illuminates the basis for the model’s predictions, and this interpretability is largely unique among data analytics methods. Variable selection occurs during model training and parameter tuning and does not need any external hyperparameter tuning routines. Finally, we apply the modified relevance vector machine to a simulated dataset of list-mode neutron collar data generated with the radiation transport code MCNP. The method can correctly identify off-normal fuel configurations, categorize the data according to four fuel defect scenarios, and rank the channels in the data according to prediction utility. For nuclear safeguards applications, it is concluded that this method has the potential to increase the sensitivity and reliability to detect missing fuel rods from a standard 17 x 17 Pressurized Water Reactor (PWR) fresh fuel assembly. Within this analysis, “off-normal” (i.e., missing fuel rods) were correctly classified in 17 simulated test scenarios with one quarter (25%) of the fresh fuel rods missing using a training data set of 58 simulated measurements.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Comparisons of Nodal diffusion and whole-core transport methods for multiple cycles of a small Light water reacto

As part of a cooperative research and development agreement between the US department of Energy's Oak Ridge National Laboratory and Holtec International subsidiary SMR LLC, core simulator models of the SMR-160, a light water-cooled and -moderated, natural circulation, small modular reactor, were developed and compared. The model results indicate pin-by-pin power histories, critical boron concentration, and other key performance indicators useful for developing multi-cycle core loading patterns. One model uses industry standard nodal methods with a simple thermal-hydraulics solution, and the other uses an advanced deterministic neutronics solver with a coupled subchannel calculation. The goals of the comparison are to provide the reactor designer with a detailed and coupled core physics reference solution, inform the reactor designer of nodal model best practices based on the code-to-code comparison, and set up a workflow for performing subsequent detailed core simulator calculations. The results of this work will be used in the development of coupled neutronic, thermal hydraulic, and fuel performance models for advanced analysis of the SMR-160. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

The current status of inelastic and capture Gamma-ray production evaluations in translated ENDF-VIII.0 GNDS files and recommended remediation actions

This report provides the status of discrete nuclear levels and inelastic gamma-ray production for neutron induced reactions in current ENDF-VIII.0 evaluations for the isotopes in the GRIN project. Different categories of issues are identified after a comparison with information from “adopted” ENSDF files. Improvement strategies are given and recommendations are described. We have included a similar, but more limited, analysis for thermal capture data where EGAF and ENSDF thermal libraries are considered. In addition, we provide a validation plan to employ different transport code simulations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Bottomonium suppression in an open quantum system using the quantum trajectories method

We solve the Lindblad equation describing the Brownian motion of a Coulombic heavy quark-antiquark pair in a strongly coupled quark-gluon plasma using the highly efficient Monte Carlo wave-function method. The Lindblad equation has been derived in the framework of pNRQCD and fully accounts for the quantum and non-Abelian nature of the system. The hydrodynamics of the plasma is realistically implemented through a 3+1D dissipative hydrodynamics code. We compute the bottomonium nuclear modification factor and compare with the most recent LHC data. The computation does not rely on any free parameter, as it depends on two transport coefficients that have been evaluated independently in lattice QCD. Our final results, which include late-time feed down of excited states, agree well with the available data from LHC 5.02 TeV PbPb collisions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Characterization of Uranium Foil Irradiations at the WSU TRIGA Reactor Using a New Reactor Model in SCALE [Slides]

A depleted uranium foil and a highly enriched uranium foil that were irradiated at the Washington State University TRIGA reactor were used to measure the relative production of 237 U and 239 U, two key isotopes of interest in nuclear forensics, as well as their decay daughters 237 Np and 239 Pu. Total fissions in the two foils were also estimated using 99 Mo measurement. A reactor model of the TRIGA was built in the SCALE neutron transport code using dimensions obtained from an existing MCNP geometry. The model was successfully used to simulate foil fissions and the reaction rate ratios in 238 U, using the laboratory-measured isotope ratios as a benchmark. The flux distribution in the foils was determined. Finally, using the continuous energy Tsunami-3d module of SCALE, a sensitivity/uncertainty analysis was performed to determine the effects on the reaction ratio caused by other reactions in the reactor and by uncertainties in the known cross-sections.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2021

Multigroup Scattering in Monte Carlo Radiation Transport Codes [Slides]

Legendre truncations to multigroup scattering distributions are not amenable to Monte Carlo sampling due to negative values. We have implemented two moment-preserving methods in MGMC that capture the shape of the truncation, are non-negative over [-1, 1], and can be efficiently sampled on CPUs and GPUs. MGMC can now simulate neutrons with anisotropic scattering mechanics and MGMC shows good agreement with LANL production codes PARTISN.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Low level coupling scheme between neutronics and thermal-hydraulics based on Anderson acceleration

The simulation of nuclear reactors is a multiphysics problem mixing, amongst other fields, neutron transport and thermal-hydraulics. The simplest and most used approach in multiphysics simulation is based on the coupling of single-physics codes in a black-box fashion. However, in order to reduce the computational time needed for such simulations, case-dependent optimizations are often required. In this paper, we aim at reducing the computational time required to solve a coupled neutronic/thermal-hydraulic steady-state problem on a simplified Pressurized Water Reactor (PWR) core. The idea is to deal simultaneously with the coupling of the energy groups of the deterministic neutronic description of the core and its thermal-hydraulic description with the Anderson acceleration. By doing so, the fission source terms are directly accelerated instead of the power map as done in most cases. The power method used to solve the k-eigenvalue problem inside the neutronic solver is thus accelerated with the Anderson acceleration. The numerical experimentations conducted in this work are performed using APOLLO3 and THEDI, and indicate that such coupling strategy improves the convergence rates in terms of number of iterations required and the total computational time. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Rotational symmetry relation for efficient response function generation in the coarse mesh transport method COMET

The coarse mesh transport code COMET is a continuous energy hybrid stochastic-deterministic neutronics solver with high fidelity and formidable computation speed in solving reactor core problems. Its method is based on the incident flux expansion theory. In this work, we take advantage of the local geometric symmetry in many reactor cores lattices (e.g., fuel lattices and reflector blocks) to develop relations among the flux response expansion coefficients for symmetric surfaces to further improve the computational efficiency of the COMET response function generation tool (method). This is done by a rigorous derivation of the transformation matrices for the angular and spatial expansion moments resulting from a rotation of a coarse mesh by an arbitrary angle. The relations for the response coefficients for the symmetric surfaces can be then written as the Kronecker product of those transformation matrices. The method is implemented into COMET and tested on two advanced high temperature reactor (AHTR) full-length single assembly benchmark problems. The COMET results using the response function library based on the symmetry relations were compared to those using the library directly generated by continuous energy Monte Carlo for all surfaces. It was found that the eigenvalues and stripe-wise fission densities using the two libraries are in statistical agreement as expected. This indicates that the new method maintains the high fidelity of the original COMET method while improving the computational efficiency in the response function generation by 270% to 400%, depending on the local geometric symmetry. This method also reduces the size of the response function library by the same magnitude (270% to 400%). (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

SERPENT 2 and TRIPOLI-4 neutronics benchmarking of the Jules Horowitz reactor

Jules Horowitz Reactor (JHR) is a future center for nuclear material testing in Europe, currently under construction at CEA Cadarache in southern France. It has been modelled extensively with CEA's HORUS3D/N neutron calculation package, that includes a Monte-Carlo neutron transport code TRIPOLI-4. In order to support the future computational analyses on the reactor core physics, VTT's Monte Carlo neutronics code SERPENT 2 has been used to create an alternative model of the JHR, taking advantage of the pre-existing Constructive Solid Geometry (CSG) structures on which the TRIPOLI-4 model was built. This work describes a benchmark, in which two JHR core configurations were analyzed using both of these models. The first configuration is a fresh fuel core at the beginning of life, in which aluminum based dummy-devices are installed in the reactor test positions. The second configuration represents the end-of-cycle equilibrium core with a heterogeneous burnup distribution, in which CALIPSO devices containing cladding material samples are installed at the in-core experimental positions, while the reflector positions are occupied by DEN-REP devices containing UO{sub 2} fuel. Neutron flux spectra and heat deposition are the parameters of interest in this study, due to their high importance in creating the desired experimental irradiation environment. The impact of delayed gamma energy deposition will be studied in detail, since the two models account for it with different methods. This work aims to establish a SERPENT 2 model of JHR, understanding its differences and agreements with TRIPOLI-4, and discussing its potential future applications. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Proton beam power limits for stationary water-cooled tungsten target with different cladding materials

The proton beam power limit for a solid-tungsten spallation target is largely determined by beam induced thermomechanical structural loads and decay heat power deposition, while its lifetime is limited by radiation damage and fatigue life of the target materials. In this paper, we studied the power limits of a stationary water-cooled solid tungsten target concept. Tantalum clad tungsten was considered as a reference case. Being a low activation material, zircaloy 2 cladding option was studied and its decay heat driven power limit was compared with the reference case. Zirconium alloys have proven operations records in spallation target and nuclear fission environments, supported by materials data obtained from post irradiation examinations. Recent study also demonstrated feasibility of diffusion bonding zirconium to tungsten using vanadium foil inter layer. Particle transport simulations code FLUKA was used to calculate energy deposition and decay heat power deposition in the target, based on the beam parameters technically feasible at the Second Target Station of the Spallation Neutron Source at Oak Ridge National Laboratory. The energy deposition data were used for flow, thermal, and structural analyses to determine the beam intensity limit on the target concept studied. The decay heat deposition data were used to calculate the transient temperature evolution in the tungsten volumes in a loss of coolant accident (LOCA) scenario to determine its beam power limit. For a 1.3 GeV proton beam, the power limit on a stationary target was 400 kW for a tantalum clad target model and 800 kW for a zircaloy 2 clad target model.

Lee, Yong Joong↗
DOE OSTI2021

EMPIRE Simulations of the July 2020 photoelectron driven cavity B-Dot experiments at the National Ignition Facility

Simulations of several of the end-irradiated cylindrical photoelectron driven cavity experiments (also known as B-Dot cavities) that were fielded during the July 1 through 2, 2020 shot series at the National Ignition Facility are presented in this report with comparisons to experimental measurements. All cavity B-Dots fielded on the second, third, fourth, fifth and seventh shots were simulated using coupled Integrated Tiger Series (ITS) Monte Carlo transport codes and the Electromagnetic Plasmas in Realistic Environments (EMPIRE) electromagnetic particle-in-cell code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

(U) A Code System for Cross-Section Uncertainty Propagation for PARTISN

A rudimentary code system has been developed for propagating neutron cross-section uncertainties to response uncertainties using the PARTISN multigroup discrete-ordinates neutron transport code. The code system uses the first-order sandwich rule with sensitivities computed using SENSMG and covariances computed using NJOY. Databases of reaction cross sections and covariances using ENDF/B-VIII.0 data in the LANL standard 30-group structure have been precalculated. The final uncertainty in a k eff example problem compares well with the result from TOFFEE, a similar code that propagates neutron cross-section uncertainties to response uncertainties using the MCNP6 Monte Carlo neutron transport code. Presently, neither code system propagates uncertainties in $\overline{v}$, $\overline{μ}$, or χ, and neither code system propagates nuclide-nuclide covariances

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Implementation and testing of the GPT sensitivities in TRIPOLI-4

Computing the sensitivity of given reactor parameters with respect to nuclear data is mandatory in key applications such as the assessment of reactor safety or core design. In recent years, the development of Monte Carlo methods capable of estimating sensitivity coefficients based on the Standard or Generalized Perturbation Theory (SPT and GPT, respectively) has attracted intensive research efforts, especially in view of having reference golden standard simulations to be compared to faster but approximate deterministic calculations. Following our previous implementation and validation of the SPT functionalities in the Monte Carlo TRIPOLI-4, devoted to the sensitivity of the k eigenvalue, in this work we detail the development of GPT methods for ratios of reaction rates, and validate them against a few relevant benchmark results. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Simulation of Peach Bottom turbine trip transients TT1, TT2 and TT3 with SIMULATE-3K

The Peach Bottom 2 Turbine Trip transients tests TT1, TT2 and TT3 are simulated with SIMULATE-3K. This work provides the initial basis for evaluation of the capabilities and accuracy of coupled 3D neutronics/thermal hydraulics best estimate calculations for support of safety analysis. Selected results, obtained with SIMULATE-3K, are compared to available measurements and benchmark results. The comparisons are analyzed and the application of SIMULATE-3K for best estimate analysis approach within the safety analysis methodology is discussed. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Simulation of Peach Bottom turbine trip transient with RETRAN-3D using BXFGEN

The Peach Bottom 2 Turbine Trip transient is used to evaluate the capabilities and accuracy of coupled 3-D neutronics/thermal hydraulics best estimate calculations for support of safety analysis. Selected results, obtained with RETRAN-3D, are compared to available measurements and benchmark results. The comparisons are analyzed and the application of RETRAN-3D for a best estimate analysis approach within the safety analysis methodology is discussed. Overall, both the thermal-hydraulic and neutronic results show a good comparison to the experimental data, and the model is considered to accurately predict test data. In particular, the initial acoustic response of the main steam line and the subsequent main pressurization stages are well-predicted. During the main pressurization stage, the timing and magnitude of the in-core local pressurization and associated void compression are properly predicted based upon the agreement between local power range monitors (LPRM) predictions upwards, but the shift is slightly underpredicted. This indicates the methodology is well-posed for simulating this class of transient response.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗