Search NASA⌕ Search

SEARCH · Search NASA

Results for “nuclear shell model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Coupled Calculations in Helical Steam Generator: Validation on Legacy Data

The NEAMS program aims to develop an integrated multi-physics simulation capability “pellet-to-plant” for the design and analysis of future generations of nuclear power plants. In particular, the Reactor Product Line code suite's multi-resolution hierarchy is being designed to ultimately span the full range of length and time scales present in relevant reactor design and safety analyses, as well as scale from desktop to petaflop computing platforms. Flow-induced vibration (FIV) is widespread problem in energy systems because they rely on fluid movement for energy conversion. Vibrating structures may be damaged as fatigue or wear occurs. Given the importance of reliable components in the nuclear industry, flow-induced vibration has long been a major concern in safety and operation of nuclear reactors. In particular, nuclear fuel rods and steam generators have been known to suffer from flow-induced vibration and related failures. Advanced reactors, such as integral Pressurized Water Reactors (PWRs) considered for Small Modular Reactors (SMR), often rely on innovative component designs to meet cost and safety targets. One component that is the subject of advanced designs is the steam generator, some designs of which forego the usual shell-and-tube architecture in order to fit within the primary vessel. In addition to being more cost- and space-efficient, such steam generators need to be more reliable, since failure of the primary vessel represents a potential loss of coolant and a safety concern. A significant amount of data exists on flow-induced vibration in shell-and-tube heat exchangers, and heuristic methods are available to predict their occurrence based on a set of given assumptions. In contrast, advanced designs have far less data available. Advanced modeling and simulation based on coupled structural and fluid simulations have the potential to predict flow-induced vibration in a variety of designs, reducing the need for expensive experimental programs, especially at the design stage. Over the past five years, the Reactor Product Line has developed the integrated multi-physics code suite SHARP. The goal of developing such a tool is to perform multi-physics neutronics, thermal/fluid, and structural mechanics modeling of the components inside the full reactor core or portions of it with a user-specified fidelity. In particular SHARP contains high-fidelity single-physics codes Diablo for structural mechanics and Nek5000 for fluid mechanics calculations. Both codes are state-of-the-art, highly scalable tools that have been extensively validated. These tools form a strong basis on which to build a flow-induced vibration modeling capability. In this report we discuss one-way coupled calculations performed with Nek5000 and Diablo aimed at simulating available FIV experiments in helical steam generators in the turbulent buffeting regime. In this regime one-way coupling is judged sufficient because the pressure loads do not cause substantial displacements. It is also the most common source of vibration in helical steam generators at the low flows expected in integral PWRs. The legacy data is obtained from two datasets developed at Argonne and B&W.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Global properties of nuclei at finite-temperature within the covariant energy density functional theory

In stellar environments nuclei appear at finite temperatures, becoming extremely hot in core-collapse supernovae and neutron-star mergers. However, due to theoretical and computational complexity, most model calculations of nuclear properties are performed at zero temperature, while those existing at finite temperatures are limited only to selected regions of the nuclide chart. In this study we perform the global calculation of nuclear properties for even-even 8 ≤ Z ≤ 104 nuclei at temperatures in range 0 ≤ T ≤ 2 MeV. Calculations are based on the finite-temperature relativistic Hartree-Bogoliubov model supplemented by the Bonche-Levit-Vautherin vapor subtraction procedure. We find that near the neutron-drip line the continuum states have significant contribution already at moderate temperature T ≈ 1 MeV, thus emphasizing the necessity of the vapor subtraction procedure. Results include neutron emission lifetimes, quadrupole deformations, neutron-skin thickness, proton and neutron pairing gaps, entropy and excitation energy. Up to the temperature T ≈ 1 MeV, the nuclear landscape is influenced only moderately by the finite-temperature effects, mainly by reducing the pairing correlations. Here, as the temperature increases further, the effects on nuclear structures become pronounced, reducing both the deformations and the shell effects.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Two-phonon wobbling in 135 Pr

The second-phonon (n ω = 2) wobbling band has been established in the nucleus 135 Pr. Conclusive evidence for its wobbling nature comes from the ΔI = 1, E2 character of the transitions between the new band and the previously identified transverse wobbler band (n ω = 1) in this nucleus. Theoretical calculations in the framework of the quasiparticle triaxial rotor and triaxial projected shell models are found to be in good agreement with the experimental results.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

First-excited state g factor of Te 136 by the recoil in vacuum method

The g factor of the first 2 + state of radioactive 136 Te with two valence protons and two valence neutrons beyond double-magic 132 Sn has been measured by the recoil in vacuum (RIV) method. The lifetime of this state is an order of magnitude longer than the lifetimes of excited states recently measured by the RIV method in Sn and Te isotopes, requiring a new evaluation of the free-ion hyperfine interactions and methodology used to determine the g factor. In this paper, the calibration data are reported and the analysis procedures are described in detail. The resultant g factor has a similar magnitude to the g factors of other nuclei with an equal number of valence protons and neutrons in the major shell. However, an unexpected trend is found in the g factors of the N = 84 isotones, which decrease from 136 Te to 144 Nd. Finally, shell model calculations with interactions derived from the CD Bonn potential show good agreement with the g factors and E2 transition rates of 2 + states around 132 Sn, confirming earlier indications that 132 Sn is a good doubly magic core.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2018

Theory of Carbon Dioxide Absorption by Hanford Waste - 18208

Hanford nuclear waste has a very high pH, and Hanford tanks are ventilated with air. Air has CO{sub 2(g)} in it, which is an acid, so it has some affinity for the waste. The absorption of CO{sub 2} by the waste is important because CO{sub 2} absorption consumes hydroxide at the liquid-air interface, and the hydroxide is important for inhibiting tank corrosion. NaOH is periodically added to some tanks in order to mitigate the surface acidification. The NaOH added then becomes part of the waste and must be treated as waste, increasing the cost of the mission. The Hobbs Equation, developed at the Savannah River site, is used at Hanford to predict CO{sub 2} absorption from the tank ventilation system. In this study, the Hobbs Equation functional form is theoretically evaluated from thermodynamic calculations of CO{sub 2} absorption, and diffusion rates in the waste surface. This study theoretically justifies the empirical equation reported by Hobbs. Thermodynamic calculations performed using the Environmental Simulation Program (ESP) software on a typical double-shell tank waste composition indicate that essentially 100 % of CO{sub 2} that contacts the waste will react with the waste whenever the pH is above 10. In the 1990's, a substantial amount of data was collected on the concentration of CO{sub 2} in the off-gas of ventilated double-shell tanks at Hanford. This historic data was mined in order to determine the quantity of CO{sub 2} absorbed by the waste. The results of this data mining activity determined that tanks below approximately pH 10 (i.e., tank 241-AN-107) will show limited CO{sub 2} absorption, consistent with both the Hobbs model and the theory developed here. Two other tank with higher pH showed contrasting CO{sub 2} absorption. One showed substantial CO{sub 2} absorption, and the other showed minimal CO{sub 2} absorption despite the high waste pH in one tank, presumably because there was minimal contact between the waste surface and most of the air flow. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Analog-antianalog isospin mixing in 47 K 𝛽−decay

Here, we have measured the isospin mixing of the 𝐼 𝜋 =1/2 + , 𝐸 𝑥 = 2.599 MeV state in nearly doubly magic 47 Ca with the isobaric analog 1/2 + state of 47 K . Using the TRIUMF atom trap for 𝛽 decay, we have measured a nonzero asymmetry of the progeny 47 Ca with respect to the initial 47 K spin polarization, which together with the 𝛽 asymmetry implies a nonzero ratio of Fermi to Gamow-Teller matrix elements 𝑦 = 0.098 ± 0.037 for the 1/2 + → 1/2 + transition. Interpreting 𝑦 as mixing between this state and the isobaric analog state implies a Coulomb matrix element magnitude 101 ± 37 keV. This relatively large matrix element supports a model from the literature of analog-antianalog isospin mixing, which predicts large matrix elements in cases involving excess neutrons over protons occupying more than one major shell. The result supports pursuing a search for time-reversal odd, parity-even, isovector interactions using a correlation in 47 K 𝛽 decay.

beta decay↗
DOE OSTI2017

Unexpected distribution of ν 1 f 7 / 2 strength in Ca 49

Here, the calcium isotopes have emerged as a critical testing ground for new microscopically derived shell-model interactions, and a great deal of experimental and theoretical focus has been directed toward this region. We investigate the relative spectroscopic strengths associated with 1f 7/2 neutron hole states in 47,49 Ca following one-neutron knockout reactions from 48,50 Ca. The observed reduction of strength populating the 7/2 – 1 state in 49 Ca, as compared to 47 Ca, is inconsistent with shell-model calculations using both phenomenological interactions such as GXPF1, and interactions derived from microscopically based two- and three-nucleon forces. The result suggests a fragmentation of the l = 3 strength to higher-lying states as suggested by the microscopic calculations, but the observed magnitude of the reduction is not reproduced in any shell-model description.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Impact of areal-density asymmetries on the loss of confinement and ignition threshold in inertial confinement fusion capsules

In inertial confinement fusion implosion experiments, variations in the shell areal density reduce the shell’s inertia to confine the core pressure. Distorted capsules with large areal-density modulations decompress faster than uniform capsules in the disassembly phase. A simple 3-D analytic hot-spot model is derived to include the effects of low-mode areal-density modulations in the ignition criterion. The generalized 3-D ignition criterion for low modes is shown to depend on both the harmonic mean and the arithmetic mean of the areal density. Here, the “thin spots” in the shell are shown to dominate the loss of confinement as reflected by the harmonic mean definition of areal densities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Generator-coordinate reference states for spectra and 0 ν β β decay in the in-medium similarity renormalization group

We use a reference state based on symmetry-restored states from deformed mean-field or generator-coordinate-method (GCM) calculations in conjunction with the in-medium similarity- renormalization group (IMSRG) to compute spectra and matrix elements for neutrinoless double- beta (0$v$ββ) decay. Because the decay involves ground states from two nuclei, we use evolved operators from the IMSRG in one nucleus in a subsequent GCM calculation in the other. We benchmark the resulting IMSRG+GCM method against complete shell-model diagonalization for both the energies of low-lying states in 48 Ca and 48 Ti and the 0$v$ββ matrix element for the decay of 48 Ca, all in a single valence shell. Our approach produces better spectra than either the IMSRG with a spherical-mean-field reference or GCM calculations with unevolved operators. For the 0$v$ββ matrix element the improvement is slight.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Comprehensive Analysis of the Neutrino Process in Core-collapsing Supernovae

We investigate the neutrino flavor change effects due to neutrino self-interaction and shock wave propagation, as well as the matter effects on the neutrino process in core-collapsing supernovae (CCSNe). For the hydrodynamics, we use two models: a simple thermal bomb model and a specified hydrodynamics model for SN1987A. For the presupernova model, we take an updated model, adjusted to explain SN1987A, which employs recent developments in the (n, γ) reaction rates for nuclei near the stability line (A ~ 100). As for the neutrino luminosity, we adopt two different models: equivalent neutrino luminosity and nonequivalent luminosity models. The latter is taken from a synthetic analysis of CCSN simulation data, which quantitatively presented the results obtained by various neutrino transport models. Relevant neutrino-induced reaction rates are calculated using a shell model for light nuclei and a quasiparticle random phase approximation model for heavy nuclei. For each model, we present abundances of the light nuclei ( 7 Li, 7 Be, 11 B, and 11 C) and the heavy nuclei ( 92 Nb, 98 Tc, 138 La, and 180 Ta) produced by the neutrino process. The light nuclei abundances turn out to be sensitive to the Mikheyev–Smirnov–Wolfenstein (MSW) region around O-Ne-Mg layer while the heavy nuclei are mainly produced prior to the MSW region. Through detailed analyses, we find that neutrino self-interaction becomes a key ingredient, in addition to the MSW effect, for understanding the neutrino process and the relevant nuclear abundances. The normal mass hierarchy is shown to be more compatible with the meteorite data. The main nuclear reactions for each nucleus are also investigated in detail.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Impact of hohlraum cooling on ignition metrics for inertial fusion implosions

This paper extends the evaluation of ignition metrics to include the impact of hohlraum cooling before peak implosion velocity in radiation driven implosions. First, we provide an extension of the results for the key hot spot stagnation quantities from the 2018 paper. The modified analytic expressions presented here match the Hydra results for these National Ignition Facility scale implosions both with and without hohlraum cooling before peak velocity if the effective ablation pressure P abl (effective) = P abl (t pv - 0.5 ns) is used in the analytic formulas, where t pv is the time of peak implosion velocity. Second, we provide an analysis that enables a comparison of the Hydra radiation hydrodynamics code calculations utilized here with the predictions of the analytic piston model of an ICF implosion, which focused on sensitivity to time duration of the hohlraum cooling phase before peak velocity (often called the “coast time”) and the shell radius at peak velocity R pv . Third, we provide a set of ignition metrics that are valid across a wide range of capsule designs valid for implosions both with and without hohlraum cooling before peak implosion velocity is reached.

42 ENGINEERING↗
DOE OSTI2018

Demonstration of Fully Coupled Calculations in Helical Steam Generator: Toward Predictions of Fluid Elastic Instability

The NEAMS program aims to develop an integrated multi-physics simulation capability “pellet-to-plant” for the design and analysis of future generations of nuclear power plants. In particular, the Reactor Product Line code suite's multi-resolution hierarchy is being designed to ultimately span the full range of length and time scales present in relevant reactor design and safety analyses, as well as scale from desktop to petaflop computing platforms. Flow-induced vibration (FIV) is widespread problem in energy systems because they rely on fluid movement for energy conversion. Vibrating structures may be damaged as fatigue or wear occurs. Given the importance of reliable components in the nuclear industry, flow-induced vibration has long been a major concern in safety and operation of nuclear reactors. In particular, nuclear fuel rods and steam generators have been known to suffer from flow-induced vibration and related failures. Advanced reactors, such as integral Pressurized Water Reactors (PWRs) considered for Small Modular Reactors (SMR), often rely on innovative component designs to meet cost and safety targets. One component that is the subject of advanced designs is the steam generator, some designs of which forego the usual shell-and-tube architecture in order to fit within the primary vessel. In addition to being more cost- and space-efficient, such steam generators need to be more reliable, since failure of the primary vessel represents a potential loss of coolant and a safety concern. A significant amount of data exists on flow-induced vibration in shell-and-tube heat exchangers, and heuristic methods are available to predict their occurrence based on a set of given assumptions. In contrast, advanced designs have far less data available. Advanced modeling and simulation based on coupled structural and fluid simulations have the potential to predict flow-induced vibration in a variety of designs, reducing the need for expensive experimental programs, especially at the design stage. Over the past seven years, the Reactor Product Line has developed the integrated multi-physics code suite SHARP. The goal of developing such a tool is to perform multi-physics neutronics, thermal/fluid, and structural mechanics modeling of the components inside the full reactor core or portions of it with a user-specified fidelity. In particular SHARP contains high-fidelity single-physics codes Diablo for structural mechanics and Nek5000 for fluid mechanics calculations. Both codes are state-of-the-art, highly scalable tools that have been extensively validated. These tools form a strong basis on which to build a flow-induced vibration modeling capability. In the last two years a significant effort was devoted to 1) develop and validate a fully coupled capability in SHARP for Fluid-Structure calculations (FSI) and 2) validate the use of SHARP for one-way coupled calculations in helical steam generators on legacy data available as well as novel experiments. In this report we discuss an extension of previous one-way coupled calculations performed with Nek5000 and Diablo aimed at simulating available FIV experiments in helical steam generators in the turbulent buffeting regime. It was demonstrated that in the buffeting regime one-way coupling was judged sufficient because the pressure loads do not cause substantial displacements. However, higher speeds led to a reduction in accuracy and the need for a fully coupled capability for velocity or approaching the critical velocity and fluid elastic instability. In this report we demonstrate such fully coupled capability on the Argonne test case discussed in previous reports.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Prediction of Redox Potentials for U, Np, Pu, and Am in Aqueous Solution

The redox properties of the actinides in aqueous solution are important for fuel production/reprocessing and understanding the environmental impact of nuclear waste. The redox potentials for U, Np, Pu, and Am in oxidation states from 0 up to VII (as appropriate) in aqueous solutions have been predicted at the density functional theory level with the B3LYP functional, Stuttgart small core pseudopotential basis sets for the actinides, and explicit (30H 2 O molecules)/implicit treatment of the aqueous solvent using the self-consistent reaction field COSMO and SMD approaches for the implicit solvation. The predictions of the structural parameters of clusters incorporating first and second solvation shells are consistent with the available experimental data., Our results are typically within 0.2 V of the available experimental data using two explicit solvation shells with an implicit solvent model. The use of the PW91 functional substantially improved the prediction of the Pu(VI/V) redox couple. The redox couples for An(VI/IV) and An(V/IV) which involve the addition of protons and removal of the actinyl oxygens led to slightly larger differences from experiment. Here, the An(IV/0) and An(III/0) couples were reliably predicted with our approach. Predictions of the unknown An(II/I) redox potentials were negative, consistent with expectations, and predictions for unknown An(VII/VI), An(III/II), and An(II/0) redox couples improve prior estimates.

Actinides↗
DOE OSTI2017

Wait for it: post-supernova winds driven by delayed radioactive decays

In most astrophysical situations, the radioactive decay of {sup 56}Ni to {sup 56}Co occurs via electron capture with a fixed half-life of 6.1 days. However, this decay rate is significantly slowed when the nuclei are fully ionized because K-shell electrons are unavailable for capture. In this paper, we explore the effect of these delayed decays on white dwarfs (WDs) that may survive Type Ia and Type Iax supernovae (SNe Ia and SNe Iax). The energy released by the delayed radioactive decays of {sup 56}Ni and {sup 56}Co drives a persistent wind from the surviving WD’s surface that contributes to the late-time appearance of these SNe after emission from the bulk of the SN ejecta has faded. We use the stellar evolution code MESA to calculate the hydrodynamic evolution and resulting light curves of these winds. Our post-SN Ia models conflict with late-time observations of SN 2011fe, but uncertainties in our initial conditions prevent us from ruling out the existence of surviving WD donors. Much better agreement with observations is achieved with our models of post-SN Iax bound remnants, providing evidence that these explosions are due to deflagrations in accreting WDs that fail to completely unbind the WDs. Future radiative transfer calculations and wind models utilizing simulations of explosions for more accurate initial conditions will extend our study of radioactively powered winds from post-SN surviving WDs and enable their use as powerful discriminants among the various SN Ia and SN Iax progenitor scenarios.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Properties of Medium Nuclei from First Principles

During this project we have extended the Auxiliary Field Diffusion Monte Carlo (AFDMC) to calculate several properties of light/medium nuclei. In particular, we have included and tested several models of three-body forces and the associated approximations required, and the comparison with other exact calculations for simpler problems demonstrated the accuracy of the AFDMC method. In this project we have also developed new subroutines to solve for the ground-state of open shell nuclei. This case requires very complex wave function with respect to the case of closed-shell nuclei. We have used the AFDMC code to calculate many properties of nuclei up to the Oxygen, including energies, radii, spatial and momentum distributions and others.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Ab initio based examination of the kinetics and thermodynamics of oxygen in Fe-Cr alloys

The behavior of oxygen in Fe-Cr binary alloys is important to understand for a variety of applications including fuel cell interconnects and nuclear energy systems. In this work, we performed an ab initio investigation of the binding of oxygen in the entire compositional range of Fe-Cr at minima and saddles. This database was subsequently used to parametrize concentration-dependent local cluster expansions which were utilized to examine the kinetics of oxygen transport in the alloy and the thermodynamics of the system. The behavior of oxygen in the alloy system was investigated using our cluster expansion model and a convex hull diagram was recorded. The most favorable composition for the incorporation of oxygen was found to occur at 70% Cr. The dependence of the oxygen solution energy on the Cr occupation of the nearest neighbor shells was carefully characterized. It was determined that the third nearest neighbor Cr occupation was unfavorable and that, generally speaking, the first and second nearest neighbor shell Cr occupation was favored with a few exceptions. Kinetic Monte Carlo models at dilute oxygen levels were performed using two models [a kinetically resolved activation (kRA) model and a cluster expansion model for the saddle points]. The models qualitatively agreed with the observed trends in the literature which reported a decrease in diffusivity at dilute Cr levels, but the models predicted different behaviors beyond the dilute Cr limit. Finally, using grand canonical Monte Carlo, we further examined the thermodynamics of oxygen incorporation into the alloy. Furthermore, our results offer new insight into the behavior of oxygen in Fe-Cr alloys that have ramifications on the early stages of the corrosive behavior of the material.

36 MATERIALS SCIENCE↗
NASA NTRS1978

Barium and neodymium isotopic anomalies in the Allende meteorite

The discovery of Ba and Nd isotopic anomalies in two inclusions from the Allende meteorite is reported. The inclusions are Ca-Al-rich objects typical of the type considered as high-temperature condensation products in the solar nebula and contain distinctive Mg and O isotopic anomalies of the FUN (mass Fractionation, Unknown Nuclear processes) type. Mass-spectrometry results are discussed which show that inclusion C1 has anomalies in Ba at masses 134 and 136, while inclusion EK1-4-1 exhibits large marked negative anomalies at 130, 132, 134, and 136, as well as a positive anomaly at 137. It is also found that inclusion EK1-4-1 shows marked negative anomalies in Nd at masses 142, 146, 148, and 150, in addition to a positive anomaly at 145. These isotopic shifts are attributed to addition of r-process nuclei rather than mass fractionation. It is suggested that an onion-shell supernova explosion followed by injection into the solar nebula is the most likely generic model that may explain the observations.

Mcculloch, M. T.↗
DOE OSTI2018

Structure of the exotic He 9 nucleus from the no-core shell model with continuum

Here, the exotic 9 He nucleus, which presents one of the most extreme neutron-to-proton ratios, belongs to the N = 7 isotonic chain famous for the phenomenon of ground-state parity inversion with decreasing number of protons. Consequently, it would be expected to have an unnatural (positive) parity ground state similar to 11 Be and 10 Li. Despite many experimental and theoretical investigations, its structure remains uncertain. Apart from the fact that it is unbound, other properties including the spin and parity of its ground state, and the very existence of additional low-lying resonances are still a matter of debate

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗