Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “NEUTRON 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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Vacancy defect control of colossal thermopower in FeSb2

Abstract Iron diantimonide is a material with the highest known thermoelectric power. By combining scanning transmission electron microscopic study with electronic transport neutron, X-ray scattering, and first principle calculation, we identify atomic defects that control colossal thermopower magnitude and nanoprecipitate clusters with Sb vacancy ordering, which induce additional phonon scattering and substantially reduce thermal conductivity. Defects are found to cause rather weak but important monoclinic distortion of the unit cell P n n m → P m . The absence of Sb along [010] for high defect concentration forms conducting path due to Fe d orbital overlap. The connection between atomic defect anisotropy and colossal thermopower in FeSb 2 paves the way for the understanding and tailoring of giant thermopower in related materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Uncovering Where Compensating Errors Could Hide in ENDF/B-VIII.0

Unconstrained physics spaces between two or more nuclear data observables in a library occur when their values can be simultaneously adjusted without violating the uncertainties in either differential information or simulations of relevant integral experiments. Differential data are often too imprecise to fully bound all nuclear data observables of interest for application simulations. Integral data are simulated with combinations of nuclear data so that an error in one observable may be hidden by a counterbalancing error in another. In this manner compensating errors may lurk within nuclear data libraries and these errors have the potential to undermine the predictive power of neutron transport simulations, particularly in situations where there is no conclusive validation experiment that resembles the application of interest. The EUCLID project (Experiments Underpinned by Computational Learning for Improvements in Nuclear Data) developed a preliminary workflow to identify these unconstrained physics spaces by bringing together results from a large collection of integral experiments with their simulated counter-parts as well as differential information that have a one-to-one correspondence to nuclear data. This wealth of information is processed by machine learning tools for subsequent refinement by human experts. Here, we show how the EUCLID work-flow is executed by applying it first to 239 Pu and then to 9 Be nuclear data in ENDF/B-VIII.0.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Optical design for the single crystal neutron diffractometer Pioneer

Pioneer is a single-crystal neutron diffractometer optimized for small-volume samples and weak signals at the Second Target Station at Oak Ridge National Laboratory. This paper presents the preliminary optical design progress, focusing on the rationale behind key design choices. It covers the T 0 and bandwidth disk choppers, guide and beam control system, incident-beam polarizer, scattering beam collimators, and additional strategies. The chopper locations are selected to maximize neutron transport while taking advantage of standardized shielding structures. To accommodate the maintenance shield, operational shutter, and polarizing V-cavity, the guide design includes significant gaps. When these optical components are moved out of the beam path, oversized collimators, rather than guides, will be translated in. Pioneer will utilize slit packages to control beam size and divergence and a translatable polarizing V-cavity. Absorbing panels are strategically placed near the end station to minimize background. An oscillating radial collimator, operating in a shift mode, will be used with the vertical cylindrical detector, while a fixed multi-cone collimator will be used with the bottom flat detector. In conclusion, these collimators will enable the detection of weak signals when complex sample environments are used.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Multilevel-in-Space-and-Energy CMFD in VERA

For full-core modeling in the Virtual Environment for Reactor Analysis (VERA), the three-dimensional multigroup eigenvalue neutron transport problem is solved by MPACT. To improve the efficiency of MPACT, advancements have been made in the transport accelerator. Multilevel-in-energy and multilevel-in-space coarse mesh finite difference (CMFD) solvers were developed to improve the efficiency of the CMFD accelerator. In this paper a new multilevel-in-space-and-energy CMFD solver is developed with coarsening in both space and energy on every level. Several different strategies are investigated for coarsening groups in energy. Modified V-cycle and multiple-cycle algorithms are evaluated for solving the multilevel equations. The performance of these solvers is compared for typical full-core reactor physics problems.

42 ENGINEERING↗
DOE OSTI2022

Multiphysics analysis system for heat pipe cooled micro-reactors employing PRAGMA-OpenFOAM-ANLHTP

A multiphysics analysis system for neutronics/thermo-mechanical/heat-pipe analysis of heat pipe cooled micro-reactors was developed using the PRAGMA code as the neutronics engine. PRAGMA, which used to be a GPU-based continuous-energy MC code for power reactor applications, now has an extended geometry package to handle geometries with unstructured meshes generated by the ANSYS Design-Modeler and Meshing. The NVIDIA ray tracing engine OptiX was exploited for efficient neutron transport on unstructured geometry. On the multiphysics side, the open-source CFD tool OpenFOAM and one-dimensional heat pipe analysis code ANLHTP were adopted. The manager-worker system based on the MPI dynamic process management (DPM) model enables efficient coupling of codes employing different parallelization schemes. With all the features, the multiphysics analysis of the one-sixth symmetrical MegaPower 2D core was performed and it demonstrated the benefit of the tight integration of the three-way coupling system and one-to-one geometry coupling strategy. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Enhancing Monte Carlo Workflows for Nuclear Reactor Analysis with Metamodel-Driven Modeling

Monte Carlo codes are essential components of many reactor physics simulation workflows as high-fidelity continuous-energy neutron transport solvers. Among Monte Carlo radiation transport codes, MCNP is particularly notable due to its diverse simulation capabilities, large user base, and long validation history. Despite being a powerful simulation tool, MCNP provides limited capabilities to allow automated execution, model transformation, or support for user-defined logic and abstractions that limit its compatibility with modern workflows. Here, to better integrate MCNP into a modern scientific workflow, we have developed an intuitive yet full-featured MCNP Application Program Interface (API) in Python, named MCNPy, which provides a specialized set of classes for MCNP input development. Moreover, to guarantee that our reading, writing, and modeling capabilities remain self-consistent (and to render the huge scope of the MCNP API manageable), we have adopted a strategy of model-driven software development in which a generalized model of the MCNP input format has been created. From this generalized model, or “metamodel,” problem-specific implementations such as an engine for input validation or a codebase for programmatic operations may be automatically generated. Since MCNPy primarily acts as a Python front-end to the underlying Java API that directly interfaces with the metamodel, it is intrinsically linked to the metamodel and thus remains maintainable. With MCNPy, users can programmatically read, write, and modify any syntactically valid MCNP input file regardless of its origin. These capabilities allow users to automate complicated tasks like design optimization and model translation for nuclear systems. As examples, this work demonstrates the use of MCNPy to find the critical radius of a plutonium sphere and to translate a 9000+ line MCNP input file into a corresponding OpenMC model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Demonstration of TOFFEE: A Response Uncertainty Quantification Tool

A key characteristic in neutron transport is nuclear data. Cross-section uncertainty is not used in MCNP6.3 to propagate response uncertainty without external analysis. Here, the TOol For Fast Error Estimation (TOFFEE) is a Python-based code developed to automate the propagation of cross-section uncertainty for MCNP evaluations. TOFFEE implements the sandwich rule to calculate the uncertainty from cross sections with sensitivity coefficients from MCNP6.3 and ENDF/B covariance data. In this paper, TOFFEE has been tested with benchmark experiments, and it has been compared to the uncertainty quantification capabilities of Sampler and TSUNAMI, within SCALE, to verify the application’s capabilities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A Polarization Tensor for Neutron Optics

A rank 2 dimension 3 tensor formulation for polarized neutron transport in a magnetic field is derived. Using this formulation, the rotating magnetic field problem is revisited to produce a single 3 × 3 vector-transformation matrix that describes the polarization transport. That result is then used to fashion a new strategy for a numerical approximation of polarization transport in an arbitrary magnetic field using time-dependent perturbation theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Flat-band-induced itinerant ferromagnetism in RbCo 2 Se 2

ACo 2 Se 2 (A=K,Rb,Cs) is a homologue of the iron-based superconductor, AFe 2 Se 2 . From a comprehensive study of RbCo 2 Se 2 via measurements of magnetization, transport, neutron diffraction, angle-resolved photoemission spectroscopy, and first-principle calculations, we identify a ferromagnetic order accompanied by an orbital-dependent spin-splitting of the electronic dispersions. However, we identify the ordered moment to be dominated by a d x 2 –y 2 flatband near the Fermi level, which exhibits the largest spin splitting across the ferromagnetic transition, suggesting an itinerant origin of the ferromagnetism. In the broader context of the iron-based superconductors, we find this d x 2 –y 2 flatband to be a common feature in the band structures of both iron chalcogenides and iron pnictides, accessible via heavy electron doping.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Growth and characterization of large (Y,La)TiO 3 and (Y,Ca)TiO 3 single crystals

The Mott-insulating rare-earth titanates ( R TiO 3 , with R being a rare-earth ion) are an important class of materials that encompasses interesting spin-orbital phases as well as ferromagnet-antiferromagnet and insulator-metal transitions. The growth of these materials has been plagued by difficulties related to overoxidation, which arises from a strong tendency of Ti 3 + to oxidize to Ti 4 + . We describe our efforts to grow sizable single crystals of YTiO 3 , Y 1 - x La x TiO 3 ( x ≤ 0.25 ), and Y 1 - y Ca y TiO 3 ( y ≤ 0.35 ) with the optical traveling-solvent floating-zone technique. Finally, we present sample characterization via chemical composition analysis, magnetometry, charge transport, neutron scattering, x-ray absorption spectroscopy, and x-ray magnetic circular dichroism to understand macroscopic physical property variations associated with overoxidation. Furthermore, we demonstrate a good signal-to-noise ratio in inelastic magnetic neutron scattering measurements of spin-wave excitations. A superconducting impurity phase, found to appear in Ca-doped samples at high doping levels, is identified as TiO.

36 MATERIALS SCIENCE↗
DOE OSTI2024

SOURCES4D

SOURCES is a code for computing neutron source rates and spectra from spontaneous fission (including delayed neutrons) and (alpha,n) reactions in homogeneous materials and (alpha,n) reactions in single-interface and two-interface geometries. SOURCES is a Los Alamos National Laboratory (LANL) code that is written in FORTRAN and distributed through the Radiation Safety Information Computation Center (RSICC). LANL’s last release of SOURCES to RSICC was SOURCES4C in 2002. This disclosure covers the latest version of SOURCES, SOURCES4D. This version adds sensitivity capabilities for (alpha,n) sources in homogeneous materials. Specifically, SOURCES4D writes new output that can be used to calculate, in post-processing, first and second derivatives of the (alpha,n) source rate density and spectrum with respect to nuclide densities in a homogeneous material and first derivatives of the (alpha,n) source rate density and spectrum with respect to nuclide stopping powers and (alpha,n) cross sections (nuclear data). These derivatives are useful for uncertainty quantification, predictive modeling, and other applications in neutron transport problems.

Favorite, Jeffrey A.↗
DOE OSTI2025

All Particle In Fission Reactor - Energy Deposition

This product includes software developed by Members of the Geant4 Collaboration (http://cern.ch/geant4). The basic principle of ALFRED consists of a k-eigenvalue module in which all generated particles are tracked and all deposited energy is accounted for. An eigenvalue module updates the neutron source after each run based on the neutrons emitted at each fission in the previous run. As a result, the source distribution converges to the fundamental mode of the steady-state eigenvalue problem of the associated critical reactor. ALFRED leverages the High Precision neutron transport package.

Ferney, PaulA. [Idaho National Laboratory (INL), I↗
DOE OSTI2024

pymcnp v1

This is a software package that makes it easier to handle MCNP, a neutron transport simulation toolkit, simulations. It helps with reading and writing output from MCNP and modifying the input files programmatically. Furthermore, it also allows the visualization of input geometry and of the results. Pymcnp provides access to these features using the Python programming language.

Persaud, Arun↗
DOE OSTI2022

Coupling coarse-mesh CFD with fine-mesh CFD for modeling molten-salt reactors in the Virtual Test Bed (VTB)

The Nuclear Energy Advanced Modeling and Simulation (NEAMS) program aims at developing a simulation tool kit to accelerate the development and deployment of nuclear power technologies. NEAMS multiphysics tools have been designed to provide numerical simulation support for the design and licensing of GEN IV reactors. Pronghorn is NEAMS's coarse-mesh computational fluid dynamics (CFD) tool, which is designed to run 3D core transients in GEN IV reactors at a reduced computational cost. To increase their accuracy, coarse-mesh CFD simulations require calibrated closure coefficients. One way of computing these coefficients is via the Nek5000, NEAMS's high-fidelity CFD tool. This article discusses our current research lines in informing Pronghorn closure coefficients via Nek5000 to enable multiphysics simulations of the core cavity of the molten-salt fast reactor. We present an application in which Pronghorn mixing length turbulent viscosity has been calibrated from Nek5000 simulations. The resulting Pronghorn thermal-hydraulics model is then coupled to Griffin, the NEAMS neutron transport solver, to solve for the steady-state configuration of the molten-salt fast reactor. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Characterization of Uranium Foil Irradiations at the WSU TRIGA Reactor Using a New Reactor Model in SCALE [Slides]

A depleted uranium foil and a highly enriched uranium foil that were irradiated at the Washington State University TRIGA reactor were used to measure the relative production of 237 U and 239 U, two key isotopes of interest in nuclear forensics, as well as their decay daughters 237 Np and 239 Pu. Total fissions in the two foils were also estimated using 99 Mo measurement. A reactor model of the TRIGA was built in the SCALE neutron transport code using dimensions obtained from an existing MCNP geometry. The model was successfully used to simulate foil fissions and the reaction rate ratios in 238 U, using the laboratory-measured isotope ratios as a benchmark. The flux distribution in the foils was determined. Finally, using the continuous energy Tsunami-3d module of SCALE, a sensitivity/uncertainty analysis was performed to determine the effects on the reaction ratio caused by other reactions in the reactor and by uncertainties in the known cross-sections.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2020

Characterization of Uranium Foil Irradiations at the WSU TRIGA Reactor Using a New Reactor Model in Scale

A new reactor model of the Washington State University TRIGA was developed in the SCALE neutron transport code, and its fidelity was verified by comparison to MCNP and available data for several reactor parameters. The model was used to characterize irradiations designed to produce the short-lived actinides 237 U and 239 U, two key isotopes for nuclear forensics. These short-lived actinides, their decay daughters 237Np and 239Pu, and total fissions (via 99 Mo) were measured in irradiated foils at Los Alamos National Laboratory and other labs with good agreement among parent/daughter pairs and among labs. The laboratory-measured isotope ratios were used as a benchmark for the model determination of reaction rates. The flux distribution in the foils was also determined. Finally, using the continuous energy TSUNAMI-3D module of SCALE, a sensitivity/uncertainty analysis was performed to determine the effects on foil reactions caused by other system-wide reactions and by uncertainties in the known cross-sections.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Integration of Kokkos into MonteRay [Slides]

Monte Carlo Neutron Transport codes are C++-based and simulate interaction of nuclear particles with materials. MonteRay is a library for accelerating Monte Carlo ray-casting tallies with GPUs. Kokkos is commonly used to write performance portable applications for HPC platforms. A goal is to implement Kokkos into the Expected Path Length file that uses Cuda as a backend. The methodology and results are summarized.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

High Flux Isotope Reactor Low-Enriched Uranium Low Density Silicide Fuel Design Parameters

High Flux Isotope Reactor (HFIR) highly enriched uranium (HEU) to low-enriched uranium (LEU) conversion activities are ongoing as part of the Department of Energy (DOE) National Nuclear Security Administration (NNSA)’s nuclear nonproliferation mission. Design activities studying the conversion of HFIR from HEU to LEU fuel explored different fuel design features and shapes with a low density uranium-silicide dispersion (U 3 Si 2 -Al) fuel, which has a uranium density of 4.8 gU/cm 3 . The goal of these studies is to generate several HFIR LEU fuel designs of varying fuel fabrication complexity that meet the current HEU performance metrics and safety requirements. The documented designs will serve as references for fuel fabrication and qualification activities. Recent advancements in modeling and simulation tools enable quick prototyping of fuel designs. Shift, a Monte Carlo neutron transport and depletion tool optimized for high-performance computing (HPC) architectures, is used for efficient fuel cycle and performance metrics calculations. The HFIR Steady State Heat Transfer Code (HSSHTC) is used to vet the thermal safety margin. Also, a new automation tool that connects all fuel design analysis steps, named Python HFIR Analysis and Measurement Engine (PHAME), has been developed to expedite the design study in an efficient and reproducible manner. Leveraging these tools, several candidate fuel designs were selected for varying fabrication complexity. This report provides design feature details for four selected HFIR LEU low density U 3 Si 2 -Al fuel designs and their corresponding performance and safety metrics. Nominal, best-estimate design parameters and irradiation conditions, including fission rate densities, power densities, heat fluxes, and cumulative fission densities are provided for candidate fuel designs relevant to framing irradiation experiments to support fuel qualification efforts. Simulations show that the low density U 3 Si 2 -Al, with design features to enhance safety, can meet HEU core performance metrics and safety requirements if the reactor power is increased from 85 MW (HEU) to 95 MW (LEU) and if the active fuel length is increased from 50.80 cm (HEU) to 55.88 cm (LEU).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗