SEARCH · Engineering Papers

Results for “Neutron transport equation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Implementation of Monte Carlo Transport in the General Relativistic SpEC Code

Neutrino transport and neutrino-matter interactions are known to play an important role in the evolution of neutron star mergers and of their post-merger remnants. Neutrinos cool remnants, drive post-merger winds, and deposit energy in the low-density polar regions where relativistic jets may eventually form. Neutrinos also modify the composition of the ejected material, impacting the outcome of nucleosynthesis in merger outflows and the properties of the optical/infrared transients that they power (kilonovae). So far, merger simulations have largely relied on approximate treatments of the neutrinos (leakage, moments) that simplify the equations of radiation transport in a way that makes simulations more affordable but also introduces unquantifiable errors in the results. To improve on these methods, we recently published a first simulation of neutron star mergers using a low-cost Monte Carlo algorithm for neutrino radiation transport. Our transport code limits costs in optically thick regions by placing a hard ceiling on the value of the absorption opacity of the fluid, yet all approximations made within the code are designed to vanish in the limit of infinite numerical resolution. We provide here an in-depth description of this algorithm, of its implementation in the SpEC merger code, and of the expected impact of our approximations in optically thick regions. We argue that the last is a subdominant source of error at the accuracy reached by current simulations and for the interactions currently included in our code. We also provide tests of the most important features of this code.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

On the novel 3-D neutron transport kinetic tRAPID algorithm and its validation

The Real-time Analysis for Particle-transport and In-situ Detection (RAPID) Code System, based on the Multi-stage Response-function Transport (MRT) methodology, allows for real-time simulation of nuclear systems based on 3-D continuous-energy particle transport. RAPID's steady-state (criticality) neutron transport algorithm is based on the Fission Matrix (FM) method, and has been extensively verified and validated against computational benchmarks and experiments. This paper introduces the novel 3-D time-dependent transport algorithm that has been implemented into the code, tRAPID, and its validation using the JSI TRIGA Mark-II reactor. tRAPID accurately and efficiently calculates neutron kinetics parameters (such as β{sub eff}, l{sub eff} , Λ, α{sub Rossi}) and 3-D time-dependent neutron fission source distribution and neutron importances for both prompt and delayed neutrons. tRAPID is used to simulate a rod insertion experiment performed at the JSI TRIGA Mark-II reactor, during which signals from four fission chambers at four different locations in the core were collected. The results demonstrate how tRAPID is capable of calculating detailed and accurate results with only a minimal use computational resources and time.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

Characterization of the inner edge of the neutron star crust

The poorly known crustal equation of state plays a critical role in many observational phenomena associated with a neutron star. Here using semiclassical Monte Carlo simulations, we explore the possible configurations of the inner edge of the neutron star crust for a variety of baryon densities and proton fractions. Applying the Kirkwood-Buff theory to these two-component systems, we observe how the isothermal compressibility reaches a maximum when isolated nonsymmetric (or quasispherical) clusters are formed in an extremely dilute neutron gas. To determine the neutron fraction, we suggest a geometrical model based on the behavior of the proton-neutron pair correlation function. Accordingly, the equation of state of the inner crust is calculated, illustrating that the nuclear energy in β equilibrium follows a power-law behavior with baryon density. As a possible astrophysical outcome of this study, our results could help refine the mass-radius relation. Finally, our results pave the way towards further investigations of the impact of the proton-neutron pair correlation function on transport properties within the neutron star crust.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Impact of Magnetic-field-driven Anisotropies on the Equation of State Probed in Neutron Star Mergers

Binary neutron star mergers can produce extreme magnetic fields, some of which can lead to strong magnetar-like remnants. While strong magnetic fields have been shown to affect the dynamics of outflows and angular momentum transport in the remnant, they can also crucially alter the properties of nuclear matter probed in the merger. In this work, we provide a first assessment of the latter, determining the strength of the pressure anisotropy caused by Landau-level quantization and the anomalous magnetic moment. To this end, we perform the first numerical relativity simulation with a magnetic polarization tensor and a magnetic-field-dependent equation of state using a new algorithm we present here, which also incorporates a mean-field dynamo model to control the magnetic field strength present in the merger remnant. Our results show that—in the most optimistic case—corrections to the anisotropy can be in excess of 10% and are potentially largest in the outer layers of the remnant. This work paves the way for a systematic investigation of these effects.

General relativity↗

DOE OSTI2023

Neutrino transport in general relativistic neutron star merger simulations

Abstract Numerical simulations of neutron star–neutron star and neutron star–black hole binaries play an important role in our ability to model gravitational-wave and electromagnetic signals powered by these systems. These simulations have to take into account a wide range of physical processes including general relativity, magnetohydrodynamics, and neutrino radiation transport. The latter is particularly important in order to understand the properties of the matter ejected by many mergers, the optical/infrared signals powered by nuclear reactions in the ejecta, and the contribution of that ejecta to astrophysical nucleosynthesis. However, accurate evolutions of the neutrino transport equations that include all relevant physical processes remain beyond our current reach. In this review, I will discuss the current state of neutrino modeling in general relativistic simulations of neutron star mergers and of their post-merger remnants. I will focus on the three main types of algorithms used in simulations so far: leakage, moments, and Monte-Carlo scheme. I will review the advantages and limitations of each scheme, as well as the various neutrino–matter interactions that should be included in simulations. We will see that the quality of the treatment of neutrinos in merger simulations has greatly increased over the last decade, but also that many potentially important interactions remain difficult to take into account in simulations (pair annihilation, oscillations, inelastic scattering).

Foucart, Francois (ORCID:0000000346174738)↗

DOE OSTI2025

Single Grid Error Estimation for Neutron Transport Solvers

The method of nearby problems (MNP) is a solution verification technique that does not require the use of multiple spatial grids. To estimate spatial discretization error without requiring a high-fidelity spatial grid, an analytical curve fit is interpolated from the numerical solution. The residual between the curve fit solution and numerical solution is calculated and added as an additional source term to the governing equation. The nearby solution is estimated using the updated source term and boundary conditions to remain consistent with the curve fit interpolation. The nearby solution can be compared to the curve fit solution as a discretization error estimation while using a single spatial grid. Without the use of higher fidelity spatial grids, the MNP is able to approximate the spatial discretization error, a facet of solution verification. The application of the method of nearby problems is presented for one- and two-dimensional neutron transport problems for both fixed source and criticality problems on the spatial variable. The fixed source results demonstrate the effectiveness of nearby problems for spatial error identification using the discrete ordinates method. Criticality results are shown to identify area of high spatial error for the C5G7 problem as well as for the discrete ordinates solver. A novel approach of combining the capabilities of Monte Carlo with the discrete ordinates nearby problems is presented for one- and two-dimensional fixed source problems. In conclusion, the MNP demonstrates its effectiveness at identifying spatial error on a single structured grid with a wide variety of neutron transport problems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Impact of moment-based, energy integrated neutrino transport on microphysics and ejecta in binary neutron star mergers

We present an extensive study of the effects of neutrino transport in three-dimensional general relativistic radiation hydrodynamics (GRHD) simulations of binary neutron star (BNS) mergers using our moment-based, energy-integrated neutrino radiation transport (M1) scheme. Here, we consider a total of eight BNS configurations, while varying equation of state models, mass ratios, and grid resolutions, for a total of 16 simulations. We find that M1 neutrino transport is crucial in modeling the local absorption of neutrinos and the deposition of lepton number throughout the medium. We provide an in-depth look at the effects of neutrinos on the fluid dynamics and luminosity during the late inspiral and postmerger phases, the properties of ejecta and outflow, and the postmerger nucleosynthesis. The simulations presented in this work comprise an extensive study of the combined effect of the equation of state and M1 neutrino transport in GRHD simulations of BNS mergers, and establish that the solution provided by our M1 scheme is robust across system properties.

150 ≤ A ≤ 18959 ≤ A ≤ 8990 ≤ A ≤ 149↗

DOE OSTI2023

Progress Towards the Reality of Low-Temperature Critical Experiments with Low-Temperature TEX

Integral benchmarks for criticality safety do not exist below room temperature. Incredible support for below room temperature benchmarks is present across the international community. Validation of low-temperature neutron cross sections is necessary to verify any operation at temperatures below room temperature which is typically observed in environments far from the equator. For example, a fissile material transportation truck may routinely observe ambient temperatures down to -40°C, which is the lower temperature bound of the normal conditions of transportation defined in the United States Title 10 Code of Federal Regulations §71.71c2. Similarly, the IAEA, which contains 176 member states, has the same requirement for subcriticality for fissile material transportation. Additionally, sub-room temperature benchmarks can validate newly produced cross sections, that include novel thermal scattering laws, from North Carolina State University. A six-case experiment designed based off the TEX-HEU baseline experiments has been proposed utilizing a chilled vacuum chamber capable of cooling the HEU stack down to -50°C. The six proposed configurations span the entire neutron fission energy range and can be tested at intermediate temperatures between room temperature and -40°C. Of the six configurations, the fastest case utilizes nearly 135 kg of HEU. Understanding the system is required prior to critical experiment operations. In order to test the functionality of the system, testing with surrogate fuel materials has been conducted. The following discusses the design of the low temperature TEX (LT-TEX) experiment as well as the progress made via surrogate testing towards a reality in which low-temperature integral benchmarks exist in the International Criticality Safety Benchmark Evaluation Project (ICSBEP).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Mapping dynamical ejecta and disk masses from numerical relativity simulations of neutron star mergers

Abstract We present fitting formulae for the dynamical ejecta properties and remnant disk masses from the largest to date sample of numerical relativity simulations. The considered data include some of the latest simulations with microphysical nuclear equations of state (EOS) and neutrino transport as well as other results with polytropic EOS available in the literature. Our analysis indicates that the broad features of the dynamical ejecta and disk properties can be captured by fitting expressions, that depend on mass ratio and reduced tidal parameter. The comparative analysis of literature data shows that microphysics and neutrino absorption have a significant impact on the dynamical ejecta properties. Microphysical nuclear EOS lead to average velocities smaller than polytropic EOS, while including neutrino absorption results in larger average ejecta masses and electron fractions. Hence, microphysics and neutrino transport are necessary to obtain quantitative models of the ejecta in terms of the binary parameters.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

The Legendre Polynomial Axial Expansion Method

This work presents a new formulation of the axial expansion transport method explicitly using Legendre polynomials for arbitrarily high-order expansions. This new formulation also features an alternative method of axial leakage calculation to allow for nonextruded flat source region meshes. This alternative axial leakage is introduced alongside a balance equation requirement to ensure that neutron balance is preserved in the coarse mesh for a given axial leakage formulation, which allows for effective coarse mesh finite difference acceleration. A matrix exponential table method is derived to allow for fast computations of arbitrarily high-order matrix exponentials for this work and precludes the need for further research into matrix exponential calculations for this method. Numerical results are presented that demonstrate the stability of the axial expansion method in systems with voidlike regions, showcase the speedup from matrix exponential tables, and investigate the axial convergence of the method in terms of both expansion order and mesh size.

Herring, Nicholas↗

NASA NTRS2000

An Improved Elastic and Nonelastic Neutron Transport Algorithm for Space Radiation

A neutron transport algorithm including both elastic and nonelastic particle interaction processes for use in space radiation protection for arbitrary shield material is developed. The algorithm is based upon a multiple energy grouping and analysis of the straight-ahead Boltzmann equation by using a mean value theorem for integrals. The algorithm is then coupled to the Langley HZETRN code through a bidirectional neutron evaporation source term. Evaluation of the neutron fluence generated by the solar particle event of February 23, 1956, for an aluminum water shield-target configuration is then compared with MCNPX and LAHET Monte Carlo calculations for the same shield-target configuration. With the Monte Carlo calculation as a benchmark, the algorithm developed in this paper showed a great improvement in results over the unmodified HZETRN solution. In addition, a high-energy bidirectional neutron source based on a formula by Ranft showed even further improvement of the fluence results over previous results near the front of the water target where diffusion out the front surface is important. Effects of improved interaction cross sections are modest compared with the addition of the high-energy bidirectional source terms.

Clowdsley, Martha S.↗

DOE OSTI2026

A Monte Carlo Laplace Transform Estimator for Radiation Transport

This work formulates and implements a Laplace transform estimator in a simple Monte Carlo radiation transport code. The estimator maps flux-based quantities of interest, like reaction rates, from a desired phase-space dimension to the complex Laplace domain. This on-the-fly Monte Carlo integration technique enables the spectral analysis of arbitrary nuclear systems via the Laplace transform. A simple code tests the estimator in neutron slowing-down problems across various infinite media, and the results compare well with Ganapol’s uninverted analytical solution of the neutron slowing-down equation.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Constraining the Equation of State of Dense Neutron-rich Matter (Final Report)

Using a new Time Projection Chamber (TPC), called the SpiRIT (SAMURAI pion Reconstruction Ion Tracker) constructed and funded by DOE office of science (award # DE-SC0004835), we performed two experiments with the SAMURAI spectrometer at RIKEN, Japan to study the equation of state of neutron rich matter. As a result of the project, the SpiRIT collaboration, an international collaboration consisting of groups from US, Japan, Korea, Poland, China and Germany, has been formed to pursue the science opportunities provided by the SpiRIT TPC. After the experiments, we developed the software to analyze the SpiRIT experiments and extract constraints of symmetry energy at supra-saturation densities. As a result 270 TB of data have been obtained. Over 10 technical papers on the TPC and 4 science papers on the nuclear equation of state have been published. A total of six (2 US, 1 Japanese and 1 Korean) PhD students have graduated based on their research on the EOS research using the SpiRIT data. In addition, the data also stimulated theorists on transport models to form the “Transport Model Evaluation Project” collaboration to improve the transport model to compare calculations with the data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

ONIX: An open-source depletion code

Open Source software enables innovative, community-based software development. ONIX brings this concept to the field of depletion calculations. It is an open-source depletion software to be used for nuclear reactor simulations, for fissile material production analysis as well as for nuclear arms control applications. ONIX provides a module to solve the depletion equation using a Chebyshev Rational Approximation Method. For the generation of one-group cross sections, it includes a coupling interface for the open-source neutron transport code, OpenMC, as well as a module to read pre-computed values in a stand-alone mode. ONIX has special features to optimize nuclear data libraries, to update isomeric branching ratio during burnup, and to support automation of simulations for nuclear archaeology. In conclusion, ONIX has been validated against results from numerical and experimental benchmarks, and its results agree with other methods within expected error ranges.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1993

Transport methods and interactions for space radiations

This report presents a brief history leading to the involvement of the Langley Research Center of the National Aeronautics and Space Administration (NASA) in space-radiation physics and protection. Indeed, a relatively complete summary of technical capability as of the summer of 1990 is given. The Boltzmann equations for coupled ionic and neutronic fields are presented and inversion techniques for the Boltzmann operator are discussed. Errors generated by the straight ahead approximation are derived and are shown to be negligible for most problems of space-radiation protection. A decoupling of projectile propagation from the target fields greatly simplifies the Boltzmann equations and allows an analytic solution of the target fragment transport. Analytic and numerical methods of solving the projectile transport equations are discussed. The nuclear physics underlying the coefficients in the Boltzmann equation is discussed. A coupled-channel optical model is found as a consequence of the loose binding of nuclear matter and closure of the nuclear states in high-energy reactions. Transport solutions with the developed data base are used with laboratory experiments to validate both the transport code and the data base. Numerical benchmarks and comparison with Monte Carlo calculations are also used for code validation.

Wilson, John W.↗

NASA NTRS1974

Extension of Case's method to the two-region critical cylinder and its application as an analytic standard

The normal mode expansion technique is applied to the transformed monoenergetic integral transport equation to develop a solution for the rotationally invariant and axially infinite, critical, two-region cylinder with a finite outer reflector boundary. Isotropic scattering and identical neutron mean free paths in the core and reflector regions are assumed. The solution in terms of singular integral equations is obtained by applying a completeness theorem found for the singular eigenfunctions. Numerical results for a variety of core and reflector multiplying properties and reflector thickness are presented and compared with the results of other method. An example of this type of application is given in a study of approximations inherent in the neutronic design analysis of a small, fast-neutron-spectrum reactor concept proposed as a space power source.

Westfall, R. M.↗

NASA NTRS1993

Quasiparticle Interactions in Neutron Matter for Applications in Neutron Stars

A microscopic model for the quaisiparticle interaction in neutron matter is presented. Both particle-particle (pp) and particle-hole (ph) correlation are are included. The pp correlations are treated in semi-empirical way, while ph correlations are incorporated by solving coupled two-body equations for the particle hole interaction and the scattering amplitude on the Fermi sphere. The resulting integral equations self-consistently sum the ph reducible diagrams. Antisymmetry is kept at all stages and hence the forward-scattering sum rules are obeyed. Results for Landau parameters and transport coefficients in a density regime representing the crust of a neutron star are presented. We also estimate the S-1 gap parameter for neutron superfluidity and comment briefly on neutron-star implications.

Wambach, J.↗

DOE OSTI2021

Neutronics Calculation Advances at Los Alamos: Manhattan Project to Monte Carlo

The history and advances of neutronics calculations at Los Alamos during the Manhattan Project through the present are reviewed. Substantial improvements to neutron diffusion methods and the invention of both the Monte Carlo neutron transport methods in 1947 and deterministic discrete ordinates Sn in 1953 were all made at Los Alamos just after the Manhattan Project. We briefly summarize early simpler and more approximate neutronics methods and then describe the need to better predict neutronics behavior through consideration of theoretical equations, models and algorithms, experimental measurements, and available computing capabilities and their limitations. This paper briefly covers key advances in deterministic methods during the Manhattan Project. These capabilities, coupled with increasing postwar defense needs and the invention of electronic computing with the Electronic Numeric Integrator and Computer, known as ENIAC, and the Mathematical Analyzer Numerical Integrator and Automatic Computer Model, known as MANIAC, led to the creation of Monte Carlo and deterministic discrete ordinates neutronics transport methods. We note the important role that the scientific comradery between the Los Alamos scientists played in the process. This paper briefly covers the early methods, algorithms, computers, and electronic and women pioneers that enabled Monte Carlo to spread to all areas of science. We focus heavily on these early developments and the subsequent creation of the MCNP® code, advances in its associated nuclear data, and its applications to problems of national defense at Los Alamos.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗