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 OSTI2021

SCALE 6.2.4 Validation for Light Water Reactor Decay Heat Analysis

Energy release from the decay of radionuclides in nuclear fuel after its discharge from reactor is a critical parameter for design, safety, and licensing analyses of used nuclear fuel storage, transportation, and repository systems. Well-validated computational tools and nuclear data are essential for decay heat prediction. This paper summarizes the validation of the SCALE nuclear analysis code system version 6.2.4, used with ENDF/B-VII.1 libraries, for decay heat analysis of light water reactor used fuel. The experimental data used for validation include full-assembly decay heat measurements that cover assembly burnups of 5 to 51 GWd/tonne U, cooling times after discharge in the 2- to 27-year range, and initial fuel enrichments up to 4 wt% 235U. The comparison between calculated (C) and experimental (E) decay heat showed very good agreement, with an average C/E over all considered measurements of 1.006 (σ = 0.016) for pressurized water reactor and 0.984 (σ = 0.077) for boiling water reactor assembly measurements. The effect of using assembly-average versus axially varying modeling data on the calculated decay heat, important to thermal analyses for used fuel transportation and storage systems, is discussed.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Cask RT-100 Type B(U) Self-Certification for DOT 7A, Type A Nuclear Waste Transportation - 20098

The purpose of this Type A certification is to diversify the use of the RT-100 cask by allowing the end user to transport either Type A or Type B radioactive waste in accordance with the handling procedures for each waste category. The project's scope included the development of a matrix that contains all applicable requirements for a Type A package from the Code of Federal Regulations (CFR), Title 49 for the Department of Transportation (DOT) with a justification for each requirement. Moreover, the scope included the writing of conditions and restrictions list for the transportation of Type A nuclear waste, the writing of a compliance statement, and the design and fabrication of a Type A placard to be permanently placed on the RT-100 cask. Type A packages are self-certified by the end user using one of three methods as permitted by 49 CFR 173.461, either by performance of tests with prototypes or samples of the specimens representing the packaging, or reference to a previous satisfactory demonstration of compliance of a sufficiently similar nature, or calculations or reasoned evaluation using reliable and conservative procedures and parameters. Since the RT-100 cask was licensed by the Nuclear Regulatory Commission (NRC) as a Type B cask, Robatel relied on the method of reference to a previous satisfactory demonstration of compliance of a sufficiently similar nature for Type A certification. This method was determined to be acceptable, since Type B packages are evaluated for more demanding transport conditions than those required for Type A packages. This paper presents Robatel's process for certifying the RT-100 cask for Type A shipments, an overview of the DOT regulations related to Type A packages, and the documents that were generated for the end user. This paper is concluded with a discussion on the advantages of having the option to transport both Type A and Type B nuclear waste, and possible future projects to further improve the operation of the cask when making Type A shipments. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Enhancing Monte Carlo Workflows for Nuclear Reactor Analysis with Metamodel-Driven Modeling

Monte Carlo codes are essential components of many reactor physics simulation workflows as high-fidelity continuous-energy neutron transport solvers. Among Monte Carlo radiation transport codes, MCNP is particularly notable due to its diverse simulation capabilities, large user base, and long validation history. Despite being a powerful simulation tool, MCNP provides limited capabilities to allow automated execution, model transformation, or support for user-defined logic and abstractions that limit its compatibility with modern workflows. Here, to better integrate MCNP into a modern scientific workflow, we have developed an intuitive yet full-featured MCNP Application Program Interface (API) in Python, named MCNPy, which provides a specialized set of classes for MCNP input development. Moreover, to guarantee that our reading, writing, and modeling capabilities remain self-consistent (and to render the huge scope of the MCNP API manageable), we have adopted a strategy of model-driven software development in which a generalized model of the MCNP input format has been created. From this generalized model, or “metamodel,” problem-specific implementations such as an engine for input validation or a codebase for programmatic operations may be automatically generated. Since MCNPy primarily acts as a Python front-end to the underlying Java API that directly interfaces with the metamodel, it is intrinsically linked to the metamodel and thus remains maintainable. With MCNPy, users can programmatically read, write, and modify any syntactically valid MCNP input file regardless of its origin. These capabilities allow users to automate complicated tasks like design optimization and model translation for nuclear systems. As examples, this work demonstrates the use of MCNPy to find the critical radius of a plutonium sphere and to translate a 9000+ line MCNP input file into a corresponding OpenMC model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Validation of MPACT BWR capabilities against critical experiments

Over the years, significant validation work for the neutronics code MPACT has been performed against zero-power critical benchmarks and measured data from operating nuclear power plants. As a part of the Modeling and Analysis of Exelon BWRs for Eigenvalue and Thermal Limits Predictability project, new validation efforts relevant to boiling water reactor (BWR) core applications have been performed. This paper presents the results of the critical experiment portion of the BWR validation efforts for MPACT. The Kritz-4 experiments and IPEN/MB-01 BWR-relevant configurations are modeled with MPACT. The Kritz-4 BWR critical experiments were performed at both cold and hot conditions, which is a unique feature among other critical experiment facilities. The MPACT results with linear source (LS) Method of Characteristics (MOC) have very good agreement with the measured criticality. When using the 60-group BWR library with LS and P2 scattering, the maximum k{sub eff} error is 109 pcm, and the average cold-hot bias is only -0.11 pcm/K. The fission rate distributions of MPACT are also verified with Serpent calculations. For the IPEN facility, two configurations are considered: one with a large central void, and one with a cruciform control rod. The MPACT k{sub eff} errors are less than 50 pcm for all IPEN problems. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Uncertainty analysis for VERA problem 2 using the cell-code Condor v2.8.05

Condor is a cell-level neutronic calculation code that applies multi-group collision probabilities with heterogeneous response coupling method within generic geometry configurations. Under the Condor's code continuous development, the incorporation of up-to-date methodologies and state-of-the-art practices in reactor analysis represents a driving force. In this work, the capabilities of Condor v2.8.05 to develop an uncertainty analysis for realistic PWR-kind fuel assemblies are studied. The Total Monte Carlo approach is applied to quantify the impact of fabrication tolerances in the code's results for reactivity and power distributions, by means of randomly sampled input values using the VERA problem 2 as basis. The VERA problem 2 proposes a series of Westinghouse 2D 17 x 17-type fuel lattices, to be calculated reflected at beginning-of-life without Xe. The configurations correspond to a modern PWR. Selected neutronic parameters from Condor runs are thus analyzed in terms of the observed spread as well as the obtained distributions for the randomly perturbed cases, showing the capability of the code to handle the required input data, as well as its ability to provide valuable insights regarding uncertainty quantification.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS1994

LET spectra of fragmented 600-MeV/A iron beam on aluminum and polyethylene: A preliminary analysis

An iron beam experiment recently conducted at the Lawrence Berkeley Laboratory Bevalac by Benton et al. provides an opportunity for verifying the new Green's function computer code (GRNTRN) and assessing the related nuclear database. The iron beam with 600 MeV/A at extraction traversed a series of beam transport elements, lead foil, and several triggering devices before impacting the target. Of these, only the 2.24 g/sq cm lead foil and target are considered in the transport analysis with an assumed (inferred from calibration) beam energy of 557 MeV/A. A thin layer of CR-39 plastic nuclear track detector (PNTD) was placed in front of the target to monitor the incident flux and a stack of four PNTD's placed behind to measure the linear energy transfer (LET) distribution of the transported beam. Test data are analyzed for three separate targets: 2 g/sq cm aluminum, 5-cm polyethylene, and 8-cm polyethylene. The two-layer GRNTRN results were mapped into the detector response function for comparison with the measured LET spectra. Reasonable agreement is obtained. Future research and analysis can be improved by using a more accurate isotope set or including other important media which significantly alter the beam. The assumption that the fragmentation cross sections are too small for aluminum is inconclusive because the production of fragments by 2 g/sq cm aluminum can be more significantly affected by the fragment contribution originating from the media in front of the target than that currently considered in the analysis.

Shinn, Judy L.↗
NASA NTRSDate not supplied

Accelerating SNP Development with Radiation Source Approximations for Far-Field Environment Modeling

The successful development and deployment of space nuclear power and propulsion technology requires a high-fidelity, efficient means of modeling ex-core radiation fields. Monte Carlo particle transport codes provide enable the highest level of fidelity in radiation field modeling but are inefficient for modeling ex-core radiation without modification. This paper details the development of an approach to ex-core radiation field modeling that maintains the fidelity of a Monte Carlo transport approach without sacrificing computational efficiency. This approach is based on recent improvements to an existing Monte Carlo transport acceleration technique known as surface source banking. As the transport of many particles is required to ensure adequate uncertainty in far field transport, analytical reconstruction techniques are employed to represent surface source banks as sampleable series of distributions in particle phase from which an arbitrarily large number of particles can be generated and simulated. All steps of this approach are incorporated into an AMA model to test both the underlying mathematics of the analytical source reconstruction process, as well as the relevance of the technique as a whole towards SNP applications. These tests demonstrated that this approach is capable of accelerating far-field radiation modeling in SNP-relevant scenarios by factors of at least nine over purely eigenvalue-based scenarios without significantly sacrificing simulation accuracy.

Surface Source Banking↗
DOE OSTI2022

Energy deposition model in STREAM with photon transport

Recent efforts have been taken to implement a photon transport module for the neutron transport analysis code STREAM, considering the explicit photon energy distribution in the problem instead of using the locally deposition assumption as employed in previous STREAM versions. The current energy release model in STREAM based on the On-the-fly energy released per fission model is modified accordingly to account for the photon energy deposition computed explicitly by this module. The impacts of the modified energy release model with photon transport in STREAM on the power distribution, fuel temperature and other core parameters during depletion with feedback calculations are studied on several problems from the VERA benchmark suit. Overall, the photon transport provides a better power map for depletion and feedback calculations, especially when strong gamma emitters are present. The fuel temperature, in general, is decreased when the photon transport is employed because fission gammas have been transported out of their birthplace. The new energy deposition model in STREAM with photon transport indicates that gamma energy accounts for about 9.7% of the total energy released in a fresh VERA5 2D core and rises to 10.5% at the End of Cycle. Around 2.5% of the total energy released will be deposited in the coolant in this 2D core. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

A Primer on Nuclear "Recoil" Data

For many years the ENDF-6 format has existed to contain evaluated nuclear data. One primary driver for the format and the data is neutron transport calculations. Evaluated data can be processed by a code like NJOY into either continuous energy form (ACE format) for Monte Carlo codes like MCNP or into multi-group form (NDI tables) for deterministic codes like Partisn. The reaction cross sections are found in the MF 3 section of the ENDF-6 format. If the reaction produces one or more neutrons as outputs, then secondary neutron data must also be given in MF 4,5 or 6 sections of the format. Energy and angular distributions for the output neutrons must be given in one of several available formats. The most general formats are found in MF 6. MF 4 is for angular distributions, MF 5 is for Energy distributions, and MF 6 contains both. In recent years, more interest has developed in the other output particles (i.e., the “recoil” particles) from neutron induced reactions. This has been driven by interest in charged particle transport and in more specialized partial kermas. ( A separate total neutron kerma has been available for a long time.) Partial kermas are a breakdown of the Kinetic Energy Released into the MAterial by output particle or by neutron reaction. The sum of the partial kermas should be equal to the total kerma on a group-wise basis. The ENDF-6 format is general enough for these new data requirements. Ideally, evaluated data would exist for every output particle (including the secondary neutrons) from every neutron-induced reaction. In the MF 6 format section, data for multi-particle outputs may be entered using the LAW =1 option. This option explicitly allows energy and angular output distributions for each particle produced in the reaction. Such information should preserve the balance between partial kermas and the total kerma at the groupwise level as well as the individual particle averaged energies. For simpler 2-body reactions, LAW = 2 is available. This allows the evaluation to specify only enough data to specify the 2-body reaction fully. Such a simplification does not exist for 3 (or more) body breakups. A simplified form of LAW 2, i.e., LAW 3, also exists to generate approximate recoil output distributions in the absence of full data. However, evaluated data does not generally exist at this fine granularity for reactions involving more than 2 output particles ( e.g., the 3-body break-up reaction). When the multi-body detailed output distribution data is not available, LAW=6 may be employed in NJOY. LAW 6 produces approximate output distributions for all the particles produced in the reaction. Just like in the 2-body case, the smaller particles will generally carry off more energy. At any particular energy, the balance between the sum of the partial kermas and the total kerma will not necessarily be preserved. However, the average energy of each output particle from a reaction is preserved.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Modeling and Measurement of Axial Gas Transport in Nuclear Fuels

This document is an American Nuclear Society (ANS) transaction for the ANS annual conference. This transaction provides an overview of the Modeling and Measurement of Axial Gas Transport Laboratory Directed Research and Development project including preliminary modeling results. The goal is to provide model for axial gas transport in the Bison Fuel Performance code to support industry's desire to achieve higher burnups in existing light water reactors. A brief overview of the experimental approach and associated modeling activities is presented.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Coupling MOOSE-Wrapped MPACT to BISON

As part of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, an effort is being made to leverage codes like Virtual Environment for Reactor Analysis (VERA) and Michigan Parallel Characteristics Transport (MPACT) by coupling them with other NEAMS codes. To facilitate coupling with other Multiphysics Objected Oriented Simulation Environment (MOOSE) applications, a MOOSE-wrapped MPACT app is created. The capability of this app, named Trogdor, is demonstrated by coupling it with another MOOSE app, BISON. Several single pin problems were run to test the Trogdor and BISON coupling.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Verification of the LOTUS code with C5G7 benchmark

In this study, the verification of the neutron transport code LOTUS through the well-known C5G7 benchmark is performed. LOTUS uses the current coupled collision probability method with the expansion of the flux by orthogonal polynomials (CCCPO) for the solution of the neutron transport equation. The expansion of the flux by orthogonal polynomials allows one to avoid discretization of the calculation regions and significantly decreases the simulation time. The results of the LOTUS calculations are compared with the results of the reference OpenMC Monte Carlo simulations. The results of the comparison demonstrate almost perfect agreement with Monte-Carlo for the second order of the flux expansion. Even though the results obtained in the presented study agree well with the results of reference Monte Carlo calculations, further investigations are necessary for a better understanding of the stability and limitations of the code.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Neutronics Modeling of the pulsed plasma rocket reactor using rattlesnake

In the pulsed plasma rocket (PPR) reactor, neutron induced fission processes are utilized to implement a series of pulsed micro-explosions of very high power and eject plasma as a propellant. More specifically, projectiles (bullets) composed of moderated uranium are sent through the chamber of an unmoderated uranium barrel. By inducing rapid fission within the bullets, a plasma can be generated with appropriate delivery of neutrons. This study aims to assess the neutronics performance of the PPR reactor during normal operations using the MOOSE-based Rattlesnake code through the evaluation of the impact of the movement of the fuel bullet and the rotation of control drums on the criticality of the system. The Monte Carlo (MC) code Serpent 2 was employed to generate material-based cross sections for use in Rattlesnake and the reference neutronics solution. Cubit was used to generate the mesh for the Rattlesnake model. Parametric studies were conducted to evaluate the best approaches for cross section and mesh generation to ensure accurate results from Rattlesnake. As part of the verification process, the eigenvalue results of the system at various fuel bullet positions were obtained using Rattlesnake and compared with the reference solutions. The acceptable differences show that the Rattlesnake model with appropriate cross section and mesh generation procedure is a sufficiently accurate approximation of the continuous energy (CE) MC model for micro-sized reactors like the PPR reactor. Next, the impact of the rotation of the control drums on the system criticality was evaluated using the verified model. It was shown that the high-fidelity simulation using the deterministic code Rattlesnake can produce sufficiently accurate results for the evaluation of the reactor's neutronics performance in different phases of the normal operation with acceptable computational cost.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

S/U Comparison Study with a Focus on USLs

Under a DOE Nuclear Criticality Safety Program (NCSP) task involving Analytical Methods, three Laboratories collaborated in a comparison of results obtained from Sensitivity/Uncertainty (S/U) packages relevant to validation of transport codes. The task involves Institut de Radioprotection et de Sûreté Nucléaire (IRSN), Los Alamos National Laboratory (LANL), and Oak Ridge National Laboratory (ORNL) comparing results of MORET 5/MACSENS V3.0, MCNP6.2/Whisper-1.1, and SCALE 6.2.3/TSUNAMI/USLSTATS respectively. All Monte Carlo transport code results utilize ENDF/B-VII.1. Four cases from the International Handbook of Evaluated Criticality Safety Benchmark Experiments (ICSBEP Handbook) were selected as application models: HEU-MET-FAST-013-001, HEU-SOL-THERM 001-008, PU-MET-FAST-022-001, and PU-SOL-THERM 001-001. Ultimately, comparison is made between Upper Subcritical Limits (USLs) obtained using each code package for each application case. Since differences exist in whether packages take into account margin of subcriticality (MOS), the USL may be computed using bias and bias uncertainty, also known as the calculational margin (CM) in ANSI/ANS 8.24. Application of portions of MOS to the USL for nuclear data uncertainty of and potential code margin is referred to as USL herein. In either case, additional MOS is considered for actual application cases.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

SCALE

SCALE v6.3.1 is a beta version built from existing SCALE code. This code provides modeling and simulation of nuclear reactor physics, criticality safety, radiation shielding, and associated spent fuel characterization for nuclear facilities and nuclear fuel transportation and storage package designs (please see attachment for more details).

Ade, Brian↗
DOE OSTI2020

Comparison of Modeled to Measured Spectra using MCNP and GADRAS to Benchmark and Contrast Modeling Limitations

The desire to improve nuclear material detection through portal monitors and other low resolution detectors has led to interest in benchmarking the performance of radiation transport codes. These codes can be used to generate a variety and quantity of spectra that may be cost prohibitive to measure directly. Particular characteristics of typical detection scenarios were isolated to compare the performance of radiation transport codes to laboratory experiments. These benchmark experiments were performed to validate simulations in a number of key areas. These experiments included high areal density configurations in both one- and three-dimensional configurations. In addition, off-angle measurements were conducted to investigate the impact of detector response assumptions on complex geometries, for example where an unknown source is not physically collocated with an apparent hot spot. Furthermore, by examining both simple shielding and backscatter configurations, isolation of the scatter emanating from the object and environmental background scatter contributions to the spectra could be elucidated. A series of experiments were performed in which the order of shielding layers and the thickness were evaluated to examine the ability of the simulations to address these variables. Experiments with increasing thickness of polyethylene around a neutron source also allows for examinations of the ability of the transport simulations to properly model n,γ reactions. Finally, depleted uranium (DU) measurements have more complicated spectra than the cobalt-60 and consequently allow for investigations into the ability of transport simulations to model self-shielding. In each case, the results of the experiments were compared to MCNP simulations to judge the performance of the code, as was done previously with GADRAS 18.7.9 simulations [7]. All of the benchmark measurements were taken with a liquid nitrogen cooled 140% high purity germanium (HPGe) detector with a bismuth collimator that had the front tin filter removed. The measurement location and detector configuration were subsequently used for all experimental configurations. The remainder of this report details the experimental measurements, Section 2 and the MCNP model, Section 3. Comparisons between the experimental results and MCNP-generated spectra are performed, Section 4. Section 5 compares the results obtained with MCNP with the previously reported results obtained with GADRAS [2]. Finally, Section 6 presents the conclusions.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2021

Validation of Light Water Reactor Ex-Core Calculations with VERA

The Consortium for Advanced Simulation of Light Water Reactors (CASL) Virtual Environment for Reactor Applications (VERA) is a reactor simulation software. It offers unique capabilities by combining high-fidelity in-core radiation transport with temperature feedback by using MPACT (a deterministic neutron transport code) and COBRA-TF (a thermal-hydraulic code) with follow-on, fixed-source transport calculations using the Shift Monte Carlo code to calculate ex-core quantities of interest. In these coupled calculations, MPACT provides Shift with the fission source for follow-on ex-core calculations. These ex-core simulations can be set up to calculate detector responses, as well as the flux and fluence in ex-core regions of interest, such as the reactor pressure vessel, nozzle, and irradiated capsules. A Watts Bar Nuclear Plant Unit 1 (WBN1) ex-core model was developed, as described in this paper, and this model was used to perform coupon calculations. The results for the coupon flux calculations show close agreement with the reference values for cycle 1 produced by the two-dimensional Discrete Ordinates Transport (DORT) code and presented in a BWXT Services Inc. report. However, differences in the results (10%) seen in cycles 2 and 3 and the reasons for these differences are discussed in this paper. The VERA WBN1 model was also used to perform a vessel fluence calculation for cycle 1. Additionally, a collaboration between CASL and Duke Energy led to the first code-to-code validation of VERA for reactor ex-core applications that used a model for the Shearon Harris reactor. Results from this collaboration show excellent agreement between VERA and the Monte Carlo N-Particle Transport Code for the detector response calculations. The work performed under this collaboration is also detailed in this paper.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗