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.

SourceResearch connectionFrom yearThrough year

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

At least 1,189 records · Page 66

All sources NASA NTRS DOE OSTI

DOE OSTI2012

MPACT Fast Neutron Multiplicity System Design Concepts

This report documents work performed by Idaho National Laboratory and the University of Michigan in fiscal year (FY) 2012 to examine design parameters related to the use of fast-neutron multiplicity counting for assaying plutonium for materials protection, accountancy, and control purposes. This project seeks to develop a new type of neutron-measurement-based plutonium assay instrument suited for assaying advanced fuel cycle materials. Some current-concept advanced fuels contain high concentrations of plutonium; some of these concept fuels also contain other fissionable actinides besides plutonium. Because of these attributes the neutron emission rates of these new fuels may be much higher, and more difficult to interpret, than measurements made of plutonium-only materials. Fast neutron multiplicity analysis is one approach for assaying these advanced nuclear fuels. Studies have been performed to assess the conceptual performance capabilities of a fast-neutron multiplicity counter for assaying plutonium. Comparisons have been made to evaluate the potential improvements and benefits of fast-neutron multiplicity analyses versus traditional thermal-neutron counting systems. Fast-neutron instrumentation, using for example an array of liquid scintillators such as EJ-309, have the potential to either a) significantly reduce assay measurement times versus traditional approaches, for comparable measurement precision values, b) significantly improve assay precision values, for measurement durations comparable to current-generation technology, or c) moderating improve both measurement precision and measurement durations versus current-generation technology. Using the MCNPX-PoliMi Monte Carlo simulation code, studies have been performed to assess the doubles-detection efficiency for a variety of counter layouts of cylindrical liquid scintillator detector cells over one, two, and three rows. Ignoring other considerations, the best detector design is the one with the most detecting volume. However, operational limitations guide a) the maximum acceptable size of each detector cell (due to PSD performance and maximum-acceptable per-channel data throughput rates, limited by pulse pile-up and the processing rate of the electronics components of the system) and b) the affordability of a system due to the number of total channels of data to be collected and processed. As a first estimate, it appears that a system comprised of two rows of detectors 5" Ø $\times$ 3" would yield a working prototype system with excellent performance capabilities for assaying Pu-containing items and capable of handling high signal rates likely when measuring items with Pu and other actinides. However, it is still likely that gamma-ray shielding will be needed to reduce the total signal rate in the detectors. As a first step prior to working with these larger-sized detectors, it may be practical to perform scoping studies using small detectors, such as already-on-hand 3" Ø $\times$ 3" detectors.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2014

Stochastic Particle Real Time Analyzer (SPARTA) Validation and Verification Suite

This report presents the test cases used to verify, validate and demonstrate the features and capabilities of the first release of the 3D Direct Simulation Monte Carlo (DSMC) code SPARTA (Stochastic Real Time Particle Analyzer). The test cases included in this report exercise the most critical capabilities of the code like the accurate representation of physical phenomena (molecular advection and collisions, energy conservation, etc.) and implementation of numerical methods (grid adaptation, load balancing, etc.). Several test cases of simple flow examples are shown to demonstrate that the code can reproduce phenomena predicted by analytical solutions and theory. A number of additional test cases are presented to illustrate the ability of SPARTA to model flow around complicated shapes. In these cases, the results are compared to other well-established codes or theoretical predictions. This compilation of test cases is not exhaustive, and it is anticipated that more cases will be added in the future.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2014

Stochastic Particle Real Time Analyzer (SPARTA) Validation and Verification Suite

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2016

Semi-Analytic Benchmarks for MCNP6

This series of slides shows results for benchmarks for code verifications in MCNP6. For additional related information see report LA-UR--16-28568.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2017

A Pulsed Sphere Tutorial

Here I attempt to explain what physically happens when we pulse an object with neutrons, specifically what we expect the time dependent behavior of the neutron population to look like. Emphasis is on the time dependent emission of both prompt and delayed neutrons. I also describe how the TART Monte Carlo transport code models this situation; see the appendix for a complete description of the model used by TART. I will also show that, as we expect, MCNP and MERCURY, produce similar results using the same delayed neutron model (again, see the appendix).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2017

Modeling blur in various detector geometries relevant to DARHT

Monte Carlo transport codes have been used to model the detector blur and energy deposition in various detector geometries for MeV radiography which are applicable to DARHT. Segmented scintillator detectors, where low Z scintillators combined with a high-Z metal matrix, are designed such that the resolution increases with increasing metal fraction. The combination of various types of metal intensification screens and storage phosphor imaging plates has also been studied. A storage phosphor coated directly onto a metal intensification screen has superior performance over a commercial plate. Stacks of storage phosphor plates and tantalum intensification screens show an increase in energy deposited and detective quantum efficiency with increasing plate number, at the expense of resolution which degrades with plate number. Excellent agreement between simulation and experimental modulation transfer functions for select detector geometries was found.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2014

LDRD Report : Analysis of Defect Clustering in Semiconductors using Kinetic Monte Carlo Methods

The transient degradation of semiconductor device performance under irradiation has long been an issue of concern. Neutron irradiation can instigate the formation of quasi-stable defect structures, thereby introducing new energy levels into the bandgap that alter carrier lifetimes and give rise to such phenomena as gain degradation in bipolar junction transistors. Typically, the initial defect formation phase is followed by a recovery phase in which defect-defect or defect-dopant interactions modify the characteristics of the damaged structure. A kinetic Monte Carlo (KMC) code has been developed to model both thermal and carrier injection annealing of initial defect structures in semiconductor materials. Following the development of a set of verification tests, the code is employed to investigate annealing in electron-irradiated, p-type silicon as well as the recovery of base current in silicon transistors bombarded with neutrons at the Los Alamos LANSCE "Blue Room" facility. The results reveal that KMC calculations agree well with experiment once adjustments are made to significant defect parameters within the appropriate uncertainty bounds.

36 MATERIALS SCIENCE↗

DOE OSTI2018

Verification of the prompt radiation shielding design of the ESS Accelerator (Stage 2)

The purpose of this CRADA is to perform additional analyses using MARS Monte-Carlo calculations code for verification of the prompt radiation shielding design of the ESS Accelerator. Stage 2 includes the following activities: 1.Activation assessment of the soil surrounding the wall of the ESS accelerator that maymigrate with the groundwater 2.Shielding assessment of the HVAC exhaust duct The existing ESS accelerator model (with possible minor modifications) shall be used in calculations. This work builds on the analysis originally performed under agreement CRADA FRA-2015-0302.

43 PARTICLE ACCELERATORS↗

DOE OSTI2019

Shortening the Design and Certification Cycle for Additively Manufactured Materials by Improved Mesoscale Simulations and Validation Experiments: Fiscal Year 2019 Status Report.

This report outlines the fiscal year (FY) 2019 status of an ongoing multi-year effort to develop a general, microstructurally-aware, continuum-level model for representing the dynamic response of material with complex microstructures. This work has focused on accurately representing the response of both conventionally wrought processed and additively manufactured (AM) 304L stainless steel (SS) as a test case. Additive manufacturing, or 3D printing, is an emerging technology capable of enabling shortened design and certification cycles for stockpile components through rapid prototyping. However, there is not an understanding of how the complex and unique microstructures of AM materials affect their mechanical response at high strain rates. To achieve our project goal, an upscaling technique was developed to bridge the gap between the microstructural and continuum scales to represent AM microstructures on a Finite Element (FE) mesh. This process involves the simulations of the additive process using the Sandia developed kinetic Monte Carlo (KMC) code SPPARKS. These SPPARKS microstructures are characterized using clustering algorithms from machine learning and used to populate the quadrature points of a FE mesh. Additionally, a spall kinetic model (SKM) was developed to more accurately represent the dynamic failure of AM materials. Validation experiments were performed using both pulsed power machines and projectile launchers. These experiments have provided equation of state (EOS) and flow strength measurements of both wrought and AM 304L SS to above Mbar pressures. In some experiments, multi-point interferometry was used to quantify the variation is observed material response of the AM 304L SS. Analysis of these experiments is ongoing, but preliminary comparisons of our upscaling technique and SKM to experimental data were performed as a validation exercise. Moving forward, this project will advance and further validate our computational framework, using advanced theory and additional high-fidelity experiments.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Patchy snapshots of nuclear chain reactions

Stochastic fluctuations of the neutron population within a nuclear reactor are typically prevented by operating the core at a sufficiently high power. This regime, where the evolution of the neutron density is essentially deterministic, is key for automatic protection and safety systems to safely detect unwanted power excursions during an accident, and to rapidly initiate a reactor trip procedure in case it is needed. Recent works, supported by numerical simulations, have however reported that, for large reactors, the branching nature of the fission reactions might induce strongly non-Poissonian patterns in the neutron spatial distribution, and that stochastic fluctuations might still persist at reactor powers close to operating conditions (startup phase). An international program conducted by LANL, IRSN and CEA was therefore setup to experimentally detect and characterize such fluctuations and correlations. An experiment took place in 2017 at the Reactor Critical Facility (RCF) of the Rensselaer Polytechnic Institute (USA). In this paper we will report the main findings of this experimental program, supported by stochastic models and by the development of a dedicated high-fidelity Monte Carlo simulation code. We will in particular describe and explain the strong patchiness in neutron power distributions measured at the RCF, as well as a peculiar ‘blinking’ behavior of the core, and discuss the consequences of these findings on nuclear safety.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2018

Transient Thermal Analysis of Calorimeters Used in Characterization of the ACRR Radiation Environments

Silicon calorimeters have been used for active radiation dosimetry in the central cavity of the Annular Core Research Reactor (ACRR) for over a decade. Recently, there has been interest in using other materials for calorimetry to accurately measure the prompt gamma-ray energy deposition in the mixed neutron and gamma-ray environment. The calorimeters used in the ACRR use a thermocouple (TC) to measure the change in temperature of specific materials in the radiation environment. The temperature change is related to the instantaneous dose received by the material in a pulse-transient operation. SOLIDWORKS Simulation and ANSYS Mechanical were used to model the calorimeter and analyze the thermal behavior under pulse-transient conditions. This report compares the results from modeling to experimental results for selected calorimeter materials and radiation environments. These materials include bismuth, tin, zirconium, and silicon. Calorimeters assembled with each material were irradiated in the ACRR central cavity in the free- field, LB44, CdPoly, and PLG radiation environments. The neutronics code Monte-Carlo N- Particle (MCNP) was used to calculate the neutron and gamma-ray response of the calorimeter materials at the experimental locations in the central cavity. Different response tallies were used and found to give different results for the gamma-ray energy deposition. It was determined that performing the neutron/gamma-ray/electron transport in MCNP using the *F8 electron tally gave the overall best agreement with the experimental results. The *F8 tally, however, is much more computationally intensive than the neutron/gamma-ray transport calculations. Also, this report contains parametric analyses that examine the ways to improve the current design of the calorimeters. One finding from the parametric analysis was that the TC should be placed closer to the outer radius of the disks to obtain a measurement closer to the maximum temperature of the disk. Also, the parametric analysis showed that the most dominant mechanism of heat loss in the calorimeters is conduction through the alumina posts. In future designs, the conduction should be minimized to reduce the effect of heat loss on the measurements.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2020

ADETS User Manual

ADETS (Automated DEpletion Transport System) calculates coupled neutronic/isotopic results for nuclear systems and produces a large number of criticality and burnup results based on various material feed/removal specifications, power(s), and time intervals. ADETS is a fully automated tool that links the LANL MCNP Monte Carlo transport code with the SCALE system, specifically ORIGEN-ARP (radioactive decay and burnup code). Additionally, ADETS can compute the dose at various distances of the reactor using a point gamma source, whose characteristics are computed by ORIGEN-ARP. In addition, ADETS is fully coupled with the Uncertanty Quantification, Probabilistic Risk Assessment, Parameter Optimization and Data Analysis code RAVEN (developed at the Idaho National Laboratory as well). RAVEN is released with ADETS but MCNP and SCALE (ORIGEN-ARP) are not included.

42 ENGINEERING↗

DOE OSTI2020

Verification of the CP2020 Library

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

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

YAHFC: A Code Framework to Model Nuclear Reactions and Estimate Correlated Uncertainties

Reaction modeling is a key ingredient in designing experiments and interpreting their results, and is an essential component in the process of evaluating nuclear data and assembling nuclear data libraries used in nuclear technology applications. Typically, experimental data are available only for a handful of reaction channels and theory models are used to fill in the gaps. In addition, theory is often called upon as the arbitrator between discrepant data. Most importantly, theory and modeling are required for an accurate determination of uncertainties in the evaluated data and the correlations between the multiple channels. A fast, accurate, and flexible modeling capability has been developed at LLNL with the code system YAHFC (Yet Another Hauser-Feshbach Code). YAHFC is a Monte Carlo, Hauser-Feshbach code framework, making full use of dynamic memory allocation, derived types, and parallel computing. YAHFC can generate events to simulate experiments and is guiding experiments designed to measure inelastic neutron scattering from actinide targets. YAHFC is also being used to analyze decays from surrogate experiments, thereby enabling the inference of reaction cross sections inaccessible by direct measurement. Finally, by modeling nuclear reactions with constraints from experimental data, YAHFC can deliver complete nuclear data libraries, with evaluated uncertainties, using the modernized Generalized Nuclear Data Structure (GNDS).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

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 OSTI2021

Application of Shift Ex-Core Calculations to a Detailed MSR Model

This report documents the collaborative research conducted between Argonne National Laboratory (ANL) and Oak Ridge National Laboratory (ORNL) on applying the Monte Carlo neutronics code Shift to calculate the radiation dose rates for a detailed molten salt reactor (MSR) model based on the Molten Salt Breeder Reactor (MSBR). In this work scope, the ANL and ORNL teams worked together and applied Shift to calculate the radiological environments of this detailed MSR facility modeled with explicit geometries. The radiological conditions within the MSR facility were modeled for when the reactor is at two different operation modes: normal full power and drained state. The FW-CADIS hybrid method in Shift was applied successfully to calculate the ex-core neutron and gamma dose rates for the MSBR at full power. Dose rate maps showed that in the current MSBR numerical model, potential pathways exist for neutrons and gammas to stream through the 8-foot concrete shield to reach the top of the MSBR reactor cell. Neutron and gamma dose rates within the drain cell were also calculated for the MSBR at the drained state by integrating the source terms obtained from an ORIGEN-S depletion calculation into the Shift simulation. The results indicate that in the drained state, the delayed gammas from the depleted fuel salt are the main contributors to the dose rates, which suggests that a steel liner in the tank model must be added for calculating the dose rates in the drain cell.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2022

Fluid-Kinetic Coupling: Advanced Discretizations for Simulations on Emerging Heterogeneous Architectures (LDRD FY20-0643)

Plasma physics simulations are vital for a host of Sandia mission concerns, for fundamental science, and for clean energy in the form of fusion power. Sandia's most mature plasma physics simulation capabilities come in the form of particle-in-cell (PIC) models and magnetohydrodynamics (MHD) models. MHD models for a plasma work well in denser plasma regimes when there is enough material that the plasma approximates a fluid. PIC models, on the other hand, work well in lower-density regimes, in which there is not too much to simulate; error in PIC scales as the square root of the number of particles, making high-accuracy simulations expensive. Real-world applications, however, almost always involve a transition region between the high-density regimes where MHD is appropriate, and the low-density regimes for PIC. In such a transition region, a direct discretization of Vlasov is appropriate. Such discretizations come with their own computational costs, however; the phase-space mesh for Vlasov can involve up to six dimensions (seven if time is included), and to apply appropriate homogeneous boundary conditions in velocity space requires meshing a substantial padding region to ensure that the distribution remains sufficiently close to zero at the velocity boundaries. Moreover, for collisional plasmas, the right-hand side of the Vlasov equation is a collision operator, which is non-local in velocity space, and which may dominate the cost of the Vlasov solver. The present LDRD project endeavors to develop modern, foundational tools for the development of continuum-kinetic Vlasov solvers, using the discontinuous Petrov-Galerkin (DPG) methodology, for discretization of Vlasov, and machine-learning (ML) models to enable efficient evaluation of collision operators. DPG affords several key advantages. First, it has a built-in, robust error indicator, allowing us to adapt the mesh in a very natural way, enabling a coarse velocity-space mesh near the homogeneous boundaries, and a fine mesh where the solution has fine features. Second, it is an inherently high-order, high-intensity method, requiring extra local computations to determine so-called optimal test functions, which makes it particularly suited to modern hardware in which floating-point throughput is increasing at a faster rate than memory bandwidth. Finally, DPG is a residual-minimizing method, which enables high-accuracy computation: in typical cases, the method delivers something very close to the $L^2$ projection of the exact solution. Meanwhile, the ML-based collision model we adopt affords a cost structure that scales as the square root of a standard direct evaluation. Moreover, we design our model to conserve mass, momentum, and energy by construction, and our approach to training is highly flexible, in that it can incorporate not only synthetic data from direct-simulation Monte Carlo (DSMC) codes, but also experimental data. We have developed two DPG formulations for Vlasov-Poisson: a time-marching, backward-Euler discretization and a space-time discretization. We have conducted a number of numerical experiments to verify the approach in a 1D1V setting. In this report, we detail these formulations and experiments. We also summarize some new theoretical results developed as part of this project (published as papers previously): some new analysis of DPG for the convection-reaction problem (of which the Vlasov equation is an instance), a new exponential integrator for DPG, and some numerical exploration of various DPG-based time-marching approaches to the heat equation. As part of this work, we have contributed extensively to the Camellia open-source library; we also describe the new capabilities and their usage. We have also developed a well-documented methodology for single-species collision operators, which we applied to argon and demonstrated with numerical experiments. We summarize those results here, as well as describing at a high level a design extending the methodology to multi-species operators. We have released a new open-source library, MLC, under a BSD license; we include a summary of its capabilities as well.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Modeling, Performance Assessment, and Nodal Data Analysis of TRISO-Fueled Systems with Shift

This technical report documents several enhancements to the Shift Monte Carlo (MC) code under the US Department of Energy (DOE) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program in fiscal year (FY) 2022. Performance enhancements were added to Shift specifically for tristructural isotropic (TRISO)–fueled reactor systems and guided based on performance analysis in FY 2021. For the pebble performance model developed in previous studies, the runtime improved by ~ 91× compared to the original model and ~ 2× compared to the user-optimized model. Compared to Serpent, Shift is ~ 3× slower if Serpent delta-tracking is enabled but ~ 2× faster when delta-tracking is disabled. The multigroup cross section generation was improved through simplifying tally input definitions, porting several post-processing tally operations from Python scripts into the Shift code base, and accounting for production reactions in the scattering multiplicity. Progress was also made on two emerging capabilities: (1) the development of Titan (a Shift reactor physics user interface) and (2) initial investigation into path-length tallies for computing multigroup scattering matrices.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗