Search NASA⌕ Search

SEARCH · Search NASA

Results for “mcnp”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Benchmarking GEANT4 and PHITS for 14.8-MeV neutron transport in polyethylene and graphite materials

The feasibility of GEANT4 and PHITS with ENDF/B-VIII.0, JEFF-3.3 and JENDL-4.0 evaluated nuclear data libraries in 14.8-MeV neutron transport simulation are verified by comparing with the MCNP calculations and the experimental data. Here the measured data are from neutron integral experiments including graphite and polyethylene samples. For 6 cm-thick polyethylene and 20 cm-thick graphite samples, the results from GEANT4 and PHITS with three evaluated nuclear data libraries agree well with the MCNP calculations. For a 2 cm-thick graphite sample, the simulated results from GEANT4 and PHITS agree with the MCNP calculations, and they all underestimate the experiment data in below 3 MeV and overestimate them in the energy range of 5–7 MeV. The difference between the simulation results and the experimental data depends on the evaluated nuclear data libraries. In overall, GEANT4 and PHITS are performing reasonable jobs for 14.8-MeV neutron transport simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

PyTrac

SAND2025-00635O PyTrac is a software tool that analyzes and visualizes PTRAC event files generated by MCNP 6.3. It converts these files into a graph network that makes it easier to interpret individual histories. The software includes command line tools for viewing the graph data structure in both 2D and 3D formats. PyTrac also integrates with MCNP to run simulations and manage data files. This provides a streamlined approach to analyzing and understanding the complex data generated by MCNP simulations. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Nowack, Aaron [Sandia National Lab. (SNL-CA), Live↗
DOE OSTI2024

RadSim: Particle Transport Wrapper

RadSim is being developed to provide the capability to simulate radiation source emissions, interpolate results from radiation transport tools into a common format to prepare incident flux, and model radiation detector response. This software package includes wrapper API code which is utilized in the development of the second (transport) task of the RadSim project. This software package only contains tools to interface RadSim with radiation transport tools (GEANT4 and MCNP) and does not include any part of the transport source code itself. The users will need prior access to the transport tools (GEANT4 and MCNP), and can use this software package to interface with the transport tool of their choice. This software package takes the geometry file for a radiation transport problem and converts it into an input deck for the transport tools (GEANT4 and MCNP). Additionally, it also contains filters that convert the output of the transport tools to a common flux format.

Cheung, Hoi Sing↗
DOE OSTI2020

Demonstration of NEAMS Multiphysics Tools for Fast Reactor Applications

The SHARP toolkit is a high-fidelity reactor simulation tool developed under the U.S. Department of Energy, Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) Campaign. SHARP toolkit is comprised of the neutronics module PROTEUS thermal hydraulics module Nek5000, and structural mechanics module Diablo. During FY17 and FY18, the PROTEUS and Nek5000 components of SHARP were applied to solve challenging sodium-cooled fast reactor (SFR) problems. In particular, selected hot channel factors (HCF) for a prototype metal-fueled SFR design (the AFR-100) were analyzed in high fidelity, and the “SHARP zooming capability” for SFRs was developed and demonstrated to reduce computational expense for full core problems in cases where detailed data is needed in selected fuel assemblies. After the previous success applying SHARP to challenging SFR problems, the focus in FY19 and FY20 expanded to additional fast reactor applications including lead cooled fast reactors (LFR) and sodium cooled fast reactors (SFR). The specific technical tasks were (1) assessment of hot channel factors for LFR, for which no data currently exists, and (2) demonstration of zooming capability in assemblies of the Versatile Test Reactor (VTR). First-of-a-kind hot channel factor (HCF) estimation for LFR with high fidelity codes (PROTEUS/Nek5000) was successfully demonstrated in this study which began in FY19 and continued in FY20. Selected HCF were computed and compared with SFR data (AFR-100, EBR-II). The findings confirm that different reactor types, design parameters and uncertainties lead to different HCFs. Careful estimation of HCF for a specific design is necessary to obtain appropriate HCFs. In addition to improvement in HCF accuracy, high fidelity tools generate data to help the designer better understand the mechanism of the impact from these uncertainties. For example, the impact of cladding thickness manufacturing tolerance resulted in non-intuitive effects in the corner pins of the LFR assembly. This procedure of computing HCF using high fidelity models shows promise and flexibility for being repeated for any arbitrary reactor of choice. Along with the application on SFR and LFR, the capability of the tools has also been matured to deal with different reactor types and designs. Progress was made towards extending the previously demonstrated SHARP zooming capability to non-fueled SFR assemblies. In particular, in FY19 a gamma transport capability was implemented in both high fidelity PROTEUS solvers in order to accurately account for heat deposition caused by gamma particles, which accounts for ~10% of total core power. Neutronics verification cases were carried out for a candidate Versatile Test Reactor (VTR) design using the new gamma transport capability in PROTEUS. Comparisons were made with continuous energy MCNP calculations and shown to agree well. The models for the full core design with heterogeneous control and fuel assemblies is in progress for PROTEUS-SN and completed with MCNP. The MCNP power distributions were transferred to Nek5000 to perform thermal hydraulic calculations of the control and fuel assembly.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Analysis of Approximations in Modeling of BWR Bundle Void Distributions

In boiling water reactors, complex heterogeneous bundle designs, control blades adjacent to the corner of bundles, and the presence of boiling can lead to complex internal void distributions. A few approximations exist to model these void distributions. They could be modeled using a 1D axial solver in which each axial node is assumed to be at an average void, or each pin cell could be modeled with its own void concentration. In the latter case, the void could be discretized in pin-centered or coolant-centered channels. The goal of this project was to quantify the effect of using the different approximations for modeling internal void distributions on neutronics calculations. Using 3D void distributions calculated with CTF, Monte Carlo Neutral Particle (MCNP) transport code models were created for GE-9 and GE-14 lattices. For each model, the internal void distribution from CTF at a given axial node was selected, and a lattice calculation was carried out with MCNP. Comparisons between models using a lattice-averaged void, or using a void distribution in coolant-centered channels, showed large differences in reactivity which in some cases were well above 1,000 pcm, and it also showed differences in normalized fission rates greater than 20%. It was also found that using a lattice average void can lead to a significant difference in the estimation of the worth of a control blade. The differences found when comparing results from models using pin-centered and coolant-centered channels were up to 200 pcm in reactivity and up to 1.4% in the normalized fission rates. In addition to these two sets of comparisons, MCNP models were set up so that each subchannel had a saturated liquid component around the fuel pins and a saturated vapor component in the center to approximate annular flow. In comparison to the models using coolant-centered subchannels, up to 1–3% differences in normalized fission rates could be found.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Initial Investigation of Coupled Electron-Photon Transport with SCEPTRE and RAPTURE (Rev.1.1)

The SCEPTRE1 code models coupled photon and electron transport in material media using discrete ordinates methods. The code can use different solvers in 1D and 3D in parallel, promising users the ability to run problems that are significant both in physics and computational scale. One capability that is perhaps unique at this time in the NNSA complex is SCEPTRE’s coupled photon– electron transport physics in a discrete ordinates approach. For example, incident x–rays can generate a significant electron population, whose transport can be simulated within the same code. In this context, fully–coupled means that photons can produce electrons and electrons can produce photons; uncoupled or no–coupling means photons do not produce electrons. Partial–coupling, an intermediate option in which photons can produce electrons but electrons do not produce photons, is not considered in the present work. The MCNP code can treat this physics with a Monte Carlo approach. MCNP should be viewed as complementary to SCEPTRE in both computational method (Monte Carlo versus discrete ordinates) as well as data libraries. SCEPTRE and MCNP should provide important checks on each other, especially from complementary modeling frameworks of a S n (discrete ordinates) mesh based method versus Monte Carlo particle based method. In order to learn more about SCEPTRE, and how it might inform engineering work at the Laboratory, we used it to carry out a set of initial studies. This was a series of 1D simulations of photon beams hitting Aluminum slabs. Our studies mainly used SCEPTRE directly, and made little use of the closely associated RAPTURE code (more on this below). The goals of these studies were first, to familiarize ourselves with operating the code, and second, to begin to assess the importance of coupled electron–photon transport. To support these studies, we began to develop a set of scripts to initialize inputs, run tools, process outputs, and plot results. An initial version of these scripts has been committed to the Lightning workflow repository. This memo relates the results of those studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Monte Carlo Analyses of the FOEHN Experiment

The FOEHN critical experiment, which has been carried out at Cadarache (France) in the reactor EOLE in the early 70s, was designed to verify RHF design analyses. This study focuses on the impact of the FOEHN experiment uncertainties on the calculation of the energy deposition. The latter has been obtained by the Monte Carlo codes Serpent and MCNP. In the calculation of the energy deposition, the first code relies on KERMA factors whereas the latter code relies on Q-values. The two Monte Carlo codes share the same geometry and material specifications. In addition to the design and experiment parameters uncertainties, nuclear data also impacts the calculation of the energy deposition. In this study, the MCNP simulations use three different nuclear data library sets: ENDF/B versions VII.0, VII.1, and VIII.0. Finally, the ksens card of MCNP, using perturbation theory, has been used to identify the cross sections and isotopes with largest impact on the effective multiplication factor. By showing that the use of different software and computational methodologies and design parameters uncertainties do not significantly impact the core power distribution, this work demonstrate the appropriateness of the tools and methods employed for the conversion of involute-plate reactors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Assessment of the Impact of Realistic Sensor Physics and the Integration of Ex-Core Sensors on Reactor Power Synthesis

In the work documented in this report, a weighting function–based core power synthesis method was applied to multiple Monte Carlo N-Particle (MCNP) reactor models, which are informed based on simulated self-powered neutron detector (SPND) responses. The weighting function method used has been coined the point-based iterative (PBI) method. The goal of this application is to assess the impact of considering realistic sensor physics in the generation of the simulated SPND outputs as well as to consider how the synthesis is impacted based on the inclusion of ex-core detectors in the model. The NuScale small modular reactor (SMR) and Westinghouse AP1000 pressurized water reactor (PWR) are the models that served as the testbeds for the assessment of realistic sensor physics; this was achieved by using Geant4 SPND models in comparison with analytical models, such that the effect of electron transport in realistic SPND geometries in the Geant4 model can be understood in terms of synthesis error and convergence time. The comparison was considered for fuel burnup–induced perturbations, for a range of sensor string densities and synthesized power distribution axial fidelities. The Texas A&M Testing, Research, Isotopes, General Atomics Reactor (TAMU TRIGA) reactor MCNP model was used to assess the impact of ex-core sensors; this was done by performing synthesis with and without the ex-core detectors and by quantifying the synthesis error and number of iterations associated with Gaussian-type perturbations in many locations in the core. The TAMU TRIGA model was particularly pertinent for this study because of the interest in future experimental tests with SPNDs in this reactor, as well as the ease of modifying the MCNP model to include ex-core detectors with heterogeneously described response functions. Results from the comparison between the Geant4 and analytical SPND models indicate that similar average and maximum synthesis errors were obtained for burnup-induced perturbations in both the NuScale SMR and the AP1000. This was true for a range of sensor string densities and axial fidelities. However, there were marked differences between both the Geant4 and analytically informed models in terms of the iterations required to converge on the synthesized power distribution. Namely, the Geant4-informed models tended to lead to fewer iterations, except for a few sensor–core configurations that had particularly numerous iterations. Results from the ex-core sensor assessment with the TAMU TRIGA model indicate that the inclusion of ex-core sensors drastically reduces the synthesis error of Gaussian-type perturbations close to the edge of the core, and it slightly reduces synthesis errors for perturbations closer to the center of the core. This was achieved with a minimal increase in computational cost—that is, the number of iterations required for convergence. The errors were identified to be in the same location as the perturbation in the core, indicating that the methodology remains robust for unperturbed regions of the core. A secondary result from this study with the TAMU TRIGA was yielded by analysis of the neutron flux levels in the in-core and ex-core sensor locations of the core; these flux levels indicate that SPNDs could be used as both in-core and ex-core sensors, so long as the emitter material is sensitive to thermal neutrons. The results from these studies provide a quantitative understanding of the importance of considering realistic sensor physics and including ex-core sensors to perform accurate and timely power distribution synthesis of a reactor core.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Methodology to Generate the Decay Gamma Source for a Second Target Station Target Wedge

This report details the specifics of applying the Position-Averaged Method methodology for the calculation of the rotating lasagna-style target wedge activation and the resulting decay gamma source terms. The methodology uses MCNP® Code Version 6.2.0 with the RNUCS patch coupled with CINDER2008 from the AARE V1.0 package. The process begins with model preparations that include ensuring impurities are included in the materials in the regions of interest along with segmentation of the geometry in order to capture the decay gamma source gradients. This report also details the physics models and other MCNP®, CINDER2008, and ADVANTG 3.2.0 related options used to calculate the decay gamma source terms. ADVANTG 3.2.0 is used to generate a weight window set to aid in convergence of the neutron fluxes in the initial MCNP® protons-on-target particle transport calculation to calculate the neutron fluxes and spallation products in the regions of interest. The neutron fluxes and spallation products are then used in the CINDER2008 transmutation calculation, and the decay gamma sources are generated by the AARE_GAMMA_SOURCE_SCRIPT as a part of the AARE V1.0 package. These decay gamma sources include the energy distributions as well as the necessary distributions to sample the location of the source using the cell rejection method. The gamma sources that are produced using the methodology may be used for a wide variety of analysis involving the target wedge. This report includes the material specifications, physics options, proton source term description, geometrical configuration, and post-processing details.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Parallel Programming in MCNP6

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

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

JMOCUP Physics Depletion Calculations for the As-Run AGR-5/6/7 TRISO Particle Experiment in ATR Northeast Flux Trap

This ECAR documents the as-run Jim Sterbentz’s MCNP ORIGEN Coupled Utility Program (JMOCUP) physics depletion calculation and the calculated results for the AGR-5/6/7 irradiation experiment in the Advanced Test Reactor (ATR). AGR 5/6/7 was irradiated for nine power cycles in the northeast flux trap. The depletion calculations were performed to provide input data for a variety of other engineering analyses supporting the AGR-5/6/7 experiment along with post-irradiation characterization of the tri-structural isotropic (TRISO) particle fuel compacts. Detailed full-core MCNP models and Oak Ridge Isotope Generation (ORIGEN2) radionuclide generation models were specifically developed as part of the JMOCUP Monte Carlo depletion calculations. The MCNP ORIGEN2 computer codes were coupled using the well-established JMOCUP utility modules to couple the two codes and run the depletion calculations. The physics calculations were performed in support of the Advanced Gas Reactor (AGR) program.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Predicting deviations in ex-core detector response in Krsko NPP during an earthquake with Monte-Carlo neutron transport methods

To study the response of ex-core detectors during an earthquake, an approach using Monte Carlo code MCNP was developed. As a first step to predict the response of the ex-core detectors a detailed MCNP model of the Krsko reactor core was created. A detailed MCNP core model was used to calculate 3D power distributions inside the core. Calculated power distributions were verified by comparison to the CORD-2 calculations, which is currently used for core design calculation verification of the Krsko nuclear power plant. Due to the outside the core position of the neutron detectors, the hybrid code ADVANTG is used to generate weight windows to speed up neutron transport outside the reactor core. To be able to use ADVANTG, fixed neutron source had to be reconstructed from the criticality core calculation. It was found that the power range ex-core detector movement has a negligible effect on the value of the thermal neutron flux at the active part of the detector. However, the 5 cm radial movement of the intermediate-range detector leads to 7 % - 8 % change in thermal neutron flux within the active intermediate-range detector region. The analysis continued with the evaluation of the effect of core barrel movement on the ex-core detector response. It was determined that the 2 mm core barrel radial oscillation can lead to ∼ 4% change in thermal neutron flux within the active detector region. It was evaluated that reactor core movement inside the baffle can lead up to ∼ 1.9% and ∼ 1.3% deviation in power-range and intermediate-range detector signal respectively. Analysis showed that the mechanical movement of ex-core neutron detectors could not explain the fluctuations in the ex-core detector signal. However, combined core barrel and reactor core inside baffle oscillations could be a probable reason for the observed fluctuations in the ex-core detector signal during an earthquake. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Cross-Section Comparison for Pu-238 Production in the Advanced Test Reactor at Idaho National Laboratory

Qualification of Advanced Test Reactor (ATR) positions for Pu-238 production has been ongoing at Idaho National Laboratory (INL). The ATR qualifications have stretched over multiple years during which new techniques have been developed and made available for ATR experiment neutronic analysis. As part of the transition to newer codes, new cross-section libraries have been evaluated for use in the Pu-238 production experiment analysis. A comparative study was done using the MCNP ORGIEN Activation Analysis (MOAA) tool between ENDF/B-VII.0 and ENDF/B-VIII.0 cross sections to capture the impact of the change in cross-sections on the analysis needed to qualify Pu-238 production targets. All comparisons were done assuming the ATR GEN-I targets were located in the south flux trap of the ATR. Finally, an overview of how this qualification and potential irradiation fits into Pu-238 is discussed.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2022

Nondestructive and destructive assay for forensics characterization of weapons-grade plutonium produced in LEU irradiated in a thermal neutron spectrum

A post-irradiation examination (PIE) of an enriched uranium sample was performed to advance a nuclear forensics methodology, based on intra-elemental nuclide ratios. The PIE was carried out on a few milligrams of low enriched uranium dioxide (LEUO 2 ) sample, irradiated to a burnup of 1 GWd/MTU. The PIE of the LEUO 2 sample provided the concentration of selected plutonium and fission product nuclides ( 239 Pu, 240 Pu, 241 Pu, 91 Y, 95 Zr, 95 Nb, 103 Ru, 133 Cs, 134 Cs, 135 Cs, 137 Cs, 136 Ba, 138 Ba, 140 Ba, 140 La, 141 Ce, 144 Ce, 149 Sm, 150 Sm, 152 Sm, 153 Eu, and 154 Eu). The PIE supported validation of the MCNP (Monte Carlo N-Particle neutronics simulation code) predicted concentration of these nuclides in the irradiated LEUO 2 sample. This study showed that most of the MCNP predicted nuclide concentrations were within 15% and several within 10% compared to the PIE results. Furthermore, gamma and mass spectrometry of nuclides in irradiated LEUO 2 enabled the determination of uranium burnup and time since irradiation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Comparison of active interrogation methods for source location in a scattering and absorbing medium, consisting of PGNAA, and an AmBe quasi-forward biased directional source

Source location of Special Nuclear Material (SNM) encompassing 95% 235 U and 239 Pu is identified by utilizing a directional source from Patent No: US20190013109A1 using Prompt Gamma Neutron Activation Analysis (PGNAA) and neutron spectroscopy simulated with Monte-Carlo N-Particle transport 6.2 (MCNP). BC-408, HPGe, LaBr 3 detector arrays were used to identify the location of the SNM using total counts incident on each detector, and PGNAA photopeaks from HPGe and LaBr 3 detector arrays in a polyethylene shield. The conducted simulations varied the volume and location of the SNM in the MCNP input files to observe how the source location method behaved. PGNAA photopeaks used for source identification include 61 keV from fission, 2.223 MeV prompt gamma from hydrogen, 511 keV annihilation, and a single and double escape peaks from the prompt gamma interaction from hydrogen. The capabilities of each detector systems to acquire well resolved photopeaks with a 1% relative error or less, and total relative error for F4 and F8 tallies were less than 0.015% relative error. Source predictions of the SNM with uneven amounts of polyethylene shielding between the source and detectors was observed to overpredict and give invalid source location predictions. Source locations of the SNM with even amounts of polyethylene material between the source and each detector were found to be valid. With a 1 Ci 241 Am source activity, it was determined that 1630 s were needed to obtain the results for each detector system with the quasi-forward directional AmBe source. Coupling source and material identification together would increase acquisition time but would only require one system to determine.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Benchmarking Monte Carlo codes for the modelling of low-energy neutron production target reactions

The increasing adoption of accelerator-based neutron sources (ABNS) for applications including neutron capture therapy (NCT) research has highlighted the need for accurate simulation tools. Precise modelling of the neutron production target is crucial to ensure that simulated predictions of neutron beam characteristics used for subsequent beam shaping assembly design are reliable. This work presents a comprehensive benchmarking of four widely-used Monte Carlo codes - Geant4, PHITS, FLUKA (CERN), and MCNP - for modelling low-energy neutron production target reactions. Using their recommended physics models and cross-section libraries, we evaluate each code’s performance in simulating four beam-target reactions: 7 Li(p,n) 7 Be, 9 Be(p,n) 9 B, 9 Be(d,n) 10 B, and C(d,n)N. Predictions of neutron yield, angular distributions, and energy spectra are compared against available thick target experimental data. Results show varying levels of agreement between the codes depending on the reaction type, energy range, and beam characteristics. Geant4, MCNP and PHITS are the overall best performing codes for the simulation of total neutron yield and yield in the forward direction across most reactions. Across energies where experimental benchmarks exist, inter-code discrepancies in total and forward-directed yield are typically 10 to 30%, with larger deviations at near-threshold incident ion energies. PHITS provides the best overall reproduction of experimental spectra, particularly for the 9 Be(p,n) 9 B reaction. Additionally, PHITS demonstrates superior computational performance for most reactions. These findings provide valuable guidance for ABNS design, highlighting the strengths and limitations of each code for the simulation of low-energy neutron production reactions.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

An unstructured mesh based neutronics optimization workflow

We have developed a fully automated workflow to optimize the neutronics performance of the Second Target Station (STS) at the Oak Ridge National Laboratory’s Spallation Neutron Source. The optimization workflow starts with the parametrized solid CAD engineering models and converts them into the unstructured mesh (UM) models for the neutronics calculations with MCNP6.2. Calculations are executed and their results are loaded into the Dakota optimization toolkit. Dakota analyzes the results and proposes new geometry parameters for the next design iteration. The cycle repeats until the optimal parameters are found. The automated CAD to MCNP conversion, the use of high-fidelity UM models, and the use of modern optimizer are the key elements that advance the entire optimization workflow in comparison with the original workflow. The original workflow was based on a simplified constructive solid geometry (CSG) modeling with MCNPX, mcnp_pstudy tool, and an in-house optimizer. Herein to demonstrate the new workflow, we present a case of neutronics optimization of the moderator–reflector assembly (MRA). Apart from the MRA, the workflow can optimize other major STS components, such as the spallation target, neutron beamlines, radiation shielding, and various accelerator components. Importantly, the new workflow opens the door to the advanced multi-physics multi-parameter optimization and has the potential for use in other nuclear physics and accelerator applications.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Advanced Modeling and Simulation Methods for Evaluation of Thermal Neutron Scattering Materials

With the rise of interest in thermal neutron scattering data for advanced reactor, criticality safety, and shielding applications, new experimental data are required for evaluation of new materials or for re-evaluation (or validations) of previously evaluated materials. New experimental data are evaluated in a three-step process: (1) computing the phonon characteristics, (2) computing the dynamic structure factor (DSF) from the data, and (3) using the experimental setup to simulate the experimental data. All three steps have challenges, ranging from the need for a sufficiently general material simulation code—a processing code that can compute the corresponding DSF—to having a detailed layout of the instrument/beamline/facility where the data were measured. Whereas phonon characteristics of materials can be computed using various methods (molecular dynamics, density functional theory, etc.), a high-fidelity computation of the DSF and the simulation of the experiment based on the DSF is vital to the accuracy of the evaluation. The latter two steps can be achieved by using the two corresponding code systems developed by instrument scientists at the Spallation Neutron Source (SNS) at Oak Ridge National Laboratory: (1) OCLIMAX, a program that calculates the dynamic structure factor from DFT and MD simulation results, and (2) MCViNE, a Monte Carlo neutron ray-tracing program designed to simulate neutron scattering experiments. Recently, polyethylene and yttrium hydride were measured at the Wide Angular-Range Chopper (ARCS) and SEQUOIA instrument stations of the SNS. These experiments are simulated using the density functional theory code, the Cambridge Serial Total Energy Package (CASTEP), to compute its phonon characteristics (eigenvalues/vectors and PDOS), which is then processed using OCLIMAX to yield the DSF, and finally the data at each instrument station are simulated by the MCViNE for comparison to the measured data for evaluation. For comparison to conventional evaluation methods, the scattering data processed from OCLIMAX are compared against those processed from the LEAPR module of NJOY, and the results from MCViNE simulations are compared against previously used simplified beamline models implemented in the Monte Carlo N-Particle (MCNP) code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗