Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo Particle Transport”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Efficiency Studies of Fast Neutron Tracking Using MCNP

Fast neutron identification and spectroscopy is of great interest to nuclear physics experiments. Using the neutron elastic scattering, the fast neutron momentum can be measured. Wang and Morris introduced the theoretical concept that the initial fast neutron momentum can be derived from up to three consecutive elastic collisions between the neutron and the target, including the information of two consecutive recoil ion tracks and the vertex position of the third collision or two consecutive elastic collisions with the timing information. Here, we also include the additional possibility of measuring the deposited energies from the recoil ions. In this paper, we simulate the neutron elastic scattering using the Monte Carlo N-Particle Transport Code (MCNP) and study the corresponding neutron detection and tracking efficiency. The corresponding efficiency and the scattering distances are simulated with different target materials, especially natural silicon (92.23%28Si, 4.67%29Si, and 3.1%30Si) and helium-4 (4He). The timing of collision and the recoil ion energy are also investigated, which are important characters for the detector design. We also calculate the ion traveling range for different energies using the software, “The Stopping and Range of Ions in Matter (SRIM)”, showing that the ion track can be most conveniently observed in 4He unless sub-micron spatial resolution can be obtained in silicon.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Impact of the delayed energy deposition in TREAT on transient dynamics

Estimating the energy released during a reactivity insertion accident (RIA) is of prime importance. TREAT is an experimental reactor designed to simulate these transients. Typically, models of these transients assume instantaneous energy deposition. However, a delayed energy deposition component exists. This work leverages the advanced capabilities of ALFRED, a high-fidelity Monte Carlo code for particle transport, to study the impact of this delayed component on RIA transient dynamics in TREAT. The study reveals that the delayed component can affect the energy released by up to $3\%$ for natural pulse transients. Understanding the impact of the delayed energy deposition component will enable the development of more accurate models, which will assist in designing experiments in TREAT and evaluating their safety.

73 - NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

ECAR-7300 Rev 1 Verification and Validation of MCNP6.2 for MARVEL Neutronic Analysis

This report documents the verification and validation (V&V) efforts of the Monte Carlo N-Particle transport code (MCNP) version 6.2 on the Sawtooth supercomputer for the Microreactor Applications Research Validation and Evaluation (MARVEL) microreactor required for the preliminary documented safety analysis. This document records V&V for a safety, hazards, analysis, and design software used for design and analysis of safety class structures, system and components (SSC)s.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
DOE OSTI2025

COG11.3 Abstract

COG is a high-resolution code for the Monte Carlo simulation of coupled particle transport in arbitrary3-D geometry. COG will transport neutrons, protons, deuterons, alpha particles with energies up to hundreds of GeV, and photons with energy ranges limited by the available cross section sets and physics models. Electrons can be transported via the EGS5 electron transport kernel, electrons can also be transported. The COG code is a significant upgrade from earlier Monte Carlo transport codes and has been written specifically to makes it more versatile, accurate, and easy to use. COG has provisions for calculating deep penetration (shielding) problems, criticality problems, and neutron activation problems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Forward Modeling of Gamma Reaction History Signatures From Anticipated Deuterium-Tritium Filled MagLIF Implosions on Sandia’s Z-Machine

Nuclear reaction history measurements provide a bang time and burn width of Inertial Confinement Fusion (ICF) implosions and are essential for understanding implosion performance to constrain ICF capsule design. When fusion fuel contains Deuterium (D) and Tritium (T) gas, reaction history is informed by measuring the 16.75 MeV gamma rays generated from the D(T,γ) 5 He reaction. Such DT based reaction history measurements have not been made on the Magnetized Laser Inertial Fusion (MagLIF) platform on Sandia’s Z-Machine due to the lack of Tritium being used. The recent development of ICF implosions with tritiated fuel will open the possibility of measuring the gamma reaction history on the Z-Machine. A forward model of the Gamma Reaction History diagnostic on Z (GRH-Z) has been developed using the MCNP6.3 (Monte-Carlo N-Particle) radiation transport code. The model included the Z-Machine geometry of interest to characterize the impact of neutron induced gamma rays on the DT reaction history measurements. In addition, the impulse response functions of the GRH-Z diagnostic to understand the temporal response of the detector, and the minimum yields required to make a reaction history measurement were calculated. This approach also predicted that with T 2 gas doping of MagLIF implosions a reaction history may be made for high performance shots >8e12-2.4e13 depending on the chosen threshold for the detector, with a maximum signal to background ratio of 25%. It was found that for long duration ICF implosions that additional collimation will be needed to prevent the neutron induced gamma rays from modifying the shape of the measured DT reaction history curve.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Calculations versus measurements for residual dose rates from SNS spent structures

Here, residual dose rate measurements were conducted on target vessel #13 and proton beam window #5 after extraction from their service locations. These measurements are used to verify calculation methods of radionuclide inventory assessment that are typically performed for nuclear waste characterization and transportation of these structures. Neutronics analyses for predicting residual dose rates are carried out using the transport code MCNPX and the transmutation code CINDER90. For transport analyses a complex and rigorous geometry model of the structures and their surroundings are applied. The neutronics analyses are carried out using the Bertini and CEM high energy physics models for simulating particles interactions above the table-based cross section range. Finally, obtained calculational results are analyzed and compared to the measured dose rates and overall show good agreement within 25%, which shows applicability of the methods used in analyses.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Validation of the single-event method for low-energy electron transport via stopping power calculations with $\mathrm{MCNP}$

Monte Carlo simulations of low-energy ( <50 keV) electron transport in matter are essential for a broad range of application fields. Several Monte Carlo codes have developed specialized treatments for this case, but a comprehensive validation of low-energy electron transport for general-purpose simulations remains lacking in the literature. One approach to accomplish this validation is calculation of stopping power using low-energy electron transport physics, as stopping power is a fundamental radiation transport quantity which must be simulated accurately for nearly any application. Here in this work, we use the Monte Carlo N-Particle (MCNP) radiation transport code with the single-event method for electron transport to calculate stopping powers of low-energy electrons (50 eV to 30 keV) in 41 elemental solids, 14 compound solids, and five rare gas solids, comparing simulation results to published semi-empirical stopping power calculations from optical measurements. In general, the simulations give good agreement (typically within ±10%) with semi-empirical stopping power values at higher energies: 300 eV and above for most elemental solids, 1 keV and above for compound solids, and 400 eV and above for rare gas solids. Agreement between MCNP and semi-empirical values is worse below these energies. The most significant source of error is the EPRDATA14 cross section data, which does not account for changes in electronic structure due to solid-state bonding, particularly in compound materials. The simplistic model of atomic excitation used to generate the EPRDATA14 cross sections is another key source of error. Additionally, the breakdown of the continuous slowing-down approximation introduces significant uncertainty at low energies, although this is a limitation of the calculation method and not of the simulation procedure. Accounting for these and other uncertainty sources, the single event method in MCNP is robust and able to give good accuracy for a variety of low-energy electron transport problems through diverse kinds of materials.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2022

Verification of Upcoming MCNP Features For Estimating Nuclear Data Sensitivities in Fixed Source Simulations [Abstract]

Predictive simulation codes, like the Monte Carlo N-Particle (MCNP) transport code, are used throughout the nuclear community. These simulations are based on nuclear data. Maximizing the accuracy and precision of nuclear data maximizes the accuracy and precision of the overall simulation. This is imperative to applications that rely on simulations. For example, improving nuclear data for special nuclear material improves simulation accuracy in stockpile stewardship applications, which results in larger safety margins and decreased operational costs. The improvement and validation of nuclear data is completed through integral benchmark experiments. Past benchmarks have primarily been limited to focus on the effective multiplication factor ($\kappa$ eff ); broadening the purview of benchmarks beyond $\kappa$ eff -dependent nuclear data addresses nuclear data deficiencies. Different response types depend on different areas of nuclear data. This dependence is quantified as nuclear data sensitivity: the change in response due to perturbation of a contributing parameter. The larger the nuclear data sensitivity of a response, the more the experiment is influenced by the uncertainties of the nuclear data. The optimization of nuclear data sensitivities in future benchmarks would result in more detailed validation of lesser studied areas of nuclear data. Currently, direct sensitivity capabilities are not easily found for all experiment types and parameters. An MCNP tool to directly estimate the cross section sensitivities of tallied values is under development. Additionally, updates have been made to the perturbation feature of MCNP, which can be used in a less direct approach to estimating sensitivities. This work verifies these features to estimate nuclear data sensitivities in fixed source simulations of a 4.5-kg sphere of alpha- phase weapons-grade plutonium surrounded by differing amounts of copper and polyethylene. Integrated estimates made using MCNP’s tools were found to statistically agree with integrated estimates made from manual perturbation of nuclear data proving the validity of the MCNP tools.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
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

Carrier gas attenuation in gamma assay of radioactive xenon samples

Gamma attenuation attributable to xenon and carbon dioxide carrier gas in radioactive xenon assay samples has been calculated using Monte Carlo n-Particle (MCNP) transport code for the 0.65 mL flame sealed ampoule (bean), and 15 mL and 50 mL Schlenk tube geometries. The attenuation was correlated as a % loss per standard mL of carrier gas. In this work, high Purity Ge gamma spectroscopy experiments were performed to measure against the MCNP results. The experimental results from the 133 Xe in xenon carrier gas in a bean geometry showed 0.97 ± 0.4% loss per standard mL xenon, which was in agreement with the MCNP predictions.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2023

Utilization of ACE nuclear data file toolkit ACEtk to calculate relative sensitivity coefficients of point-kinetics parameters

Sensitivity and uncertainty methods are quintessential for nuclear criticality safety and experiment design. This type of analysis relies on calculations of sensitivity coefficients; sensitivity coefficients of the effective neutron multiplication factor with respect to some nuclear data are predominantly calculated and used. As a part of the Laboratory Directed Research & Development project EUCLID (Experiments Underpinned by Computational Learning for Improvements in nuclear Data) at Los Alamos National Laboratory, sensitivity coefficients of many radiation detector measurement responses with respect to nuclear data were investigated. Specifically, this paper outlines a method to calculate point-kinetics parameters relative sensitivity coefficients with respect to nuclear data. Point-kinetics parameters such as the prompt neutron decay constant, effective delayed neutron fraction, and neutron generation time are especially important to experimenters and reactor operators designing systems with dynamic neutron populations. This method couples capabilities of the ACE (A Compact ENDF) nuclear data file toolkit, ACEtk, with the ability to load cross sections into the radiation transport code Monte Carlo N-Particle (MCNP). In conclusion, key aspects of optimizing this method for a particular application and sensitivity profiles of the Jezebel criticality experiment are examined and discussed.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2023

Preliminary verification of the MCNP perturbation and fixed-source tally sensitivity tools

Integral benchmark experiments are vital in the adjustment and validation of the nuclear data that govern predictive simulations across the nuclear community. The nuclear data sensitivity capabilities of the Monte Carlo N-Particle (MCNP ®) transport code are currently limited; however, expanding sensitivity capabilities will allow benchmark experiments to be designed to resolve compensating errors and adjust nuclear data where previously prohibitively difficult. This paper provides details of a preliminary verification for the use of (i.) the recently revised perturbation and (ii.) developmental fixed-source sensitivity tools within MCNP to calculate sensitivities of tallied responses (such as current integrated over a surface, F1, and flux averaged over a cell, F4) to nuclear data in fixed-source simulations. Energy-binned and energy-integrated sensitivities calculated with these tools are compared against sensitivities calculated using a central-difference approximation. Here, the verification is completed for four configurations of a benchmarked system using a 4.5-kg plutonium sphere surrounded by varying amounts of copper and/or polyethylene. The results show that sensitivities calculated with the perturbation and fixed-source sensitivity tools agree with the central-difference-based approach.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Multi-modal tomographic imaging system for poolside characterization of nuclear test fuels: Design considerations and studies

Testing and qualification of advanced nuclear fuels involves an iterative process of prototyping, in-pile irradiation testing, and in-situ or ex-situ examination. Fuel restructuring and fission product migration during burnup are among the most important aspects of fuel evolution that affect several important performance characteristics such as heat removal, accident tolerance, and fission product retention. Pool-side non-destructive characterization techniques provide fuel developers with tools to understand fuel evolution at different points of burnup. A design for a compact, submersible, and multi-modal gamma-ray tomography system for imaging irradiated nuclear fuel is presented here. Detector selection, collimator geometry and fabrication, mechanical design, imaging protocol and acquisition speed are discussed. Modeling calculations show that sub-millimeter resolution can be achieved in both transmission computed tomography images as well as in emission computed tomography images in a matter of hours. Several design compromises and fabrication challenges are discussed which should be taken into consideration for future submersible gamma-ray tomography instruments.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Optimization of the Second Target Station cold source moderators using an automated workflow

The Second Target Station (STS) at the US Department of Energy’s Oak Ridge National Laboratory is designed to become the world’s highest peak-brightness spallation source of cold neutrons. Successful completion of the STS, which is currently in the preliminary design phase, will provide transformative new capabilities to examine novel materials for future technologies. At STS, neutrons will be generated by spallation reactions in a solid tungsten target. They will be moderated and thermalized in two cold (20 K) para-hydrogen moderators. Careful optimization of these moderators is essential to the project’s success. To find optimal moderator designs, an advanced optimization workflow integrates high-fidelity neutronics calculations using the Monte Carlo N-Particle (MCNP) transport code MCNP6.2 with state-of-the-art optimization algorithms in the Dakota optimization toolkit. For each design iteration, a parametrized solid CAD geometry is generated in Creo and automatically converted into an unstructured mesh geometry by Attila 4MC for the neutronics calculation with MCNP. Iterations repeat until optimal designs are found. Herein this paper presents the results of a sensitivity and optimization study for the cylindrical and tube moderators. Both moderators can be optimized for maximum peak brightness, maximum time-integrated brightness, or any combination between these extremes. Maximum peak brightness is achieved by using smaller optimal dimensions of the moderators, whereas maximum time-integrated brightness is achieved by using larger dimensions. A Pareto front details the designs that optimally balance both brightness metrics. The Pareto front can be found in only 40–110 iterations with 4–5 design parameters when using the efficient global and Pareto-set optimization algorithms in Dakota. Additionally, important engineering constraints can be taken into account, such as the coupling between the cylindrical moderator radius and aluminum vessel wall thicknesses required to ensure structural integrity of the vessels. This interaction has a significant impact on the resulting optimal designs. Our new, highly efficient, fully automated optimization workflow will be used to optimize additional STS components in the future and can be adopted for design and optimization studies at other experimental neutron and accelerator facilities.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

An NCrystal extension for thermal neutron scattering in MCNP

A new extension for MCNP has been developed that incorporates the NCrystal thermal scattering package in the form of a Fortran to C++ wrapper, which adds the capability to describe thermal neutron scattering in the NCrystal formalism instead of the ACE tabular format. This can be used to described poly- or single-crystalline materials as well as liquids. In conclusion, the functionality of this extension to MCNP is illustrated using a layered crystal monochromator, leakage of neutrons from a moderating sphere of material, and by simulating criticality benchmarks.

MCNP↗
DOE OSTI2024

Shielding development for the spallation neutron source VENUS instrument

VENUS is a neutron imaging instrument that will use a broad range of neutron wavelengths, from epithermal to cold energies, and will include enhanced contrast mechanisms. It will offer novel energy-selective imaging techniques directly connecting complex engineering materials and systems’ structures, properties, and functions to reveal practical and fundamental answers about their real-world performance. Here, the instrument will be built at SNS beamline 10, facing the decoupled poisoned hydrogen moderator. The driving cost for the instrument is the beamline and instrument cave shielding. Final analyses were performed to evaluate the thickness and composition of shielding materials for the instrument cave and beamline to meet radiation safety criteria for the instrument to start up in 2024 after completing the SNS proton power upgrade.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Determining the 9 $\mathrm{Be}$($n, γ$) 10 $\mathrm{Be}$ integral cross section at fission neutron energies

The 9 B neutron capture cross section has significant implications for Be materials in the nuclear industry as well as the α process in stellar nucleosynthesis. While the cross section is well constrained at thermal neutron energies, there is a lack of experimental data at higher neutron energies, and the evaluated nuclear data libraries can differ by up to two orders of magnitude. We calculate the 9 Be(n, γ) 10 Be integral cross section at fission neutron energies in an effort to resolve disagreements amongst the nuclear data libraries. Foil irradiation experiments were performed using the Flattop critical assembly at the National Criticality Experiments Research Center with either the highly enriched U or Pu cores, with target foil stacks placed at multiple locations to exploit different neutron energy profiles. Accelerator mass spectrometry was used to measure the 10 Be/ 9 Be ratio in irradiated Be foils, while all other activation products were quantified through gamma spectrometry. The experiments were simulated using the Monte Carlo N-Particle radiation transport code and combined with experimental results to determine the total neutron fluence, while the staysl-pnnl suite and fispact-ii code were used to validate the model and assess the systematic uncertainty. The new 9 Be(n, γ) 10 Be integral cross sections calculated in this work are 26.5 ± 2.2µb at 0.59 ± 0.07 MeV, 24 ± 3 µb at 0.98 ± 0.14 MeV, 21.7 ± 1.3 µb at 1.26 ± 0.11 MeV, 21.8 ± 1.4 µb at 1.32 ± 0.11 MeV, and 18.6 ± 1.1 µb at 1.46 ± 0.13 MeV. These results do not agree with integral cross sections from any of the nuclear data library evaluations. Discrepancies between the new integral cross sections reported here and the nuclear data libraries suggest a more complex cross-section structure in the MeV range which allows for more resonance contributions, and more work is needed to further constrain the evaluated cross sections.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2020

Processing MCNP Elemental Edit Outputs

The Los Alamos National Laboratory's (LANL) Monte Carlo N-Particle® (MCNP®) transport code version 6 has the capability for tracking particles on unstructured mesh (UM) geometry models. The MCNP UM feature was originally designed for the UM models created by the Abaqus/CAE software product. The MCNP version 6.2.0 and later recognizes the UM models read from Abaqus INP files or MCNPUM files converted from Abaqus INP files. Other finite element analysis software packages may generate UM models and then convert these models into Abaqus INP formats.

97 MATHEMATICS AND COMPUTING↗