SEARCH · Engineering Papers
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.
New cross section measurements for reactions induced by neutrons with energies greater than 70 MeV are used to calculate the production rates for cobalt isotopes in meteorites and these new calculations are compared to previous estimates. Additional information is contained in the original extended abstract.
The 93 Nb (𝑡, 3 He ) reaction at 115 MeV/nucleon was studied to demonstrate that nuclear level densities and 𝛾-ray strength functions can be extracted from charge-exchange reactions at intermediate energies using the Oslo technique. The matrix of excitation energy in 93 Zr, reconstructed from the (𝑡, 3 He ) reaction, versus the energy of 𝛾 rays emitted by the excited 93 Zr nuclei, was obtained in an experiment with the S800 Spectrograph operated in coincidence with the GRETINA 𝛾 -ray detector. The extracted level density and 𝛾 -ray strength function obtained by applying the Oslo method to this matrix were used to estimate the 92 Zr(𝑛,𝛾) 93 Zr cross section by combining the new results with other experimental data and theoretical calculations for 𝐸1 and 𝑀1 strength functions at higher energies. Good agreement with direct measurements of the 92 Zr(𝑛,𝛾) 93 Zr cross section was found. The contribution from the upbend in the extracted 𝛾-ray strength function was important to achieve the consistency, as the neutron-capture cross section without this contribution is significantly below the direct measurements otherwise. Since charge-exchange reactions at intermediate energies have long been used for extracting Gamow-Teller strengths, the successful demonstration of the charge-exchange Oslo method enables experiments in which (𝑛,𝛾) cross sections and Gamow-Teller strengths can be measured simultaneously, which is of benefit for astrophysical studies.
Uranium silicide (U 3 Si 2 ) is a candidate material for the high-density nuclear fuel in commercial light water reactors [1], [2]. Its higher uranium density, 11.3 g-U/cm3, compared to that of uranium dioxide (UO 2 ), 9.7 g-U/cm3, can improve the performance of a nuclear reactor while using low enriched uranium (LEU) and diversify the choice of cladding materials [1]–[3]. It also has a higher thermal conductivity than UO 2 , which can reduce the thermal stress on the material caused by a temperature gradient across the fuel pellet and provide a larger margin for some postulated accidents [1], [2], [4]. Furthermore, compared to U3Si, another high-density fuel candidate, it has better resistance to in-pile swelling due to less irradiation-induced rapid amorphization [1], [3]. Corresponding to its importance in nuclear engineering, many previous studies have reported the properties of U3Si2. Experiments showed that U 3 Si 2 is a paramagnetic (PM) metal, where a slight linear increase in magnetic susceptibility was measured with increasing temperature [5], [6]. In addition, thermodynamic quantities such as thermal expansion coefficient, heat capacity, and thermal conductivity were experimentally determined over a wide temperature range [1], [7], [8]. In several computational studies, ab initio atomistic simulations based on density functional theory (DFT) were performed to calculate various properties including elastic constants, electronic density of states (DOS), and phonon dispersion curves [9]–[12]. Nevertheless, thermal neutron scattering cross sections, which are critical to the prediction of the parameters in reactor physics that are ultimately related to reactor criticality, have not yet been evaluated for U3Si2. The scattering cross section can be calculated from the phonon DOS, or the energy spectrum of lattice vibrations, of the crystalline system [13], [14]. However, there is also no experimental data available for the phonon DOS of U 3 Si 2 . While some computational studies reported the phonon DOS and/or dispersion curves from ab initio simulations [9]–[12], the accuracy cannot be guaranteed because it is unclear whether the spin-polarization behavior of PM U 3 Si 2 was properly described. In the present study, the thermal neutron scattering cross section for U 3 Si 2 is evaluated for the first time by calculating the phonon DOS for U3Si2 from ab initio lattice dynamics (AILD) simulations based on DFT. First, U 3 Si 2 is modeled based on the experimental structure, and AILD simulations are performed on the modeled U3Si2 to optimize the structure. Next, AILD simulations are performed for supercells with atomic displacement to calculate Hellmann-Feynman forces. Based on the calculated forces, partial phonon DOSs for U and Si are obtained, and the thermal neutron scattering law (TSL) for U 3 Si 2 is finally evaluated. To verify the accuracy of the calculations in the present study, the calculation results are compared with experimental data on the structure and heat capacity of U3Si2 [1], [7], [8], [15].
The materials which compose the ICRP Voxel phantoms used in the computation of conversion coefficients involve neutron interaction cross-sections that have resonances at specific energies. Depending on the energy bin structure used in the computations, these cross-section resonances may occur at energies that fall between energies at which dose coefficients are computed, thus their effects may not be completely accounted for in the reported coefficients. Here In this study, a highly refined energy grid that closely follows the resonance structure in the phantom material cross-sections was identified and used to calculate dose coefficients. Both the equivalent organ/tissue doses for male and female voxel phantoms were computed as well as their summation to obtain the effective dose coefficients. The used refined energy grid tracks very closely the cross-sections in the vicinity of the resonances. The resulting refined energy grid coefficients are compared to coefficients for the coarser energy grid used in ICRP Publication 116. Additionally, reference spectra have been folded with both the fine and coarse sets of conversion coefficients. The resulting total effective doses for these reference spectra are used to assess the adequacy of the dose coefficients calculated on the original ICRP 116 energy grid. The dose coefficients were similarly computed for the local skin dose on the trunk of the body using the ICRU Report 95 phantom. The overall impact of the resonances on the organ/tissue equivalent dose, the effective dose, and the local skin dose are presented and discussed. In general, it was found that resonances can impact neutron dose coefficients, but in most cases the wide range of neutron energies encountered minimized this effect. The impact of resonances was further limited when computing effective dose due to organ/tissue summing and sex-averaging. For the neutron fields studied here, the impact was below 5%.
A best set of neutron-capture cross sections has been evaluated for the most important s-process isotopes. With this data base, s-process studies have been carried out using the traditional model which assumes a steady neutron flux and an exponential distribution of neutron irradiations. The calculated sigma-N curve is in excellent agreement with the empirical sigma-N-values of pure s-process nuclei. Simultaneously, good agreement is found between the difference of solar and s-process abundances and the abundances of pure r-process nuclei. The abundance pattern of the iron group elements where s-process results complement the abundances obtained from explosive nuclear burning is discussed. The results obtained from the traditional s-process model such as seed abundances, mean neutron irradiations, or neutron densities are compared to recent stellar model calculations which assume the He-burning shells of red giant stars as the site for the s-process.
In this work, the double-differential cross section of the deuteron breakup reaction 2 H(n, 2n) has been studied experimentally with a neutron beam energy of 15 MeV. Special attention has been devoted to estimation of background condition and multiple scattering effect in the scattering sample. Experimental data have been compared with models based on phase-space approximation used in the ENDF/B-VIII.0 data library and in the MCNP neutron transport code, as well as with rigorous model based on Faddeev equations used for cross section evaluations in JENDL data library. It was found that experimental data are better reproduced by Faddeev model, however, the model overestimates data in the low-energy region of the neutron spectrum (<4 MeV)
Current and future accelerator neutrino oscillation experiments need neutrino interaction models with smaller systematic uncertainties to resolve much of delta CP phase space. Final state interactions (FSI) and scattering off of correlated nuclei (2p2h) are poorly understood processes that currently contribute large uncertainties to leading models. These processes have proven difficult to study because they often produce relatively low energy nucleons. Protons up to about 100 MeV are below the detection threshold of some accelerator neutrino detectors, and neutrons are usually discounted as undetectable. This poster presents a measurement of the multi-neutron antineutrino cross section at low available energy using the MINERvA detector at Fermilab. This interaction channel is particularly sensitive to FSI and 2p2h interactions. A sideband-driven background constraint that greatly reduces uncertainties on the result will be presented. The measured cross section is compared to GENIE v3 models with different FSI treatments and the SuSA model's 2p2h predictions.
The primary goal of this study is to develop a national facility concept for transmuting long-lived fission products (LLFP) to substantially reduce the disposal impact by minimizing the need for a geologic-timescale repository. As a charter for this study, the national transmutation facility was required to reduce the radiotoxicity and decay heat of LLFP isotopes by at least 90% relative to their values at discharge from a commercial LWR, while consuming less than 10% of the reactor's energy. The identified LLFP isotopes are Se-79, Zr-93, Tc-99, I-129, Sn-126, and Cs-135, whose radiotoxicity is about 99% of the total radiotoxicity of all fission products at 1,000 years. Approximately ~72 kg of LLFPs is discharged every year from a 1,000 MWe commercial or advanced nuclear reactor. First, LLFP transmutation options with non-neutron beams (photons and protons) were explored. The study concluded that LLFP transmutation is feasible with high-energy, high-intensity photons or protons, but impractical on an engineering scale due to low transmutation rates and the high energy requirements to produce the desired photon or proton beams. As alternatives, LLFP transmutation options with neutrons from fission, fusion, and spallation reactions were additionally explored. The transmutation options using advanced critical reactors are attractive only for selective LLFP isotopes because the production rates of several LLFP isotopes (Zr-93, Sn-126, and Cs-135) from fission reactions are larger than the transmutation rates. The transmutation options with only spallation neutrons are favorable to transmute all LLFP isotopes, but as a tradeoff, the net transmutation rates are reduced. The national transmutation facility concept was developed following an exploration of transmutation options using various incident particles. The proposed national LLFP transmutation comprises a dedicated molten-salt reactor (MSR), a proton accelerator, and a spallation neutron-based transmuter. The MSR power was set at 300 MWt and 120 MWe, with the thermal power approximately 10% of that of a commercial 1,000 MWe PWR. The electricity generated by the MSR powers the accelerator and transmuter. The accelerator produces 1 GeV, 30 mA protons, which are introduced into the spallation neutron-based transmuter. The spallation neutron-based transmuter consists of a central spallation target and LLFP target pins merged in a heavy water tank. The six LLFP isotopes are separated into two groups. Tc-99, I-129, and Se-79, having larger neutron cross sections, belong to group A, while Zr-93, Sn126, and Cs-135, having smaller neutron cross sections, belong to group B. Then, for effective transmutation, LLFPs in groups A and B are transmuted in the dedicated MSR and in a spallation neutron-based transmuter, respectively. The estimated capital cost of the national transmutation facility is approximately $\$$3.1B, and its annual O&M cost is expected to be ~$\$$182M. Radiotoxicity and decay heat of LLFPs were calculated and compared with those of the original LLFPs. It was assumed that the targets were made with elementwise LLFP rather than isotopic LLFP, owing to the potentially high cost of isotopic separation from used nuclear fuels. The decay heat of LLFPs can be reduced by more than 90% using a single national transmutation facility. However, radiotoxicity decreases by 79–84%, which does not meet the transmutation performance requirement, primarily because Cs-135 is produced rather than depleted. Thus, to meet the design requirement, Cs-135 should be separated from other Cs isotopes and irradiated in a spallation neutron-based transmuter. Then, radiotoxicity decreases by ~92%.
Solid metal hydrides have long been considered viable moderators for nuclear reactor designs due to their moderating ratios, hydrogen densities, high dissociation temperatures, and mechanical properties. CaH 2 is an orthorhombic saline-hydride that has recently been investigated and shows promise for use as a moderator in microreactors. At present, there is no Thermal Scatting Law (TSL) evaluation for CaH 2 in the ENDF/B-VIII.0 database. An evaluation does exist in the JEFF-3.3 database, performed by Serot, however it was limited by the implemented methods and thus both drastically over-predicts the incoherent elastic contribution and completely ignores the coherent elastic contribution to scattering from the metal ions. The thermal neutron scattering cross sections of CaH 2 are evaluated in this work to establish accurate data for use in reactor design. The present evaluation corrects the inaccuracies of the JEFF-3.3 data, and yields three distinct libraries: Ca in CaH 2 , H 1 in CaH 2 , and H 2 in CaH 2 .
Carbon materials are commonly found in both nuclear reactors and experimental systems. Various carbon structures occur in nuclear applications ranging from crystalline and nuclear graphite to the amorphous carbon seen in next-generation advanced reactor designs. Amorphous carbon is based on a randomized graphite-like structure and offers the unique ability to disperse impurities throughout the bulk composition. A graphite-like amorphous carbon system was modeled using the classical molecular dynamics (MD) code LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator). An improved version of the temperature-dependent Adaptive Intermolecular Reactive Empirical Bond Order (AIREBO) potential was used to model the carbon-carbon atomic interactions for the temperature at 300 K along with densities 1.60, 1.70, 1.85, and 2.23 g/cm{sup 3}. From the normalized velocity autocorrelation function (VACF), the phonon density of state (DOS) was then calculated as the Fourier transform of the normalized VACF. This DOS was then used as the primary input for the evaluation of the thermal scattering law (TSL, i.e. S(α,β)) and associated neutron thermal scattering cross sections. The TSL was analyzed using the Full Law Analysis Scattering System Hub (FLASSH). The amorphous structure results in shifts of the phonon DOS to lower energy modes than typically displayed for ideal crystalline graphite. This impact on the DOS is directly reflected in the TSL. Furthermore, the typical optical peak at 0.25 eV for the ideal graphite disappears for amorphous carbon, in good agreement with the expected structure. (authors)
The capability to utilize neural networks to predict macroscopic and microscopic cross section parametric spaces has been developed for the Griffin reactor physics application. The LibTorch interface enables Griffin's MOOSE-based materials to interact with LibTorch-trained models, allowing for the evaluation of complex macroscopic or microscopic cross section spaces, which are then used to evaluate the neutronic properties of the Griffin finite element model. This study benchmarks traditional ISOXML-formatted tabulation libraries against neural network-based models for 279 nuclides on 20,160 grid points for zero-dimensional and two-dimensional reactor models. Benchmark metrics include the fundamental mode eigenvalue, fission and absorption rates, and various temperature coefficients of reactivity (isothermal, fuel, and moderator). From the perspective of storage space, the complete set of LibTorch models uses 11 MB on disk, compared to the 10 GB for the ISOXML multigroup library that covers the same grid space. For the two-dimensional performance case considered in Griffin, the Torch model uses 97% less RAM than the reference ISOXML dataset while runtime increases by a factor of 3 when using the LibTorch model compared to the ISOXML dataset with multi-linear interpolation. The LibTorch model consistently yields errors within 0.01% for most analyzed quantities except for the temperature coefficients of reactivity where the maximum discrepancies are up to 0.3 $\frac{pcm}{K}$. Due to the neural network attempting to best predict quantities with no regard for a positive or negative bias for any given quantity, predictions may experience random fluctuations, resulting in both positive and negative errors. Future work will entail both depletion and coupled transient analysis to determine the predictive capabilities of Griffin with neural network-based cross sections.
Total thermal neutron cross section measurements serve as the primary means of validation for thermal neutron scattering kernels, an important quantity for neutron transport calculations. In an effort to improve the quality of thermal neutron scattering kernels, researchers at Rensselaer Polytechnic Institute (RPI) designed and constructed a polyethylene based cold moderation system to enhance neutron flux below 10 meV when coupled with the Enhanced Thermal Target (ETT) at the RPI Gaerttner LINAC. The final design yielded an increase in sub-thermal neutron flux (below 10 meV) by a factor of 4.5 for a moderator temperature of 37.5 K relative to the ETT alone. A further increase to a factor of 6 is expected after a minor geometry modification and decrease in polyethylene temperature to 25 K. This novel capability will be used to conduct total thermal neutron cross section measurements from 0.0005–10 eV for different materials including moderator materials.
Carbon materials are commonly found in both nuclear reactors and experimental systems. Various carbon structures occur in nuclear applications ranging from crystalline and nuclear graphite to the amorphous carbon seen in next-generation advanced reactor designs. Amorphous carbon is based on a randomized graphite-like structure and offers the unique ability to disperse impurities throughout the bulk composition. A graphite-like amorphous carbon system was modeled using the classical molecular dynamics (MD) code LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator). An improved version of the temperature-dependent Adaptive Intermolecular Reactive Empirical Bond Order (AIREBO) potential was used to model the carbon-carbon atomic interactions for the temperature at 300 K along with densities 1.60, 1.70, 1.85, and 2.23 g/cm 3 . From the normalized velocity autocorrelation function (VACF), the phonon density of state (DOS) was then calculated as the Fourier transform of the normalized VACF. This DOS was then used as the primary input for the evaluation of the thermal scattering law (TSL, i.e. S(α,β)) and associated neutron thermal scattering cross sections. The TSL was analyzed using the Full Law Analysis Scattering System Hub (FLASSH). The amorphous structure results in shifts of the phonon DOS to lower energy modes than typically displayed for ideal crystalline graphite. This impact on the DOS is directly reflected in the TSL. Furthermore, the typical features and the optical graphitic peak at 0.25 eV for the ideal graphite DOS disappear for graphite-like amorphous carbon, which shows good agreement with the expected structure.
Recent parameterization of absorption cross sections for any system of charged ion collisions, including proton-nucleus collisions, is extended for neutron-nucleus collisions valid from approx. 1 MeV to a few GeV, thus providing a comprehensive picture of absorption cross sections for any system of collision pairs (charged or uncharged). The parameters are associated with the physics of the problem. At lower energies, optical potential at the surface is important, and the Pauli operator plays an increasingly important role at intermediate energies. The agreement between the calculated and experimental data is better than earlier published results.
Tungsten neutron cross section, resonance parameters, and scattering matrixes for certain nuclear reaction concepts
Reaction cross sections of Si/n,p/ and Si/n,alpha/ at 14.5 MeV neutron energy
Neutrinos are among a handful of particles in the Standard Model of particle physics whose properties have not yet been fully measured. Neutrino oscillations, discovered in the late 1990s, still have yet to be fully characterized. Early evidence suggests that neutrinos and antineutrinos may undergo neutrino oscillations differently. So, multiple collaborations are building and operating kTon-scale particle detectors to try to measure this CP violation in the neutrino sector of the Standard Model. Measuring neutrino and antineutrino energies is essential to characterizing oscillation phenomena. But neutrons interfere with (anti)neutrino energy measurements because they are hard to detect in (anti)neutrino detectors. Furthermore, there are very few measurements of neutrons produced by (anti)neutrinos in the energy range of interest to accelerator-based oscillation experiments. The MINERvA collaboration initially tried to fill this gap by measuring the rate of neutron-induced activity i n plastic. The data deviated significantly from leading model predictions, but it could not be determined which models needed to be tuned. MINERvA’s first neutron counting result also cannot be compared to future models because it was not suitable for unsmearing. This thesis first expands upon MINERvA’s first foray into neutron detection at higher antineutrino energy. Neutron detection efficiency is enhanced by advances in analysis algorithms. Then, an unsmeared cross section for producing two or more neutrons is measured. Existing antineutrino interaction models again over-predict the rate of neutron production by antineutrinos. Preliminary neutron production rate measurements in MINERvA’s iron and lead targets suggest that models are also over-predicting neutron production in heavier nuclei.
We present a detailed account of neutron capture experiments of astrophysical relevance on 78,80,84,86 Kr (n, γ) reactions at the border between weak and main s process. The experiments were performed with quasi-Maxwellian neutrons from the Liquid-Lithium Target (LiLiT) and the mA-proton beam at 1.93 MeV (2-3 kW) of the Soreq Applied Research Accelerator Facility (SARAF). The setup yields high-intensity approximate to 40 keV quasi-Maxwellian neutrons (3-5 x10 10 n/s) closely reproducing the conditions of s-process stellar nucleosynthesis. A sample of 100 mg of atmospheric, pre-nuclear-age Kr gas contained in a Ti spherical shell was activated in the LiLiT neutron field. The abundances of long-lived Kr isotopes ( 81,85g Kr) were measured by atom counting via atom trap trace analysis (ATTA) at Argonne National Laboratory and low-level counting (LLC) at University of Bern. This work is the first measurement of a nuclear cross section using atom counting via ATTA. The activities of short-lived Kr isotopes ( 79,85m,87 Kr) were measured by gamma-decay counting with a high-purity germanium detector. Maxwellian-averaged cross sections for s-process thermal energies are extracted. By comparison to reference values, our nucleosynthesis network calculations show that the experimental cross sections have a strong impact on calculated abundances of krypton and neighboring nuclides, in some cases improving agreement between theory and observations.