Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo 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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

ParaStell: parametric modeling and neutronics support for stellarator fusion power plants

The three-dimensional variation inherent to stellarator geometries and fusion sources motivates three-dimensional modeling to obtain accurate results from computational modeling in support of design and analysis of first wall, blanket, and shield (FWBS) systems. Manually constructing stellarator fusion power plant geometries in computer-aided design (CAD) and defining the corresponding fusion source can be cumbersome and challenging. The open-source parametric modeling toolset ParaStell has been developed to automate construction of such geometries in low-fidelity. Low-fidelity modeling is useful during the conceptual phase of engineering design as a means of rapidly exploring the design space of a given device. The modeling capability of ParaStell includes in-vessel components and magnets, for any given stellarator configuration, using a parametric definition and plasma equilibrium data. Furthermore, the toolset automates the generation of detailed, tetrahedral neutron source definitions and DAGMC geometries for use in neutronics modeling. ParaStell assists rapid design iteration, parametric study, and design optimization of stellarator fusion cores. As a demonstration of the design iteration capability, the effect of the three-dimensional parameter space on tritium breeding and magnet shielding is investigated, using the WISTELL-D configuration as a design basis. Blanket and shield thicknesses are varied in three dimensions, using the space available between the plasma edge and magnet coils as a constraint. The corresponding effects on tritium breeding ratio and magnet heating are tallied using the open-source Monte Carlo particle transport code OpenMC. The inclusion of additional and higher-fidelity modeling capabilities is planned for ParaStell’s future, as well as its implementation in machine-driven optimization.

Moreno, Connor↗
DOE OSTI2021

A Broad Grid of 2D Kilonova Emission Models

Depending upon the properties of their compact remnants and the physics included in the models, simulations of neutron star mergers can produce a broad range of ejecta properties. The characteristics of this ejecta, in turn, define the kilonova emission. To explore the effect of ejecta properties, we present a grid of 2-component 2D axisymmetric kilonova simulations that vary mass, velocity, morphology, and composition. The masses and velocities of each component vary, respectively, from 0.001 to 0.1 M ⊙ and 0.05 to 0.3c, covering much of the range of results from the neutron star merger literature. The set of 900 models is constrained to have a toroidal low electron fraction (Y e ) ejecta with a robust r-process composition and either a spherical or lobed high-Y e ejecta with two possible compositions. We simulate these models with the Monte Carlo radiative transfer code SuperNu using a full suite of lanthanide and 4th row element opacities. We examine the trends of these models with parameter variation, show how it can be used with statistical tools, and compare the model light curves and spectra to those of AT2017gfo, the electromagnetic counterpart of GW170817.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

A Comparison of Three Neodymium Atomic Data Sets for Kilonova Modeling

We examine the impact of input neodymium (Nd) atomic data on the light curves and spectra of kilonovae (KNe), probing the sensitivity of kilonova observables to the atomic physics of this important lanthanide element. We use the SuperNu Monte Carlo radiative transfer code, simulating a simple semianalytic 1D kilonova (KN) with a pure Nd atmosphere, fixing the radiative transfer method while using input atomic data generated by three different codes: the LANL suite of atomic physics codes, HULLAC, and Autostructure. We see that the choice of atomic data significantly shapes the resulting light curves and spectra. Peak bolometric luminosities differ by a ratio of nearly 1.5 between HULLAC/Autostructure and LANL data sets. Moreover, we observe significant near- to mid-IR differences in the structure of the spectra. We specifically attribute these differences to the choice of atomic data for neutral Nd I. Many of the results here have been adapted from a presentation at “Radiative Transfer and Atomic Physics of Kilonovae” in Stockholm, 2023. We additionally present a LANL data set with energies calibrated to available values in the NIST Atomic Spectra Database, and demonstrate that this calibration also significantly affects IR spectral structure at late time. The substantial differences in KN observables that arise from tuning the atomic data of just one lanthanide element highlight the special attention that must be paid to atomic physics uncertainties when modeling KNe, from AT2017gfo to beyond.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Experimental Steps toward a Density Law for Chlorine-Crediting Criticality Models of Aqueous Plutonium solutions

Criticality accident prevention is an essential safety consideration for all operations involving fissionable and fissile materials. Predictive criticality calculations using advanced neutron transport codes, such as Monte Carlo N-Particle Transport (MCNP), are invaluable tools for designing and implementing operational limits. While indispensable, these tools are limited by the quality and accuracy of the inputs that the user provides to define the modeled system. Parameters such as atomic composition, shape, and material density must be accurately defined to obtain a meaningful result. In the case of material property details for a plutonium model, accurate material characterization data is sometimes sparse. Thus limitations in calculation accuracy can be reduced by improving our knowledge of material properties in the targeted fissile systems. One of the biggest remaining challenges to accurately defining fissile systems is the description of aqueous fissile solutions. Even simple properties, such as density, are not well-known for fissile solutions relevant to nuclear energy and security. Described here are initial efforts undertaken to improve criticality calculation inputs for fissile plutonium chloride solutions in water. This effort is focused on experimentally determining accurate solution characteristics for ternary plutonium chloride/hydrochloric acid/water systems, by measurement of water activity and solution density. The effect of inputting experimental densities for these solutions into MCNP criticality calculations is compared to the traditional approach of modeling an idealized (and fictitious) plutonium-water mixture. Expansion of these efforts to a working density law for aqueous plutonium chloride solutions is also discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Overview on Shielding Analyses for the VENUS Instrument at SNS

VENUS, a world‐class versatile neutron imaging instrument, is under construction and is expected to be completed and ready to start commissioning in 2023. The range of cold to epithermal neutrons at SNS will give users of VENUS access to novel imaging methods, as well as to significantly improved existing methods. The instrument is being built on beam line 10 at Spallation Neutron Source (SNS) First Target Station (FTS) facing a decoupled poisoned hydrogen moderator. Instrument design, which includes optics, Front-End components and instrument enclosure components and shape started over 10 years ago and during this time had significant changes. All changes were supported by neutronics analyses to provide adequate shielding for both Front-End and instrument enclosure.Final optics design will provide the field of view (FOV) at the detector position (where the image is formed) to be as high as 0.20 by 0.20 m. Both the beam, which contains a large fraction of high-energy neutrons, and the desire for a large footprint on the detector are challenges for shielding design and budget, because the driving cost for the instrument is the beam line and enclosure shielding. For cost reason VENUS baseline design is optimized with respect to both the instrument cave footprint and its wall thickness, and Front-End shielding is tailored along the beam. Analyses are performed with the Monte Carlo particle transport code MCNPX version 2.7.0 to make choices on materials, thicknesses and configurations of the shielding. Numerous calculations are performed to meet space and coast constrain and to satisfy instrument physics needs and to comply with radiation protection requirements.

Gallmeier, Franz X.↗
DOE OSTI2023

Local Power Impact Experiment Design for a New Fuel Type for use in the Advanced Test Reactor

The Advanced Test Reactor (ATR), and complimentary zero-power ATR Critical (ATRC) reactor, located at Idaho National Labs (INL), are undergoing conversion from Highly Enriched Uranium (HEU) to Low Enriched Uranium (LEU). Both have a variety of testing locations that can receive large variations in flux due to its unique serpentine design, consisting of five lobes (see Figure 1). Initial criticality and power distribution throughout the core are controlled by core-external outer shim control cylinders (OSCCs). Distinct test loops allow for testing at specific temperatures, pressures, and irradiation conditions. The ATR is one of the key nuclear engineering research and testing facilities within the DOE National Laboratory Complex, and the ATRC supports its operation [1]. Currently, the Office of Material Management and Minimization (M3) within the National Nuclear Security Administration of the DOE is working to convert the remaining research reactors, including the ATR, from 93% HEU fuel to 19.75% LEU fuel (LEU) to support non-proliferation [2]. Extensive materials testing at INL and internationally has demonstrated that a high-density uranium molybdenum (U 10Mo) alloy can meet the performance requirements of the remaining high powered research reactors. However, there are many technical challenges to address before the conversion to LEU can be successful, including the accurate characterization of the reactor core physics with LEU fuel. Reactor physics safety evaluations currently use Monte Carlo for the 21st Century (MC21), a continuous-energy Monte Carlo radiation transport code [3]. Existing MC21 models of the ATR and ATRC cores have a validation basis for use in neutronics analyses with HEU fuel. The models are used to support safety analyses that include comparisons to the safety requirements for the reactors. However, the use of the LOWE element in the ATR and ATRC is not currently covered by the current model validation basis. To deploy the new fuel type, extensive computational reactor physics support is necessary to support the use of LOWE in the ATR and ATRC. Therefore, LOWE requires a rigorous validation basis, aligned with that of HEU fuel, that takes advantage of the existing software tools and processes currently used for the ATR and ATRC. The experiment to validate of the MC21 models for determining power, the Power Impact Validation Experiment, will consist of two flux runs in the ATRC, one with fully HEU loading and one with a single LOWE element. Both flux runs will be instrumented with 20 sets of azimuthal fission wires and 3 sets of axial fission wires, as shown in Figure 4. Standard flux run methodology will be used [4]. Power Impact Validation Experiment data will be compared against MC21 calculated data, both for absolute fission rate accuracy and to determine the relative change in fission rates between the two runs. The results of the Power Impact Validation Experiment and subsequent evaluations will provide the validation basis for MC21 for use with LOWE elements. Key features of the Power Impact Validation Experiment include: (1) Two flux runs to allow for LOWE perturbed measurements to be compared to already validated measurements taken from a full core of HEU fuel, (2) Optimization of instrumentation to balance analytical needs with practical considerations (e.g., limited time window to count beta particles from fission products), and (3) Standard ATRC core loading, including both driver positions and flux traps, to minimize cost while remaining representative of typical ATR core loading.

42 ENGINEERING↗
DOE OSTI2024

Overview of Shielding Analyses at Oak Ridge National Laboratory’s Spallation Neutron Source Second Target Station

The Neutronics Group of the Second Target Station project at Oak Ridge National Laboratory is responsible for all neutronics analyses related to design and construction. This paper provides an overview of four neutronics analyses performed for the Second Target Station project, which is representative of all the ongoing work but especially tasks related to shielding. These analyses extend from the proton accelerator through the target monolith and bunker to the end of a neutron beamline. Each example highlights the tools and methods the Neutronics Group uses for analysis. The primary tool is the Monte Carlo radiation transport code MCNP 6.2, but this is augmented by other codes that supplement the input of MCNP. The other codes specifically highlighted in this paper include Attila4MC, ADVANTG, and AARE.

Miller, Thomas↗
DOE OSTI2025

Investigation of the approach used in the unresolved resonance region

The typical practice used in Monte Carlo neutron transport codes in the unresolved resonance region (URR) is to take advantage of the probability table (PT) approach. Cross sections are sampled from PTs and used as needed, along with generated average cross sections. The PTs are generated based on cross-section calculations performed using the single-level Breit–Wigner approximation. Although the approach used in the URR seems plausible, a detailed examination was needed to understand the benchmark results obtained using specific tests.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

A Bottom-Up Cost Estimation Tool for Nuclear Microreactors

The rising interest in nuclear microreactors has highlighted the need for comprehensive technoeconomic assessments. However, the scarcity of publicly available designs and cost data has posed significant challenges. To address this issue, the Microreactor Optimization Using Simulation and Economics (MOUSE) tool is developed. MOUSE is a tool that integrates nuclear microreactor design with reactor economics. The design calculations encompass core simulations using the OpenMC Monte Carlo Particle Transport Code [romano2015], along with simplified balance of plant calculations. On the economic side, MOUSE provides detailed bottom-up cost estimates, calculating both the total capital cost and the levelized cost of energy for first-of-a-kind and nth-of-a-kind microreactors. The cost estimation correlations are developed using data from the MARVEL project and additional literature sources. MOUSE has released as an open-source tool on GitHub (MOUSE Tool). By combining design calculations with cost equations, MOUSE enables users to evaluate the impact of various technological consideration, advanced moderators, design changes, material/fuel changes, and geometry modifications—as well as economic parameters like interest rates and construction duration. This comprehensive framework can guide stakeholders towards technological solutions that enhance microreactor competitiveness. Additionally, powered by the WATTS toolkit [romano2022], MOUSE supports optimization studies, parametric analyses, and uncertainty calculations/propagation. Currently, preconceptual designs of three microreactor types are included in MOUSE: a liquid metal thermal microreactor (LTMR), gas cooled TRISO-fueled microreactor (GCMR) and heat-pipe TRISO fueled microreactor (HPMR). To showcase its ability, MOUSE was used to conduct detailed bottom-up cost estimates for the first of a kind (FOAK) and Nth of a kind (NOAK) of the following microreactors • A 20MWt LTMR that is built on the ongoing MARVEL demonstration at Idaho National Laboratory (INL) • A 15 MWt GCMR that was designed to be more representative of the typical commercial microreactor • A 7 MWt HPMR that was built on previous work (Choi 2024) The The reader should note that these three designs and corresponding cost estimates are examples to demonstrate the MOUSE capability. The designs are pre-conceptual, the reactor designs were not optimized, and the cost estimates were developed with incomplete information. Additionally, stakeholders might be interested in a variety of designs that may differ from the examples provided in this report. The MOUSE tool can also be used to study how design choices affect economics. To demonstrate its capability, MOUSE was used conduct parametric studies such as examining the economic impact of the reflector's material and thickness, the moderator's booster material and dimensions, fuel composition and enrichment, core size, and power level. Several insights were gained from these parametric studies.

Hanna, Botros↗
DOE OSTI2024

Celeritas: Accelerating Geant4 with GPUs

Celeritas [1] is a new Monte Carlo (MC) detector simulation code designed for computationally intensive applications (specifically, High Lumi- nosity Large Hadron Collider (HL-LHC) simulation) on high-performance heterogeneous architectures. In the past two years Celeritas has advanced from prototyping a GPU-based single physics model in infinite medium to implementing a full set of electromagnetic (EM) physics processes in complex geometries. The current release of Celeritas, version 0.3, has incorporated full device-based navigation, an event loop in the presence of magnetic fields, and detector hit scoring. New functionality incorporates a scheduler to offload electromagnetic physics to the GPU within a Geant4-driven simulation, enabling integration of Celeritas into high energy physics (HEP) experimental frameworks such as CMSSW. On the Summit supercomputer, Celeritas performs EM physics between 6 and 32 faster using the machine’s Nvidia GPUs compared to using only CPUs. When running a multithreaded Geant4 ATLAS test beam application with full hadronic physics, using Celeritas to accelerate the EM physics results in an overall simulation speedup of 1.8–2.3× on GPU and 1.2× on CPU.

Johnson, Seth R.↗
DOE OSTI2021

The PENELOPE Physics Models and Transport Mechanics. Implementation into Geant4

A translation of the PENELOPE physics subroutines to C++, designed as an extension of the GEANT4 toolkit, is presented. The Fortran code system PENELOPE performs Monte Carlo simulation of coupled electron-photon transport in arbitrary materials for a wide energy range, nominally from 50 eV up to 1 GeV. PENELOPE implements the most reliable interaction models that are currently available, limited only by the required generality of the code. In addition, the transport of electrons and positrons is simulated by means of an elaborate class II scheme in which hard interactions (involving deflection angles or energy transfers larger than pre-defined cutoffs) are simulated from the associated restricted differential cross sections. After a brief description of the interaction models adopted for photons and electrons/positrons, we describe the details of the class-II algorithm used for tracking electrons and positrons. The C++ classes are adapted to the specific code structure of GEANT4. They provide a complete description of the interactions and transport mechanics of electrons/positrons and photons in arbitrary materials, which can be activated from the G4ProcessManager to produce simulation results equivalent to those from the original PENELOPE programs. The combined code, named PENG4, benefits from the multi-threading capabilities and advanced geometry and statistical tools of GEANT4.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Cross Section Generation Capability in Griffin

The Griffin code is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based reactor multiphysics analysis application jointly developed by Idaho National Laboratory and Argonne National Laboratory. The code includes a variety of steady-state solvers for fixed-source, k-eigenvalue, adjoint, and subcritical multiplication, as well as transient solvers for point-kinetics, improved quasi-static, and spatial dynamics. The code reads multigroup cross sections in the ISOXML format generated from external deterministic or Monte Carlo cross section generation codes. The implementation of the cross section generation capability in Griffin was initiated last year by plugging in the cross section application programming interface (CSAPI) and reviewing the methodologies for treating particulate fuels. The focus this year was on improving the CSAPI integration and implementing advanced self-shielding methods for applications to advanced reactor problems with TRISO fuels. First, the process for cross section library generation was updated to accurately and rigorously produce isotopic cross section data. Second, the equivalent Dancoff factor cell method performing slowing down calculations on the fly for the resonance treatment was implemented in CSAPI to improve the accuracy of effective multigroup cross sections in the resonance energy range. Third, the iterative local spatial self-shielding method was implemented under the calculation framework of the equivalent Dancoff factor cell method to accurately deal with the double heterogeneity effect of particulate fuel. The updated CSAPI with the advanced self-shielding methods, together with the cross section libraries generated based on the improved process, were tested for pin-cell, unit-cell, and fuel assembly problems with various resonance self-shielding conditions based on very high temperature reactor, high temperature test reactor, and Empire benchmark cores, indicating that the updated CSAPI in Griffin is able to produce multigroup cross sections accurately and efficiently. We also showed that the methodology worked well for pebble bed fuel from HTR-10, but the capability still needs to be fully integrated into CSAPI. In the future, further benchmark tests will be performed for various thermal reactor core problems, including particulate fuel-based pebble bed reactors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Cross Section Generation Capability in Griffin

The Griffin code is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based reactor multiphysics analysis application jointly developed by Idaho National Laboratory and Argonne National Laboratory. The code includes a variety of steady-state solvers for fixed-source, k-eigenvalue, adjoint, and subcritical multiplication, as well as transient solvers for point-kinetics, improved quasi-static, and spatial dynamics. The code reads multigroup cross sections in the ISOXML format generated from external deterministic or Monte Carlo cross section generation codes. The implementation of the cross section generation capability in Griffin was initiated last year by plugging in the cross section application programming interface (CSAPI) and reviewing the methodologies for treating particulate fuels. The focus this year was on improving the CSAPI integration and implementing advanced self-shielding methods for applications to advanced reactor problems with TRISO fuels. First, the process for cross section library generation was updated to accurately and rigorously produce isotopic cross section data. Second, the on-the-fly slowing down method for the resonance treatment was implemented in CSAPI to improve the accuracy of effective multigroup cross sections in the resonance energy range. Among various on-the-fly slowing down methods, the equivalent Dancoff factor cell method was employed. Third, the iterative local spatial self-shielding method was implemented under the calculation framework of the equivalent Dancoff factor cell method to accurately deal with the double heterogeneity effect of particulate fuel. The updated CSAPI with the advanced self-shielding methods, together with the cross section libraries generated based on the improved process, were tested for the very high temperature reactor (VHTR), high temperature test reactor (HTTR), and Empire benchmark problems with various resonance self-shielding conditions, indicating that the updated CSAPI in Griffin is able to produce multigroup cross sections accurately and efficiently. We also show that the methodology works well for pebble bed fuel from HTR-10, but the capability still needs to be fully integrated into CSAPI. In the future, further benchmark tests will be performed for various thermal reactor core problems, including particulate fuel-based pebble bed reactors.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

MCNP ® Code Version 6.3.1 Theory & User Manual

This document acts as a repository of knowledge for the Monte Carlo N-Particle (MCNP) transport computer code. It is maintained alongside the source code and attempts to introduce new users and re-familiarize experienced users with the theory and practices of using the MCNP code for the wide range of particle transport analyses that it is appropriate for. The latest version of the MCNP code, version 6.3.1, provides the Monte Carlo particle transport community with the latest feature developments and bug fixes in the MCNP code. The MCNP code version 6.0 and later is also known as the MCNP6 code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

MontePy: a Python library for reading, editing, and writing MCNP input files.

The Monte Carlo N-Particle (MCNP) radiation transport code is a highly capable and accurate code with a long legacy. MCNP uses the Monte Carlo simulation process to simulate the path of particles (e.g., neutrons, photons, charged particles, etc.), and their interaction with materials. It is widely used in nuclear engineering, high-energy physics, and other fields. Its origins in the mid-twentieth century predate many modern software conventions. MCNP users provide an input file to MCNP, which it then uses to create an internal representation of the simulation problem. These input files originally had to be stored as punchcard decks, and the user manual still uses the terminology of cards and decks, despite moving beyond punchcards. MCNP predates nearly all modern human readable markup or data serialization languages, such as the extensible Markup Language (XML), the Standard Generalized Markup Language (SGML), YAML (YAML Ain’t Markup Language), and Javascript Object Notation (JSON). Due to this, MCNP uses an entirely custom defined syntax language for its input, making off-the-shelf libraries for XML, YAML, and JSON impossible to use for scripting various operations on MCNP input files (Kulesza et al., 2022).

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2022

MCNP ® Code Version 6.3.0 Theory & User Manual

This document acts as a repository of knowledge for the Monte Carlo N-Particle (MCNP) transport computer code. It is maintained alongside the source code and attempts to introduce new users and re-familiarize experienced users with the theory and practices of using the MCNP code for the wide range of particle transport analyses that it is appropriate for. The latest version of the MCNP code, version 6.3.0, provides the Monte Carlo particle transport community with the latest feature developments and bug fixes in the MCNP code. The MCNP code version 6.0 and later is also known as the MCNP6 code.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2026

MCNP® Code Version 6.3.2 Theory & User Manual (Rev. 1)

This document acts as a repository of knowledge for the Monte Carlo N-Particle (MCNP) transport computer code. It is maintained alongside the source code and attempts to introduce new users and re-familiarize experienced users with the theory and practices of using the MCNP code for the wide range of particle transport analyses that it is appropriate for. The latest version of the MCNP code, version 6.3.2, provides the Monte Carlo particle transport community with the latest feature developments and bug fixes in the MCNP code. The MCNP code version 6.0 and later is also known as the MCNP6 code.

42 ENGINEERING↗
DOE OSTI2023

Parallel simulation via SPPARKS of on-lattice kinetic and Metropolis Monte Carlo models for materials processing

Abstract SPPARKS is an open-source parallel simulation code for developing and running various kinds of on-lattice Monte Carlo models at the atomic or meso scales. It can be used to study the properties of solid-state materials as well as model their dynamic evolution during processing. The modular nature of the code allows new models and diagnostic computations to be added without modification to its core functionality, including its parallel algorithms. A variety of models for microstructural evolution (grain growth), solid-state diffusion, thin film deposition, and additive manufacturing (AM) processes are included in the code. SPPARKS can also be used to implement grid-based algorithms such as phase field or cellular automata models, to run either in tandem with a Monte Carlo method or independently. For very large systems such as AM applications, the Stitch I/O library is included, which enables only a small portion of a huge system to be resident in memory. In this paper we describe SPPARKS and its parallel algorithms and performance, explain how new Monte Carlo models can be added, and highlight a variety of applications which have been developed within the code.

36 MATERIALS SCIENCE↗