Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo N-Particle Transport 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.

84 records · Page 5

Particle Swarm Optimization Algorithm for Critical Experiment Design

Nuclear criticality experiments are used to validate nuclear cross section data used by simulation software. This is typically achieved by designing a critical system with a high sensitivity to a certain material’s cross section. Once the experiment has been carried out, a high fidelity model of the system is developed into a benchmark. When this benchmark model is simulated by a transport code, some of the difference between the experimental and computational effective neutron multiplication factor can be attributed to inaccurate nuclear data. Nuclear data evaluators then can make adjustments accordingly to improve cross section data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Verification and validation testing and tools: comparison between MCNP code versions and nuclear data libraries [Slides]

This presentation discusses the primary goal of software testing which is to test the code for correctness. It also discusses the results for individual suites and the role of validation and verification also referred to in the presentation as V&V. In summation, the V&V framework enables easy comparison between calculations performed with different code versions and/or nuclear data libraries. This entire framework will be distributed with the upcoming MCNP6.3 release. V&V test suites shown and several that were not (Criticality, LAQGSM, Lockwood) will be distributed in the new framework.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Neutronics Analysis of Cold Critical KRUSTY Experiments using MCNP and Serpent

Kilowatt Reactor Using Stirling Technology (KRUSTY) is a prototype designed as a proof of concept for NASA’s Kilo power program, which was funded in order to create a small reactor for various space applications such as providing power for a colony on Mars or the moon and eventually as a possible option for powering rockets on deep space missions. KRUSTY was tested at the National Criticality Experiments Research Center (NCERC) located at the Nevada National Security Site (NNSS) from November 2017, through March 2018. The experimental campaign was conducted in four phases: com ponent critical experiments, cold critical experiments, warm critical runs, and high temperature demonstration. The component critical and cold critical experiments were used to determine the worth of the beryllium oxide (BeO) reflector rings and the boron carbide (B4C) control rod disks under different conditions. These conditions were altered by the addition and subtraction of other components in the assembly. The component critical phase was performed at atmospheric pressure and had four solid support rods spanning the center core; whereas the cold critical phase had the center core placed in a vacuum chamber and eight sodium-filled heat pipes spanning the center core region. This summary discusses the static neutronic analysis completed on the cold critical phase of the campaign using MCNP ® and Serpent. These two codes were used to assess any potential differences or biases in the simulation results. Understanding these potential biases is a key step toward the final goal of being able to fully analyze the transient experiments performed during phases 3 and 4 of the campaign.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Easy_PERT: a Python tool for writing PERT cards and parsing PERT card results [Slides]

This presentation begins by providing an overview of the PERT card. The PERT card uses differential operator method to compute first- and second-order tally variations due to density, composition, and reaction cross-sections. It is possible to have multiple PERT cards in one MCNP input deck to study tally variations for several sets of nuclides, reactions, and energy ranges. Furthermore, the METHOD option tells MCNP to calculate either the perturbed tally (METHOD=-1, -2, -3) or the change in the unperturbed tally (METHOD=1, 2, 3). In summation, a powerful use-case for the MCNP code PERT card is that it facilitates calculating tally sensitivities to nuclear data. Writing PERT card entries and parsing output MCTAL files is tedious and error prone. however, Easy_PERT makes use of existing tools (Faust and MCNPTools) to handle writing PERT card entries and parsing the output MCTAL files. The PERT card is early in the development process and planned upcoming capabilities include calculating sensitivities and combining MCTAL files from separate runs into one JSON file.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

MCNP6.3 Code and Nuclear Data Installation Guide [Slides]

This presentation discusses the MCNP6.3 release package and examines how it offers more MCNP6-specific executables and options available than ever before. It also discusses the installer and that it streamlines the MCNP6.3 installation and automatically downloads all production nuclear data and generates xsdir files for the user. Additionally, the nd_manager is also discussed regarding how it allows the LANL Nuclear Data and Monte Carlo teams to work and distribute their products independently. It also allows updates, fixes (errata), and newly released data to be made available in a much faster timescale than ever before.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗

Verifying MCNP Models of the TEX High 240 Plutonium Benchmark

Computational modeling programs are invaluable tools that allow us to understand systems, safely develop new processes, and make reliable predictions about future designs. However, the effectiveness of these codes is limited by the degree to which their parameters match the real world. In the field of nuclear engineering, cross section data is one of these vital parameters. Accurate cross section data on important fissile and fissionable isotopes promotes the design of safer and more efficient fabrication, transportation, storage, and stockpiling of nuclear fuel. Unfortunately, there are knowledge gaps in data on key isotopes. In 2011, a multinational meeting hosted by the US Department of Energy Nuclear Criticality Safety Program ranked the priority of certain cross section data needs. In response, Lawrence Livermore National Lab (LLNL) designed the Thermal and Epithermal eXperiment (TEX) series of benchmark experiments. Benchmark experiments are used to validate current cross section data. They validate data by comparing the results of an actual experiment to the predicted results from a computational model. The data a benchmark applies to depends on the isotope and energy range the experiment’s neutron multiplication factor ( k eff ) is most sensitive to. The development and testing of the TEX High 240 Plutonium Benchmark will help validate 240 Pu cross section data. The configuration and materials of this benchmark are designed to be most sensitive to 240 Pu's intermediate energy range (from 0.625 ev to 100 keV ). MCNP® models of the assembly have been developed by LLNL and the results have been written in the final design report. In order for the discrepancies between benchmark models and experiments to be attributed to cross section inaccuracies, the accuracy of the models needs to be verified. The goal of this project is to verify of the results of LLNL's modeling by creating a new set of MCNP models and comparing the results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Coincident Capture through Post-processing PTRAC [Slides]

This presentation discusses the new PTRAC capabilities and workflows. The PTRAC capability in MCNP6.3 has seen a massive overhaul since MCNP6.2. The new HDF5 file format allows for both MPI- and thread-based parallelism. MCNPTools has been updated to handle the new HDF5 PTRAC format and is now open sourced on GitHub. Built-in capabilities, such as the pulse-height tally coincident capture special treatment, can largely be replicated through separate postprocessing scripts that leverage both PTRAC and MCNPTools. This allows for greater flexibility in user-specified and controlled detector response functionality, ultimately using the MCNP code for what it is best at (i.e., particle transport).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

MCNP ® Particle Transport Code: Its History and Future [Slides]

The history, development and application uses of Monte Carlo particle transport codes at LANL are extensive. Recent and ongoing investments in MCNP modernization are enabling new code developments using modern programming tools and languages. New and challenging applications continue to drive new research and developments beyond the rich set of established MCNP capabilities. The LANL MCNP radiation transport code has been the international gold standard for particle transport applications for over three decades.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Second Target Station High-Fidelity Target Activation Comparison

The development of the Second Target Station (STS) target system at the Spallation Neutron Source (SNS) at Oak Ridge National Laboratory (ORNL) is well underway. The target system at STS consists of a rotating target disk that contains 21 segments of tungsten clad in tantalum clad in steel. A key aspect of the design of the target system is to account for the delayed heating and material damage caused by the delayed dose from decaying radionuclides. These radionuclides are a product of either spallation reactions or transmutation of the nuclei in the target system. These radionuclides build up in the target system components over the lifetime of the facility, and the radiation that is emitted can deposit energy in the components causing significant component heating and material damage. Monte Carlo N-Particle (MCNP) Version 6.2 transports the various particle species and calculates the spallation products and neutron fluxes throughout the target system. These spallation products and neutron fluxes along with the material definition of each component are relayed to the CINDER2008 transmutation code to calculate the radionuclide inventories and the corresponding decay gamma emission spectra. MCNP6.2 coupled with CINDER2008 is the computational method-of-choice for the analysis discussed in the following sections of this report. The analysis focuses on validating major assumptions in calculating the radionuclide inventory in the STS target system: all of the target segments are fresh, unirradiated material when the protons are incident on the segment, the average of the 21 segments of the target is sufficient to represent a single segment, and that averaging the proton pulse structure over time does not significantly affect the radionuclide inventory. The position-averaged, high-fidelity, and single-tally computational methods are used to validate the assumptions and provide a point of comparison to evaluate how the assumptions impact the radionuclide inventories. A more detailed explanation of the three computational methods is provided in Section 2. The position-averaged and single-tally methods are less computationally expensive when compared with the high-fidelity method where 54,000 individual calculations are needed to calculate 1 hr of STS operation. Section 3 details the comparison of the three methods to show that the assumptions made in the position-averaged method do not significantly impact the radionuclide inventory after 1 hr of operation. The discussions and results in this report are for 1 hr of operation. Due to the computational cost associated with calculating the transmutation and activation using the high-fidelity method, only 1 hr of operation has been calculated. The discrepancies observed after 1 hr of operation are not extrapolated out to longer operational times, and this report does not address how the discrepancies between the computational methods may manifest for longer operational periods.

43 PARTICLE ACCELERATORS↗

Improvements in MCNP6.3 (Rev. 1) [Slides]

This presentation outlines improvements in the Monte Carlo N-Particle Code (MCNP) version 6.3. MCNP is an internationally renowned Monte Carlo particle transport code developed by Los Alamos National Laboratory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

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↗

Hybrid Particle Swarm-Interpolation Algorithm for Improved Critical Experiment Design

The typical goal of designing a critical experiment is twofold: a system that is both critical and optimized for some other value. This value could be an energy-integrated sensitivity, percent fissions in a certain energy range, or some other value that can be calculated by a transport code. By simulating different combinations of reflector, moderator, and fuel thicknesses a designer can identify such a desirable configuration. The domain of all possible combinations of these thicknesses is referred to as the experiment search space. As more dimensions are added, the search space increases in size exponentially. For a three-dimensional case, which includes three thickness values between zero and ten centimeters to the nearest tenth of a millimeter, a total of 1,000 3 , configurations exists. Rather than check each configuration individually, which would be extremely computationally expensive, it has been shown to be more efficient to use a conventional particle swarm optimization (PSO) algorithm to identify critical and optimal configurations. This work presents the theory and implementation of a novel hybrid PSO interpolation algorithm to perform these optimizations faster than a conventional PSO algorithm. To demonstrate this, an example optimization will be carried out by the conventional and hybrid PSO algorithms and their performances will be compared.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗