Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MCNP6.2”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Coupled beam-target-moderator optimization for the Second Target Station

To support the design of the Second Target Station, that aims to provide the world’s highest peak brightness of cold neutrons, studies that optimize the dimensions of the target and moderators are invaluable. In this work, we investigate the influence of the target dimensions and beam profile on the performance and optimal size of the moderators. We perform optimization runs with detailed MCNP6.2 simulations using high-fidelity unstructured mesh geometries generated from parametrized CAD models. We demonstrate that small changes in target height do not influence the moderator performance if the beam dimensions are chosen adequately. We quantify the effect of beam footprint and target width on the moderator performance and show that the optimal moderator dimensions are insensitive to limited changes in target and beam profile.

Dakota↗
DOE OSTI2022

Dose conversion factors for absorbed dose in a mobile phone to absorbed dose in critical organs in an anthropomorphic phantom for emergency dosimetry applications: OSL and TL experimental results, and Monte Carlo simulations

Here, an anthropomorphic phantom was equipped with three smartphones and irradiated in five separate exposure geometries in air using a 137 Cs source. The phantom was fitted internally with thermoluminescence dosimeters (TLD-700 chips) to measure the absorbed dose to the individual organs, while the absorbed doses in the phones were evaluated using optically stimulated luminescence (OSL) from the phones’ surface mount resistors (SMRs), along with reference dosimeters, namely Luxel OSL films and TLD-700 chips. The results indicate a strong dependence on the phone location with respect to the source due to shielding and scattering from the body. Dose conversion factors were calculated for several plausible exposure scenarios in order to correlate the absorbed doses measured in the phones to the average body dose received by the phantom. The measured conversion factors were compared with calculated factors using MCNP6.2 Monte Carlo simulations and an ADAM voxel phantom. Experiments were conducted to assess how the phone orientation and angle impacts the dose received by the phone components. The overall results indicate that the doses measured can vary significantly depending upon exposure geometry (including phone position and angle). However, in rotational irradiation symmetry, the phone doses lead to an over-estimate of the average whole-body dose by an average of ~10%, independent of the phone locations tested in this study. This result, if confirmed in future studies, may be acceptable for triage.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2025

The MCNP ® 6 code: A decade of progress

After several years of effort involved in merging the Los Alamos National Laboratory MCNP5 and MCNPX codes, in 2013 the first production release of version 6 of the Monte Carlo N-Particle ® , or MCNP ® , code MCNP6.1 was distributed publicly. Since then, three significant releases have been issued: MCNP6.1.1beta in 2014, MCNP6.2 in 2018, and MCNP6.3 in 2023. While each release always contains new features, code enhancements, and bug fixes, each version has had a different primary focus, ranging from improved calculational efficiency to new powerful utilities and tools, to software modernization of the code base. With all that has been learned over the first decade of the MCNP6 code, continuous progress is being made toward a modernized, general-purpose Monte Carlo radiation transport code that remains a trusted resource for the global community of practitioners. This paper describes these first 10+ years of the MCNP6 code and its continually improving data libraries, and gives some insight into how the next decade is expected to unfold.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Decay Dose Shielding Analysis with Hybrid Unstructured Mesh/Constructive Solid Geometry Monte Carlo Calculation and ADVANTG Acceleration

Here, the target segments of the Oak Ridge National Laboratory Second Target Station (STS) neutron production facility become highly activated due to spallation reactions or nuclei transmutation by primary protons and emitted neutrons. Once the target segments are removed from their location within the core vessel, decay dose rates must be accurately quantified to determine the shielding configurations of remote-handling tools and transport casks and to aid in planning maintenance activities. For this analysis, we utilized a hybrid unstructured mesh (UM)/constructive solid geometry approach for calculating spallation products and neutron fluxes, activation calculations using the AARE package that includes the CINDER2008 activation code to calculate the decay photon source at different cooling times, and the ADVANTG code to accelerate the final decay photon transport calculation. Both Type 316 stainless steel (SS-316) and lead were investigated as candidates for shielding materials. The decay photon transport calculation through the thick SS-316 or lead shields exhibited between 25 and 30 orders-of-magnitude attenuations in the radial direction, depending on the shield. Such a difficult shielding calculation required advanced variance reduction. ADVANTG has some missing features, which limits its usability in spallation neutron source applications. It does not support volumetric sources created for MCNP6.2 UM capability. An approximate source was created for this problem. Not only was this approximate source needed for running the ADVANTG calculation to generate the weight windows, but also it was essential to develop source biasing (SB) parameters that were crucial for dramatically accelerating the decay photon transport in this problem. With this approximate source, the analysis was completed in a very reasonable computational time, and the design of the STS remote-handling equipment was finalized. This paper compares the efficiency of Monte Carlo simulations with different weight window and SB parameters calculated using different approximate ADVANTG calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Development and validation of fully open-source R2S shutdown dose rate capabilities in OpenMC *

We present the first fully open-source capabilities for shutdown dose rate (SDR) calculations of fusion energy facilities based on the Rigorous 2-Step (R2S) methodology. These capabilities have been implemented in the OpenMC Monte Carlo particle transport code, building on its existing capabilities while also leveraging new features that have been added to the code to support SDR calculations, such as decay photon source generation. Each of the individual physics components in the R2S workflow—neutron transport, activation, decay photon source generation, and photon transport—have been verified through code-to-code comparisons with MCNP6.2 and FISPACT-II 4.0. These comparisons generally demonstrate excellent agreement between codes for each of the physics components. The full cell-based R2S workflow was validated by performing a simulation of the first experimental campaign from the Frascati Neutron Generator (FNG) ITER dose rate benchmark problem from the Shielding INtegral Benchmark Archive and Database (SINBAD). For short cooling times, the dose calculated by OpenMC agrees with the experimental measurements within the stated experimental uncertainties. For longer cooling times, an overprediction of the shutdown dose was observed relative to experiment, which is consistent with previous studies in the literature. Altogether, these features constitute a combination of capabilities in a single, open-source codebase to provide the fusion community with a readily-accessible option for SDR calculations and a platform for rapidly analyzing the performance of fusion technology.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Historical review and proof-of-concept future method demonstration of adaptive mesh refinement in nuclear engineering for increased fidelity and computational efficiency

As the nuclear industry's use of computational tool increases, the need for increased fidelity and computational efficiency is well known. While most approaches to increased fidelity rely on applying a fine mesh over the problem domain, a more efficient method is to apply an adaptive mesh refinement (AMR) algorithm to the mesh definition. In the field of nuclear engineering, AMR has previously been used in conjunction with deterministic methods, including: S{sub N} transport methods, Lattice Boltzmann Methods, and COMSOL. The future of AMR in nuclear engineering is to couple it to a Monte Carlo code with the goal of reducing calculation time. A proof-of-concept example yielded positive results for using the gradient of the flux as a refinement criteria. The refinement criteria was varied from 0.01 to 0.10, which yielded a recommended range of 0.01 to 0.04, and the number of refinement iterations was varied from 0 to 7, with diminishing returns seen after 5 iterations. After the success of the proof-of-concept exercise, work began on creating a full program coupling MCNP6.2 and the AMR algorithm in the deal.II library. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

A Verification of Flux Sensitivity Estimates Using the MCNP Tally Perturbation Tool

Nuclear data is commonly used in applications such as nuclear nonproliferation, safeguards, and criticality safety. More specifically, nuclear data is used in predictive simulation codes like the Monte-Carlo N-Particle (MCNP ® ) transport code, Serpent, and similar radiation transport codes. The improvement of nuclear data enables more precise and accurate simulations, which result in higher fidelity designs and reduced operational/procedural costs. Therefore, the improvement of nuclear data is of paramount importance across the nuclear community. Nuclear data is improved and validated through integral benchmark experiments. The design of benchmark experiments is an extensive process; therefore, these experiments are often optimized on multiple characteristics, including sensitivity to the nuclear data, during the design process. Sensitivity is a measure of how much a quantity changes due to changes in independent variables such as experimental configuration. An experimental design that has a larger sensitivity to the nuclear data of interest will have a larger impact on the accuracy and precision of the validated data. Past integral benchmark experiments have primarily used the effective multiplication factor ($k_{eff}$) as the predominant measured quantity; however, experiments designed with other quantities in mind would be able to optimize on validating different areas of the nuclear data. A primary goal of the EUCLID project is to design, constrain, and reduce compensating errors in experiments focused on quantities other than $k_{eff}$ to better validate nuclear data across the board. Currently, there is a capability in MCNP to easily calculate the sensitivity of $k_{eff}$ to specific nuclear data of numerous reactions types and isotopes (KSEN card); however, the sensitivity of other quantities must be estimated in more strenuous manners. For example, the perturbation feature (PERT card) of MCNP can be used to estimate first-order sensitivities of some response in fixed source simulations. A recent announcement revealed that the first- and second-order perturbation features in previous releases of MCNP contained a bug. It was identified that particles were being scored into the wrong energy bin. The bug is in the most recent public release (MCNP6.2); however, a patch has been added to the most up to date version (MCNP6.2.2) that has not been released publicly. A direct comparison of the PERT card results for an F4 (neutron flux averaged over a cell) tally before and after the patch are shown in figure 1. All simulations used in the sensitivity estimates in this report were performed with MCNP6.2.2. This work verifies the patched MCNP perturbation tool by comparing first order sensitivities made using the PERT card to estimates made using manual perturbation of the compact ENDF (ACE) files.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Verification and Validation Tests of the Gamma Library of the ARC Software Package (Rev.1)

As part of the verification and validation (V&V) efforts of the Argonne Reactor Code (ARC) software system to support the Versatile Test Reactor (VTR) project of the U.S. Department of Energy (DOE), work has been performed to verify and validate the gamma (or photon) library of the multigroup cross section generation code MC 2 -3 of the ARC system, which is the essential part of the coupled neutron and gamma heating calculations to determine the heat generation rate or power distribution in the core. This V&V study was carried out in three steps: 1) review of the procedures and utility programs for generating the MC 2 -3 gamma library, 2) verification tests of cross sections and KERMA factors by comparing the total heat generation rate of coupled neutron and gamma calculation with that obtained with reaction Q values under the assumption of local gamma energy deposition and by comparing the prompt heat generation rate with that calculated with the MCNP6.2 code, and 3) validation tests of cross sections and KERMA factors by analyzing the ZPPR-15D gamma dose data measured with thermoluminescent dosimeters (TLD).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

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↗
DOE OSTI2022

MCNP6.3 Workshop [Slides]

This workshop presentation revisits several new features of MCNP6.2 and tools. Additionally, it offers a preview of the new features and improvements included in the upcoming MCNP6.3 code release that are most relevant to the nuclear radiation protection and shielding communities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Open-source release of CGMF 1.1 and Integration into the MCNP6.3 ® Code [Slides]

As a result of a multi-year NA-22 project, CGMF was integrated into MCNP6.2 and publicly released. CGMF was open-sourced and publicly released and MCNP6.3 was updated to include the latest version and is in the process of being publicly released. Current and future plans include global optimization and uncertainty quantification within CGMF, model parameter fitting such that CGMF may be used in ENDF/B evaluations, and improving both standalone and MCNP-integrated CGM (non-fission) simulations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Summary of LANL Critical Benchmark Comparison Study and Revisions for Cases Involving HEU, LEU, MIX, and Pu

This report documents results obtained for revisions made to cases involving Highly Enriched Uranium (HEU), Intermediate Enriched Uranium (IEU), a mixture of Pu and Uranium (MIX), as well as Pu cases. A previous summary of revisions for HEU an Pu cases was reported and additional investigations into four cases originally presented therein uncovered further revisions which led to better agreement with other transport codes, those cases are updated in this report. The summary of all cases reported in Reference 2 is updated in this report. In addition, a previous summary of revisions for LEU and MIX was reported, a summary of those revisions in reproduced in this report for a comprehensive summary of changes to benchmarks beginning in fiscal year 2020 to current date. The report focuses on the changes made to LANL benchmarks modeled with MCNP6 using ENDF/B-VII.1 nuclear data that appeared to have discrepant results when compared with results of other codes. Feedback was used to pinpoint review of benchmark input files and to revise them when necessary. This report documents the results of review and revision of specific benchmarks highlighted as possibly discrepant in the comparison study. In addition, there is an effort tied to this work involving collaboration between LANL XCP and NCS Divisions in the development of a shared review/revision procedure and use of a new benchmark repository. LANL has a benchmark library of critical experiments from the International Criticality Safety Benchmark Evaluation Project (ICSBEP) Handbook modeled for use with MCNP. This collection is now over 1100 benchmarks, referred to as the Whisper-1.1 library because it is used with the sensitivity/uncertainty package, Whisper, which supports nuclear criticality safety validation and is released with MCNP6.2. The collection, originally created several decades ago, is a combination of smaller collections, which has been revised and expanded, by various groups at LANL over the years. The original authors are no longer at the laboratory and little formal documentation of review and revision of these benchmarks exists today. A branch of the benchmark collection was already the subject of a formal review undertaken by the LANL NCS Division and expanded to include XCP Division.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

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↗
DOE OSTI2024

IER 480: TEX-Pu Benchmark (PU-MET-THERM-004) to Test Polyethylene and Lucite Thermal Scattering Laws [Slides]

This presentation finds that PMT-004 has four new benchmark cases highly sensitive to PE TSL (2 cases) and PMMA TSL (2 Cases). PE cases were well predicted using MCNP6.2 and ENDF/B-VIII.0. PMMA cases overpredicted by approximately 0.6-0.7% in k eff at 20°C. Accepted into 2023 version of the ICSBEP Handbook. Temperature had a large impact on reactivity of the critical configuration. Implications for validation work for thermal cases- need to adjust TSL data to correct temperature as it can have hundreds of pcm effects for a few °C. Future thermal experiments should try and measure reactivity at multiple temperatures to aid in data testing and benchmark adjustment.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Estimate of Gamma Dose Rates from Arrays of Fermi-1 Blanket Elements During the MEDE Process

The Enrico Fermi Atomic Power Plant Unit 1 (Fermi-1) was a sodium-cooled fast breeder reactor located in Monroe County, Michigan. The reactor was powered by a core of enriched uranium metal alloy driver fuel, which was enveloped by an axial and radial blanket material consisting of depleted uranium metal alloyed with 3 wt.% molybdenum. There are 406 axial and 559 radial irradiated sodium-bonded Fermi-1 blanket assemblies in storage at INL, totaling 34 metric tons of heavy metal. Disposal of the Fermi-1 blanket material directly into a geological repository is prohibited due to the reactive characteristic of its bond sodium. A melt drain evaporate (MEDE) process can effectively remove bond sodium from Fermi-1 blanket material. Consequently, planning is underway to apply a MEDE process to treat the 34 metric tons of heavy metal of irradiated Fermi-1 blanket material. Given the irradiated Fermi-1 blanket material’s relatively low power history and 50+ years of decay, it can be treated in a shielded glovebox. To assess the requisite shielding, the dose rates of the Fermi-1 blanket elements in various process configurations are needed. The method to perform the dose rate calculations is to first generate an average source term and associated photon source spectra for the Fermi-1 blanket material in SCALE and then use the associated photon spectra to calculate dose rates using MCNP6.2 in seven representative geometries. Dose rates in rem/h were calculated on contact (1 cm from outer geometry surface), 30 cm away from outer geometry surface, and 1 m from outer geometry surface at axial heights spanning the length of the blanket material within the geometry. The maximum average dose rate for the single Fermi-1 radial blanket element is ~0.60 rem/h, for the Fermi-1 radial blanket assembly ~1.10 rem/h, for the single Fermi-1 axial blanket element ~0.5 rem/h, for the Fermi-1 axial blanket assembly ~1.01 rem/h, for the MEDE can ~ 0.98 rem/h, for the Fermi Storage Canister ~ 0.83 rem/h, and for the MEDE cans in the DOE Standard Canister ~0.51 rem/h. Based on these maximum average values, the bounding dose rates are assumed to be 3× the average dose rates calculated for the axial blanket material and 7.6× the average dose rates calculated for the radial blanket material, based on the distribution of Fermi-1 blanket material assembly burnup.

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

A Python Script to Read MCNP6.3 Surface-Source Files

This report provides a Python script to read an MCNP ® surface source file created with the SSW card with SYM = 0 (i.e., the default symmetry treatment). For background: the general format of an MCNP surface-source file is described in; however, that document did not provide coding and/or a tool to interrogate such files. The current format will not be given in this document other than through the record-read statements necessary for the script to function. The reader capability in this report is augmented with the ability to directly write a couple demonstrative outputs: 1. A comma-separated value (CSV) file containing particle phase-space state information and 2. A Matplotlib histogram of the energy distribution of the particles. This report also describes accompanying verification work that shows the script performing as required with MCNP6.2, MCNP6.3, and (expected) MCNP6.4 surface-source files. However, users of the enclosed script must still verify that the script is behaving correctly for their own work.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

An update to δ-ray production in MCNP6

As originally introduced in MCNP6.2, the δ-ray production subroutine had a number of flaws related to treatment of particle or nuclear spin contributions to the production cross section. A recent update by the author to the MCNP6 code corrects these flaws. This report documents those changes for the benefit of MCNP6 users who use the delta ray production capability. Additionally, this report aims to provide a more comprehensive documentation of the delta ray production, building on the prior documentation by C. Anderson and coworkers. On this latter point, the goal is to save future MCNP developers from a visit to the library to dig out ancient texts from 1952, or at least to provide useful points of reference should such a visit be necessary.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Coupled Target-Beam-Moderator Optimization for the Second Target Station

This report describes the results for a coupled target-beam-moderator optimization analysis for the Second Target Station (STS) at ORNL's Spallation Neutron Source. This study is a continuation of the optimization analysis for the moderators in the preliminary design of STS performed in 2022. In the 2022 analysis the dimensions of the moderators are parameterized, while the target and the proton beam profile are kept constant. In this analysis the target height and the proton beam profile are added as parameters. This allows to study the coupled effects of changing target, moderator and beam dimensions. Similar to the 2022 analysis, this work is performed with an automated optimization workflow that uses the optimization toolbox DAKOTA, parameterized geometries in CREO and SpaceClaim, the unstructured mesh generation in Attila4MC, and the particle transport code MCNP6.2©. This workflow enables an efficient optimization using high-fidelity geometries. The main conclusions of this analysis are the following: • Coupled beam-target-moderator optimization provides a few additional percent performance gain over stand-alone moderator optimization. • The moderator performance is not very sensitive to the target height (between ≈60 and ≈80 mm) as long as the beam profile is chosen adequately. • The moderator performance is sensitive to the choice of beam spatial standard deviations, even when the footprint is kept constant. • The optimal moderator radius is the same for a beam footprint of 30 cm 2 , 62.5 cm 2 , and 90 cm 2 . Also the slope of the super-gaussian beam profile does not significantly impact the optimal moderator radius. • The optimal parameters and sensitivities are very similar to the 2022 optimization analysis. These results only indicate a a difference in the optimal radius of the cylindrical moderator, however, this has been corrected in the final design moderator optimization. The main purpose of this report is to document the simulations, results and lessons learned. The most impactful results are summarized in. We also note that the target geometry used in this work is not the final design.

43 PARTICLE ACCELERATORS↗