Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Neutron physics”

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 289 records · Page 16

Level density within a micro-macroscopic approach

Statistical level density $ρ$($E, A$) is derived for nucleonic system with a given energy $E$, particle number $A$ and other integrals of motion in the micro-macroscopic approximation beyond the standard saddle-point method of the Fermi gas model. This level density reaches the two limits; the well-known Fermi gas grand-canonical ensemble limit for a large entropy $S$ related to large excitation energies, and the finite micro-canonical limit for a small combinatorical entropy $S$ at low excitation energies. In conclusion, the inverse level density parameter $K$ as function of the particle number $A$ in the semiclassical periodic orbit theory, taking into account the extended Thomas-Fermi and Strutinsky shell corrections, is calculated and compared with experimental data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Simultaneous extraction of the weak radius and the weak mixing angle from parity-violating electron scattering on 12 C

We study the impact of nuclear structure uncertainties on a measurement of the weak charge of 12 C at the future Mainz Energy Recovering Superconducting Accelerator (MESA) facility in Mainz. Information from a large variety of nuclear models, accurately calibrated to the ground-state properties of selected nuclei, suggest that a 0.3% precision measurement of the parity-violating asymmetry at forward angles will not be compromised by nuclear structure effects, thereby allowing a world-leading determination of the weak charge of 12 C . Furthermore, we show that a combination of measurements of the parity-violating asymmetry at forward and backward angles for the same electron-beam energy can be used to extract information on the nuclear weak charge distribution. We conclude that a 0.34% precision on the weak radius of 12 C may be achieved by performing a 3% precision measurement of the parity-violating asymmetry at backward angles.

Electroweak interaction↗

Improvements to Contributions from Neutron Inelastic Scattering for Next-Event Estimators in MCNP ® Software

An improvement to the handling of contributions from neutron inelastic scattering to next-event estimators has been implemented in the MCNP6 ® software for release with version 6.3.1. The kinematic equations that govern the outgoing energy of inelastic neutron scattering contributions to next-event estimators have two roots. Historically, the implementation in the MCNP coding has only used the upper root to the quadratic equation. This includes all versions predating version 6.3.1 including all versions of MCNP5 and MCNPX software. However, a review of the neutron next-event estimator physics has shown that this does not reproduce the track-length estimator results at low energies. Several examples are presented that test single neutron inelastic scattering reaction types, Level Scattering (Law 3), Tabulated Energy Angle (Law 61), and Kalbach-Mann distribution (Law 44). The test problems compare the track-length estimator (f4 tally) in the MCNP software, with the existing implementation of the neutron next-event estimator (f5 tally), and the modified changes to the neutron next-event estimator implementation. As the MCNP implementation ignores the lower root, the unmodified neutron next-event estimators will generally underestimate the lower energy contribution. However, a second issue with the Kalbach-Mann distribution (Law 44) implementation allows contributions to backward scattering in the center-of-mass frame that is not kinematically possible, thus overestimating backward scattering contributions. A third issue with the way the MCNP implementation handles floating point comparison for scattering directly ahead or directly backward in the center-of-mass frame generally leads to underestimation (except for backwards scattering for Law 44).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DRiFT - Release 2.1.0: Organic Scintillators and Gas Detectors

DRiFT (a Detector Response Function Toolkit) is LANL-developed software that post-processes output from the extensively validated radiation transport code, MCNP, and generates realistic nuclear instrumentation response. DRiFT is designed to be flexible, enabling users to specify detector type and many experimental settings, as well as accommodating the addition of their own desired features. The focus of this release is on organic scintillator, gas detector, and associated capabilities, while semiconductor features are still under development. Organic scintillators are widely used in the areas of nuclear safeguards and nuclear non-proliferation efforts. DRiFT has several diagnostic and detector physics features relevant to detailed scintillator simulations including: tracking source particle information, scintillation light production, the effects of PMT quantum efficiency and gain, and digitizer settings. Users can select responses from many scintillator and PMT types supported natively by DRiFT, or add their own by following the instructions in this document. DRiFT also has several diagnostic and detector physics features relevant to neutron gas detector simulations including: gas detector wall effects, inactive areas at the end of detector tubes, and effects of a preamplifier. We acknowledge that DRiFT is under active development, bug reports and general questions and comments should be directed to Madison Andrews, madison@lanl.gov. This manual is divided into four parts: I) An overview of DRiFT, including how to obtain and install the executable, II) A description of the detector physics related to scintillators available, III) a description of more general DRiFT features the user may find useful, and IV) a description of the test suite and examples made available with the code release.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Informing nuclear physics via machine learning methods with differential and integral experiments

Information from differential nuclear-physics experiments and theory is often too uncertain to accurately define nuclear-physics observables such as cross sections or energy spectra. Integral experimental data, representing the applications of these observables, are often more precise but depend simultaneously on too many of them to unambiguously identify issues in the observable with human expert analysis alone. Here, we explore how we can leverage physics knowledge gained from differential experimental data, nuclear theory, integral experiments, and neutron-transport calculations to better understand nuclear-physics observables in the context of the application area represented by integral experiments. We support this task with machine-learning methods to discern trends in a large amount of convoluted data. Differential and integral information was used in an analysis augmented by the random forest and the Shapley additive explanations metric. We chose as an application area one that is represented by criticality measurements and pulsed-sphere neutron-leakage spectra. We show one representative example ( 241 Pu fission observables) where the combination of differential and integral information allowed to resolve issues in data representing these observables. As a starting point, the machine learning (ML) algorithms highlighted several observables as leading potentially to bias in simulating integral experiments. Differential information, paired with sensitivity to integral quantities, allowed us then to pinpoint one specific observable ( 241 Pu fission cross section) as the main driver of bias. The comparison to integral experiments, on the other hand, allowed us to indicate a likely reliable experiment among several discrepant ones for this observables. In other cases (e.g., 239 Pu observables), we were not able to resolve the confounding introduced by integral experiments but instead highlighted the need for targeted new experiments and theory developments to better constrain the nuclear-physics space for the application area represented by integral experiments. We were able to combine information from differential experimental data, nuclear-physics theory, integral experiments, and neutron-transport simulations of the latter experiments with the help of the random forest algorithm and expert judgment. This combination of knowledge allows to improve our description of nuclear-physics observables as applied to a particular application area.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Sensitivity of Neutron-Rich Nuclear Isomer Behavior to Uncertainties in Direct Transitions

Nuclear isomers are populated in the rapid neutron capture process (r process) of nucleosynthesis. The r process may cover a wide range of temperatures, potentially starting from several tens of GK (several MeV) and then cooling as material is ejected from the event. As the r-process environment cools, isomers can freeze out of thermal equilibrium or be directly populated as astrophysically metastable isomers (astromers). Astromers can undergo reactions and decays at rates very different from the ground state, so they may need to be treated independently in nucleosythesis simulations. Two key behaviors of astromers—ground state ↔ isomer transition rates and thermalization temperatures—are determined by direct transition rates between pairs of nuclear states. We perform a sensitivity study to constrain the effects of unknown transitions on astromer behavior. Detailed balance ensures that ground → isomer and isomer → ground transitions are symmetric, so unknown transitions are equally impactful in both directions. We also introduce a categorization of astromers that describes their potential effects in hot environments. We provide a table of neutron-rich isomers that includes the astromer type, thermalization temperature, and key unmeasured transition rates.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Scattering Calorimeter FY24 Deliverable Report

A simulation-based method has been developed to prototype new detector designs for nuclear data measurements utilizing neutron scattering. This method uses representative physics inputs for signal and background generation, full detector resolution smearing benchmarked by experimental data, and a neutron beam timing simulation to produce analyzable output like a physical measurement. A test case has been studied using a hybrid time-of-flight calorimeter detector for scattering cross-section measurements on 239 Pu with 1-5 MeV incident monoenergetic neutrons. Data analysis methods have been developed to perform event-level particle reconstruction and reaction channel discrimination. This analysis has been used to estimate the capability of the test detector to perform simultaneous scattering and fission cross section measurements, as well as its ability to provide neutron spectra and particle angular information.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The calculation of light element impurity (α, n) yield curves in a PuO 2 matrix and associated specific yield coefficients: Influence of the reaction cross sections

Most of the Pu separated from irradiated commercial nuclear fuel is stored as PuO 2 . The primary quantitative nondestructive measurement technique used to verify the amount of Pu in storage containers is passive neutron correlation counting. An important physical property of the oxide material is the ratio, α, of the rate of (α, n) neutrons produced inside the item to the rate of neutrons produced by spontaneous fission. This ratio influences the precision of the correlated counting method and affects the interpretation of the data because of how it changes both the primary total neutron production rate and the rate of induced fission events taking place inside the item. In addition to the main O(α, n) contribution, additional contributions come from α-particle interactions with light element impurities that are inevitably present. In this work, we calculate specific (α, n) yield coefficients, expressed in units of neutrons per second per gram of α-emitting nuclide per part per million by mass of the specified impurity element distributed in a pure PuO 2 matrix, for some key α-emitting actinides commonly present in reprocessed Pu ( 238–242 Pu+ 241 Am). These coefficients are directly applicable to nuclear safeguards verification work in which the α ratio is often calculated from the Pu-isotopic composition and chemical information obtained by other means. Further, they also provide a convenient up-to-date reference set against which values generated by other methods can be compared. Results are presented for impurities with atomic number from 3 to 17 inclusive, plus K and Fe. In most cases, these coefficients are not expected to change by more than 5%–10% at any time in the future. However, as new data become available, changes as large as 20% may be needed for some targets (e.g., F). The present yield calculations are limited by the general shortage of quality experimental total (α, n) reaction cross section data, which, together with unexplained variation between determinations, means that an objective and coherent evaluation is not possible. The situation is even less satisfactory for the partial differential cross section needed to calculate neutron spectra.

(𝛼, n) reactions↗

Multiplicity counting using organic scintillators to distinguish neutron sources: An advanced teaching laboratory

In this advanced instructional laboratory, students explore complex detection systems and nondestructive assay techniques used in the field of nuclear physics. After setting up and calibrating a neutron detection system, students carry out timing and energy deposition analyses of radiation signals. Through the timing of prompt fission neutron signals, multiplicity counting is used to carry out a special nuclear material (SNM) nondestructive assay. Our experimental setup is comprised of eight trans-stilbene organic scintillation detectors in a well-counter configuration, and measurements are taken on a spontaneous fission source as well as two (α,n) sources. By comparing each source's measured multiplicity distribution, the resulting measurements of the (α,n) sources can be distinguished from that of the spontaneous fission source. Such comparisons prevent the spoofing, i.e., intentional imitation, of a fission source by an (α,n) neutron source. This instructional laboratory is designed for nuclear engineering and physics students interested in organic scintillators, neutron sources, and nonproliferation radiation measurement techniques.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Measurements of effective actinides neutron capture cross sections in a cadmium-filtered epithermal neutron spectrum

In this work, a total of 20 actinide samples were irradiated in the Advanced Test Reactor (ATR) at Idaho National Laboratory (INL): 233 U (1), 235 U (1), 236 U (1), 238 U (1), 237 Np (2), 239 Pu (1), 240 Pu (1), 242 Pu (2), 244 Pu (2), 241 Am (2), 243 Am (2), 244 Cm (2), and 248 Cm (2). The samples were loaded in a basket containing a 1.1-mm thick cadmium filter. The presence of the filter drastically changes the neutron spectrum in the samples and accentuates the reaction rates in the resonance region. Isotopic ratios were determined using Multi-Collector Inductively Coupled Plasma Mass Spectrometry (MC-ICP-MS) that allowed inferring effective capture cross sections. Detailed best-estimate MCNP models of the experiment using ENDF/B-VII.0 nuclear data were developed to calculate the effective neutron cross sections and, subsequently, compare them with the cross sections inferred from experimental results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Development of HELIOS/BIPR/PARCS/MCNP6 Computation Route for WWER RPV Neutron Fluence Analysis and Validation Against Ex-Vessel Detector Measurement Data

Here, a computational route was developed for precise calculation of fast neutron fluence on a WWER-type reactor pressure vessel (RPV). The method is based on the transfer of neutronics data from HELIOS-2 lattice calculations and nodal diffusion neutronics data (power, density, and temperature) from BIPR7.1 and PARCS 3.36/PATHS core calculations into a three-dimensional (pinwise axially distributed) fixed neutron source for modeling of transport of fast neutrons from the reactor core to the outer surface of the RPV using MCNP6.2. Validation of the proposed computational method was carried out based on comparative analysis of MCNP6.2-predicted and neutron dosimetry–measured reaction rates [ 54 Fe(n,p) 54 Мn, 93 Nb(n,n') 93 mNb, and 58 Ni(n,р) 58 Со] on the outer surface of the Armenian Nuclear Power Plant (ANPP) Unit 2 RPV. Validation revealed that the MCNP6.2-predicted fast neutron fluence results are very sensitive to the ENDF-B neutron data. Particularly, MCNP6.2 with ENDF/B-VIII.0 significantly underpredicts (20% to 30%) fast neutron fluence while using ENDF/B-VII.1 data overpredicts it. Adding revised beta-released evaluations of 54 Fe, 56 Fe, 57 Fe, and 16 O from the International Nuclear Data Evaluation Network (INDEN) to ENDF/B-VIII.0 allows one to obtain reasonable agreement with measurement results for all types of measured reaction rates.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Fission and the r -process nucleosynthesis of translead nuclei in neutron star mergers

In this report we study the impact of fission on the production and destruction of translead nuclei during the r-process nucleosynthesis occurring in neutron-star mergers. Abundance patterns and rates of nuclear energy production are obtained for different ejecta conditions using three sets of stellar reaction rates, one of which is based on microscopic and consistent calculations of nuclear masses, fission barriers, and collective inertias. We show that the accumulation of fissioning material during the r process can strongly affect the free neutron abundance after the r-process freeze-out. This leads to a significant impact on the abundances of heavy nuclei that undergo α decay or spontaneous fission, affecting the radioactive energy production by the ejecta at timescales relevant for kilonova emission.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A multiphysics model of the versatile test reactor based on the MOOSE framework

The traditional modeling approach for sodium fast reactor cores relies on separate physics models, where the fuel performance, thermal–hydraulics, and neutronics calculations required to predict the core physics characteristics for nominal conditions are decoupled by relying on user-imposed boundary conditions. Here, this paper aims at evaluating the impact of multiphysics simulations for predicting the core characteristics of the Versatile Test Reactor, which is being designed as a 300-MWt sodium-cooled fast reactor. The purpose of the Versatile Test Reactor is to accelerate the testing of advanced nuclear materials in the United States. The proposed multiphysics model relies on the Griffin reactor physics code, the SAM thermal–hydraulic system code, the BISON fuel performance code, as well as generic Multiphysics Object-Oriented Simulation Environment capabilities implemented in the open-source tensor mechanics module. For k eff calculations, the introduction of a tight coupling between the neutronics, thermo-mechanical and thermal–hydraulics models induces a change of around 543 pcm in the eigenvalue, compared to the traditional standalone neutronics calculation where approximate temperature profiles are used. The multiphysics model is then employed for quantifying the impact of the thermal conductivity uncertainties on some of the key figures of merit, such as the fuel centerline temperature, assembly powers, and keff for nominal core conditions. As anticipated, uncertainties on fuel thermal conductivity mostly impact the fuel centerline temperature, and to a lesser extend the k eff .

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Future Physics Perspectives on the Equation of State from Heavy Ion Collisions to Neutron Stars

With the computational power and algorithmic improvements available today, the ongoing STAR/RHIC and HADES/GSI experiments, the future FAIR and NICA facilities becoming operational, and the new precise measurements from NICER and LIGO/VIRGO, the high-energy nuclear physics and astrophysics communities are in the unique position to set very stringent constraints on the equation of state of strongly interacting matter. Here, we review the state-of-the-art of different approaches used in the description of hot and ultradense baryonic matter in and out of equilibrium, and discuss the regions in the phase diagram where heavy-ion collisions and neutron star mergers can overlap. Future perspectives are discussed to help define a comprehensive, multi-disciplinary strategy to map out the phase diagram of strongly interacting matter from heavy ion collisions to neutron stars.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Estimating β eff from static, late-time Feynman neutron noise measurements

The effective delayed neutron fraction, β eff , is a kinematic parameter describing the contribution of delayed neutrons to the effective multiplication factor, $\kappa$ eff , and is of great interest to the nuclear criticality safety community and beyond. In fact, $\kappa$ eff presently must be inferred (as opposed to directly measured) since variables like β eff must be obtained from simulation, reference tables, or a dynamic or pulsed measurement where the assembly of fissile material is perturbed and not in steady state. Several recent works have proposed neutron noise techniques to estimate β eff from static measurements (such as benchmark experiments). This work experimentally investigates an approach based on looking at late-time (up to thousands of seconds) correlated neutrons in Feynman histograms. This work utilizes measurements of weapons grade plutonium and highly enriched uranium to construct late-time Feynman histograms and investigate a proposed methodology of estimating β eff from a fit of the distribution. Data is shown from both He-3 detectors and an organic scintillator array. This work is preliminary in nature, demonstrating some experimental feasibility and proposing improvements for future investigations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

SERPENT 2 and TRIPOLI-4 neutronics benchmarking of the Jules Horowitz reactor

Jules Horowitz Reactor (JHR) is a future center for nuclear material testing in Europe, currently under construction at CEA Cadarache in southern France. It has been modelled extensively with CEA's HORUS3D/N neutron calculation package, that includes a Monte-Carlo neutron transport code TRIPOLI-4. In order to support the future computational analyses on the reactor core physics, VTT's Monte Carlo neutronics code SERPENT 2 has been used to create an alternative model of the JHR, taking advantage of the pre-existing Constructive Solid Geometry (CSG) structures on which the TRIPOLI-4 model was built. This work describes a benchmark, in which two JHR core configurations were analyzed using both of these models. The first configuration is a fresh fuel core at the beginning of life, in which aluminum based dummy-devices are installed in the reactor test positions. The second configuration represents the end-of-cycle equilibrium core with a heterogeneous burnup distribution, in which CALIPSO devices containing cladding material samples are installed at the in-core experimental positions, while the reflector positions are occupied by DEN-REP devices containing UO{sub 2} fuel. Neutron flux spectra and heat deposition are the parameters of interest in this study, due to their high importance in creating the desired experimental irradiation environment. The impact of delayed gamma energy deposition will be studied in detail, since the two models account for it with different methods. This work aims to establish a SERPENT 2 model of JHR, understanding its differences and agreements with TRIPOLI-4, and discussing its potential future applications. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Correlating isothermal compressibility to nucleon fluctuations in the inner crust of neutron stars

The question of how and which physical observables or thermodynamic parameters can best predict the onset of a possible phase transition in the inner crust of neutron stars remains largely unresolved. Here, using semiclassical Monte Carlo simulations, we investigate the isothermal compressibility and density fluctuations in a region of relevance to the dynamics of the inner crust. We show that the isothermal compressibility serves as a robust observable to characterize the transition from the nonuniform crust to the uniform core for proton fractions over 0.2. Moreover, we show explicitly how the two-component isothermal compressibility, computed using the Kirkwood-Buff theory, is directly connected to the fluctuations in the number density, recorded in the grand canonical ensemble by monitoring the number of particles in a small volume located at the center of the simulation box. That is, we compute mean-square particle fluctuations and compare them against the isothermal compressibility for different proton fractions. Although our results show that the mean-square particle fluctuations are proportional to the isothermal compressibility, the lack of a perfect correlation is attributed to the relatively small number of particles included in the simulations. The nonunity slope observed in the dimensionless isothermal compressibility—total nucleon fluctuation variance relationship suggests that the inner crust of neutron stars is composed of anisotropic and inhomogeneous matter.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗