SEARCH · Search NASA

Results for “Small Core”

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.

SourceResearch connectionFrom yearThrough year

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

At least 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2025

A small core in Vesta inferred from Dawn’s observations

Vesta’s large-scale interior structure had previously been constrained primarily using the gravity and shape data from the Dawn mission. However, these data alone still allow a wide range of possibilities for the differentiation state of the body. The moment of inertia is arguably the most diagnostic parameter related to the radial density distribution of a planetary body, making it crucial for assessing the body’s state of internal differentiation. Determining the moment of inertia requires additional measurements of the amplitudes of small rotational motions, such as precession and nutation. Here we report an updated estimate of the moment of inertia of Vesta inferred from Dawn’s Doppler tracking via the Deep Space Network and onboard imaging data. The recovered value for Vesta’s normalized polar moment of inertia is $\overline{C}$/MR 2 = 0.4208 ± 0.0047 (where M is the mass of Vesta and R is the reference radius), which is only 6.6% lower than the homogeneous value of 0.4505. This value, combined with the gravity field and global shape, suggests that Vesta’s interior has limited density stratification beneath its howardite–eucrite–diogenite-dominated crust. We propose two possible origin scenarios that are consistent with the observed constraints. In the first scenario, Vesta’s interior did not undergo full differentiation due to late accretion. In the second scenario, Vesta originated as an impact remnant of a larger differentiated body re-accreted with non-chondritic bulk composition produced from a catastrophic impact. Vesta did not experience complete differentiation in either scenario, suggesting that its current state reflects a complex interplay between its accretion timing, thermal evolution, redistribution of 26 Al bearing melt and/or impact processes.

CNEOS 2014-01-08 bolide↗

DOE OSTI2024

Liquid Structure of Iron and Iron–Nitrogen–Carbon Alloys Within the Cores of Small Terrestrial Bodies

Abstract Nitrogen has been proposed to be stored within planetary cores, but its effects on the structure and density of molten Fe–alloys have not been explored experimentally. Using energy‐dispersive X‐ray diffraction, we determined the structure of Fe–N(–C) liquids at core conditions (1–7 GPa and 1700–1900°C) within a Paris‐Edinburgh press. Variation of N up to 7 wt.% and C up to 1.5 wt.% results in near‐linear changes in Fe–Fe atom distances and structure factor with increasing light element content. We did not observe a significant pressure‐driven structural transition in Fe–N(–C) liquids. We model the expansion of the Fe–Fe bonds using a modified Birch‐Murnaghan equation of state. With this model, we demonstrate that N or C contamination could lead to an overestimation of the Fe–Fe distances of pure Fe. We observe that the incorporation of 1 wt.% N or C into Fe results in a change in Fe–Fe distances that is twice as significant as the effect of 1 GPa. By approximating the change in volume, we infer that N and C incorporated in liquid iron could contribute to the density deficit observed in the cores of terrestrial bodies.

Pease, Allison [Department of Earth and Environmen↗

DOE OSTI2025

Advancements in NEAMS Tool Capabilities for Multiphysics Simulation of Fast Reactor Core Bowing and Identification of Validation Test Data

Under the U.S. Department of Energy Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) Program, an integrated multiphysics approach is being developed to model the core bowing phenomena important to liquid metal-cooled fast reactors. Core bowing is an important passive safety mechanism in liquid metal-cooled fast reactors and involves Multiphysics effects including radiation transport, fluid flow, heat transfer, and mechanical response to temperature and flux gradients. This report summarizes recent progress on developing a multiphysics, MOOSE-based workflow to predict core bowing and associated reactivity feedback. Last year, thermal fluids and mechanics were coupled on a multi-assembly benchmark problem based on ABR-1000 design. This year, the reactor physics code Griffin was assessed for readiness of core bowing calculations. Preliminary integration of Griffin’s ring-heterogeneous model with thermal fluids and thermal mechanics solvers was performed. Specifically, thermal-mechanics and reactor physics were coupled for single- and multi-assembly problems, and reactor physics and subchannel methods were coupled for a single assembly model. Finally, the workflow of all three physics was preliminarily demonstrated on a single assembly model. Caveats and future development needed have been identified. To supplement the multiphysics demonstration, verification and assessment efforts of thermos-mechanical capabilities for modeling thermo-mechanical core bowing behavior were continued by analyzing IAEA Verification Problem 5 which includes radiation swelling and creep. Additionally, a small core reactor physics benchmark defined by Japan Atomic Energy Agency (JAEA) was performed to assess neutronics models for estimating reactivity feedback. Finally, Fast Flux Test Facility (FFTF) validation test data for core bowing phenomena has been identified and summarized, with a recommended path forward for validation once this capability is mature.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Robust Implicit Adaptive Low Rank Time-Stepping Methods for Matrix Differential Equations

In this work, we develop implicit rank-adaptive schemes for time-dependent matrix differential equations. The dynamic low rank approximation (DLRA) is a well-known technique to capture the dynamic low rank structure based on Dirac–Frenkel time-dependent variational principle. In recent years, it has attracted a lot of attention due to its wide applicability. Our schemes are inspired by the three-step procedure used in the rank adaptive version of the unconventional robust integrator (the so called BUG integrator) (Ceruti et al. in BIT Numer Math 62(4):1149–1174, 2022) for DLRA. First, a prediction (basis update) step is made computing the approximate column and row spaces at the next time level. Second, a Galerkin evolution step is invoked using an implicit solves for the small core matrix. Finally, a truncation is made according to a prescribed error threshold. Since the DLRA is evolving the differential equation projected on to the tangent space of the low rank manifold, the error estimate of the BUG integrator contains the tangent projection (modeling) error which cannot be easily controlled by mesh refinement. This can cause convergence issue for equations with cross terms. To address this issue, we propose a simple modification, consisting of merging the row and column spaces from the explicit step truncation method together with the BUG spaces in the prediction step. In addition, we propose an adaptive strategy where the BUG spaces are only computed if the residual for the solution obtained from the prediction space by explicit step truncation method, is too large. Here, we prove stability and estimate the local truncation error of the schemes under assumptions. We benchmark the schemes in several tests, such as anisotropic diffusion, solid body rotation and the combination of the two, to show robust convergence properties.

97 MATHEMATICS AND COMPUTING↗

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 OSTI2024

Prediction of Redox Potentials for Ac, Th, and Pa in Aqueous Solution

Density functional theory in conjunction with small core pseudopotentials and the associated basis sets was used to calculate potentials for multiple redox couples, covering a range of oxidation states for Ac (0 to III), Th (0 to IV), and Pa (0 to V) in aqueous solution. Solvation effects were incorporated using a supermolecule-continuum approach, with 30 water molecules representing two solvation shells, and the COSMO and SMD implicit solvation models. The calculated geometries for Ac(III), Th(IV), and Pa(V) were in reasonable agreement with the available experimental data. Using the COSMO model with the B3LYP functional, the calculated redox potentials were within ± 0.2 V from experiment for most redox couples. Several pathways were explored for the Pa(V/IV) redox couple for different forms of Pa(V) and Pa(IV). Most Pa(V/IV) redox couples have very similar potentials, ranging from 0 to -0.4 V up to a pH of 1.4. At pH = 1.4, the potentials shift to values that are more negative than -0.7 V, reflecting the growing unfavorable nature of the redox process at higher pH levels. The calculated values for An(III/II) potentials were consistent with prior estimates and the available experimental data. The predicted redox potentials for An(II/I) were highly negative, as expected. For An(I/0) potentials, Th and Pa exhibited positive values, contrasting with the negative values calculated for Ac. Furthermore, the An +m /An(0) potentials agreed better with the experimental data when using the COSMO solvation model as compared to the SMD model.

Chemical calculations↗

DOE OSTI2025

Prediction of Redox Potentials for Different Oxidation States of U, Np, Pu, and Am in Alkaline Aqueous Solution

The redox potentials for U, Np, Pu, and Am for oxidation states +III up to +VIII in alkaline aqueous solutions were predicted using density functional theory (DFT) and small-core pseudopotentials and their basis sets, with a hybrid explicit/implicit solvent model using SHE = 4.28 V. For each oxidation state, various oxo/hydroxo complexes were evaluated resulting in a variety of one-electron redox pathways. For An(VIII/VII) couples, the predicted redox potentials for the [An(VIII)O 5 (OH)] -3 /[An(VII)O 4 (OH) 2 ] -3 , or [An(VIII)O 4 (OH) 2 ] -2 /[An(VII)O 4 (OH) 2 ] -3 couples are in good agreement with existing estimates. For An(VII/VI) redox couples, all couples, particularly [An(VII)O 4 (OH) 2 ] -3 /[An(VI)O 2 (OH) 4 ] -2 , were in agreement with experimental values for U, Np, and Pu, but the results for Am showed larger differences from the estimated potentials. The An(VI/V) couples were consistent with experiment for dioxo/tetrahydroxo couples, and the An(V/IV) couples showed acceptable agreement based on actinide-specific couples, with neutral hydroxides often favored in the +IV state. The An(IV/III) couples were consistent with literature values when modeled as soluble neutral hydroxides. The use of our approach yielded calculated redox potentials that were within ±0.2 V from experimental or estimated values consistent with our prior calculations on redox potentials of actinides from Ac to Am in acidic aqueous solutions. Finally, as a result, this supports the robustness of our DFT-based methodology for predicting actinide redox potentials, offering valuable insights for actinide chemistry in aqueous solutions.

Actinides↗

DOE OSTI2025

Prediction of Redox Potentials for the Late Actinides Cm to Lr Using Electronic Structure Methods

Our previously developed computational method for calculating the aqueous redox potentials of the early actinides has been extended to the later elements in the actinide series: Cm, Bk, Cf, Es, Fm, Md, No, and Lr in multiple oxidation states. These calculations were performed using density functional theory with small-core pseudopotentials and their associated basis sets. Solvation effects were considered via a supermolecule-continuum approach, with 30 water molecules representing two solvation shells. Both the COSMO and SMD implicit solvation models were utilized. The structural parameters and hydration numbers for Cm(III), Bk(III), Bk(IV), and Cf(III) are in reasonable agreement with the available experimental data. For redox processes involving atomic cations in solution, the B3LYP/COSMO approach predicted redox potentials to within ±0.2 V of experiment for most redox couples, consistent with our prior work. Inclusion of spin-orbit corrections in specific redox pairs, especially those with the later actinides in high oxidation states, yields improved results relative to calculations including only scalar-relativistic corrections. The An +m /An(0) redox potentials were calculated using a Born-Haber cycle incorporating sublimation, ionization, and hydration energies. Due to a lack of experimental data, three sets of ionization energies were used for the Born-Haber cycle. The calculated An(III/0) potentials showed better agreement with experimental data when using the COSMO solvation model and the test set comprising the NIST recommended ionization energies. Furthermore, the Md(II/0) potential was better described with the SMD model, whereas No(II/0) was not well described by all methods. Finally, the computational approach was able to predict redox potentials that for most cases agreed with the current available experimental or estimated data.

Actinides↗

DOE OSTI2026

WEC Irradiation Can Thermal Testing Report

With renewed interest for the next generation of small modular reactors and space-based nuclear reactors, metal hydrides have been a focus of research for moderating materials to reduce the need for highly enriched fuel to maintain small core sizes. Yttrium hydride (YH) is an attractive candidate for these metal hydrides due to its ability to retain hydrogen at elevated temperatures as well as the low neutron absorption cross section of yttrium metal. Given the necessity for reactors to operate at elevated temperatures for extended periods of time, understanding hydrogen retention in moderator channels in simulated operating thermal conditions is paramount in predicting reactor performance and the long-term stability of YH. Additionally, coating strategies on the moderator encapsulating will be explored to ascertain if it aids in retention of hydrogen within the capsule and moderator.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

Planet Formation by Gas-assisted Accretion of Small Solids

We compute the accretion efficiency of small solids, with radii 1 cm ≤ R s ≤ 10 m, on planets embedded in gaseous disks. Planets have masses 3 ≤ M p ≤ 20 Earth masses (M ⊕ ) and orbit within 10 au of a solar mass star. Disk thermodynamics is modeled via 3D radiation-hydrodynamics calculations that typically resolve the planetary envelopes. Both icy and rocky solids are considered, explicitly modeling their thermodynamic evolution. The maximum efficiencies of 1 ≤ R s ≤ 100 cm particles are generally ≲10%, whereas 10 m solids tend to accrete efficiently or be segregated beyond the planet’s orbit. A simplified approach is applied to compute the accretion efficiency of small cores, with masses M p ≤ 1 M ⊕ and without envelopes, for which efficiencies are approximately proportional to $M^{2/3}_{p}$. The mass flux of solids, estimated from unperturbed drag-induced drift velocities, provides typical accretion rates dM p /dt ≲ 10 -5 M ⊕ yr -1 . In representative disk models with an initial gas-to-dust mass ratio of 70–100 and total mass of 0.05–0.06 M ⊙ , the solids’ accretion falls below 10 -6 M ⊕ yr -1 after 1–1.5 Myr. The derived accretion rates, as functions of time and planet mass, are applied to formation calculations that compute dust opacity self-consistently with the delivery of solids to the envelope. Assuming dust-to-solid coagulation times of ≈0.3 Myr and disk lifetimes of ≈3.5 Myr, heavy-element inventories in the range 3–7 M ⊕ require that ≈90–150 M ⊕ of solids cross the planet’s orbit. The formation calculations encompass a variety of outcomes, from planets a few times M ⊕ , predominantly composed of heavy elements, to giant planets. The peak luminosities during the epoch of the solids’ accretion range from ≈10 -7 to ≈10 -6 L ⊙ .

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

$S^5$: Tidal Disruption in Crater 2 and Formation of Diffuse Dwarf Galaxies in the Local Group

We present results of a spectroscopic campaign around the diffuse dwarf galaxy Crater 2 (Cra2) and its tidal tails as part of the Southern Stellar Stream Spectroscopic Survey ($S^5$). Cra2 is a Milky Way dwarf spheroidal satellite with extremely cold kinematics, but a huge size similar to the Small Magellanic Cloud, which may be difficult to explain within collisionless cold dark matter. We identify 143 Cra2 members, of which 114 belong to the galaxy's main body and 29 are deemed part of its stellar stream. We confirm that Cra2 is dynamically cold (central velocity dispersion $2.51^{+0.33}_{-0.30}\,{\rm km\,s^{-1}}$) and also discover a $\approx$7$Ï$ velocity gradient consistent with its tidal debris track. We separately estimate the stream velocity dispersion to be $5.74^{+0.98}_{-0.83}\,{\rm km\,s^{-1}}$. We develop a suite of $N$-body simulations with both cuspy and cored density profiles on a realistic Cra2 orbit to compare with $S^5$ observations. We find that the velocity dispersion ratio between Cra2 stream and galaxy ($2.30^{+0.41}_{-0.35}$) is difficult to reconcile with a cuspy halo with fiducial concentration and an initial mass predicted by standard stellar mass$-$halo mass relationships. Instead, either a cored halo with relatively small core radius or a low-concentration cuspy model can reproduce this ratio. Despite tidal mass loss, Cra2 is metal-poor ($\langle \rm[Fe/H]\rangle=-2.16\pm0.04$) compared to the stellar mass$-$metallicity relation for its luminosity. Other diffuse dwarf galaxies similar to Cra2 in the Local Group (Antlia 2 and Andromeda 19) also challenge galaxy formation models. Finally, we discuss possible formation scenarios for Cra2, including ram-pressure stripping of a gas-rich progenitor combined with tides.

Limberg, Guilherme [Chicago U., KICP; Chicago U.] ↗

DOE OSTI2024

Strength, plasticity, and spin transition of Fe-N compounds in planetary cores

Elastic and plastic properties of Fe-light element alloys and compounds are needed to determine the compositions and dynamics of planetary cores. Elastic strength and plastic deformation mechanisms and their relationship to electronic properties of ε-Fe 7 N 3 and γ’-Fe 4 N mixture were investigated by x-ray diffraction and x-ray emission spectroscopy in the diamond anvil cell from 1 bar up to 60 GPa. X-ray diffraction shows that ε-Fe 7 N 3 reaches a pressure of 15–20 GPa before undergoing bulk plasticity at a differential stress of 4.4–10.4 GPa. ε-Fe 7 N 3 is stronger than γ’-Fe 4 N and hcp-Fe which achieve a flow stress of 1.5–3.6 GPa at 10–15 GPa and 2–3 GPa at ~20 GPa, respectively. X-ray emission spectroscopy shows that a decrease in electronic spin moment begins before and completes after plastic flow onset for each nitride, suggesting that pressure-driven changes in electronic arrangement do not trigger a plastic response although they may modify the strength and plastic behavior of Fe-N compounds. Plastic deformation in ε-Fe 7 N 3 and hcp-Fe results in a preferred orientation of (0001) normal to maximum compression, while γ’-Fe 4 N develops a maximum in the (110). Furthermore, these observations may be combined with measurements of elasticity to model seismic properties of cores of small planetary bodies such as Mars, Mercury, and the Moon.

58 GEOSCIENCES↗

DOE OSTI2025

SPC-71188 Rev 0 MARVEL Control Drum Rotary Shaft Seal Procurement Specification

Idaho National Laboratory (INL) is developing a microreactor to produce electrical power utilizing a small nuclear core under the Microreactor Application Research Validation and Evaluation (MARVEL) program. This procurement specification defines the requirements for fabrication of the MARVEL Control Drum Rotary Shaft Seal.