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

The MP{sub N} method: a new angular discretization method based on piecewise polynomial interfaces fluxes

In transport calculations, it is well known how S{sub N} method is extremely inefficient in problems where the particle physics is dominated by streaming. The ray-effect eventually produced by the insufficient angular discretization, appears to be extremely persistent with respect to the refinement of the angular quadrature. The MP{sub N} method, that relies on continuous angular representation, offers a robust remedy to such an issue. MP{sub N} is based on the decomposition of the unit sphere into solid angles and on a piecewise continuous definition of interface fluxes, which are expanded in polynomials in each solid angle. This allows propagating more than one angular degree of freedom simultaneously while maintaining unaltered the block-diagonal pattern of the displacement plus removal operator. The method is therefore well suited for the flux resolution by means of a conventional sweep algorithm. Furthermore, unlike the S{sub N} method, MP{sub N} does not rely on discrete directions and, thus, on an angular quadrature formula, but rather constructs a set of linear equations solving for the angular moments of the flux for all discrete solid angles within the sweep. MP{sub N} shows an error convergence rate higher than S{sub N} at the expense of an increased size of the coefficient matrices, so of the computational cost. Although MP{sub N} is not free from ray-effect, the latter is effectively mitigated and less persistent with respect to the increase of the angular refinement order. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Atomic dynamics in 𝑀⁢Cr⁢𝑋 2 (𝑀=Ag, Cu; 𝑋 = S, Se) across magnetic and superionic transitions

Here, a systematic study of atomic dynamics and thermal properties of the family of layered chalcogenide compounds 𝑀⁢Cr⁢𝑋 2 (𝑀= Ag, Cu; 𝑋 = S, Se) was performed, including neutron and x-ray scattering, thermal characterization, and first-principles simulations. In all compounds, we observe a breakdown of specific phonon modes across the superionic phase transition, for phonons whose eigenvectors exhibit large contributions of mobile ions. In particular, the nondispersive portions of transverse acoustic (TA) branches at short-wavelengths and the low-energy optical phonons with large contributions from Ag + or Cu + become severely damped in the superionic phase. However, well-defined quasiparticles persist in the superionic state for long-wavelength TA phonons. In the case of AgCrS 2 , the coupling of lattice dynamics with its antiferromagnetic transition was also investigated. The magnetic ordering couples with the monoclinic–rhombohedral structural transition, and the Cr 3+ spin arrangement strongly affects the phonon dispersions. We qualitatively reproduce the magnetic and nuclear components of the INS measurement for antiferromagnetic AgCrS 2 by combining models of spin-waves and spin-polarized first-principles phonon simulations. Quasielastic magnetic fluctuations persist in the paramagnetic phase up to high temperature, but are clearly distinguished from the nuclear component through their momentum dependence. Finally, we report measurements of the thermal properties of the selenide compounds and find good agreement with our DFT simulations.

36 MATERIALS SCIENCE↗
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 OSTI2020

Steady detonation propagation in thin channels with strong confinement

We examine asymptotically the dynamics of two-dimensional, steady detonation wave propagation and failure for a strongly confined high explosive (HE), in which the width of the explosive is small relative to the reaction zone length. An energy balance equation is derived, which shows how the longitudinal acceleration of subsonic flow behind the detonation shock is influenced both by chemical reaction and by the effects of HE boundary streamline deflection, specifically via the induced rate of change of mass flux through the detonation wave. The latter serves to either counteract or reinforce the acceleration of longitudinal flow, depending on the sign of the gradient of the boundary streamline deflection at the detonation shock. The analysis is valid for general equations of state and chemical reaction rates in the HE. The asymptotically derived form of the energy equation represents an eigenvalue problem for the determination of the steady detonation propagation speed, solved via a shooting method. We explore specific results for ideal and stiffened equations of state, along with a pressure-dependent reaction rate for which changes in the pressure exponent and reaction order are also studied. We consider the influences of both straight and curved HE boundary streamline shapes. The asymptotic analysis reveals significant physical insights into how detonation propagation and failure are affected by strong confinement.

42 ENGINEERING↗
DOE OSTI2020

Sensitive neutron transverse polarization analysis using a 3 He spin filter

Here we report an experimental implementation for neutron transverse polarization analysis that is capable of detecting a small angular change ($\ll$10 –3 rad) in neutron spin orientation. This approach is demonstrated for monochromatic beams, and we show that it could be extended to polychromatic neutron beams. Our approach employs a 3 He spin filter inside a solenoid with an analyzing direction perpendicular to the incident neutron polarization direction. The method was tested with polarized neutron beams and a spin rotator placed inside a μ -metal shield just upstream of the analyzer. No cryogenic superconducting shields or additional neutron spin manipulations are needed. With a counting detector, we experimentally show that the angular resolution δ θ = 1 / ( P n A N ) rad is only determined by the counting statistics for the total counts N and the product of the neutron polarization P n and the analyzing power A . With a high-flux neutron beam, 10 –6 rad angular sensitivity is feasible within a day. This simple, classical-quantum-limited transverse polarization analysis scheme may reduce the overall complexity of experimental implementation for applications requiring sensitive neutron polarimetry and improve the precision in fundamental science studies and polarized neutron imaging.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Resonance Scattering Treatment with the Windowed Multipole Formalism

A new method for directly sampling the resonance upscattering effect is presented. Alternatives have relied on inefficient rejection sampling techniques or large tabular storage of relative velocities. None of these approaches, which require pointwise energy data, are particularly well suited to the windowed multipole cross-section representation. The new method, called multipole analytic resonance scattering, overcomes these limitations by inverse transform sampling from the target relative velocity distribution where the cross section is expressed in the multipole formalism. The closed-form relative speed distribution contains a novel special function we deem the incomplete Faddeeva function, and we present the first results on its efficient numerical evaluation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Collinear antiferromagnetic order in a quasi-two-dimensional triangular lattice compound DyNiAl 4 Ge 2

The two-dimensional magnetic triangular lattice, with rare-earth local moments coupled through long-range interactions and hosting strong spin-orbit coupling, serves as an ideal platform for exploring novel quantum phenomena. Here, we report the physical properties of DyNiAl 4 Ge 2 single crystals through comprehensive structural, magnetic, heat capacity, electrical transport, and neutron powder diffraction measurements. DyNiAl 4 Ge 2 crystallizes in a trigonal lattice with the space group $R\bar{3}m$, where Dy atoms form two-dimensional triangular lattice layers within the 𝑎𝑏 plane and stack along the 𝑐 axis. Magnetic susceptibility measurement reveals a magnetic phase transition at 𝑇 𝑁 = 9.0 K, further confirmed by the 𝜆-shaped peak in the heat capacity curve. The obtained temperature-magnetic field phase diagram is relatively simple compared to other isostructural compounds, suggesting the absence of magnetic frustration. This is further evidenced by neutron powder diffraction, which shows a collinear antiferromagnetic ground state with a propagation vector 𝒌 = (0, 0, 1.5), implying the absence of magnetic frustration in DyNiAl 4 Ge 2 . In this magnetic structure, Dy moments are aligned ferromagnetically within each triangular layer and arranged antiferromagnetically between adjacent layers along the 𝑐 axis. In conclusion, our results establish DyNiAl 4 Ge 2 as a model system for understanding the interplay between RKKY interaction and geometrical frustration in this family of compounds.

Wu, HengHeng [Chinese Academy of Sciences (CAS), S↗
DOE OSTI2021

Tuning magnetism and band topology through antisite defects in Sb-doped MnBi 4 ⁢Te 7

The fine control of magnetism and electronic structure in a magnetic topological insulator is crucial in order to realize various novel magnetic topological states including axion insulators, magnetic Weyl semimetals, Chern insulators, etc. Through crystal growth, transport, thermodynamic, neutron diffraction measurements, we show that under Sb doping the newly discovered intrinsic antiferromagnetic (AFM) topological insulator MnBi 4 Te 7 evolves from AFM to ferromagnetic (FM) and then ferrimagnetic. We attribute this to the formation of Mn(Bi,Sb) antisites upon doping, which results in additional Mn sublattices that modify the delicate interlayer magnetic interactions and cause the dominant Mn sublattice to go from AFM to FM. We further investigate the effect of antisites on the band topology using the first-principles calculations. Without considering antisites, the series evolves from AFM topological insulator (x=0) to FM axion insulators. In the exaggerated case of 16.7% of periodic antisites, the band topology is modified and a type-I magnetic Weyl semimetal phase can be realized at intermediate dopings. Therefore, this doping series provides a fruitful platform with continuously tunable magnetism and topology for investigating emergent phenomena, including quantum anomalous Hall effect, Fermi arc states, etc.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Interplay of itinerant electrons and Ising moments in a hybrid honeycomb quantum magnet TmNi 3 Al 9

The interplay between itinerant electrons and local magnetic moments in quantum materials brings about rich and fascinating phenomena and stimulates various developments in the theoretical framework. In this work, thermodynamic, electric transport, and neutron diffraction measurements were performed on a newly synthesized honeycomb lattice magnet TmNi 3 Al 9 . Based on the experimental data, a magnetic-field–temperature phase diagram was constructed, exhibiting three essentially different magnetic regions. Below T N =2.97±0.02KTm 3+ , moments order antiferromagnetically in zero field. We found that the Tm 3+ ions form a pseudodoublet ground state with the Ising-like moments lying normal to the two-dimensional honeycomb layers. Application of a magnetic field along the easy axis gradually suppresses the antiferromagnetic order in favor of an induced ferromagnetic state above the critical field B c =0.92±0.05T. In the vicinity of B c , a strong enhancement of the quantum spin fluctuations was observed. The quantum Ising nature of the local moments and the coupling to itinerant electrons are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Realization of an intrinsic ferromagnetic topological state in MnBi 8 Te 13

Novel magnetic topological materials pave the way for studying the interplay between band topology and magnetism. However, an intrinsically ferromagnetic topological material with only topological bands at the charge neutrality energy has so far remained elusive. Using rational design, we synthesized MnBi 8 Te 13 , a natural heterostructure with [MnBi 2 Te 4 ] and [Bi 2 Te 3 ] layers. Thermodynamic, transport, and neutron diffraction measurements show that despite the adjacent [MnBi 2 Te 4 ] being 44.1 Å apart, MnBi 8 Te 13 manifests long-range ferromagnetism below 10.5 K with strong coupling between magnetism and charge carriers. First-principles calculations and angle-resolved photoemission spectroscopy measurements reveal it is an axion insulator with sizable surface hybridization gaps. Our calculations further demonstrate the hybridization gap persists in the two-dimensional limit with a nontrivial Chern number. Therefore, as an intrinsic ferromagnetic axion insulator with clean low-energy band structures, MnBi 8 Te 13 serves as an ideal system to investigate rich emergent phenomena, including the quantized anomalous Hall effect and quantized magnetoelectric effect.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Competing magnetic phases and fluctuation-driven scalar spin chirality in the kagome metal YMn 6 Sn 6

Identification, understanding, and manipulation of novel magnetic textures are essential for the discovery of new quantum materials for future spin-based electronic devices. In particular, materials that manifest a large response to external stimuli such as a magnetic field are subject to intense investigation. Here, we study the kagome-net magnet YMn 6 Sn 6 by magnetometry, transport, and neutron diffraction measurements combined with first-principles calculations. We identify a number of nontrivial magnetic phases, explain their microscopic nature, and demonstrate that one of them hosts a large topological Hall effect (THE). We propose a previously unidentified fluctuation-driven mechanism, which leads to the THE at elevated temperatures. This interesting physics comes from parametrically frustrated interplanar exchange interactions that trigger strong magnetic fluctuations. Our results pave a path to chiral spin textures, promising for novel spintronics.

36 MATERIALS SCIENCE↗
DOE OSTI2021

OpenMC

OpenMC is an open source, community-developed framework for performing Monte Carlo particle transport calculations.It is capable of simulating neutron and photon transport below ~50 MeV and also includes a model for bremsstrahlung production of photons by electrons/positrons. Fixed source, k-eigenvalue, and subcritical multiplication problems can be solved. Problem geometries can be modeled either using constructive solid geometry or a CAD representation. A flexible and efficient tally system enables a wide variety of physical quantities to be tallied and analyzed. OpenMC can run in parallel using a hybrid MPI and OpenMP programming model and has been extensively tested on leadership class supercomputers. A built-in Bateman equation solver enables the modeling of material composition changes due to irradiation. One of the unique features of OpenMC is its rich, extensible Python and C/C++ programming interfaces that enable programming pre- and post-processing, multigroup cross section generation, workflow automation, multiphysics coupling, and the visualization of geometry and tally results. In addition to the core Monte Carlo transport solver and associated APIs, OpenMC includes a Python-based nuclear data interface that enables power users to inspect, modify,and perform various types of nuclear data processing on ENDF, ACE, and OpenMC¿s native HDF5 files

ROMANO, PaulK.↗
DOE OSTI2026

MCNPy

SAND2026-20425O MCNPy runs and analyzes simulations from MCNP, a software that models radiation transport of neutrons and gamma rays. MCNPy uses Python to start MCNP, retrieve event data files, and convert them into graph structures for detailed analysis. It offers visualization tools, including 2D views of particle histories, making complex simulation data easier to interpret for researchers and engineers. 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 OSTI2026

ndi

The Nuclear Data Interface (NDI) is an application programming interface (API) that allows access to standard nuclear data parameters while hiding the underlying details of the data libraries and their storage. It allows access to multigroup transport data (neutron and gamma), thermonuclear burn data, dosimetry data, production/depletion chain data, radiochemistry data, and secondary neutron multiplicity data. The name NDI refers to both the code and data formats supported by the code.

Saller, Thomas↗
DOE OSTI2022

Demonstration of MOOSE-based Griffin reactor physics, code for heterogeneous lead-cooled fast reactor analysis

The MOOSE-based reactor physics code Griffin was assessed on a heterogeneous pin-resolved model of a prototype lead-cooled fast reactor assembly. This model was developed in preparation for future use in MOOSE-based multiphysics calculations for computing hot channel factors. Heterogeneous multigroup cross sections were prepared using the fast reactor multi-group cross section processing code MC{sup 2}-3 using a two-step method. Griffin simulations were performed using the DFEM-SN solver on 576 cores on Argonne's LCRC cluster. Diffusion-based acceleration methods were applied (NDA and CMFD). Reference solutions were generated with continuous energy MCNP and the hybrid MOC/finite element solver PROTEUS-MOC for code-to-code comparison. Space-angle convergence studies were conducted to observe convergence in k-eigenvalue and axial pin power distributions. The fully resolved Griffin calculation was within 68 pcm of the MCNP eigenvalue and exhibited max 1.2% relative error in the axial pin power distribution. Griffin produced nearly identical results to PROTEUS-MOC when using the same 9-group multigroup cross-section set. Griffin demonstrated favorable scaling in wall-clock time and memory usage when using diffusion-based acceleration methods. Griffin is capable of simulating the pin-resolved heterogeneous LFR assembly with good accuracy and performance, and is suitable for future use in coupled high-fidelity hot channel factor simulations. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Validation of MPACT BWR capabilities against critical experiments

Over the years, significant validation work for the neutronics code MPACT has been performed against zero-power critical benchmarks and measured data from operating nuclear power plants. As a part of the Modeling and Analysis of Exelon BWRs for Eigenvalue and Thermal Limits Predictability project, new validation efforts relevant to boiling water reactor (BWR) core applications have been performed. This paper presents the results of the critical experiment portion of the BWR validation efforts for MPACT. The Kritz-4 experiments and IPEN/MB-01 BWR-relevant configurations are modeled with MPACT. The Kritz-4 BWR critical experiments were performed at both cold and hot conditions, which is a unique feature among other critical experiment facilities. The MPACT results with linear source (LS) Method of Characteristics (MOC) have very good agreement with the measured criticality. When using the 60-group BWR library with LS and P2 scattering, the maximum k{sub eff} error is 109 pcm, and the average cold-hot bias is only -0.11 pcm/K. The fission rate distributions of MPACT are also verified with Serpent calculations. For the IPEN facility, two configurations are considered: one with a large central void, and one with a cruciform control rod. The MPACT k{sub eff} errors are less than 50 pcm for all IPEN problems. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Stereolithographic geometry model of the IBR-2M experimental facility

The IBR-2M is a fast research reactor that operates in supercritical condition for ∼ 800 μs every 200 ms. Two reflector parts in nickel rotate in opposite directions generating 1.8 GWth peak power when they align with the fuel zone changing the reactor status from deep subcritical to supercritical. The reactor core uses high-enriched plutonium fuel and is cooled by sodium. This reactor has been modeled by MCNP and SERPENT computer programs. The MCNP model uses combinatorial geometry, whereas the SERPENT model employs the Stereolithographic (STL) geometry representation that can be used by 3D printers. The STL geometry was constructed using the CUBIT computer program. The CUBIT program was also used for a three-dimensional visualization of the Monte Carlo models. SERPENT and MCNP models use the same geometry, material specifications, and nuclear data. The latter are based on the ENDF/B-7.0 library. SERPENT and MCNP using same geometry and same material specifications produce similar k{sub eff} values within 120 pcm.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Statistical uncertainty of fission matrix eigenvalues using perturbation theory

Eigenvalue search of high dominance ratio systems may be slow to converge. The fission matrix element is defined by its element (FM){sub ij}, which are the probability for a neutron born in cell i to create a fission in cell j for a spatial mesh of n{sub i}*n{sub j} elements. Fission matrices are used in Monte Carlo criticality simulations to enhance computing speed, but also to find higher order eigenvalues. However, few studies have been made on the link between statistical uncertainties of fission matrix elements and eigenvalues uncertainties. Thus, dominance ratio statistical uncertainties remain unknown. This paper uses a new generalized perturbation theory (GPT) method to estimate sensitivities of eigenvalues to fission matrix elements and then to calculate dominance ratio uncertainties.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗