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.
Integral test of inelastic scattering cross sections using measured neutron spectra from thick shells of Ta, W, Mo, and Be
Analysis of interplanetary cargo and piloted missions requires calculations of the performances and masses of subsystems to be integrated in a final design. In a preliminary and scoping stage the designer needs to evaluate options iteratively by using fast computer simulations. The Oak Ridge National Laboratory (ORNL) has been involved in the development of models and calculational procedures for the analysis (neutronic and thermal hydraulic) of power sources for nuclear electric propulsion. The nuclear modules will be integrated into the whole simulation of the nuclear electric propulsion system. The vehicles use either a Brayton direct-conversion cycle, using the heated helium from a NERVA-type reactor, or a potassium Rankine cycle, with the working fluid heated on the secondary side of a heat exchanger and lithium on the primary side coming from a fast reactor. Given a set of input conditions, the codes calculate composition. dimensions, volumes, and masses of the core, reflector, control system, pressure vessel, neutron and gamma shields, as well as the thermal hydraulic conditions of the coolant, clad and fuel. Input conditions are power, core life, pressure and temperature of the coolant at the inlet of the core, either the temperature of the coolant at the outlet of the core or the coolant mass flow and the fluences and integrated doses at the cargo area. Using state-of-the-art neutron cross sections and transport codes, a database was created for the neutronic performance of both reactor designs. The free parameters of the models are the moderator/fuel mass ratio for the NERVA reactor and the enrichment and the pitch of the lattice for the fast reactor. Reactivity and energy balance equations are simultaneously solved to find the reactor design. Thermalhydraulic conditions are calculated by solving the one-dimensional versions of the equations of conservation of mass, energy, and momentum with compressible flow.
This presentation gives an overview of thermal cross section measurements at the RPI LINAC. Work discussed includes novel cold moderator capability for use in extending neutron flux down to 0.0005 eV and total thermal neutron cross section measurements performed for polyethylene, polystyrene, Plexiglas (Lucite), and yttrium hydride over the entire thermal region (0.0005 – 1 eV).
Background: Modeling nuclear reaction networks for nuclear science applications and for simulations of astrophysical environments relies on cross section data for a vast number of reactions, many of which have never been measured. Cross sections for neutron-induced reactions on unstable nuclei are particularly scarce, since they are the most difficult to measure. Consequently, we must rely on theoretical predictions or indirect measurements to obtain the requisite reaction data. For compound nuclear reactions, the surrogate reaction method can be used to determine many cross sections of interest. Purpose: Earlier work has demonstrated that cross sections for neutron-induced fission and radiative neutron capture can be determined from a combination of surrogate reaction data and theory. For the fission case, it was shown that the Weisskopf-Ewing approximation, which significantly simplifies the implementation of the surrogate method, can be employed. Capture cross sections cannot be obtained, and require a detailed description of the surrogate reaction process. Here, we examine the validity of the Weisskopf-Ewing approximation for determining unknown (n, n') and (n, 2n) cross sections from surrogate data. Methods: Using statistical reaction calculations with realistic parametrizations, we investigate first whether the assumptions underlying the Weisskopf-Ewing approximation are valid for (n, n') and (n, 2n) reactions on representative target nuclei. We then produce simulated surrogate reaction data and assess the impact of applying the Weisskopf-Ewing approximation when extracting (n, n') and (n, 2n) cross sections in situations where the approximation is not strictly justified. Results: We find that peak cross sections can be estimated using the Weisskopf-Ewing approximation, but the shape of the (n, n') and (n, 2n) cross sections, especially for low neutron energies, cannot be reliably determined without accounting for the angular-momentum differences between the neutron-induced and surrogate reaction. Conclusions: To obtain reliable (n, n') and (n, 2n) cross sections from surrogate reaction data, a detailed description of the surrogate reaction mechanisms is required. To do so for the compound-nucleus energies and decay channels relevant to these reactions, it becomes necessary to extend current modeling capabilities.
Effective interactions for nucleon-nucleus elastic scattering from first principles require the use of the same nucleon-nucleon interaction in the structure and reaction calculations and a consistent treatment of the relevant operators at each order. Systematic investigations of the effect of truncation uncertainties of chiral nucleon-nucleon ( N N ) forces have been carried out for scattering observables in the two- and three-nucleons system as well as for bound-state properties of light nuclei. Here we extend this type of study to proton and neutron elastic scattering for 16 O and 12 C. Using the frameworks of the spectator expansion of multiple scattering theory as well as the no-core shell model, we employ one specific chiral interaction from the LENPIC collaboration and consistently calculate the leading-order effective nucleon-nucleus interaction up to the third chiral order, from which we extract elastic scattering observables. Then we apply pointwise as well as correlated uncertainty quantification for the estimation of the chiral truncation error. We calculate and analyze proton elastic scattering observables for 16 O and 12 C as well as neutron elastic scattering observables for 12 C between 65 and 185 MeV projectile kinetic energy. We find in all cases qualitatively similar results for the chiral truncation uncertainties as in few-body systems and assess them using similar diagnostic tools. The order-by-order convergence of the scattering observables for 16 O and 12 C is very reasonable around 100 MeV, while for higher energies the chiral expansion parameter becomes too large for convergence. Comparing proton and neutron differential cross sections reveals that their description is comparable up to around 100 MeV. Here, we also find a nearly perfect correlation between the differential cross section for neutron scattering and the N N Wolfenstein amplitudes for small momentum transfers. The diagnostic tools for studying order-by-order convergence in observables in few-body systems can be employed for observables in nucleon-nucleus scattering with only minor modifications provided the momentum scale in the problems is not too large. We also find that the chiral N N interaction on which our study is based gives a very good description of differential cross sections for 16 O and 12 C as low as 65-MeV projectile energy, particularly in the forward direction. In addition, the very forward direction of the neutron differential cross section mirrors the behavior of the N N interaction amazingly well.
Neutron cross section matrices for fission and scattering data are required for each material, temperature, and enrichment level to calculate the neutron transport equation accurately. Here, this information can be a limiting factor when using the multigroup discrete ordinates (S N ) method when the number of energy groups is large. Machine Learning (ML) can be used to replace the need for the cross section matrices by reproducing the function that maps the scalar flux to the scattering and fission sources. Through the use of autoencoders and Deep Jointly-Informed Neural Networks (DJINN), the data storage requirements are reduced by 94% of the original data for a 618 group problem. This is accomplished while preserving the scalar flux, maintaining generality, and decreasing wall clock times.
The development of nuclear reaction models for the production of evaluated nuclear data has traditionally been performed by comparing measured cross sections with predictions from reaction model codes whose physical input parameters are adjusted to obtain the best agreement between measured and modeled results. To more directly probe reaction model inputs, this work introduces a forward modeling approach to experimental reaction cross-section determination, where the most important physical input parameters to reaction model calculations are obtained via 𝜒 2 minimization between measured and calculated observables. This was demonstrated using data collected by the Gamma Energy Neutron Energy Spectrometer for Inelastic Scattering (GENESIS) at the 88-inch cyclotron at Lawrence Berkeley National Laboratory, a detection array consisting of organic liquid scintillators and high-purity germanium (HPGe) detectors. Using a broad-spectrum neutron beam and a 99.98%-enriched 56 Fe target, GENESIS was used to perform a simultaneous measurement of 56 Fe 𝛾-ray production cross sections and secondary neutron energy and angle distributions. The results of the forward modeling approach to the determination of energy-differential 𝛾-ray production cross sections for the yrast 4 + → 2 + and 6 + → 4 + transitions, as well as eight other off-yrast transitions, were compared against those obtained using conventional techniques, and the results are in good agreement. In addition to discrete 𝛾-ray yield total scattered neutron energy-angular distributions as a function of incident neutron energy were also obtained using forward modeling and found to agree with evaluated data, with the exception of elastic scattering at small angles. The fitted reaction model parameters obtained through forward modeling were also used to calculate the cross section for the unobserved (𝑛, 2𝑛) reaction; excellent agreement with the current evaluation was obtained, providing a validation of the predictive capabilities of the forward model approach. This work bridges the gap between nuclear data experiment and evaluation by providing a new means for extracting inelastic neutron-scattering cross sections and neutron-induced 𝛾-ray production data while directly probing reaction model physics.
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.
One nuclear electric propulsion (NEP) reactor systems under consideration is a hydride moderated thermal spectrum reactor fueled by high assay low enriched uranium (HALEU). While such a reactor is expected to yield the lightest HALEU reactor design, its development challenges grow exponentially with increasing mission demands, most notably power output, specific weight of the overall system (which may require operation at temperatures exceeding 1200 K), service lifetime, and human-rated reliability. Two of the greatest cost drivers are full-powered nuclear demonstrations and extensive material development campaigns, so it is important to consider options that can minimize the need for or complexity of such tasks. This paper discusses a structured framework being developed for assessing how NEP design choices, such as materials selection, neutronic features, and heat-removal technologies, can translate into project risk and how project performance goals can be traded with development cost.Reactors operating at high temperatures often require cutting-edge heat transfer technologies and creep-resistant materials. Use of new materials in high temperature reactors brings additional complication beyond those common to any new space materials development campaign. For example, such materials may not possess necessary neutronic cross-sectional or neutronic irradiation data. Similarly, use of new materials may significantly influence core neutronics; in some extreme cases, neutronic reactivity feed-back of certain new materials can vary during their service life as radiation damage impacts the scattering cross-section. In an affordable development approach, high fidelity modeling and simulation tools are used to identify and characterize potential ‘knees-in-the-curves’ in the relationship that exists between the mission characteristics and the project risk. Of particular significance is use of modern uncertainty management and variance reduction methods to perform gap analyses that feed into phenomena identification and ranking tables (PIRT) commonly used to communicate nuclear readiness levels. Model-based measurements techniques are used to design sub-scale experiments as a substitute to minimize orcompletely eliminate the need for nuclear demonstrations.This paper will describe the approach and present preliminary results. It will lay the groundwork for developing a set of metrics that can be broadly characterized as system nuclear readiness levels and advancement degree of difficulty for nuclear systems. Equally importantly, a goal of this paper is to initiate a dialogue among stakeholders on what is the sufficient level of maturity that is required for launching a demonstration unit.
The program reviewed is a study of the feasibility, design, and optimization of the GCATR. The program is designed to take advantage of initial results and to continue work carried out on the Gas Core Breeder Reactor. The program complements NASA's program of developing UF6 fueled cavity reactors for power, nuclear pumped lasers, and other advanced technology applications. The program comprises: (1) General Studies--Parametric survey calculations performed to examine the effects of reactor spectrum and flux level on the actinide transmutation for GCATR conditions. The sensitivity of the results to neutron cross sections are to be assessed. Specifically, the parametric calculations of the actinide transmutation are to include the mass, isotope composition, fission and capture rates, reactivity effects, and neutron activity of recycled actinides. (2) GCATR Design Studies--This task is a major thrust of the proposed research program. Several subtasks are considered: optimization criteria studies of the blanket and fuel reprocessing, the actinide insertion and recirculation system, and the system integration. A brief review of the background of the GCATR and ongoing research is presented.
Accurate neutron cross section data are a vital input to the simulation of nuclear systems for a wide range of applications from energy production to national security. The evaluation of experimental data is a key step in producing accurate cross sections. There is a widely recognized lack of reproducibility in the evaluation process due to its artisanal nature and therefore there is a call for improvement within the nuclear data community. This can be realized by automating/standardizing viable parts of the process, namely, parameter estimation by fitting theoretical models to experimental data. This automation effort could greatly benefit from a synthetic data resource. This work leverages problem-specific physics, Monte Carlo sampling, and a general methodology for data synthesis to generate unlimited, labelled experimental cross-section data that is statistically indistinguishable to the observed data. Heuristic and, where applicable, rigorous statistical comparisons to observed data support this claim. The demonstration is based on/limited to transmission measurements at Rensselaer Polytechnic Institute (RPI) and energy-differential cross sections in the resolved resonance region (RRR). An open-source software is published alongside this article that executes the complete methodology to produce high-utility synthetic datasets. The goal of this work is to provide an approach and corresponding tool that will allow the evaluation community to begin exploring more data-driven, ML-based solutions to long-standing challenges in the field.
For the accurate assessment of the heat generation rate in fast reactors, the gamma library of MC2 -3 and the MC2 -3 + GAMSOR procedure employing the coupled neutron and gamma heating calculation has been thoroughly verified against Monte Carlo results and validated using the ZPPR-15D gamma dose measurement data. NJOY outputs are post-processed in a consistent way with the NJOY procedure to avoid any missing data or double counting of data. Prompt heating for 379 out of 391 isotopes in the gamma library was verified against MCNP6.2 to within 1% relative error in total heating for most isotopes. For both a simple one-dimensional slab problem representing a sodium cooled fast reactor and the Experimental Breeder Reactor II (EBR-II) benchmark problem, the root-mean-square values of assembly power error were less than 0.5% for fuel assemblies, ~1% in blankets and ~1 to ~3% in reflectors compared to MCNP6.2 results. The most plausible cause for the 3% error in a reflector assembly is believed to be the error in the multigroup neutron cross sections for the reflector assembly. For validation, gamma doses measured with thermoluminiscent dosimeters (TLDs) in the ZPPR-15D experiment were calculated using GAMSOR. Due to the uncertainty in the TLD measurement with regards to the energy deposition of photons and neutrons, the validation data leads to a 12.7% uncertainty on the experimental measurement. With this uncertainty bound, the calculated doses all fell within one standard deviation of the measured value. Combined with the accurate calculation of reaction rate distributions and neutron spectrum measurements, these results indicate good agreement for the neutron and gamma heating calculations that were performed.
Radioactivity induced by high-energy protons and secondary neutrons (from nuclear interactions) in various samples returned from different locations in Skylab was measured directly by gamma-ray spectroscopy measurements of decay gamma rays from the samples. Incident fluxes were derived from the activation measurements, using known nuclear cross-section. Neutron and proton flux values were found to range from 0.2 to 5 particles/sq cm-sec, depending on the energy range and location in Skylab. The thermal neutron flux was less than 0.07 neutrons/sq cm-sec. The results are useful for data analysis and planning of future high-energy astronomy experiments.
The space radiation transport code, HZETRN, has been used extensively for research, vehicle design optimization, risk analysis, and related applications. One of the simplifying features of the HZETRN transport formalism is the straight-ahead approximation, wherein all particles are assumed to travel along a common axis. This reduces the governing equation to one spatial dimension allowing enormous simplification and highly efficient computational procedures to be implemented. Despite the physical simplifications, the HZETRN code is widely used for space applications and has been found to agree well with fully 3D Monte Carlo simulations in many circumstances. Recent work has focused on the development of 3D transport corrections for neutrons and light ions (Z < 2) for which the straight-ahead approximation is known to be less accurate. Within the development of 3D corrections, well-defined convergence criteria have been considered, allowing approximation errors at each stage in model development to be quantified. The present level of development assumes the neutron cross sections have an isotropic component treated within N explicit angular directions and a forward component represented by the straight-ahead approximation. The N = 1 solution refers to the straight-ahead treatment, while N = 2 represents the bi-directional model in current use for engineering design. The figure below shows neutrons, protons, and alphas for various values of N at locations in an aluminum sphere exposed to a solar particle event (SPE) spectrum. The neutron fluence converges quickly in simple geometry with N > 14 directions. The improved code, 3DHZETRN, transports neutrons, light ions, and heavy ions under space-like boundary conditions through general geometry while maintaining a high degree of computational efficiency. A brief overview of the 3D transport formalism for neutrons and light ions is given, and extensive benchmarking results with the Monte Carlo codes Geant4, FLUKA, and PHITS are provided for a variety of boundary conditions and geometries. Improvements provided by the 3D corrections are made clear in the comparisons. Developments needed to connect 3DHZETRN to vehicle design and optimization studies will be discussed. Future theoretical development will relax the forward plus isotropic interaction assumption to more general angular dependence.
Second order optical model solutions to the elastic scattering amplitude were used to evaluate total, absorption, and abrasion cross sections for neutron nucleus scattering. Improved agreement with experimental data for total and absorption cross sections is found when compared with first order (coherent approximation) solutions, especially below several hundred MeV. At higher energies, the first and second order solutions are similar. There are also large differences in abrasion cross section calculations; these differences indicate a crucial role for cluster knockout in the abrasion step.
The Critical Unresolved Region Integral Experiment (CURIE) critical experiment was performed at the National Criticality Experiments Research Center (NCERC) at the Device Assembly Facility (DAF) at the Nevada Nuclear Security Sites (NNSS) in 2020. The objective of CURIE was to improve the quality of integral nuclear data in the uranium-235 ( 235 U) unresolved resonance region (URR) by performing benchmark integral experiments that were sensitive to the URR energy ranges. The CURIE experiment was evaluated for the International Criticality Safety Benchmark Evaluation Project (ICSBEP) handbook and the benchmark evaluation was accepted in 2022. The URR is a region within the intermediate neutron energy range (the intermediate energy ranges from 0.7 eV to 100 keV). The observed resonance structure in neutron cross sections is due to discrete energy levels in the nucleus and are characterized by resonance parameters. In the URR region, the resonance parameters are only partially resolved as the resolution of the experiment techniques becomes comparable to the average width of the resonances themselves, and the resonances are so close to one another that the structure cannot be determined empirically. The ENDF/B-VIII.0 and JEFF-3.3 nuclear data libraries define the URR as beginning at 2.25 keV and continuing until 25 keV. There are minimal intermediate neutron energy benchmarks available in the ICSBEP benchmark handbook, and aside from CURIE there are none that are highly sensitive in the URR energy region. There is a current need for additional experiments sensitive to the URR energy region. There is a new proposed subgroup for the Organization for Economic Co-operation and Development and Nuclear Energy Agency (OECD NEA) Working Party on International Nuclear Data Evaluation Cooperation (WPEC), so any new models or information will need experimental validation and testing. The NEA Working Party on Nuclear Criticality Safety (WPNCS) recent experimental needs and priority list includes intermediate energy 235 U and 238 U experiments, as does the recent Integral Experiments to Address Nuclear Criticality Safety Needs Meeting (May 2023). A variant on the original CURIE experiment called the Multiple Critical Unresolved Region Integral Experiment (MCURIE) is proposed to provide further investigation of the intermediate and URR energy region for uranium. MCURIE will utilize existing fuels at NCERC but use alternative moderators and reflectors to modify the neutron absorption, scattering, and fission spectra of the experiments, allowing for precise targeting of nuclear data sensitivities. The overall goal of MCURIE is to develop a framework for designing and performing integral benchmark experiments with high sensitivities in the intermediate and URR energy regions for nuclear data validation.
Detection of special nuclear materials (SNMs) is of vital importance in the prevention of nuclear terrorism and to secure states’ national security. Neutron detection is a particularly useful tool to identify SNM, and neutron-sensitive scintillators have many promising properties, such as ease of use, good time resolution, and high detection efficiency. In this work, we develop highly stable, self-oriented, ultrafast 1D ZnO:Li (and codoped with Al, Ga, and In) nanorods (NRs) as thermal neutron-sensitive scintillators. Lithium-6 has high thermal neutron cross section for the (n, α) reaction in ZnO:Li scintillators which have a vertical nano array design greatly increasing the effective surface area and scintillation efficiency. Cost-effective low-temperature (95 °C) hydrothermal growth is used to obtain highly crystalline ZnO:Li nano scintillators by combining nuclear range data and electron transport mechanisms. Among the studies using low-temperature hydrothermal synthesis and a relatively low annealing temperature (≈350 °C) along with optimized NRs (length ≈ 5–8 μm, mean diameter ≈ 700 nm) for thermal neutron detection, this study reports the shortest scintillation decay time (≈ 470 ps) so far to the best of our knowledge. In conclusion, this nano array scintillator combines the advantages of a low-cost growth technique with environmentally friendly and widely available materials.
Neutron fluences were measured from 435 MeV/nucleon Nb ions stopping in a Nb target and 272 MeV/nucleon Nb ions stopping in targets of Nb and Al for neutrons above 20 MeV and at laboratory angles between 3 degrees and 80 degrees. The resultant spectra were integrated over angles to produce neutron energy distributions and over energy to produce neutron angular distributions. The total neutron yields for each system were obtained by integrating over the angular distributions. The angular distributions from all three systems are peaked forward, and the energy distributions from all three systems show an appreciable yield of neutrons with velocities greater than the beam velocity. Comparison of the total neutron yields from the two Nb + Nb systems suggests that the average neutron multiplicity decreases with decreasing projectile energy. Comparison of the total yields from the two 272 MeV/nucleon systems suggests that the total yields show the same dependence on projectile and target mass number as do total inclusive neutron cross sections. The data are compared with Boltzmann-Uehling-Uhlenbeck model calculations.