Application of thermochemical modeling to assessment/evaluation of nuclear fuel behavior
Not Available
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.
Not Available
A large fraction of the rare-earth elements observed in the solar system is produced in the astrophysical rapid neutron capture process (r-process). However, current stellar models cannot completely explain the relative abundance of these elements partially because of nuclear physics uncertainties. To address this problem, a β-decay spectroscopy experiment was performed at RI Beam Factory (RIBF) at RIKEN, aimed at studying a wide range of very neutron-rich nuclei with Z∼60 that are progenitors of the rare-earth elements with mass number A∼460. The experiment provides a test of nuclear models as well as experimental inputs for r-process calculations. This contribution presents the experimental setup and some preliminary results of the experiment.
We present results of modeled cross sections for neutron- and proton-induced reactions leading to the final product nucleus 58 Co. In each case the gamma-cascade branching ratios given in the ENSDF database circa 2014 predict modeled nuclear cross sections leading to the ground and first excited metastable state that are incompatible with measured cross sections found in the NNDC experimental cross section database EXFOR. We show that exploring the uncertainty in the mixing ratio used to calculate the gamma-cascade branching ratios for the 53.15 keV 2 nd excited state leads to changes in the predicted partial cross sections by amounts that give good agreement with measured data.
The optical model of composite particle scattering is considered and compared to the accuracies of other models. A nonrelativistic Schroedinger equation with two-body potentials is used for the scattering of a single particle by an energy-dependent local potential. The potential for the elastic channel is composed of matrix elements of a single scattering operator taken between the ground states of the projectile and the target; the coherent amplitude is considered as dominating the scattering in the forward direction. A multiple scattering series is analytically explored and formally summed by the solution of an equivalent Schroedinger equation. Cross sections of nuclear scattering are then determined for He-4 and C-12 nuclei at 3.6 GeV/nucleus and O-16 projectiles at 2.1 GeV/nucleus, and the optical model approximations are found to be consistently lower and more accurate than approximations made by use of Glauber's theory.
The general-purpose Monte Carlo radiation-transport code MCNP6{sup TM} came from the integration of MCNP5{sup TM} and MCNPX{sup TM}. The most recent release of MCNP, version 6.2.0, is scheduled for release in early 2017, which includes a variety of bug fixes and new features to the code. MCNP traditionally models nuclear reactions using Monte Carlo sampling techniques on measured and modeled cross-sectional data contained within ACE (A Compact Evaluated Nuclear Data File) libraries. In addition to reaction cross-sections, ACE libraries contain information relating to the production of secondary particles, such as secondary neutrons and photons. For example, an incident neutron on a nucleus can produce a nuclear reaction in the form (n, M{sub n}n'), where n is the incident neutron and M{sub n} is the number of secondary neutrons n' produced. Using this notation, different reactions can be described based on their neutron multiplicity, such as neutron capture (M{sub n} = 0), elastic and inelastic scattering (M{sub n} = 1), and reactions with multiple secondary neutrons (M{sub n} > 1). Accounting for secondary γ-rays emitted from the residual nuclei, this reaction can further be generalized as n, Mnn'M to describe a reaction producing M{sub n} neutrons and M -rays. These values M{sub n} and M are referred to as the neutron and multiplicities, respectively. ACE libraries contain the cross-sectional data for each reaction, and their corresponding statistically averaged multiplicities for neutrons and photons, M-bar{sub n} and M-bar{sub γ}. However, there are several limitations to reaction sampling using ACE libraries. First, the data contained within ACE libraries is limited by the accuracy of the model and/or experiment, which out of necessity can greatly simplify the true physics behind a reaction. For example, the value M{sub γ} in a (n, M{sub n}n'M{sub γ}γ) reaction is determined by the statistically averaged multiplicity M-bar{sub γ} instead of a distribution. Consequently, MCNP, unable to produce a continuous distribution, will utilizes a binary sampling around M-bar{sub γ}. While this sampling method provides a statistically average M{sub γ} over a large simulation, it is an inaccurate representation of an individual reaction and its corresponding secondary particles. Secondary neutrons are affected in a similar manner, such that different kinds of secondary neutrons can be sampled from the same interaction (eg., an inelastic neutron with an [n, 2n] neutron). Lastly, the random sampling technique for secondary particles remains completely independent of the sampled reaction, eliminating the ability to correlate secondary particles for a given reaction. To supplement these data libraries, additional codes and models describing certain nuclear reactions are employed. These codes can fill in the information missing from ACE, provide a better model to the true physics distribution, and can provide true correlated secondary particles, which has many useful applications in research and industry.
The ground state half-lives of 69 Ge, 73 Se, 83 Sr, 63 Zn, and the half-life of the 1/2 - isomer in 85 Sr have been measured with high precision using the photoactivation technique at an unconventional bremsstrahlung facility that features a repurposed medical electron linear accelerator. The γ -ray activity was counted over about 6 half-lives with a high-purity germanium detector, enclosed into an ultra low-background lead shield. The measured half-lives are: T 1/2 ( 69 Ge) = 38.82 ± 0.07 (stat) ± 0.06 (sys) h; T 1/2 ( 73 Se) = 7.18 ± 0.02 (stat) ± 0.004 (sys) h; T 1/2 ( 83 Sr) = 31.87 ± 1.16 (stat) ± 0.42 (sys) h; T 1/2 ( 85m Sr) = 68.24 ± 0.84 (stat) ± 0.11 (sys) min; T 1/2 ( 63 Zn) = 38.71 ± 0.25 (stat) ± 0.10 (sys) min. These high-precision half-life measurements will contribute to a more accurate determination of corresponding ground-state photoneutron reaction rates, which are part of a broader effort of constraining statistical nuclear models needed to calculate stellar nuclear reaction rates relevant for the astrophysical p -process nucleosynthesis.
The precise measurement of neutrino properties is among the highest priorities in fundamental particle physics, involving many experiments worldwide. Since the experiments rely on the interactions of neutrinos with bound nucleons inside atomic nuclei, the planned advances in the scope and precision of these experiments require a commensurate effort in the understanding and modeling of the hadronic and nuclear physics of these interactions, which is incorporated as a nuclear model in neutrino event generators. Lastly, this model is essential to every phase of experimental analyses and its theoretical uncertainties play an important role in interpreting every result.
Fission neutron multiplicity data are imperative for nuclear models that support criticality safety, nuclear nonproliferation, and stockpile stewardship. This work presents the first experimental measurement to determine the average neutron multiplicity ($\bar{ν}$) and the neutron multiplicity distribution using a surrogate reaction. 241 Pu has a short half-life of 14.3 yr. To circumvent experimental difficulties, inelastic scattering of 53.96-MeV α particles on 242 Pu served as a surrogate for 241 Pu (n , f). The present text reports and discusses the average neutron multiplicity, in addition to the neutron multiplicity moments and distribution for equivalent neutron energies up to 20 MeV. Finally, the $\bar{ν}$ measured in this experiment agrees with previous measurements until 12 MeV, where discrepancies arise due to preequilibrium neutron emission.
Changes in the nuclear charge radii of lithium isotopes were determined using a combination of precise isotope shift measurements and theoretical atomic structure calculations. We discuss the choice of the reference isotope for absolute charge radii determinations in the lithium isotopic chain and report a new value for the charge radius of {sup 6}Li, based on the analysis of the world scattering data. A summary of the lithium nuclear charge radii obtained in this way is presented. Additionally, new calculations in fermionic molecular dynamics for the lithium isotopes were performed. We summarize the status of the lithium nuclear charge radii, magnetic dipole and electric quadrupole moments from experimental investigations and compare them to the results of various microscopic and three-body nuclear models.
Purpose: Monte Carlo (MC) codes are increasingly used in the ion beam therapy community due to their detailed description of radiation transport and interaction with matter. The suitability of a MC code demands accurate and reliable physical models for the transport and the interaction of all components of the mixed radiation field. This contribution will address an overview of the recent developments in the FLUKA code oriented to its application in ion beam therapy. Methods: FLUKA is a general purpose MC code which allows the calculations of particle transport and interactions with matter, covering an extended range of applications. The user can manage the code through a graphic interface (FLAIR) developed using the Python programming language. Results: This contribution will present recent refinements in the description of the ionization processes and comparisons between FLUKA results and experimental data of ion beam therapy facilities. Moreover, several validations of the largely improved FLUKA nuclear models for imaging application to treatment monitoring will be shown. The complex calculation of prompt gamma ray emission compares favorably with experimental data and can be considered adequate for the intended applications. New features in the modeling of proton induced nuclear interactions also provide reliable cross section predictions for the production of radionuclides. Of great interest for the community are the developments introduced in FLAIR. The most recent efforts concern the capability of importing computed-tomography images in order to build automatically patient geometries and the implementation of different types of existing positron-emission-tomography scanner devices for imaging applications. Conclusion: The FLUA code has been already chosen as reference MC code in many ion beam therapy centers, and is being continuously improved in order to match the needs of ion beam therapy applications. Parts of this work have been supported by the European FP7 project ENVISION (grant agreement no. 241851)
Developing automated methods for data collection and analysis that can facilitate nuclear nonproliferation assessment is an important research area with significant consequences for the effective global deployment of nuclear energy. Facility modeling that can integrate and interpret observations collected from monitored facilities in order to ascertain their functional details will be a critical element of these methods. Although improvements are continually sought, existing facility modeling tools can characterize all aspects of reactor operations and the majority of nuclear fuel cycle processing steps, and include algorithms for data processing and interpretation. Assessing nonproliferation status is challenging because observations can come from many sources, including local and remote sensors that monitor facility operations, as well as open sources that provide specific business information about the monitored facilities, and can be of many different types. Although many current facility models are capable of analyzing large amounts of information, they have not been integrated in an analyst-friendly manner. This paper addresses some of these facility modeling capabilities and illustrates how they could be integrated and utilized for nonproliferation analysis. The inverse problem of inferring facility conditions based on collected observations is described, along with a proposed architecture and computer framework for utilizing facility modeling tools. After considering a representative sampling of key facility modeling capabilities, the proposed integration framework is illustrated with several examples.
The Sampler sequence was introduced into the SCALE nuclear modeling and simulation suite in SCALE 6.2 to perform uncertainty quantification via random sampling of nuclear data, material number densities, and dimensions. Sampler was expanded with the introduction of a parametric capability in SCALE 6.2.2. This paper discusses input for the Sampler parametric sequence and presents two case studies of analyses performed using the sequence. These case studies include preconceptual design of a package for transporting high assay low-enriched uranium (HALEU) oxide and scoping calculations to support subcritical limit development for a future update of the ANSI/ANS-8.1 (ANS-8.1) standard. The parametric capability within Sampler provides many benefits to analysts. For instance, parametric sweeps are frequently used to identify optimum parameter values as part of safety analysis or system design, but such sweeps can require substantial engineering time or may rely on custom-written scripts or scripts such as Write One, Run Many (or WORM) developed outside of any software quality assurance program. With the parametric capabilities in Sampler, however, a large number of inputs can be generated automatically without recourse to scripting by individual analysts. The parametric capability can also be used in lieu of the CSAS5S search sequence to identify optimum parameters more simply with straightforward inputs and outputs. Sampler can also be used to calculate input parameters from engineering specifications. For example, diameters can be converted to radii, or masses can be used to calculate number densities. Overall, the Sampler parametric capability provides a robust feature within SCALE, eliminating the need for user-developed scripting.
Deposition models were built to evaluate contaminant deposition on spent nuclear fuel (SNF) canisters. The primary contaminant of concern is chloride, which is dispersed in the atmosphere and then deposits onto the canisters. During dry storage, the primary degradation process is likely to be Chloride Induced Stress Corrosion Cracking (CISCC) at the heat-affected zones of the canister welds. It is known that stainless steel canisters are susceptible to CISCC; however, the rate of chloride deposition onto the canisters is poorly known, based on sparse field data from a small number of sites. The models presented in this study could be useful for determining the rate of deposition on the canisters and the likelihood of CISCC to help with SNF canister ageing management. The deposition models were developed with the commercial computational fluid dynamics (CFD) code STAR-CCM+. Various deposition mechanisms were considered and incorporated into the models, and a sensitivity study was conducted to determine the most important mechanisms for deposition within a SNF storage system. The models included both a vertical and horizontal configuration storage system: NAC International’s Modular, Advanced Generation, Nuclear All-purpose STORage System (MAGNASTOR®) and a NUHOMS® horizontal storage module respectively. The resulting canister deposition on the horizontal canister is visually compared with inspection data taken onsite at the Calvert Cliffs Nuclear Power Plant. These models are preliminary, and development of the models will continue. Future validation exercises are currently being planned, including the Canister Deposition Field Demonstration (CDFD) effort funded by U.S. Department of Energy office of Nuclear Energy office of Spent Fuel Waste Science and Technology. The goal of the modeling presented is to demonstrate a potential modeling technique that could be used to plan and inform SNF canister ageing management programs with predictive models for the timing and occurrence of canister CISCC.
Uncertainties in nuclear models have a major impact on simulations that aim at understanding the origin of heavy elements in the universe through the rapid neutron capture process ( r process) of nucleosynthesis. Within the framework of the nuclear density functional theory, we use results of Bayesian statistical analysis to propagate uncertainties in the parameters of energy density functionals to the predicted r -process abundance pattern, by way not only of the nuclear masses but also through the influence of the masses on β -decay and neutron capture rates. We point out the importance of the nonequilibrium end stage of the r process in determining the width of the resulting abundance pattern uncertainty bands. We additionally make the first identifications of specific parameters of Skyrme-like energy density functionals which show tentative correlations with particular aspects of the r -process abundance pattern. While previous studies have explored the reduction in the abundance pattern uncertainties due to anticipated new measurements of neutron-rich nuclei, here we point out that an even larger reduction will occur when these new measurements are used to reduce the uncertainty of model predictions of masses, which are then propagated through to the abundance pattern. We make a quantitative prediction for how large this reduction will be.
We perform atomic relativistic many-body calculations of the field isotope shifts and calculations of corresponding nuclear parameters for all stable even-even isotopes of the Yb + ion. We demonstrate that if we take nuclear parameters of the Yb isotopes from a range of the state of the art nuclear models, which all predict strong quadrupole nuclear deformation, and then calculate nonlinearity of the King plot caused by the difference in the deformation in different isotopes, the result is consistent with the nonlinearity observed in the experiment. In conclusion, the changes of nuclear rms radius between isotopes extracted from experiment are consistent with those obtained in the nuclear calculations.
Analytical studies were continued to identify the design and performance characteristics of a small-scale model of a nuclear light bulb unit cell suitable for testing in a nuclear furnace reactor. Emphasis was placed on calculating performance characteristics based on detailed radiant heat transfer analyses, on designing the test assembly for ease of insertion, connection, and withdrawal at the reactor test cell, and on determining instrumentation and test effluent handling requirements. In addition, a review of candidate test reactors for future nuclear light bulb in-reactor tests was conducted.
The International Space Station (ISS) provides the proving ground for future long duration human activities in space. Ionizing radiation measurements in ISS form the ideal tool for the experimental validation of ionizing radiation environmental models, nuclear transport code algorithms, and nuclear reaction cross sections. Indeed, prior measurements on the Space Transportation System (STS; Shuttle) have provided vital information impacting both the environmental models and the nuclear transport code development by requiring dynamic models of the Low Earth Orbit (LEO) environment. Previous studies using Computer Aided Design (CAD) models of the evolving ISS configurations with Thermo Luminescent Detector (TLD) area monitors, demonstrated that computational dosimetry requires environmental models with accurate non-isotropic as well as dynamic behavior, detailed information on rack loading, and an accurate 6 degree of freedom (DOF) description of ISS trajectory and orientation.