Search NASA⌕ Search

SEARCH · Search NASA

Results for “Monte Carlo Code”

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 1,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Modeling of Pyroprocessing Hot Cell for Process Monitoring Evaluation

In this report, we examine studies undertaken to determine the applicability of applying measurements with a neutron counter within a pyroprocessing hot cell. The studies focus on the modeling of the hot cell and neutron counter in Monte Carlo N- Particle Code (MCNP) to determine expected radiation effects and applicability of neutron counters. The neutron counter of focus was the High Dose Neutron Detector (HDND) developed at Los Alamos National Laboratory (LANL). We studied the HDND to determine what would be the detector response under normal background conditions with operating process equipment as well as the background’s potential effect on counting applications involving a U/TRU ingot produced by the pyroprocessing electrorefiner. To perform this task, we were required to develop models of the pyroprocessing process cell and unit operations and generate source terms for various different times in the process. The data for the pyroprocessing materials came from SNL’s Safeguards and Security Process Model (SSPM) as part of MPACT’s advanced integration mission. This report will discuss the development of the pyroprocessing process cell and unit operations models in MCNP as well as various simulations to determine the applicability of the HDND for process monitoring efforts in the hot cell.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2020

Enhanced Visualization Capabilities in OpenMC: Report of Activities under a DOE GAIN Voucher

This work was performed under a Cooperative Research and Development Agreement (CRADA) through the U.S. Department of Energy’s Gateway for Accelerated Innovations in Nuclear (GAIN) program in partnership with Oklo Inc., the industry participant. The broad goal of this project was to improve the visualization capabilities of OpenMC related to model development and simulation results. OpenMC is a Monte Carlo particle transport code that is developed and maintained by a broad community of contributors, of which Argonne National Laboratory plays a leading role.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Compendium of Material Composition Data for Radiation Transport Modeling

In 2011, Pacific Northwest National Laboratory (PNNL) produced a document known as the Materials Compendium, or Compendium of Material Composition Data for Radiation Transport Modeling, PNNL 15870, Rev. 1, that contains material information useable for modeling purposes for properties of 372 materials. This information is used in several modeling programs used by the radiological/nuclear community, though it is primarily tailored for the Monte-Carlo-N-Particle code produced by Los Alamos National Laboratory. This new document Revision 2 includes a complete review and update of all materials data and references, addressing discrepancies and changes in materials data or references that have occurred since the first revision, an additional 40 materials have been added, primarily newer detector materials developed since the last revision, and isotopic specificity.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2023

Cost-optimized Automated Variance Reduction for Highly Angle-dependent Radiation Transport Analyses [Dissertation]

Monte Carlo variance-reduction techniques that directly bias particle direction have historically received limited attention despite being useful for many radiation transport applications such as those in which radiation travels through large regions with a low probability of colliding. One such technique is known as DXTRAN (short for deterministic transport) in the MCNP Monte Carlo radiation transport code. Until now, effectively applying DXTRAN in calculations is largely based on empirical observations of computational performance. Optimal DXTRAN parameters are identified through manual iteration. This work develops new mathematical descriptions of the DXTRAN variance-reduction process and demonstrates a new automated variance-reduction method that applies these mathematical formulations to determine the optimal application of DXTRAN in a given problem.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Modeling of the Advanced Test Reactor Using OpenMC, Cubit and Griffin

In the pursuit of the ability to perform multiphysics simulations of the Advanced Test Reactor, it is crucial to have a fast and highly accurate deterministic model. To achieve this, a contemporary two-step method is employed. The first step involves generating homogenized cross sections using OpenMC, a cutting-edge Monte Carlo neutron transport code. OpenMC offers excellent modular capabilities, allowing for easy component integration and flexibility in incorporating new designs into the model. The second step involves deterministic transport calculations, which are performed using Griffin, a reactor multiphysics application based on the Multiphysics Object-Oriented Simulation Environment. To ensure the accurate spatial resolution and assignment of material cross sections, a Cubit-generated mesh for the Advanced Test Reactor is utilized as an intermediate step between the OpenMC and Griffin models; Griffin utilizes the mesh for its finite element solution, while OpenMC material IDs are written to the mesh file to be used in Griffin material assignments. Additionally, a Python-based script converts the cross sections generated by OpenMC into the ISOXML format required by Griffin. Preliminary comparisons indicate good agreement between the neutron multiplication factors obtained from the standalone OpenMC model and the Griffin model, with differences of less than 50 pcm in the two-dimensional geometry configuration. However, in three-dimensional calculations, an unacceptably large error is found in the Griffin solution. Future work is planned to resolve this discrepancy.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Improved Fast Reactor Capability of Griffin in FY23

Griffin is a MOOSE based reactor multiphysics analysis application jointly developed by Idaho National Laboratory and Argonne National Laboratory under the DOE-NE NEAMS program. In FY23, we enhanced capabilities required for fast reactor analysis. This effort included primarily updating the cross-section generation workflow using MC2-3 for various reactor configurations, such as homogeneous, duct-heterogeneous, ring-heterogeneous, and fully-heterogeneous geometries. In addition, we initiated the implementation of a multi-cycle depletion and shuffling capability. To support fast reactor simulation capabilities, we significantly improved the performance of the DFEM-SN-based R-Z transport solver to efficiently solve ultrafine group (over 1000 groups) transport problems. Additionally, the performance of HFEM-PN was improved by introducing red-black iteration, the cmfd acceleration technique, and various optimizations. We also completed the pin power reconstruction capability to support multiphysics simulations while identifying and addressing issues associated with SPH equivalence parameter approach. These enhanced capabilities for fast reactor core simulations, specially HFEM-PN and pin power reconstruction features, were applied to benchmark problems involving ABTR and ABR-1000. These applications showcased excellent agreement with Monte Carlo and other code solutions in terms of eigenvalue, control rod worth, and assembly and pin powers.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Radiation Dose Modeling for Niowave’s Accelerator Driven Uranium Target Assembly 3

Molybdenum-99 is a high-value radionuclide commonly used for medical purposes within the United States. The National Nuclear Security Administration (NNSA) seeks to reliably produce the radioisotope 99 Mo without the use of highly enriched uranium. NNSA’s Office of Material Management and Minimization (M3) provides funding and government laboratory expertise to private companies to expedite the production process domestically and currently funds designs that use low-enriched uranium or other 99 Mo production pathways. Several production designs are being explored across the industry, including uranium fission and photonuclear conversion of 100 Mo targets. Niowave Inc. seeks to produce 99 Mo via a high-energy electron accelerator that strikes a lead-bismuth eutectic target that ultimately produces a consistent neutron flux. The neutron flux then interacts in a subcritical reactor core configuration to produce fission in low-enriched or natural uranium targets. These fissionable targets are then processed to extract 99 Mo. The purpose of this work is to estimate the neutron and photon dose response across Niowave’s proposed facility for worker safety during operation. Owing to the size of the proposed Niowave facility and necessary shielding, unbiased Monte Carlo radiation transport is impractical, and variance reduction methods are required. This work focuses on the weight window variance reduction method to produce high confidence dose response results within a Monte Carlo radiation transport code. Specifically, an adjoint-informed weight window methodology was created to improve the dose response estimates for accelerator-driven subcritical reactor designs. This adjoint-informed methodology was implemented for Niowave’s proposed design and improved dose results at far-field locations across the facility. Acceptable dose rate contours for the proposed facility were generated across the facility and are presented in this work.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2025

Parallel Programming in MCNP6

Monte Carlo N-Particle (MCNP)1 is a general-purpose Monte Carlo particle transport code developed by Los Alamos National Laboratory (LANL). To efficiently handle long simulations, MCNP version 6 (MCNP6) supports parallel execution using two primary programming models: • Shared-memory task-based threading using OpenMP (Open Multi-Processing), and • Distributed-memory calculations using MPI (Message Passing Interface). The OpenMP and MPI programming models enable MCNP6 to scale from desktop systems to high-performance computing (HPC) clusters, allowing users to run MCNP in one of three parallel modes: • OpenMP-only, • MPI-only, and • Hybrid (MPI + OpenMP). The choice of parallelization mode depends on the underlying computer architecture and the characteristics of the simulation problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Shift Validation Plan for CASL

This document outlines the benchmarks that will be used for Shift validation for CASL. This plan is divided into five main types of problems, with each building upon the previous: small criticality benchmark problems, including small reactor configurations and critical assemblies; full-core reactor benchmarks for light water reactors; fixed source coupled neutron-photon dosimetry benchmarks; depletion/burnup benchmarks ranging from small pincell problems to full-core reactor depletion; and full-core reactor performance benchmarks. Some of these benchmarks have measured data for comparison and some only have simulated Monte Carlo radiation transport code results for comparison. Upon completion of this validation plan, we will have confidence in Shift to provide reference results for CASL benchmarking.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Modeling The Nucleosynthetic Imprint of Stellar Merger Phenomena - Final Technical Report

This DOE CAREER project developed a comprehensive, multi-physics framework for modeling stellar mergers and their observable consequences. The work successfully integrated analytical models, three-dimensional hydrodynamic simulations, stellar evolution calculations, nucleosynthesis, and radiation transport into a unified pipeline. Key scientific advances include demonstrating stellar mergers as a unifying explanation for systems such as Betelgeuse and R Coronae Borealis stars, quantifying merger-driven chemical signatures, and establishing the connection between merger physics and circumstellar environments. A major outcome of the project is the development of SuperLite, an open-source Monte Carlo radiation transport code that enables the generation of synthetic spectra for astrophysical transients. This work bridges dynamical, thermal, and radiative timescales and provides a powerful, DOE-relevant computational framework for interpreting observations of supernovae and related transient phenomena, while also contributing to workforce development through the training of postdoctoral researchers and graduate students.

Chatzopoulos, Emmanouil [Louisiana State Universit↗
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 OSTI2018

Spinning Dust Emission from Circumstellar Disks and Its Role In Excess Microwave Emission

Electric dipole emission from rapidly spinning polycyclic aromatic hydrocarbons (PAHs) is widely believed to be an origin of anomalous microwave emission (AME), but recently it has encountered a setback owing to the noncorrelation of AME with PAH abundance seen in a full-sky analysis. Microwave observations for specific regions with well-constrained PAH features would be crucial to test the spinning dust hypothesis. In this paper, we present physical modeling of microwave emission from spinning PAHs from protoplanetary disks (PPDs) around Herbig Ae/Be stars and T Tauri stars where PAH features are well observed. Guided by the presence of 10 μm silicate features in some PPDs, we also model microwave emission from spinning nanosilicates. Thermal emission from big dust grains is computed using the Monte Carlo radiative transfer code (radmc-3d). Our numerical results reflect that microwave emission from either spinning PAHs or spinning nanosilicates dominates over thermal dust at frequencies ν < 60 GHz, even in the presence of significant grain growth. Lastly, we attempt to fit millimeter–centimeter observational data with both thermal dust and spinning dust for several disks around Herbig Ae/Be stars that exhibit PAH features and find that spinning dust can successfully reproduce the observed excess microwave emission (EME). Future radio observations with ngVLA, SKA, and ALMA Band 1 would be valuable for elucidating the origin of EME and potentially open a new window for probing nanoparticles in circumstellar disks.

79 ASTRONOMY AND ASTROPHYSICS↗
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↗
NASA NTRS1984

Characterization of Surface Reflectance Variation Effects on Remote Sensing

The use of Monte Carlo radiative transfer codes to simulate the effects on remote sensing in visible and infrared wavelengths of variables which affect classification is examined. These variables include detector viewing angle, atmospheric aerosol size distribution, aerosol vertical and horizontal distribution (e.g., finite clouds), the form of the bidirectional ground reflectance function, and horizontal variability of reflectance type and reflectivity (albedo). These simulations are used to characterize the sensitivity of observables (intensity and polarization) to variations in the underlying physical parameters both to improve algorithms for the removal of atmospheric effects and to identify techniques which can improve classification accuracy. It was necessary to revise and validate the simulation codes (CTRANS, ARTRAN, and the Mie scattering code) to improve efficiency and accommodate a new operational environment, and to build the basic software tools for acquisition and off-line manipulation of simulation results. Initial calculations compare cases in which increasing amounts of aerosol are shifted into the stratosphere, maintaining a constant optical depth. In the case of moderate aerosol optical depth, the effect on the spread function is to scale it linearly as would be expected from a single scattering model. Varying the viewing angle appears to provide the same qualitative effect as modifying the vertical optical depth (for Lambertian ground reflectance).

Pearce, W. A.↗
NASA NTRS1985

The Characterization of Surface Variegation Effects on Remote Sensing

Improvements in remote sensing capabilities hinge very directly upon attaining an understanding of the physical processes contributing to the measurements. In order to devise new measurement strategies and to learn better techniques for processing remotely gathered data, it is necessary to understand and to characterize the complex radiative interactions of the atmosphere-surface system. In particular, it is important to understand the role of atmospheric structure, ground reflectance inhomogeneity and ground bidirectional reflectance type. The goals, then, are to model, analyze, and parameterize the observable effects of three dimensional atmospheric structure and composition and two dimensional variations in ground albedo and bidirectional reflectance. To achieve these goals, a Monte Carlo radiative transfer code is employed to model and analyze the effects of many of the complications which are present in nature.

Pearce, W. A.↗