Search NASA⌕ Search

SEARCH · Search NASA

Results for “neutron generation”

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 361 records · Page 20

Development of deterministic transport methods for low energy neutrons for shielding in space

Transport of low energy neutrons associated with the galactic cosmic ray cascade is analyzed in this dissertation. A benchmark quality analytical algorithm is demonstrated for use with BRYNTRN, a computer program written by the High Energy Physics Division of NASA Langley Research Center, which is used to design and analyze shielding against the radiation created by the cascade. BRYNTRN uses numerical methods to solve the integral transport equations for baryons with the straight-ahead approximation, and numerical and empirical methods to generate the interaction probabilities. The straight-ahead approximation is adequate for charged particles, but not for neutrons. As NASA Langley improves BRYNTRN to include low energy neutrons, a benchmark quality solution is needed for comparison. The neutron transport algorithm demonstrated in this dissertation uses the closed-form Green's function solution to the galactic cosmic ray cascade transport equations to generate a source of neutrons. A basis function expansion for finite heterogeneous and semi-infinite homogeneous slabs with multiple energy groups and isotropic scattering is used to generate neutron fluxes resulting from the cascade. This method, called the FN method, is used to solve the neutral particle linear Boltzmann transport equation. As a demonstration of the algorithm coded in the programs MGSLAB and MGSEMI, neutron and ion fluxes are shown for a beam of fluorine ions at 1000 MeV per nucleon incident on semi-infinite and finite aluminum slabs. Also, to demonstrate that the shielding effectiveness against the radiation from the galactic cosmic ray cascade is not directly proportional to shield thickness, a graph of transmitted total neutron scalar flux versus slab thickness is shown. A simple model based on the nuclear liquid drop assumption is used to generate cross sections for the galactic cosmic ray cascade. The ENDF/B V database is used to generate the total and scattering cross sections for neutrons in aluminum. As an external verification, the results from MGSLAB and MGSEMI were compared to ANISN/PC, a routinely used neutron transport code, showing excellent agreement. In an application to an aluminum shield, the FN method seems to generate reasonable results.

Ganapol, Barry↗

Long-term evolution of crustal neutron star magnetic fields

We have derived an analytic solution to the asymptotic behavior of dipolar magnetic fields that are generated in the crusts of neutron stars. We show that if the conductivity is due to impurity scattering, as expected for late stages of evolution, the surface field strength at the magnetic pole declines with the power law B(sub p) approximately = (t/t(sub 0))(exp -2/3). The results are shown to be qualitatively consistent with detailed numerical calculations. These latter results are consistent with some recent analyses of pulsar statistics and the magnetic fields of several binary pulsars with white dwarf companions whose ages have been determined. The dependence of the surface magnetic field on spin period of the pulsar is derived.

Urpin, V. A.↗

Neutron Tunneling: A New Mechanism to Power Explosive Phenomena in Neutron Stars, Magnetars, and Neutron Star Mergers

Neutron tunneling between neutron-rich nuclei in inhomogeneous dense matter encountered in neutron star crusts can release enormous energy on a short timescale to power explosive phenomena in neutron stars. In this work, we clarify aspects of this process that can occur in the outer regions of neutron stars when oscillations or cataclysmic events increase the ambient density. We use a time-dependent Hartree–Fock–Bogoliubov formalism to determine the rate of neutron diffusion and find that large amounts of energy can be released rapidly. The roles of nuclear binding, two-body interaction, and pairing in neutron diffusion times are investigated. We consider a one-dimensional quantum diffusion model and extend our analysis to study the impact of diffusion in three dimensions. Here, we find that these novel neutron transfer reactions can generate energy in the amount of ≃ 10 40 –10 44 erg under suitable conditions and assumptions.

79 ASTRONOMY AND ASTROPHYSICS↗

The MIGDAL experiment: Measuring a rare atomic process to aid the search for dark matter

We present the Migdal In Galactic Dark mAtter expLoration (MIGDAL) experiment aiming at the unambiguous observation and study of the so-called Migdal effect induced by fast-neutron scattering. It is hoped that this elusive atomic process can be exploited to enhance the reach of direct dark matter search experiments to lower masses, but it is still lacking experimental confirmation. Our goal is to detect the predicted atomic electron emission which is thought to accompany nuclear scattering with low, but calculable, probability, by deploying an Optical Time Projection Chamber filled with a low-pressure gas based on CF 4 . Initially, pure CF 4 will be used, and then in mixtures containing other elements employed by leading dark matter search technologies — including noble species, plus Si and Ge. High resolution track images generated by a Gas Electron Multiplier stack, together with timing information from scintillation and ionisation readout, will be used for 3D reconstruction of the characteristic event topology expected for this process — an arrangement of two tracks sharing a common vertex, with one belonging to a Migdal electron and the other to a nuclear recoil. Different energy-loss rate distributions along both tracks will be used as a powerful discrimination tool against background events. In this article we present the design of the experiment, informed by extensive particle and track simulations and detailed estimations of signal and background rates. In pure CF 4 we expect to observe 8.9 (29.3) Migdal events per calendar day of exposure to an intense D–D (D–T) neutron generator beam at the NILE facility located at the Rutherford Appleton Laboratory (UK). With our nominal assumptions, 5σ median discovery significance can be achieved in under one day with either generator.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Q1 Report for FY25 Theory and Simulation Performance Target: Development of an integrated modeling framework for fusion reactor design and assessment

This report describes the work and activities carried out towards the completion of each of the following milestones in FY25 Q1: 1. Demonstrate workflow for generating self-consistent CESOL plasma profiles + first wall and divertor loading prediction and generate the CAT plasma and neutron loading needed for further engineering analysis: $/circ$ Initially run CESOL with SOLPS and map heat flux using simple HEAT method: ▪ Generate SOLPS grid for CAT reference case, ▪ Implement and apply ‘lower’ fidelity method of using HEAT-like analytic method to map SOLPS heat flux from charged particles to first wall, and ▪ Provide initial heat flux + neutron loading for further engineering analysis. 2. Demonstrate multiphysics magnet analysis: $/circ$ Demonstrate validation of the Elmer workflow by comparing induced stresses due to EM forces generated by TF coils using multiple codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Neutronics Modeling of the pulsed plasma rocket reactor using rattlesnake

In the pulsed plasma rocket (PPR) reactor, neutron induced fission processes are utilized to implement a series of pulsed micro-explosions of very high power and eject plasma as a propellant. More specifically, projectiles (bullets) composed of moderated uranium are sent through the chamber of an unmoderated uranium barrel. By inducing rapid fission within the bullets, a plasma can be generated with appropriate delivery of neutrons. This study aims to assess the neutronics performance of the PPR reactor during normal operations using the MOOSE-based Rattlesnake code through the evaluation of the impact of the movement of the fuel bullet and the rotation of control drums on the criticality of the system. The Monte Carlo (MC) code Serpent 2 was employed to generate material-based cross sections for use in Rattlesnake and the reference neutronics solution. Cubit was used to generate the mesh for the Rattlesnake model. Parametric studies were conducted to evaluate the best approaches for cross section and mesh generation to ensure accurate results from Rattlesnake. As part of the verification process, the eigenvalue results of the system at various fuel bullet positions were obtained using Rattlesnake and compared with the reference solutions. The acceptable differences show that the Rattlesnake model with appropriate cross section and mesh generation procedure is a sufficiently accurate approximation of the continuous energy (CE) MC model for micro-sized reactors like the PPR reactor. Next, the impact of the rotation of the control drums on the system criticality was evaluated using the verified model. It was shown that the high-fidelity simulation using the deterministic code Rattlesnake can produce sufficiently accurate results for the evaluation of the reactor's neutronics performance in different phases of the normal operation with acceptable computational cost.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

ParaStell: parametric modeling and neutronics support for stellarator fusion power plants

The three-dimensional variation inherent to stellarator geometries and fusion sources motivates three-dimensional modeling to obtain accurate results from computational modeling in support of design and analysis of first wall, blanket, and shield (FWBS) systems. Manually constructing stellarator fusion power plant geometries in computer-aided design (CAD) and defining the corresponding fusion source can be cumbersome and challenging. The open-source parametric modeling toolset ParaStell has been developed to automate construction of such geometries in low-fidelity. Low-fidelity modeling is useful during the conceptual phase of engineering design as a means of rapidly exploring the design space of a given device. The modeling capability of ParaStell includes in-vessel components and magnets, for any given stellarator configuration, using a parametric definition and plasma equilibrium data. Furthermore, the toolset automates the generation of detailed, tetrahedral neutron source definitions and DAGMC geometries for use in neutronics modeling. ParaStell assists rapid design iteration, parametric study, and design optimization of stellarator fusion cores. As a demonstration of the design iteration capability, the effect of the three-dimensional parameter space on tritium breeding and magnet shielding is investigated, using the WISTELL-D configuration as a design basis. Blanket and shield thicknesses are varied in three dimensions, using the space available between the plasma edge and magnet coils as a constraint. The corresponding effects on tritium breeding ratio and magnet heating are tallied using the open-source Monte Carlo particle transport code OpenMC. The inclusion of additional and higher-fidelity modeling capabilities is planned for ParaStell’s future, as well as its implementation in machine-driven optimization.

Moreno, Connor↗

Mo-99 Concentration and processing by Solvent Extraction and Ion Exchange

Normally derived from its transient-equilibrium parent molybdenum-99 ( 99 Mo) (t½ = 66 h), the short-lived isotope technetium-99m ( 99m Tc) (t ½ = 6.0 h) is the most widely used medical isotope in diagnostic imaging today. The National Nuclear Security Administration’s (NNSA’s) Material Management and Minimization (M3) program—established under the auspices of the American Medical Isotope Production Act—has been tasked with facilitating the work of domestic 99 Mo suppliers that do not use highly enriched uranium (HEU). Superconducting electron linear accelerators that employ high-Z converter targets can generate bremsstrahlung photons and neutron fluxes that can induce photonuclear reactions and uranium (U) fission. Argonne, in collaboration with industrial partners, has been developing the process chemistry for superconducting linear accelerator (LINAC)-irradiated triuranium octoxide (U 3 O 8 ) targets to produce 99 Mo. In short, the process involves dissolving the irradiated targets in nitric acid, which simultaneously releases valuable fission products (xenon and iodine). After adjusting the acid concentration, producers process the feed using tri-n-butyl phosphate (TBP), then mix the raffinate derived from this extraction cycle with a phosphinic acid such as Cyanex 272 or di(2-ethylhexyl) phosphoric acid (HDEHP), which is selective for Mo (molybdenum liquid liquid extraction [MoLLE]). In the next step, a concentration column comprising an anion exchange platform is used to decontaminate the remaining fission products, generating a pure 99 Mo stream in sodium hydroxide (NaOH)/sodium chloride (NaCl).

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Experimental validation of a high fidelity Monte Carlo neutron transport model of the MIT graphite exponential pile

High-fidelity modeling and simulation were performed for the MIT graphite exponential pile (MGEP) using Monte Carlo neutron transport codes OpenMC and MCNP, and the results were validated by experimental data. The MGEP is being used as the test bed for the design of an autonomous control system for the pile's neutron flux distribution. The main contribution of this work is to generate the training data sets of neutron flux distributions with different locations of control rods that perturb the neutron flux profiles. First, code -to-code cross verification between OpenMC and MCNP was performed to ensure consistency of the numerical modeling within statistical uncertainties. To validate the accuracy of this high-fidelity model, a series of neutron flux measurements were conducted using a Helium-3 (He-3) neutron detector on a mobile platform that is placed inside the pile. Second, the neutron flux profiles were measured in four vertical layers of interest, and compared to the corresponding simulation results. The comparison results shows that the root mean square error is less than 2.5% in the two upper layers, and less than 4.5% in all four measured layers. Here the results validated the accuracy of the modeling and simulation. Finally, the relative change of the neutron flux profiles from moving control rods was analyzed, which identified the layer that has the best sensitivity regarding the control rods movements. Thus, this work identified and provided training data sets of both simulated and experimental neutron flux profiles in the most sensitive layer, paving the path forward to the real-time experimental demonstration of the autonomous control system.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Photoneutron Production Using an Electron Linear Accelerator for Applications in Neutron Imaging

Photoneutron production is possible using an electron linear accelerator and a target capable of generating photonuclear reactions. A short pulse neutron source can be useful for neutron imaging dynamic experiments. This study is aimed at the feasibility of photoneutron production using a 20 MeV electron linear accelerator and tungsten and depleted uranium targets of various thicknesses. MCNP6 (Monte Carlo N-Particle) code will be used to develop a computational model to estimate total neutron yield, and this will be verified at the Idaho State University’s Accelerator Center. After verification of the neutron yield and energy spectra, an additional MCNP6 model will be developed to analyze the neutron imaging processes. This study will potentially prove it is possible to conduct multi-mode imaging experiments on the anticipated Scorpius electron linear accelerator at the Nevada National Security Sites.

43 PARTICLE ACCELERATORS↗

AT2017gfo: Bayesian inference and model selection of multicomponent kilonovae and constraints on the neutron star equation of state

The joint detection of the gravitational wave GW170817, of the short γ-ray burst GRB170817A and of the kilonova AT2017gfo, generated by the the binary neutron star (NS) merger observed on 2017 August 17, is a milestone in multimessenger astronomy and provides new constraints on the NS equation of state. We perform Bayesian inference and model selection on AT2017gfo using semi-analytical, multicomponents models that also account for non-spherical ejecta. Observational data favour anisotropic geometries to spherically symmetric profiles, with a log-Bayes’ factor of ~104, and favour multicomponent models against single-component ones. The best-fitting model is an anisotropic three-component composed of dynamical ejecta plus neutrino and viscous winds. Using the dynamical ejecta parameters inferred from the best-fitting model and numerical–relativity relations connecting the ejecta properties to the binary properties, we constrain the binary mass ratio to q < 1.54 and the reduced tidal parameter to $120\lt \tilde{\Lambda }\lt 1110$. Finally, we combine the predictions from AT2017gfo with those from GW170817, constraining the radius of a NS of 1.4 M ⊙ to 12.2 ± 0.5 km (1σ level). This prediction could be further strengthened by improving kilonova models with numerical-relativity information.

79 ASTRONOMY AND ASTROPHYSICS↗

Inferring Nuclear Hamiltonians from 3 rd -Generation Gravitational-Wave Detectors [Slides]

Neutrons stars are some of the densest objects in our universe. Thus, they provide us with a laboratory to study dense nuclear matter and its properties. The neutron star’s stability at such high densities, with an immense gravitational force, is due in large part to nuclear interactions. Specifically, nuclear 3-body forces provide significant repulsive contributions and at high densities. Current and future gravitational-wave detectors can help us to better understand these nuclear interactions. My project aims to predict the constraining power of the next generation of Gravitational Wave detectors: Cosmic Explorer. I have constructed a solver for the Tolman-Oppenheimer-Volkoff equations and studied equations of state with varying strengths of the 3-body forces. Using simulated data from Phillipe Landry at CITA, I can predict to which level of uncertainty Cosmic Explorer will enable us to measure nuclear 3-body force.

47 OTHER INSTRUMENTATION↗

Improvements to the Modeling of the TREAT Reactor and Experiments

This paper summarizes the latest improvements and lessons learned from the modeling and simulation of the transient test reactor at Idaho National Laboratory using the MAM- MOTH reactor physics application. MAMMOTH is a MOOSE-based, Finite Element Method application that specializes in the analysis of the spatial dynamics behavior of nuclear reactors. Since early 2018 several transient tests have been conducted at TREAT, thus providing the opportunity to apply and benchmark modern modeling and simulation tools. MAMMOTH was used to provide predictions of the power coupling factor between the core and the experiment for various experiments. Even though the power coupling factor predictions agree very well with the experimental data, within the bounds of the experimental uncertainty, one shortcoming was the underprediction of the total energy deposited in the core and in the experiment. Determination of the sources for this discrepancy are ongoing, but several key problems have been identified and resolved, thus providing valuable insights for future research. This paper discusses several of these lessons learned. First, the heat capacity data for the TREAT fuel has some significant problems due to limitations of the measurement techniques used circa 1960s. The sensitivity of the peak power and the total energy deposition to various representations of the heat capacity is approximately 5%. Second, the effects of the biological shield and thermal column on the modeling of the core are non-negligible, since they affect the mean generation time and the effective reflection of neutrons back into the core, which is suspected to be important during the core heat up. Matching the reactor period resolves the fact that the reduced spatial domain used in the MAMMOTH model underpredicts the mean generation time. The neutron reflection from these regions is marginally improved with the use of an albedo boundary condition. Third, modeling of the control rod movement with a multi-scheme method is introduced and its current limitations are exposed. Fourth, we explore the effects of using a homogenized model with Superhomogenization equivalence and how that differs from fully heterogeneous simulations. Finally, the energy condensation effects for this graphite core are significant. Solutions with 10 and 26 energy groups show the benefits of using a finer coarse group structure.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Advanced Cross Section Library Generation using Reduced Order Models

Deterministic neutronics calculations rely on multigroup neutron cross section libraries, which consist of databases of tabulated values, used to calculate the neutron cross sections through multivariate linear interpolation. However, interpolation of the multidimensional cross section data becomes memory inefficient and time consuming as the number of tabulations increases, significantly slowing down the neutronics calculation, especially in the case of microscopic cross section libraries where every isotope (on the order of hundreds) has its own set of specific reactions and cross sections. In order to address this challenge, this work constructs efficient and robust reduced-order models (ROMs) of the multi-group cross sections to support the Griffin simulation of high-temperature gas-cooled reactors (HTGRs). The first part of the study investigates the linearity of the multigroup cross section data across isotopes, reaction types, and energy groups on pre-generated datasets for the purpose of dimensionality reduction. Secondly, a down-selection of ROM techniques is presented on representative classical machine learning (ML) techniques, including variants of linear regression, kernel-based methods, tree-based algorithms, and artificial neural networks. The selection criteria jointly consider the memory efficiency, predictive accuracy, prediction speed, and scalability in comparison to the multidimensional interpolation. Among all the ML techniques, deep neural networks (DNNs) have proven to be the best selection with sufficient accuracy, high robustness, good memory efficiency, great scalability, and superior flexibility. DNNs have been trained for all isotopes in this work and systematic Griffin testing is ongoing to ensure the feasibility of this ROM technique for predicting cross section and reducing memory requirements without a significant sacrifice in computational performance.

42 - ENGINEERING↗

Radiation effects on stress evolution and dimensional stability of large fusion energy structures

Here we assess the effects of neutron irradiation on the deformation and stress evolution of large-scale fusion energy structures. This is accomplished through non-linear finite element structural analysis of the coupled thermal and mechanical fields at the Beginning-Of-Life (BOL), at 45 dpa, and at 90 dpa. Radiation effects include volumetric swelling and the influence of radiation on the mechanical properties. The system studied here is a large section of a full inboard module of an integrated structure comprising the First Wall and Blanket (FW/B) of a Dual Cooled Lithium-Lead (DCLL) energy conversion unit in the Fusion Nuclear Science Facility (FNSF). The structural material is the ferritic/martensitic steel F82H. We analyze several radiation effects phenomena that can lead to significant impact on the mechanical design and lifetime of the structure. These include volumetric swelling effects on the assembly and disassembly of modules and the impact of spatial gradients in neutron damage and helium generation on deformation, stress evolution, and plastic strain accumulation. We show that the stress state is strongly influenced by the accumulation of swelling strains, radiation hardening and softening, and by spatial gradients in displacement damage and helium gas generation. Several key regions in the FW/B structure are identified where severe plastic strains accumulate and may be potential sites for failure, especially during cyclic reactor operations. We finally show that minimum gaps between Inboard (IB) blanket sectors must be present during assembly so as to avoid excessive stress and strain in FW/B modules critical regions and raising chances of failure. A range of structural displacements are shown, from as low as 4 mm at BOL for free side walls to as high as 46 mm at 90 dpa for fully constrained side walls.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Robust direct laser acceleration of electrons with flying-focus laser pulses

Direct laser acceleration (DLA) offers a compact source of high-charge, energetic electrons for generating secondary radiation or neutrons. While DLA in high-density plasma optimizes the energy transfer from a laser pulse to electrons, it exacerbates nonlinear propagation effects, such as filamentation, that can disrupt the acceleration process. Here, we show that superluminal flying-focus pulses (FFPs) mitigate nonlinear propagation, thereby enhancing the number of high-energy electrons and resulting x-ray yield. Three-dimensional particle-in-cell simulations show that, compared to a Gaussian pulse of equal energy (1 J) and intensity (2 × 10 20 W/cm 2 ), an FFP produces 80 × more electrons above 100 MeV, increases the electron cutoff energy by 20%, triples the high-energy x-ray yield, and improves x-ray collimation. These results illustrate the ability of spatiotemporally structured laser pulses to provide additional control in the highly nonlinear, relativistic regime of laser-plasma interactions.

Laser-produced plasmas↗

Equilibrium Core Model for Micro Pebble Bed Reactors Using OpenMC

Estimating the equilibrium state for pebble bed reactors (PBRs) presents complex challenges as it requires simultaneous consideration of changes in the pebbles’ movement as well as their fuel compositions. Whereas traditional approaches use multigroup diffusion codes for neutronics calculations of PBRs’ equilibrium state, the double-heterogeneity of PBRs complicates neutron cross-section generation. Continuous-energy Monte Carlo (MC) methods are better suited for detailed PBR analysis because of their natural handling of double-heterogeneity, but they demand substantially more computational resources. Here, this study introduces a novel method for efficiently estimating the equilibrium state in small and micro PBRs with reduced computational cost. The method is anticipated to accelerate the processes of core design and performing parametric studies for utilizing advanced fuel and structural materials. The HTR-10 reactor design was used for validating the method’s predictions and evaluating its computational efficiency. When compared to reference calculation values from the literature, criticality (k-effective) was predicted to be approximately within the margin of error of the MC transport calculation, average core power density (in megawatts per cubic meter) was predicted within 2.5% relative error, and maximum thermal flux (10 13 n/cm 2 .s −1 ) was predicted within 1.8% relative error. The calculated inventory of fission products and fuel composition in the equilibrium core were within 15% and 16.6%, respectively, when compared to reported values from the literature. The difference is attributed to variance in the considered values of the core temperature, which was found to significantly affect the depletion analyses.

Equilibrium core↗