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.
In this report, we summarize the progress made in FY23 for experiments fielded at Triangle Universities Nuclear Laboratory (TUNL) of Duke University on the neutron scattering measurements of 238 U. The preliminary results on the elastic and inelastic reaction channels will be presented in addition to details of experiment and data analysis.
Method for producing neutrons in keV energy region utilizing gamma rays from radioactive decay to induce photoneutron reaction in beryllium
The time dependences of the 2.2- and 0.51-MeV gamma-ray lines from solar flares are calculated, and the results are compared with observations of the 1972 August 4 and 7 flares. The time lag between the nuclear reactions and the formation of these two lines is caused by capture of the neutrons and subsequent deceleration of the positrons and decay of the radioactive nuclei. Our main results are that the calculation is consistent with the observed rise of the 2.2-MeV line on August 4, and it does not require different time dependences for the accelerated protons and high-energy electrons in the flare region. The above lags can explain the delayed gamma-ray emission observed on August 7. Positrons of energies greater than about 10 MeV could be detected in interplanetary space following large solar flares.
One of the major advantages of using P-B11 fusion fuel is that the reaction produces only charged particles in the form of three alpha particles and no neutrons. A fusion concept that lends itself to this fuel cycle is the Magnetically Insulated Inertial Confinement Fusion (MICF) reactor whose distinct advantage lies in the very strong magnetic field that is created when an incident particle (or laser) beam strikes the inner wall of the target pellet. This field serves to thermally insulate the hot plasma from the metal wall thereby allowing thc plasma to burn for a long time and produce a large energy magnification. If used as a propulsion device, we propose using antiprotons to drive the system which we show to be capable of producing very large specific impulse and thrust. By way of validating the confinement propenies of MICF we will address a proposed experiment in which pellets coated with P-B11 fuel at the appropriate ratio will be zapped by a beam of antiprotons that enter the target through a hole. Calculations showing the density and temperature of the generated plasma along with the strength of the magnetic field and other properties of the system will be presented and discussed.
Although the exact site for the origin of the r-process isotopes remains mysterious, most thinking has centered on matter ejected from the cores of massive stars in core-collapse supernovae [13]. In the 1970's and 1980's, however, difficulties in understanding the yields from such models led workers to consider the possibility of r-process nucleosynthesis farther out in the exploding star, in particular, in the helium burning shell [4,5]. The essential idea was that shock passage through this shell would heat and compress this material to the point that the reactions 13C(alpha; n)16O and, especially, 22Ne(alpha; n)25Mg would generate enough neutrons to capture on preexisting seed nuclei and drive an "n process" [6], which could reproduce the r-process abundances. Subsequent work showed that the required 13C and 22Ne abundances were too large compared to the amounts available in realistic models [7] and recent thinking has returned to supernova core material or matter ejected from neutron star-neutron star collisions as the more likely r-process sites.
This team is exploring a modified Z-pinch geometry as a propulsion system, imploding a liner of liquid lithium onto a pellet containing both fission and fusion fuel. The plasma resulting from the fission and fusion burn expands against a magnetic nozzle, for propulsion, or a magnetic confinement system, for terrestrial power generation. There is considerable synergy in the concept; the lithium acts as a temporary virtual cathode, and adds reaction mass for propulsion. Further, the lithium acts as a radiation shield against generated neutrons and gamma rays. Finally, the density profile of the column can be tailored using the lithium sheath. Recent theoretical and experimental developments (e.g. tailored density profile in the fuel injection, shear stabilization, and magnetic shear stabilization) have had great success in mitigating instabilities that have plagued previous fusion efforts. This paper will review the work in evaluating the pellet sizes and z-pinch conditions for optimal PuFF propulsion. Trades of pellet size and composition with z-pinch power levels and conditions for the tamper and lithium implosion are evaluated. Current models, both theoretical and computational, show that a z-pinch can ignite a small (~1 cm radius) fission-fusion target with significant yield. Comparison is made between pure fission and boosted fission targets. Performance is shown for crewed spacecraft for high speed Mars round trip missions and near interstellar robotic missions. The PuFF concept also offers a solution for terrestrial power production. PuFF can, with recycling of the effluent, achieve near 100% burnup of fission fuel, providing a very attractive power source with minimal waste. The small size of PuFF relative to today's plants enables a more distributed power network and less exposure to natural or man-made disruptions.
Generation IV reactors and future fusion reactor designs have led to more demanding materials performance requirements due to their increased operating temperatures, corrosive coolants, and increased radiation doses compared to the current light-water reactor fleet. Among the innovative nuclear technologies under development, molten salt reactors stand out for their potential to offer superior fuel utilization, intrinsic safety characteristics, and economic viability. Of the proposed Generation IV designs, the gas fast reactor operates at 450 to 850°C and the molten salt reactor operates at 565 to 850°C, with the molten salt reactor design needing molten salt corrosion resistant materials [1, 2]. These increased temperatures and more extreme corrosion environments necessitate higher material performance, such as creep strength, radiation-tolerant microstructures, corrosion resistance, and high-temperature tensile properties. Hastelloy-N, a nickel-based alloy with additions of molybdenum and chromium, has been successfully employed to contain molten fluoride salt at temperatures up to 705°C. However, Hastelloy-N becomes embrittled upon neutron irradiation, primarily due to the accumulation of helium produced by (n,a) transmutation reactions. Furthermore, the corrosive nature of molten fluoride and chloride salts presents a formidable challenge, as these salts can react with and dissolve alloying elements such as Cr, Mo, and Fe, leading to selective leaching, loss of protective oxide layers, and accelerated degradation. High entropy alloys (HEAs) and refractory high entropy alloys (RHEAs) have emerged as a prominent area of interest, due to their ability to achieve tailored chemical compositions for specific applications. Unlike conventional alloys, HEAs are characterized by having multiple principal elements in equimolar or near equimolar ratios, leading to an unconventional alloying strategy [3]. This alloying strategy is believed to promote unique properties, such as single-phase stabilization of chemically compatible elements, lattice distortion effects due to atomic radius differences, and proposed sluggish diffusion effects. For extreme-environment applications, RHEAs have garnered much research interest because of the possibility of creating relatively ductile materials that can operate in extremely high-temperature environments, beyond the operating temperatures where other Ni-based superalloys begin to lose strength [4-6]. Idaho National Laboratory (INL) initiated a joint international effort with the Czech Republic to explore the feasibility of manufacturing HEAs for high-temperature nuclear applications using advanced manufacturing. This effort was funded at INL by the United States Department of Energy's Office of Nuclear Energy under the Advanced Reactor Technologies and Advanced Materials and Manufacturing Technologies (AMMT) Program. The HEAs were specifically designed for the corrosive and irradiation environments experienced in gas-cooled fast reactors, molten salt reactors, and fusion power. These alloys have been manufactured by multiple processes to determine the impact of manufacturing processes on the performance of the alloys in corrosive and irradiation environments. Preliminary molten salt corrosion testing showed that equimolar MoNbTiV and MoNbTi alloys exhibit exceptional performance, with arc-melted variants demonstrating only minimal degradation after 1000 hours of exposure to molten chloride salt at 700°C. Conversely, Nb2TiVZr2 showed significant molten salt corrosion susceptibility and microstructural instability during high-temperature molten salt exposures, and was, therefore deemed unfit for molten salt reactor applications. The MoNbTiV, MoNbTi, and Nb2TiVZr2 alloys were further evaluated through ion irradiation experiments conducted at the Michigan Ion Beam Laboratory at the University of Michigan. The microstructural stability and the evolution of irradiation-induced defects were characterized to assess the irradiation resistance of each of these alloys.
Uranium-233 plays an important role in the Th-U fuel cycle, with substantial production in the cycle. This cycle has been proposed as an alternative to the U-Pu fuel cycle due to its reduced production of transuranic elements. An accurate measurement of the 233 U($\eta$,$\gamma$) cross section is required by the National Criticality Safety Program (NCSP) to complete the neutron-induced cross section data, where experimental capture cross section data are scarce and were measured decades ago. The most recent capture cross section data available in the literature were measured in 2007 at the n_TOF facility (CERN); in the 60s measurements were performed at Rensselaer Polytechnic Institute (RPI) and at LANL. Finally, as reported by ORNL, a new evaluation with a revised (renormalized) fission cross section is needed on 233 U. The challenge for this measurement lies in the difficulty of measuring the capture cross section data in the competing fission background, as the fission cross section is around one order of magnitude larger than the capture cross section for 233 U. The accuracy of a capture cross section measurement depends on discrimination between $\gamma$’s produced in capture and fission reactions, for which an experimental setup combining capture and fission detectors is needed. For the ($\eta$,$\gamma$) cross section measurement at LANSCE, this discrimination is achieved by combining the Detector for Advanced Neutron Capture Experiments (DANCE), to measure $\gamma$’s from capture reactions, with a Parallel Plate Avalanche Counter (PPAC) to tag the $\gamma$’s produced by fission. This method was successfully used to measure 235 U and 239 Pu capture cross sections. In these measurements, the neutron capture cross section was determined in a large fission background well above 100 keV. As part of the NCSP nuclear data effort, we have looked at past DANCE measurements on 233 U($\eta$,$\gamma$) and evaluated whether existing data is adequate to apply this technique.
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.
A framework for A(d, p)B reactions is introduced by merging the microscopic approach to computing the properties of the nucleon-target systems and the three-body n + p + A reaction formalism, thus providing a consistent link between the reaction cross sections and the underlying microscopic structure. In this first step toward a full microscopic description, we focus on the inclusion of the neutron-target microscopic properties. The properties of the neutron-target subsystem are encapsulated in the Green's function which is computed with the coupled Cluster theory using a chiral nucleon-nucleon and three-nucleon interactions. Subsequently, this many-body information is introduced in the few-body Green's Function Transfer approach to (d, p) reactions. Our benchmarks on stable targets 40,48 Ca show an excellent agreement with the data. We then proceed to make specific predictions for (d, p) on neutron rich 52,54 Ca isotopes. These predictions are directly relevant to testing the new magic numbers N = 32, 34 and are expected to be feasible in the first campaign of the projected FRIB facility.
Within material test reactor calculations, energy dependent flux and reaction rate uncertainties are typically not quantified when performing as-run analyses to determine the neutron field experienced by the experiment. When high fidelity Monte-Carlo codes are used in such analyses, straight forward methods to calculate output uncertainties are not available, instead expert opinion is used to postulate computational uncertainties. New methods to propagate uncertainties through these high fidelity simulations are available when sufficient computational power is available. A tool is developed here for sampling any part of an MCNP input from random distributions to determine output uncertainties based on those inputs. Another tool is developed to sample nuclear data cross-section in ACE format using multi-group nuclear data covariances. The Total Monte-Carlo Method and Gesellschaft für Anlagen-und Reaktorsicherheit method (GRS) are implemented and compared to one another as well as MCNP sensitivity and uncertainty calculations. The methods were applied to the Godiva critical sphere k-eigenvalue, the UAM pin-cell benchmark energy dependent flux and reaction rates, and the Advanced Test Reactor energy dependent flux within an experimental location. Furthermore, the two methods agree well, with GRS allowing for an order of magnitude speedup for reaction rate uncertainty calculations and several orders of magnitude for eigenvalue uncertainty calculations.
One means of improving performance for electrochemical ammonia production through the Li-mediated N 2 reduction reaction (Li-NRR) is by cycling the current driving the reaction between open-circuit conditions and periods of applied current density. Herein, we have investigated the dynamics of the electrode–electrolyte interface under Li-NRR conditions during current cycling using in situ time-resolved neutron reflectometry and grazing-incidence synchrotron X-ray diffraction. During cycling, measured neutron reflectivity curves indicated bilayer formation in which Li-containing species such as LiOH, Li 2 O, and small quantities of Li 3 N and metallic Li primarily appeared in a thin layer at the cathode surface, above which formed a much larger, porous, ‘solid–electrolyte interface’ (SEI) layer. Upon return to open-circuit conditions, Li-containing species quickly moved out of the thin layer, leaving a compact, stable layer of decomposition products underneath the SEI layer. This SEI layer concomitantly filled with electrolyte or dissolved, becoming indistinguishable from the electrolyte via contrast in scattering-length density (SLD). During the second current cycle, Li-containing species again preferentially deposited directly atop the cathode, with the thick SEI-like layer again appearing within a minute. This SEI layer exhibited a lower SLD more quickly than in the first cycle, which might suggest that Li-containing species become distributed within the porous SEI layer. Furthermore, these time-resolved observations of SEI and plated layers during current cycling suggest that benefits associated with return to open-circuit conditions between periods of applied current density may be related to the concomitant loss of Li-containing species from a thin layer at the cathode surface into a porous SEI layer that becomes filled with electrolyte or dissolves.
A novel methodology is developed to extract accurate skeletal reaction models for nuclear combustion. Local sensitivities of isotope mass fractions with respect to reaction rates are modeled based on the forced optimally time-dependent (f-OTD) scheme. These sensitivities are then analyzed temporally to generate skeletal models. The methodology is demonstrated by conducting skeletal reduction of constant density and temperature burning of carbon and oxygen relevant to Type Ia supernovae (SNe Ia). The 495-isotopes Torch model is chosen as the detailed reaction network. A map of maximum production of 56 Ni in SNe Ia is produced for different temperatures, densities, and proton-to-neutron ratios. The f-OTD simulations and the sensitivity analyses are then performed with initial conditions from this map. A series of skeletal models are derived and their performances are assessed by comparison against currently existing skeletal models. Previous models have been constructed intuitively by assuming the dominance of α-chain reactions. The comparison of the newly generated skeletal models against previous models is based on the predicted energy release and 44 Ti and 56 Ni abundances by each model. The consequences of ye ≠ 0.5 in the initial composition are also explored where ye is the electron fraction. The simulated results show that 56 Ni production decreases by decreasing ye as expected, and that the 43 Sc is a key isotope in proton and neutron channels toward 56 Ni production. It is shown that an f-OTD skeletal model with 150 isotopes can accurately predict the 56Ni abundance in SNe Ia for ye ≲ 0.5 initial conditions.
The experiments analyzed in this report were conducted at the Health Physics Research Reactor (HPRR), also known as the $\textit{Fast Burst Reactor}$. The reactor was designed and built at Oak Ridge National Laboratory (ORNL) in 1961. The HPRR was an unmoderated, unshielded fast reactor that used highly enriched uranium and molybdenum alloy as fuel. The reactor was initially sent to the Nevada Test Site in 1962, where it was used to evaluate radiation doses received as a result of the Hiroshima and Nagasaki bombings during World War II. A few years later, the reactor was sent back to ORNL to be part of the Dosimetry Application Research (DOSAR) facility shown in Figure 1, which included a reactor building shown on the left (west) of the picture and a control and laboratory building in the upper right corner (northeast). The critical assembly was used for numerous technical studies, including systems calibration, dosimetry, radiobiology of plants and animals, testing of radiation alarms, as well as teaching and training in radiation dosimetry and nuclear engineering. Between 1963 and 1987, the HPRR was operated for thousands of hours, achieving criticality close to 10,000 times and motivating many publications. The HPRR was decommissioned in 1987. The goal of this effort was to use historical data from operation of the HPRR to create a criticality accident alarm system (CAAS) benchmark to be included in the $\textit{International Handbook of Evaluated Criticality Safety Benchmark Experiments}$ (ICSBEP Handbook). A thorough inspection was performed of all available documentation and information available. The most promising experiments that were selected for evaluation were those described in the 1987 ORNL report entitled $\textit{Health Physics Research Reactor Reference Dosimetry}$, ORNL-6240. The report includes reference dosimetry results of the shielded and unshielded configurations of the HPRR after burst operations. Because of changes to the reactor positioning and storage systems that were made in 1985, the previous dosimetry reports became obsolete, and the newly designed experiments were needed to create the HPRR’s adjusted dosimetry data. The various results reported in ORNL-6240 include reference doses and dose equivalents from different conventions at different distances and elevations as determined using the detected neutron fluence and conversion factors. The HPRR neutron fluence was obtained through different methods, including sulfur pellet analysis and threshold detector unit data. Information about the HPRR spectrum was also obtained through Bonner sphere measurements. This benchmark is focused on a part of the measured sulfur fluences reported in Appendix H of ORNL-6240. Standard commercial sulfur pellets were placed at different distances from the HPRR centerline during burst operation and were activated due to the 32 S(n,p) 32 P reaction. The resulting 32 P activity was then measured and the information about the corresponding sulfur fluence and/or neutron dose could be extracted. Many of those measurements have 7 been performed with the HPRR in its bare configuration or with different shields (combinations of Lucite, concrete, steel). All the necessary, precise information about material and/or dimensions of the different shields was not found, so it was decided to focus only on the unshielded and steel-shielded configurations to minimize the benchmark uncertainty. A total of 31 cases (24 unshielded and 7 shielded cases at different positions) of sulfur fluence were selected before evaluation to develop the benchmark.
Accurate knowledge of the 237 Np(n, γ) 238 Np cross section at fast neutron energies is important for applied nuclear science. The presently available experimental data has large disagreements in the fast neutron region. Perform a model-independent measurement of the 237 Np(n, γ) 238 Np integral cross section using a well characterized fast neutron source and compare the result with previous measurements and current nuclear data evaluations. Provide an integral measurement that can be used as a benchmark for current evaluations. Multiple samples of 237 Np were irradiated in the Godiva-IV critical assembly. Following the irradiation, the samples placed in a γ-ray counting setup and the γ-rays emitted from the decay of 238 Np were measured over a time period of approximately 7 days. Multiple γ-ray decay branches of 238 Np were observed. The observed activity of 238 Np was used to calculate the amount of 238 Np produced during the irradiation via the 237 Np(n, γ) 238 Np reaction and an integral cross section of 342(11) mb was measured for the Godiva-IV neutron spectrum. Further, the 238 Np half-life has been measured with a result of 50.31(5) hours. The 237 Np(n, γ) 238 Np integral cross section measured in this work is in agreement with overlapping 1σ error bands to ENDF/B-VIII.0. However, the measured value is 3σ away from the calculated integral cross section using JENDL-5. This measurement offers a reliable benchmark for future 237 Np(n, γ) 238 Np cross section evaluations.
Nuclear-pumped lasing of Ar-Xe and He-Xe has been demonstrated using (U-235)F6 fission-fragment excitation. Fission fragments were created by absorption of thermal neutrons in a combination of gaseous (U-235)F6 and laser-tube wall coatings formed from UF6 chemical reaction products. At a pressure of 600 torr Ar-(3%)Xe, lasing occurred at 2.65 microns in Xe. Up to 3 torr of gaseous (U-235)F6 was added to 600 torr Ar-Xe before serious laser quenching occurred. With 3 torr of (U-235)F6 added, 38% of the energy deposition came from gaseous UF6 and the remainder from the uranium wall coating. The neutron flux at lasing threshold was found to be 4 x 10 to the 15th n/sq cm sec.
This paper describes the experimental results for an energy tuning assembly created to modify the National Ignition Facility deuterium–tritium fusion neutron source into a notional thermonuclear and prompt fission neutron spectrum, which has applications in integral measurements, nuclear data benchmarks, and radiation effects on microelectronics. The Monte Carlo neutron transport utilized MCNP5 to estimate the ETA-modified fluence using the ENDF-B/VIII.0 and IRDFF-II continuous energy nuclear data libraries, and SCALE Sampler was used to estimate the systematic nuclear data covariance using ENDF-B/VII.1 and IRDFF-II in a 252-group structure. The experiment fielded eight activation foils and a highly enriched uranium sample. This provided fifteen reaction channels that were used in a forward-fit comparison to the modeled results and to unfold the neutron spectrum using STAYSL. Gamma-ray spectrometry was performed on the activation and highly enriched uranium foils, and the reduced Χ 2 between the modeled and experimental values was 1.21. The results from the STAYSL unfold, reduced Χ 2 = 1.62, indicated that the modeled neutron spectrum was achieved and the systematic nuclear data uncertainty associated with the neutron transport and activation product cross sections was representative of the experiment. Integral cumulative fission product yield data were collected for 37 mass chains with a combination of gamma-ray spectrometry and radiochemical analysis. Fission product analysis was generally in agreement with two models using a semi-empirical fit and the General Observables of Fission code, with the exception of mass chains 88, 109, 111, 112, 113, 129, 139, 142, 144, 151, and 156.
Abstract Neutron-capture cross sections of neutron-rich nuclei are calculated using a Hauser–Feshbach model when direct experimental cross sections cannot be obtained. A number of codes to perform these calculations exist, and each makes different assumptions about the underlying nuclear physics. We investigated the systematic uncertainty associated with the choice of Hauser-Feshbach code used to calculate the neutron-capture cross section of a short-lived nucleus. The neutron-capture cross section for $$^{73}\hbox {Zn}$$ 73 Zn (n, $$\gamma $$ γ ) $$^{74}\hbox {Zn}$$ 74 Zn was calculated using three Hauser-Feshbach statistical model codes: TALYS, CoH, and EMPIRE. The calculation was first performed without any changes to the default settings in each code. Then an experimentally obtained nuclear level density (NLD) and $$\gamma $$ γ -ray strength function ( $$\gamma \hbox {SF}$$ γ SF ) were included. Finally, the nuclear structure information was made consistent across the codes. The neutron-capture cross sections obtained from the three codes are in good agreement after including the experimentally obtained NLD and $$\gamma \hbox {SF}$$ γ SF , accounting for differences in the underlying nuclear reaction models, and enforcing consistent approximations for unknown nuclear data. It is possible to use consistent inputs and nuclear physics to reduce the differences in the calculated neutron-capture cross section from different Hauser-Feshbach codes. However, ensuring the treatment of the input of experimental data and other nuclear physics are similar across multiple codes requires a careful investigation. For this reason, more complete documentation of the inputs and physics chosen is important.