Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Criticality Safety”

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 793 records · Page 44

ENDFtk: A robust tool for reading and writing ENDF-formatted nuclear data

ENDFtk is a recently developed C++ and Python interface to interact with ENDF-6 formatted nuclear data files. It provides a robust and complete interface, allowing the reading and writing of all formats currently part of the ENDF-6 formats manual, as well as some non-ENDF formats used by the NJOY processing code. It provides an interface that mimics the names in the ENDF-6 formats manual as well as an equivalent interface using human-readable attribute names. It is robust and powerful enogh for nuclear data experts to develop complex applications, while also simple enough to be used non-experts to retrieve and manipulate evaluated nuclear data. ENDFtk offers the ability to easily interrogate and manipulate data either in large-scale code projects or in simple Python scripts. Here, in this paper, a brief overview of the interface is given, as well as more substantial examples demonstrating plotting simple data, interacting with more complex data, and writing new data to files. ENDFtk is open source and available for download via GitHub (https://github.com/njoy/ENDFtk).

97 MATHEMATICS AND COMPUTING↗

Dr. Cecil Lubitz (1925–2021)

Dr. Cecil Lubitz passed away on March 5, 2021 at his home in upstate New York. A graduate of the Class of 1945 from the Naval Academy, Cecil served in the US Navy until 1947 after which he earned a MS in electrical engineering (1949) and a PhD in Physics (1960) from the University of Michigan. Following graduate school, Cecil came to work at Knolls Atomic Power Laboratory (KAPL) on applied nuclear physics in the Naval Nuclear Propulsion Program (NNPP). Cecil worked at KAPL for 54 years and had an impactful career with many substantial contributions to the NNPP, national and global nuclear data communities.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗

Measurement of the 239 Pu(n,f)/ 235 U(n,f) Cross-Section Ratio with the NIFFTE fission Time Projection Chamber

In this work, the 239 Pu(n,f)/ 235 U(n,f) cross-section ratio has been measured with the fission Time Projection Chamber (fissionTPC) from 100 keV to 100 MeV. The fissionTPC provides three-dimensional reconstruction of fission-fragment ionization profiles, allowing for a precise quantification of measurement uncertainties. The measurement was performed at the Los Alamos Neutron Science Center which provides a pulsed white source of neutrons. The data are recommended to be used as a cross-section ratio shape. A discussion of the status of the absolute normalization and comparisons to ENDF evaluations and previous measurements is included.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Rossi-alpha measurements of fast plutonium metal assemblies using organic scintillators

The reactivity and the k-effective multiplication factor (k eff ) of a fissionable assembly are quantities of widespread interest. These values can be inferred from Rossi-alpha measurements of the prompt neutron decay constant (or the inverse:α -1 ). It has been shown that 3 He-gas proportional counter-based detection systems are insensitive to α -1 of fast assemblies (much faster than tens of microseconds). Therefore, it is of interest to investigate fast detection systems such as those based on organic scintillation detectors. In this work, an array of 12 cylindrical, 5.08 cm x 5.08 cm diameter trans-stilbene organic scintillators was used to measure five subcritical assemblies. One assembly was a sphere of approximately 4.5 kg of alpha-phase, weapons-grade plutonium (k eff = 0.773, α -1 = 11.84 ns) known as the BeRP ball encased in a thin stainless-steel clad. The other assemblies used the same encased sphere and 7.62 cm of iron, nickel, copper, or tungsten reflectors (k eff = 0.884, 0.916, 0.924, 0.939, respectively and α -1 = 36.60, 41.56, 49.60, 70.32 ns, respectively). Finally, this work (1) validates Rossi-alpha measurements with organic scintillators by demonstrating good agreement between measurements and simulations; and (2), demonstrates that organic scintillator-based systems are sensitive to α -1 on the order of 10–100 ns.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Measurement of material isotopics and atom number ratio with α-particle spectroscopy for a NIFFTE fission Time Projection Chamber actinide target

In this work, we present the results of a measurement of isotopic concentrations and atomic number ratio of a double-sided actinide target using α-spectroscopy and mass spectrometry. The double-sided actinide target, with predominantly 239 Pu on one side and 235 U on the other, was used in the fission Time Projection Chamber (fissionTPC) for a measurement of the neutron-induced fission cross-section ratio between the two isotopes. The measured atomic number ratio is needed to extract an absolute measurement fission cross-section ratio. The 239 Pu/ 235 U atom number ratio was measured with a combination of mass spectrometry and α-spectroscopy with a planar silicon detector achieving uncertainties of less than 1%. Different strategies for estimating isotopic concentration from the α-spectrum are presented to demonstrate the potential of these methods for non-destructive target assay. We found that a combination of fitting spectra with constraints from mass spectrometry, and summing counts in a region of the spectrum provided the most consistent results with the lowest uncertainty.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Quantification of the 35 Cl($\textit{n, p}$) reaction channel

The evaluated 35 Cl($\textit{n, p}$) cross sections in the ENDF/B-VIII.0 nuclear data library remain uncertain because of the lack of measured data in the neutron energy range between 100-600 keV. In this energy region, the R-matrix analysis relies heavily on estimated ($\textit{n, p}$) partial widths and level density information. This resulted in significant differences of the evaluated cross sections and predicted reactivity in nuclear reactor applications with respect to previous nuclear data libraries. Here, we present a methodology to quantify the cross sections for the 35 Cl($\textit{n, p}$) reaction channel from the combination of measured total cross sections and model calculations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Trace Actinide Signatures of a Bulk Neptunium Sample

In nuclear forensic analyses, measurements of actinide elements in a sample can assist with identifying interdicted or unknown materials. While these radiochemical signatures have been extensively investigated in uranium materials, less is known about bulk neptunium samples. This paper describes the measurement of trace actinide concentrations and isotopic profiles in a 237 Np oxide sample. Uranium, plutonium, americium, and curium concentrations and isotopic profiles in the sample were determined and deemed potentially useful for distinguishing different sources of 237 Np. Furthermore, several different potential radiochronometry systems were also investigated; discordant results indicate that the Np sample was never completely purified of other actinide elements, or that subsequent contamination of the sample occurred. Few prior studies of neptunium materials have been reported, and these data suggest that trace actinide constituents could provide unique signatures to identify material out of regulatory control.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Re-evaluation of the TSL for Yttrium Hydride

Yttrium hydride (YH x ) is of interest as a high-temperature moderator material in advanced nuclear reactor systems because of its superior ability to retain hydrogen at elevated temperatures. Thermal neutron scattering laws (TSL) for hydrogen bound in yttrium hydride (H-YH 2 ) and yttrium bound in yttrium hydride (Y-YH 2 ) were previously evaluated by Naval Nuclear Laboratory using the ab initio approach and released in ENDF/B-VIII.0. In that work, density functional theory, incorporating the generalized gradient approximation (GGA) for the exchange-correlation energy, was used to simulate the face-centered cubic structure of YH 2 and calculate the interatomic Hellmann-Feynman forces for a 2×2×2 supercell containing 96 atoms. Lattice dynamics calculations using PHONON were used to determine the phonon density of states. The calculated phonon density of states for H and Y in YH 2 were then used to prepare H-YH 2 and Y-YH 2 TSL evaluations, in the incoherent approximation, using the LEAPR module of NJOY. In addition, elastic scattering was assumed to be incoherent for both H and Y. While the incoherent elastic scattering approximation is appropriate for H-YH 2 , it introduces an undesirable approximation for Y-YH 2 . In this work, we re-evaluate the TSL for Y-YH 2 using FLASSH (Full Law Analysis Scattering System Hub). Y-YH 2 is evaluated using the FLASSH generalized coherent elastic scattering capability in order to capture the Bragg peaks associated with the YH 2 crystal structure which were neglected in the prior NJOY-based evaluation due to limitations in LEAPR. An experimental approach to validate the Y-YH 2 TSL using neutron transmission measurements is discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Validation of Thermal Scattering Laws for Light Water at Elevated Temperatures with Diffusion Experiments

Light water is the most important neutron moderator in many reactor applications. Thermal neutron scattering kernels for H bound in H 2 O are represented by ENDF-format thermal scattering laws (TSLs) evaluated at discrete temperatures T . Nuclear data evaluations are conventionally validated utilizing a combination of critical benchmarks and experimental cross section data. Existing public critical benchmarks may be inadequate for testing water TSLs over the full range of T of interest in reactor applications, and experimental thermal scattering cross section data for water is sparse at elevated T . In this work, MC21 is used to simulate the decay of the equilibrium thermal neutron flux in light-water spheres of several radii. The fundamental-mode time eigenvalue $\mathcal{a}$ is calculated at 22 °C and 227 °C (at saturation pressure) for each sphere using three different H-H 2 O TSLs. Fitting $\mathcal{a}$ to a polynomial function of geometric buckling allows calculation of the thermal neutron diffusion length L . Extrapolation lengths are treated as a function of the transport mean free path and geometry. Experimental L data from 25 publications at 49 temperatures (from 10 °C to 295 °C), computed by several different time-dependent and space-dependent decay methods, is used to develop an empirical fit of L vs. T . The MC21-calculated L for the TSLs tested are within ±1% of the predicted values at 22 °C and 227 °C. This validation approach, which may be repeated at arbitrary T , constitutes an integral benchmark specific to and characterizing the detailed physics of the thermal scattering kernel applied.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Using Prompt Neutron Decay Constant Measurements to Obtain Additional Kinetics Information

Prompt neutron decay constant measurements are typically used to infer the prompt neutron lifetime in a chain-reacting system, but these measurements include much more information about the kinetics of the system. This work focuses specifically on the Rossi-α method of measuring the prompt neutron decay constant. The Rossi-α system of equations consists of five unknown parameters. With five unknowns and three equations, the system of equations solved by the Rossi-α measurement is under-determined. It has been postulated that a well posed problem can be formulated using independent measurement configurations on a single system. The external source strength, the effective delayed neutron fractions, and Diven’s parameter do not typically change for different configurations of a single experiment. If good detector geometry is chosen, the detector efficiency does not change for multiple measurements of a single experiment using Rossi-α measurements. By making these assumptions, there are only two unique variables. This knowledge can be used to solve for most of the unique parameters using two independent, unique Rossi-α measurements.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Fast Rossi-alpha Measurements of Plutonium Using Organic Scintillators

In this work, Rossi-alpha measurements were simultaneously performed with a 3 He-based detection system and an organic scintillator-based detection system. The assembly is 15 kg of plutonium (93 wt% 239 Pu) reflected by copper and moderated by lead. The goal of Rossi-alpha measurements is to estimate the prompt neutron decay constant, alpha. Simulations estimate $k_{eff}$ = 0:624 and α = 52:3 ± 2:5 ns for the measured assembly. The organic scintillator system estimated α = 47:4 ± 2:0 ns, having a 9.37% error (though the 1.09 standard deviation confidence intervals overlapped). The 3 He system estimated α = 37 μs. The known slowing down time of the 3 He system is 35-40 μs, which means the slowing down time dominates and obscures the prompt neutron decay constant. Subsequently, the organic scintillator system should be used for assemblies with alpha much less than 35 μs.

3He↗

Thermal Cross Section Measurements At The RPI LINAC

Recently, a cold moderator was designed and developed for use at the Rensselaer Polytechnic Institute Linear Accelerator (RPI LINAC). This cold moderator proved to easily and safely couple to an existing neutron producing target, while enhancing neutron flux below 0.02 eV by up to a factor of 8 by cooling polyethylene down to 29 K. This cold moderator capability allowed for significantly improved counting statistics below 0.02 eV not previously possible due to a poor signal to background ratio. Additionally, testing was performed to characterize the energy resolution of the new cold moderator system and found the system easily capable of resolving resonances in Ta-181 at 4 and 10 eV, while also clearly resolving the Bragg edges found in Be metal below 0.01 eV. Furthermore, following the design and development of a cold polyethylene moderator, a series of thermal total cross section measurements were performed for polyethylene, polystyrene, Plexiglas and yttrium hydride in the thermal region. These measurements serve to help validate thermal scattering law (TSL) evaluations in the 0.0005 – 1 eV energy range. For polyethylene and polystyrene, two sets of experiments were performed – one with the Enhanced Thermal Target (ETT) and another with the ETT plus the new cold moderator capability (ETTC). The yttrium hydride and Plexiglas measurements were only performed with the ETTC. The measurements for polyethylene help to validate the data processing methodology when using the ETTC, while extending the measured range of polyethylene down to 0.0005 eV. Two different Plexiglas, Plexiglas G and Plexiglas G-UVT, and two different concentrations of yttrium hydride, H/Y = 1.85 and 1.68, were measured. Overall, all materials had generally good agreement with their ENDF/B-VIII.0 TSL evaluations, though some discrepancies were noticed. In the case of the yttrium hydride, the high energy oscillations in the hydrogen cross section and the low energy Bragg edges in the yttrium cross section were clearly seen. These measurements represent the first total cross section measurements that encompass the entire thermal region from 0.0005 – 1 eV for polystyrene and yttrium hydride.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Neutron Capture and Transmission Measurements of 54 Fe at the RPI LINAC

54 Fe radiative capture cross section and transmission measurements were conducted at the Rensselaer Polytechnic Institute (RPI) Gaerttner Linear Accelerator (LINAC) Center using an enriched 54 Fe sample in the keV energy region. 54 Fe is a constituent of natural iron, which is present in a large variety of nuclear grade materials. Therefore, it is important to have an accurate understanding of the cross sections of 54 Fe, which can be measured experimentally. In the time-of-flight measurements conducted at the LINAC, an array of four C 6 D 6 detectors surrounded the sample and radiative capture data were collected using a digital data acquisition system. Additionally, a Li-glass detector was used to collect transmission data using an analog data acquisition system. The radiative capture yield of the 54 Fe measurements were normalized to saturated resonances observed in Au and Ta to obtain an absolute capture yield. The preliminary capture yield and preliminary transmission obtained can be compared to evaluations and existing experimental data. Some disagreements were observed in prominent d-wave capture resonances observed in 54 Fe in the low-keV neutron energy region. Both sets of experimental data along with pre-existing datasets will greatly enhance RPI’s ability to perform resonance evaluation for 54 Fe up to roughly 1 MeV.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Parameterization of Direct and Doorway Processes in R-Matrix Formalism

R-matrix formalism is extended beyond compound nuclear (CN) resonant reactions to include parameterization of direct as well as doorway processes. Direct processes in the R-matrix exterior are parameterized by a unitary matrix that introduces mixing among wave function coefficients of the incoming and outgoing wave function components at the R-matrix channel surface. Doorway processes are parameterized by separating the Hilbert space of the interior R-matrix region into its doorway and CN subspaces, from which doorway state eigenenergies, reduced width amplitudes, and the strengths of their coupling to CN levels appear as new R-matrix parameters. Parameterization of generalized as well as the conventional Reich–Moore approximation for eliminated capture channels in the presence of direct, doorway, and CN processes is presented along with a complex-valued scattering length with contributions from direct, doorway, and CN capture processes. Derivation of Brune’s alternative R-matrix parameters is extended to include doorway states. This work suggests how R-matrix formalism could be extended further by adopting the concepts from related reaction formalisms.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Bayesian Monte Carlo Evaluation of Imperfect (n, 233 U) Data and Model

Conventional nuclear data evaluation methods using generalized linear least squares make the following assumptions: prior and posterior probability distribution functions (PDFs) of all model parameters and data are normal (Gaussian); the linear approximation is sufficiently accurate to minimize the cost function (even for nonlinear models); the model (e.g., of neutron cross section) and experimental data (including covariance data) are without defect and prior PDFs of parameters and measured data are known perfectly. Neglect of covariance between model parameters and measured data in conventional evaluations contributes to imperfections. These assumptions are inherent to the generalized linear least squares minimization method commonly used for resolved resonance region neutron cross section evaluations but are often not justified due to the presence of non-normal PDFs, nonlinear models (e.g., R-matrix formalism), and inherent imperfections in data and models (e.g., imperfect covariance data). Here, these assumptions are removed in a mathematical framework of Bayes’ theorem, which is implemented using the Metropolis-Hastings Monte Carlo method. Most importantly, new parameters are introduced to parameterize discrepancies between the theoretical model and measured data to quantify judgement about discrepancies or imperfections in a reproducible manner. An evaluation of 233U in the eV region using the ENDF-B/VIII.0 library and transmission data (Guber et al.) is presented, and posterior parameters are compared to those obtained by conventional evaluation methods. This example illustrates the effects of removing the most harmful assumption: that of model-data perfection.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗