Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2016

Dual origin of pairing in nuclei

The pairing correlations of the nucleus {sup 120}Sn are calculated by solving the Nambu–Gor’kov equations, including medium polarization effects resulting from the interweaving of quasiparticles, spin and density vibrations, taking into account, within the framework of nuclear field theory (NFT), processes leading to self-energy and vertex corrections and to the induced pairing interaction. From these results one can not only demonstrate the inevitability of the dual origin of pairing in nuclei, but also extract information which can be used at profit to quantitatively disentangle the contributions to the pairing gap Δ arising from the bare and from the induced pairing interaction. The first is the strong {sup 1}S{sub 0} short-range NN potential resulting from meson exchange between nucleons moving in time reversal states within an energy range of hundreds of MeV from the Fermi energy. The second results from the exchange of vibrational modes between nucleons moving within few MeV from the Fermi energy. Short- (v{sub p}{sup bare}) and long-range (v{sub p}{sup ind}) pairing interactions contribute essentially equally to nuclear Cooper pair stability. That is to the breaking of gauge invariance in open-shell superfluid nuclei and thus to the order parameter, namely to the ground state expectation value of the pair creation operator. In other words, to the emergent property of generalized rigidity in gauge space, and associated rotational bands and Cooper pair tunneling between members of these bands.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Shell structure in 32P

The nuclear structure of 32P was investigated by detecting multi-γ coincidences along with charged light ions following the bombardment of an 18O target with a 30-MeV 16O beam from the Florida State University accelerator facility using the Clarion2-Trinity array of Clover γ spectrometers and GAGG (Gd-Al-Ga-garnet) scintillators for charged particles. A number of new states up to 12 MeV excitation and spins up to 9+ were observed by their γ decay patterns. Spins and parities were assigned by comparing measured γ angular distributions with polarizations inferred from Compton-scattering asymmetries between the Ge crystals in the Clover spectrometers. The new level scheme of 32P compares well with shell model calculations using the FSU cross-shell interaction and reasonably well with those using the sdpf-m interaction. In particular, structures with higher spin arise from promotion of more nucleons up to the f-p shell.

Allmond, James [ORNL] (ORCID:0000000165338721)↗
DOE OSTI2015

Immobilized Low-Activity Waste Glass Release Data Package for the Integrated Disposal Facility Performance Assessment

The Hanford Site in southeastern Washington State has been used extensively by the U.S. Department of Energy (DOE) to produce nuclear materials for the U.S. defense arsenal. A large inventory of radioactive and mixed waste has accumulated in 177 buried single- and double-shell tanks. Liquid waste recovered from the tanks will be pretreated to separate the low-activity fraction from the high-level and transuranic wastes. The low-activity waste (LAW) is planned to be immobilized in glass, via vitrification, and placed into a near-surface disposal system on the Hanford Site. Before the immobilized low-activity waste (ILAW) can be placed into the disposal system, DOE must approve a performance assessment (PA), which documents the assumptions and provides the quantitative demonstration of compliance with the performance objectives for the long-term protection of the public and the environment. A critical component of the PA will be to demonstrate that releases from the selected glass waste form do not result in the performance objectives being exceeded. The PA must also, for purposes of establishing limits on radionuclides that may be buried near-surface, assess impacts to a hypothetical intruder and impacts to water resources. The key information presented in this data package consists of the specific glass compositions and associated glass dissolution rate law parameters for five ILAW glasses. In addition, the results of geochemical modeling performed on several glasses that are not specifically described in this document are provided as supplemental data to help identify potential end-point alteration phases that may form during the weathering of ILAW glasses. These data (glass compositions, rate law parameters, and potential alteration phases) will be used to identify contaminant release source terms for the Integrated Disposal Facility PA. In addition, a general discussion on nuclear waste glass corrosion and a brief overview of modeling approaches and corresponding uncertainties are discussed. Table ES.1 through Table ES.4 provide the specific data to be used to support the 2017 PA.

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

Two-neutron halo structure of 31 F

Here, we apply the Gamow shell model to study 25–31 F isotopes. As both internucleon correlations and continuum coupling are properly treated therein, the structure shape of 31 F at large distance can be analyzed precisely. For this, one-nucleon densities, root-mean-square radii, and correlation densities are calculated in neutron-rich fluorine isotopes. It is then suggested that 31 F exhibits a two-neutron halo structure, built from both continuum coupling and nucleon-nucleon correlations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Shell structure of 43 S and collapse of the N = 28 shell closure

The single-particle structure of the N = 27 isotones provides insights into the shell evolution of neutron-rich nuclei from the doubly-magic 48 Ca toward the drip line. 43 S was studied employing the one-neutron knockout reaction from a radioactive 44 S beam. Using a combination of prompt and delayed γ-ray spectroscopy the level structure of 43 S was clarified. Momentum distributions were analyzed and allowed for spin and parity assignments. The deduced spectroscopic factors show that the 44 S ground-state configuration has a strong intruder component. Here, the results were confronted with shell model calculations using two effective interactions. General agreement was found between the calculations, but strong population of states originating from the removal of neutrons from the 2p 3/2 orbital in the experiment indicates that the breakdown of the N = 28 magic number is more rapid than the theoretical calculations suggest.

39 ≤ A ≤ 58↗
DOE OSTI2022

Ab Initio Prediction of the 4 $\mathrm{H}$e($d, γ$) 6 $\mathrm{L}$i Big Bang Radiative Capture

We report the rate at which helium ( 4 He ) and deuterium (d) fuse together to produce lithium-6 ( 6 Li) and a γ ray, 4 $\mathrm{H}$e($d, γ$) 6 $\mathrm{L}$i, is a critical puzzle piece in resolving the discrepancy between big bang predictions and astronomical observations for the primordial abundance of 6 Li . The accurate determination of this radiative capture rate requires the quantitative and predictive description of the fusion probability across the big bang energy window (30 keV ≲ E ≲ 400 keV), where measurements are hindered by low counting rates. We present first-principle (or, ab initio) predictions of the 4 $\mathrm{H}$e($d, γ$) 6 $\mathrm{L}$i astrophysical S factor using validated nucleon-nucleon and three-nucleon interactions derived within the framework of chiral effective field theory. By employing the ab initio no-core shell model with continuum to describe 4 He-d scattering dynamics and bound 6 Li product on an equal footing, we accurately and consistently determine the contributions of the main electromagnetic transitions driving the radiative capture process. Our results reveal an enhancement of the capture probability below 100 keV owing to previously neglected magnetic dipole (M1) transitions and reduce by an average factor of 7 the uncertainty of the thermonuclear capture rate between 0.002 and 2 GK.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Multiparticle-hole excitations in nuclei near N = Z = 20: $^{41}$K

This experimental study of high-spin structure near the N = Z = 20 region focuses on 41 K and also reports three newly observed γ transitions in 41 Ca from the same reaction. High-spin states were populated using the 26 Mg( 18 O, p2nγ) 41 K and 26 Mg( 18 O, 3nγ) 41 Ca reactions at a beam energy of 50 MeV at the Florida State University John D. Fox Superconducting Linear Accelerator Laboratory, employing the FSU high-purity germanium detector array. The level scheme of 41 K was extended to 12325 keV, possibly up to J π = 25/2 − or 27/2 + , by means of 25 new γ-ray transitions, and that of 41 Ca to 9916 keV. Linear polarization and angular-distribution measurements were used to provide spin and parity information for several states in the 41 K decay scheme. The results are compared with spsdpf cross-shell shell-model calculations using the FSU interaction. Configurations involving zero or one nucleon promoted from the sd to the fp shell reproduce the energies of many known states reasonably well, while multi-particle excitations reveal a more complex interplay of single-particle motion, collective pairing, and deformation, posing an interesting challenge for future theoretical work.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2014

X-ray relative intensities at incident photon energies across the L{sub i} (i=1–3) absorption edges of elements with 35≤Z≤92

The intensity ratios, I{sub Lk}/I{sub Lα1} (k=l,η,α{sub 2},β{sub 1},β{sub 2,15},β{sub 3},β{sub 4},β{sub 5,7},β{sub 6},β{sub 9,10},γ{sub 1,5},γ{sub 6,8},γ{sub 2,3},γ{sub 4}) and I{sub Lj}/I{sub Lα} (j=β,γ), have been evaluated at incident photon energies across the L{sub i} (i=1–3) absorption edge energies of all the elements with 35≤Z≤92. Use is made of what are currently considered to be more reliable theoretical data sets of different physical parameters, namely, the L{sub i} (i=1–3) sub-shell photoionization cross sections based on the relativistic Hartree–Fock–Slater (RHFS) model, the X-ray emission rates based on the Dirac–Fock model, and the fluorescence and Coster–Kronig yields based on the Dirac–Hartree–Slater model. In addition, the Lα{sub 1} X-ray production cross sections for different elements at various incident photon energies have been tabulated so as to facilitate the evaluation of production cross sections for different resolved L X-ray components from the tabulated intensity ratios. Further, to assist evaluation of the prominent (L{sub i}−S{sub j}) (S{sub j}=M{sub j}, N{sub j} and i=1–3, j=1–7) resonant Raman scattered (RRS) peak energies for an element at a given incident photon energy (below the L{sub i} sub-shell absorption edge), the neutral-atom electron binding energies based on the relaxed orbital RHFS calculations are also listed so as to enable identification of the RRS peaks, which can overlap with the fluorescent X-ray lines. -- Highlights: •The L X-ray relative intensities and Lα{sub 1} XRP cross sections are evaluated using physical parameters based on the IPA models. •Comparison of the intensity ratios evaluated using the DHS and DF models based photoionization cross sections is presented. •Importance of many body effects including electron exchange effects is highlighted.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Lower-length-scale modeling of chemical additions, corrosion and fission products in select MSR base salts

Molten salts play a crucial role in advancing Generation IV nuclear energy technology, with chloride-based salts like NaCl-UCl 3 garnering significant attention due to their distinctive properties. However, the corrosive nature of molten salts causes the dissolution of chromium (Cr), leading to the formation of CrCl 2 and CrCl 3 species in molten chlorides. Moreover, the radioactive decay of nuclear fuel gives rise to fission products, including Cs, Sr, and I chlorides. The first part of this report presents a comprehensive study utilizing ab initio molecular dynamics (AIMD) simulations to investigate the properties of eutectic NaCl-UCl 3 molten salt in the presence of corrosion products (CrCl 2 and CrCl 3 ) and fission product (CsCl). The study focuses on essential structural and thermophysical properties such as density, mixing energy, coordination numbers (CN), and Radial Distribution Functions (RDF) of the salts with varying compositions of corrosion products and fission product, covering a range from 0% to 13.5%. The results offer valuable insights into the behavior of corrosion and fission products in uranium-based molten salts, providing essential data that can be used as input to the MTDB-TC and MSTDB-TP property databases being developed by the NEAMS program. Due to their favorable characteristics such as low melting points, high boiling points, and low costs, MgCl 2 +NaCl+KCl (MNK) eutectic salts have recently attracted attention as high-temperature heat transfer fluids. The incorporation of LiCl into MNK salts can further reduce their melting points and increase their specific heat capacities, which is de- sirable for high-temperature heat transfer applications. The second part of this report presents the development and validation of a new shell-model potential for the MgCl 2 +NaCl+KCl+LiCl system, which captures the polarization of Cl anions. The extensive comparison with experimental data and AIMD simulations demonstrates the accuracy and reliability of the potential. Furthermore, using this potential, we elucidate the intricate network structure in MNK eutectic salts that is formed through polyhedron sharing. This research contributes to a fundamental understanding of the atomic structures and thermophysical properties of multi-component molten salts, which is crucial for future development of heat transfer fluids for applications in molten salt reactors.

36 MATERIALS SCIENCE↗