Search NASA⌕ Search

SEARCH · Search NASA

Results for “mixture density network”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Numerical Investigation of Ammonia/n-Heptane Dual-Fuel Spray Flames Using Large Eddy Simulations

Diesel engines are extensively used in heavy-duty transportation, power generation, and marine vehicles due to their superior thermal efficiency and extended high-load operability compared to spark ignition (SI) engines. However, combustion in diesel engines is generally characterized by locally rich fuel–air mixtures and high combustion temperatures, causing significant amounts of soot and NO x emissions from these engines. Utilizing carbon-free alternative fuels and enhancing fuel efficiency represent promising strategies to mitigate greenhouse gas (GHG) and other emissions in the heavy-duty transportation sector. In this context, ammonia (NH 3 ), as a hydrogen carrier, has received significant attention as a viable substitute for hydrocarbon fuels due to its carbon-free composition, relatively high energy density, and well-established infrastructure. Many previous studies have considered combustion and emission characteristics of ammonia-hydrocarbon fuel blends in engines and simplified flames. But, detailed investigations on the effects of ammonia on the performance of hydrocarbon fuels under engine conditions are lacking. In the present study, we perform large eddy simulations (LES) of the ignition and flame processes in a constant-volume combustion reactor, where n-heptane is injected in an ammonia/air ambient mixture in a diesel-like environment. A detailed and validated reaction mechanism containing 302 species and 1981 reactions is employed. The Engine Combustion Network Spray H experimental data is used to validate the spray model under both non-reacting and reacting conditions. Dual-fuel combustion is simulated using the well-stirred reactor (WSR) approach. Results are presented for two spray cases: (1) single fuel (SF) with n-heptane injected into a mixture of air and combustion products and (2) dual-fuel (DF) with the injection of n-heptane in a mixture of air, ammonia, and combustion products. It is observed that the presence of ammonia has a significant effect on the ignition and flame development processes. With ammonia addition, both the first- and second-stage ignition delay times increase, but the effect of ammonia on the second-stage ignition is significantly more prominent. In addition, the ignition kernel size and growth rate decrease noticeably. For SF spray, the main ignition is characterized by multiple ignition kernels near the spray tip, whereas for DF spray, a single relatively small ignition kernel forms and grows slowly in the downstream direction. The flame development and the final quasi-steady flame structure are also modified due to ammonia. Here, the outcome of this research would enable a better understanding of ammonia–diesel dual-fuel spray flame behavior and guide the development of associated engine combustion strategies.

ammonia↗
DOE OSTI2020

High pressure/high temperature multiphase simulations of dodecane injection to nitrogen: Application on ECN Spray-A

The present work investigates the complex phenomena associated with pressure/high temperature dodecane injection for the Engine Combustion Network (ECN) Spray-A case, employing more elaborate thermodynamic closures, to avoid well known deficiencies concerning density and speed of sound prediction using traditional cubic models. A tabulated thermodynamic approach is proposed here, based on log 10 (p)-T tables, providing very high accuracy across a large range of pressures, spanning from 0 to 2500 bar, with only a small number of interpolation points. The tabulation approach is directly extensible to any thermodynamic model, existing or to be developed in the future. Here NIST REFPROP properties are used, combined with PC-SAFT Vapor-Liquid-Equilibrium to identify the liquid in mixtures penetration, hence avoiding the use of an arbitrary threshold for mass fraction. Identified liquid and vapour penetration are compared against experimental data from the ECN database showing a good agreement, within approximately 3–8% for axial penetration of liquid, 2% for vapor axial penetration and within experimental uncertainty for radial distribution of mass fraction. Analysis of the vortex evolution indicates that driving mechanisms behind the jet break-up are vortex tilting/stretching, then baroclinic torque, leading to Rayleigh-Taylor instabilities, closely followed by vortex dilation and finally viscous effects.

42 ENGINEERING↗
DOE OSTI2022

Understanding radiation effects in friction stir welded MA956 using ion irradiation and a rate theory model

We report an outstanding challenge in the manufacturing and joining of oxide dispersion strengthened steels is retaining the nanofeatures in the alloy throughout the fabrication and welding process. MA956 was friction stir welded with two different sets of welding parameters, resulting in a medium and high heat input. After welding, 5 MeV Fe ++ ion irradiations were performed at doses ranging from 50 to 200 dpa in the temperature range of 400 to 500°C. Post-irradiation characterization was performed with scanning transmission electron microscopy and energy-dispersive x-ray spectroscopy to investigate the Y-Al-O dispersoids, voids, and dislocations. After welding, the dispersoid microstructure coarsened, resulting in fewer and larger dispersoids regardless of heat input. After irradiation, the dispersoid behavior in the welded material was sensitive to temperature, exhibiting growth behavior attributed to Ostwald coarsening at 500 °C but a mixture of nucleation and more muted growth at 400 and 450 °C, attributed to competing mechanisms of radiation-enhanced diffusion and Ostwald coarsening. Void swelling correlated to heat input; being more prevalent in the welded conditions occurring at lower doses and in higher values relative to the base material. The low values of swelling despite microstructure coarsening caused by welding demonstrate the excellent swelling resistance of MA956, even after welding with the highest swelling values of 0.5% noted in the stir zone high heat input condition at 450 °C, 200 dpa. The dislocation behavior was inconsistent: the strongest trend was that network density was higher for welded versus base material, and an increase in loop diameter with temperature was observed. A rate theory model based on the observed microstructure suggests at high temperature interstitial loss to sinks was more likely to be dominant compared to mutual annihilation via point defect recombination, because of an increase of the radiation diffusion coefficient with temperature regardless of initial welded microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Bridging adsorption behavior of confined CH 4 -CO 2 binary mixtures across scales

An accurate understanding of the competitive adsorption of CH 4 -CO 2 binary mixtures in nano-confined systems is critical for engineering CO 2 storage in shale gas reservoirs. Due to difficulties in making reliable experimental observations in nano-scale, atomistic simulations (ASs), such as the Grand Canonical Monte Carlo (GCMC) method, provide a viable approach to studying the adsorption behavior of confined fluids. ASs are, however, limited in the size of the compositional domain due to the high computational cost. This work proposes a framework that combines AS and the lattice Boltzmann (LB) method to bridge the physics of confined fluids across scales. The Peng–Robinson equation of state (PR-EoS) produces fugacity coefficients, which serve as input for conducting multi-component GCMC simulations. These GCMC simulations explore the competitive adsorption behavior of CH 4 -CO 2 in nano-slits at various composition, pressure, and channel-width conditions. Both components generate adsorption layers with high densities near the walls with CO 2 preferentially adsorbing compared to CH 4 on the organic walls of carbon sheets. At the mesoscale, a pseudopotential model represents the intermolecular forces in multi-component, multiple-relaxation-time LB simulations. The LB simulations are in good agreement with the GCMC results, allowing us to obtain values for tunable LB parameters. We then extend the use of LB to simulate adsorption behavior in complex networks with nano-sized channels. The phase behavior and fluid properties in the complex geometries of nano-channels differ from nano-slits and bulk systems. Furthermore, the bridging of physics from GCMC (microscale) to LB (mesoscale) via the macroscale PR-EoS connects the adsorption behavior of binary systems across scales.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Single-shot femtosecond laser-induced breakdown spectroscopy in low-density gases

Single-shot fs laser-induced breakdown spectroscopy (LIBS) has been used to interrogate low-density gases of H 2 ⁠, CH 4 ⁠, CO 2 ⁠, and H 2 O ranging from 0.5 × 10 15 to 9.1 × 10 15 cm −3 (15–275 mTorr). Strong spectral signatures were observed above the noise baseline for each test gas at all densities considered. Unique molecular signatures were observed for both CH 4 and CO 2 along with atomic features for each gas. The LIBS intensity response with respect to density was determined for each species using Voigt lineshape fitting or numerical integration. Finally, the collection and analysis of LIBS spectra of composite mixtures of CH 4 –CO 2 and CH 4 –H 2 demonstrates that the extraction of individual densities from gas mixtures is feasible at low densities.

Artificial neural networks↗
DOE OSTI2025

Insights from Femtosecond Transient Absorption Spectroscopy into the Structure–Function Relationship of Glyceline Deep Eutectic Solvents

This study aimed to determine the structure–function relationship (SFR) for ChCl–glycerol mixtures, a deep eutectic solvent (DES), by investigating their microscopic solvation dynamics and how it relates to their macroscopic properties across varying concentrations of ChCl. Femtosecond transient absorption (fs-TA) spectroscopy revealed two distinct solvation dynamics time constants: τ 1 , governed by glycerol–glycerol interactions, and τ 2 , dominated by the choline response. The τ 2 minimum at 25–30 mol % ChCl closely aligned with the eutectic composition (~33.33 mol % ChCl), where the glycerol network was the most organized and the choline ions exhibited the fastest relaxation. The viscosity decreased sharply up to ~25 mol % ChCl and then plateaued, while the conductivity increased monotonically with ChCl concentration, reflecting enhanced ionic mobility. The density decreased with both increasing ChCl concentration and temperature, indicating disrupted hydrogen bonding and reduced molecular packing. The polarity, measured using betaine-30 (B30) and the E T (30) polarity scale, increased steeply up to approximately 25 mol % ChCl before reaching a plateau. These findings identified the eutectic composition as the optimal concentration range for balancing stability, fluidity, conductivity, and enhanced dynamics within the glycerol system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Characterization of partial wetting by CMAS droplets using multiphase many-body dissipative particle dynamics and data-driven discovery based on PINNs

The molten sand that is a mixture of calcia, magnesia, alumina and silicate, known as CMAS, is characterized by its high viscosity, density and surface tension. The unique properties of CMAS make it a challenging material to deal with in high-temperature applications, requiring innovative solutions and materials to prevent its buildup and damage to critical equipment. Here, we use multiphase many-body dissipative particle dynamics simulations to study the wetting dynamics of highly viscous molten CMAS droplets. The simulations are performed in three dimensions, with varying initial droplet sizes and equilibrium contact angles. We propose a parametric ordinary differential equation (ODE) that captures the spreading radius behaviour of the CMAS droplets. The ODE parameters are then identified based on the physics-informed neural network (PINN) framework. Subsequently, the closed-form dependency of parameter values found by the PINN on the initial radii and contact angles are given using symbolic regression. Finally, we employ Bayesian PINNs (B-PINNs) to assess and quantify the uncertainty associated with the discovered parameters. In brief, this study provides insight into spreading dynamics of CMAS droplets by fusing simple parametric ODE modelling and state-of-the-art machine-learning techniques.

Mechanics↗
NASA NTRS2016

Generalized Fluid System Simulation Program, Version 6.0

The Generalized Fluid System Simulation Program (GFSSP) is a general purpose computer program for analyzing steady state and time-dependent flow rates, pressures, temperatures, and concentrations in a complex flow network. The program is capable of modeling real fluids with phase changes, compressibility, mixture thermodynamics, conjugate heat transfer between solid and fluid, fluid transients, pumps, compressors, and external body forces such as gravity and centrifugal. The thermofluid system to be analyzed is discretized into nodes, branches, and conductors. The scalar properties such as pressure, temperature, and concentrations are calculated at nodes. Mass flow rates and heat transfer rates are computed in branches and conductors. The graphical user interface allows users to build their models using the 'point, drag, and click' method; the users can also run their models and post-process the results in the same environment. Two thermodynamic property programs (GASP/WASP and GASPAK) provide required thermodynamic and thermophysical properties for 36 fluids: helium, methane, neon, nitrogen, carbon monoxide, oxygen, argon, carbon dioxide, fluorine, hydrogen, parahydrogen, water, kerosene (RP-1), isobutene, butane, deuterium, ethane, ethylene, hydrogen sulfide, krypton, propane, xenon, R-11, R-12, R-22, R-32, R-123, R-124, R-125, R-134A, R-152A, nitrogen trifluoride, ammonia, hydrogen peroxide, and air. The program also provides the options of using any incompressible fluid with constant density and viscosity or ideal gas. The users can also supply property tables for fluids that are not in the library. Twenty-four different resistance/source options are provided for modeling momentum sources or sinks in the branches. These options include pipe flow, flow through a restriction, noncircular duct, pipe flow with entrance and/or exit losses, thin sharp orifice, thick orifice, square edge reduction, square edge expansion, rotating annular duct, rotating radial duct, labyrinth seal, parallel plates, common fittings and valves, pump characteristics, pump power, valve with a given loss coefficient, Joule-Thompson device, control valve, heat exchanger core, parallel tube, and compressible orifice. The program has the provision of including additional resistance options through User Subroutines. GFSSP employs a finite volume formulation of mass, momentum, and energy conservation equations in conjunction with the thermodynamic equations of state for real fluids as well as energy conservation equations for the solid. The system of equations describing the fluid network is solved by a hybrid numerical method that is a combination of the Newton-Raphson and successive substitution methods. The application and verification of the code has been demonstrated through 30 example problems.

Majumdar, A. K.↗
NASA NTRS1993

Intelligent neuroprocessors for in-situ launch vehicle propulsion systems health management

Efficacy of existing on-board propulsion systems health management systems (HMS) are severely impacted by computational limitations (e.g., low sampling rates); paradigmatic limitations (e.g., low-fidelity logic/parameter redlining only, false alarms due to noisy/corrupted sensor signatures, preprogrammed diagnostics only); and telemetry bandwidth limitations on space/ground interactions. Ultra-compact/light, adaptive neural networks with massively parallel, asynchronous, fast reconfigurable and fault-tolerant information processing properties have already demonstrated significant potential for inflight diagnostic analyses and resource allocation with reduced ground dependence. In particular, they can automatically exploit correlation effects across multiple sensor streams (plume analyzer, flow meters, vibration detectors, etc.) so as to detect anomaly signatures that cannot be determined from the exploitation of single sensor. Furthermore, neural networks have already demonstrated the potential for impacting real-time fault recovery in vehicle subsystems by adaptively regulating combustion mixture/power subsystems and optimizing resource utilization under degraded conditions. A class of high-performance neuroprocessors, developed at JPL, that have demonstrated potential for next-generation HMS for a family of space transportation vehicles envisioned for the next few decades, including HLLV, NLS, and space shuttle is presented. Of fundamental interest are intelligent neuroprocessors for real-time plume analysis, optimizing combustion mixture-ratio, and feedback to hydraulic, pneumatic control systems. This class includes concurrently asynchronous reprogrammable, nonvolatile, analog neural processors with high speed, high bandwidth electronic/optical I/O interfaced, with special emphasis on NASA's unique requirements in terms of performance, reliability, ultra-high density ultra-compactness, ultra-light weight devices, radiation hardened devices, power stringency, and long life terms.

Gulati, S.↗
NASA NTRS2012

Particle Filtering for Model-Based Anomaly Detection in Sensor Networks

A novel technique has been developed for anomaly detection of rocket engine test stand (RETS) data. The objective was to develop a system that postprocesses a csv file containing the sensor readings and activities (time-series) from a rocket engine test, and detects any anomalies that might have occurred during the test. The output consists of the names of the sensors that show anomalous behavior, and the start and end time of each anomaly. In order to reduce the involvement of domain experts significantly, several data-driven approaches have been proposed where models are automatically acquired from the data, thus bypassing the cost and effort of building system models. Many supervised learning methods can efficiently learn operational and fault models, given large amounts of both nominal and fault data. However, for domains such as RETS data, the amount of anomalous data that is actually available is relatively small, making most supervised learning methods rather ineffective, and in general met with limited success in anomaly detection. The fundamental problem with existing approaches is that they assume that the data are iid, i.e., independent and identically distributed, which is violated in typical RETS data. None of these techniques naturally exploit the temporal information inherent in time series data from the sensor networks. There are correlations among the sensor readings, not only at the same time, but also across time. However, these approaches have not explicitly identified and exploited such correlations. Given these limitations of model-free methods, there has been renewed interest in model-based methods, specifically graphical methods that explicitly reason temporally. The Gaussian Mixture Model (GMM) in a Linear Dynamic System approach assumes that the multi-dimensional test data is a mixture of multi-variate Gaussians, and fits a given number of Gaussian clusters with the help of the wellknown Expectation Maximization (EM) algorithm. The parameters thus learned are used for calculating the joint distribution of the observations. However, this GMM assumption is essentially an approximation and signals the potential viability of non-parametric density estimators. This is the key idea underlying the new approach.

Solano, Wanda↗
NASA NTRS1998

A Generalized Fluid System Simulation Program to Model Flow Distribution in Fluid Networks

This paper describes a general purpose computer program for analyzing steady state and transient flow in a complex network. The program is capable of modeling phase changes, compressibility, mixture thermodynamics and external body forces such as gravity and centrifugal. The program's preprocessor allows the user to interactively develop a fluid network simulation consisting of nodes and branches. Mass, energy and specie conservation equations are solved at the nodes; the momentum conservation equations are solved in the branches. The program contains subroutines for computing "real fluid" thermodynamic and thermophysical properties for 33 fluids. The fluids are: helium, methane, neon, nitrogen, carbon monoxide, oxygen, argon, carbon dioxide, fluorine, hydrogen, parahydrogen, water, kerosene (RP-1), isobutane, butane, deuterium, ethane, ethylene, hydrogen sulfide, krypton, propane, xenon, R-11, R-12, R-22, R-32, R-123, R-124, R-125, R-134A, R-152A, nitrogen trifluoride and ammonia. The program also provides the options of using any incompressible fluid with constant density and viscosity or ideal gas. Seventeen different resistance/source options are provided for modeling momentum sources or sinks in the branches. These options include: pipe flow, flow through a restriction, non-circular duct, pipe flow with entrance and/or exit losses, thin sharp orifice, thick orifice, square edge reduction, square edge expansion, rotating annular duct, rotating radial duct, labyrinth seal, parallel plates, common fittings and valves, pump characteristics, pump power, valve with a given loss coefficient, and a Joule-Thompson device. The system of equations describing the fluid network is solved by a hybrid numerical method that is a combination of the Newton-Raphson and successive substitution methods. This paper also illustrates the application and verification of the code by comparison with Hardy Cross method for steady state flow and analytical solution for unsteady flow.

Majumdar, Alok↗
DOE OSTI2025

The Triple Catalytic Action of Tertiary Nitrogen Catalysts in Recyclable Epoxy-Anhydride Thermosets

The thermosetting polymer matrix in fiber reinforced composites is an important component for energy related applications, such as the lightweighting of vehicles or their use in wind and waterpower turbine blades, due to their ability to provide superior adhesion, stiffness, and applicability to a wide range of manufacturing processes. Despite these benefits, today's thermosets are widely considered to be unrecyclable; thus, there is a large interest in redesigning these materials to be inherently recyclable so that energy intensive production of fibers and monomers can be circumvented, bolstering composite manufacture supply chains. Polyester covalent adaptable networks (PECANs) are one such promising alternative to the incumbent, nonrecyclable epoxy-amine thermosets. PECANs can be formed from the ring-opening co-polymerization (ROCOP) of epoxy-anhydride monomer mixtures and subsequent curing at mild temperatures to exhibit similar performance to conventional epoxies while also possessing unique dynamic chemistries along the ester-hydroxyl backbone that are capable of transesterification and thus reprocessability. While significant advancements have been made in formulating these materials for improved mechanical properties or optimizing solvolysis and reprocessing strategies, less attention has been placed on the impact of the residing amine catalyst used to generate the polyester network. In this work, we evaluated the triple-catalytic efficacy of 12 tertiary amines that act as a curing (bulk ROCOP), a transesterification (internal bond exchange), and a deconstruction (methanolysis) catalyst for PECAN thermosets. Specifically, we first distinguish between chain-growth and step-growth polymerization mechanisms for epoxy-amine and epoxy-anhydride mechanisms. We also utilized density functional theory (DFT) to estimate the basicity (pKb) of each catalyst. Of the tested catalysts, the ROCOP of the studied PECAN network can be completed between 95 and 247 min (at 80 degrees C), with variable gelation phenomena. Additionally, the stress relaxation (transesterification metric) efficiency of the tested PECAN networks with alternative embedded catalysts ranged from 95% to 15% reduction in stress after 5 h at 200 degrees C, and the depolymerization efficacy ranged from 2.5% to 9.8% deconstruction after 36 h at 130 degrees C. Overall, the nitrogen-based moieties were demonstrated to influence polymerization kinetics, catalyze the dynamic transesterification exchange mechanism, and aid in the solvolysis of the thermosets at end-of-life.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Tuning Cation (Dis)Order in Cr‐Based Li‐Excess Oxide Cathode Materials to Improve Li+ Transport Properties

Abstract Li‐excess disordered rocksalts (DRXs) hold promise as next‐generation cathodes for Li‐ion batteries due to their high capacity and energy density, along with the potential to eliminate the need for Co and Ni. However, due to their low Li + diffusivity, DRXs need to be pulverized into nanoparticles to achieve high performance. Herein, a new strategy for overcoming this limitation is demonstrated, involving the design of as‐synthesized partially disordered oxides with a structure that lies between ordered layered and fully disordered, exhibiting a varying degree of local (dis)order. This unique structure activates new Li + diffusion channels, improving percolation and transport properties. This strategy allows a large content of Li + to be accessed in material with large, micron‐sized particles through highly reversible Cr 3+/6+ and O redox, yielding a first discharge capacity of 286 mAh g −1 (881 Wh kg −1 ). The Li + percolation network is further improved by substituting Ti with a mixture of multiple metals, which appears to locally decrease the migration barrier through lattice distortion. Proper tuning of the chemical composition, especially the content of metals with empty d orbitals, is established as a crucial factor for controlling the degree of disorder and mitigating voltage fade and hysteresis growth upon cycling.

Moździerz, Maciej↗
DOE OSTI2020

Nanoscale Mg–B via Surfactant Ball Milling of MgB 2 : Morphology, Composition, and Improved Hydrogen Storage Properties

Metal borides have attracted the attention of researchers due to their useful physical properties and unique ability to form high hydrogen-capacity metal borohydrides. We demonstrate improved hydrogen storage properties of a nanoscale Mg–B material made by surfactant ball milling MgB 2 in a mixture of heptane, oleic acid, and oleylamine. Transmission electron microscopy data show that Mg–B nanoplatelets are produced with sizes ranging from 5 to 50 nm, which agglomerate upon ethanol washing to produce an agglomerated nanoscale Mg–B material of micron-sized particles with some surfactant still remaining. X-ray diffraction measurements reveal a two-component material where 32% of the solid is a strained crystalline solid maintaining the hexagonal structure with the remainder being amorphous. Fourier transform infrared shows that the oleate binds in a “bridge-bonding” fashion preferentially to magnesium rather than boron, which is confirmed by density functional theory calculations. The Mg–B nanoscale material is deficient in boron relative to bulk MgB 2 with a Mg–B ratio of ~1:0.75. The nanoscale MgB 0.75 material has a disrupted B–B ring network as indicated by X-ray absorption measurements. Hydrogenation experiments at 700 bar and 280 °C show that it partially hydrogenates at temperatures 100 °C below the threshold for bulk MgB 2 hydrogenation. In addition, upon heating to 200 °C, the H–H bond-breaking ability increases ~10-fold according to hydrogen–deuterium exchange experiments due to desorption of oleate at the surface. This behavior would make the nanoscale Mg–B material useful as an additive where rapid H–H bond breaking is needed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

KCl-UCl 3 molten salts investigated by Ab Initio Molecular Dynamics (AIMD) simulations

Ab Initio Molecular Dynamics (AIMD) simulations are performed on molten KCl-UCl3 salt mixtures to determine energies, heat capacities, and densities. The density-dependent energy correction (DFT-dDsC), Grimme et al.'s DFT-D3, and Langreth & Lundqvist (vdW-cx) models are used for dispersion forces and combined with the Perdew-Burke-Ernzerhof (PBE) exchange-correlation potential with a Hubbard U parameter for the 5f electrons of uranium. After validating predictions for the end-member systems to literature data, KCl-UCl 3 mixtures are studied at select temperatures. Densities and energies both deviate from ideal solution behavior, with the maximum deviation occurring around 36% UCl 3 for mixing energies and slightly lower (29% UCl3) for densities. Compared to the NaCl-UCl 3 system, which was previously investigated using the same simulation methodologies, the KCl-UCl 3 density and mixing energy deviations from ideal solution behavior are larger by almost a factor of two. No deviation from ideal solution behavior for heat capacity was observed. The AIMD predictions for mixing energies and densities agree qualitatively with experimental data, though the spread in data obtained from the various dispersion force models utilized, measurements, and empirical estimates makes strong conclusions difficult. The dependence of thermodynamic and thermophysical properties on composition is correlated with the local chemistry of the solution phase, in particular, the tendency of UCl 3 to form network structures.

36 MATERIALS SCIENCE↗
DOE OSTI2024

KCl-UCl 3 molten salts investigated by Ab Initio Molecular Dynamics (AIMD) simulations: A comparative study with three dispersion models

Ab Initio Molecular Dynamics (AIMD) simulations are performed on molten KCl-UCl 3 salt mixtures to determine energies, heat capacities, and densities. The density-dependent energy correction (DFT-dDsC), Grimme et al.’s DFT-D3, and Langreth & Lundqvist (vdW-cx) models are used for dispersion forces and combined with the Perdew-Burke-Ernzerhof (PBE) exchange-correlation potential with a Hubbard U param eter for the 5$f$ electrons of uranium. After validating predictions for the end-member systems to literature data, KCl-UCl 3 mixtures are studied at select temperatures. Densities and energies both deviate from ideal solution behavior, with the maximum deviation occurring around 36% UCl 3 for mixing energies and slightly lower (29% UCl 3 ) for densities. Compared to the NaCl-UCl 3 system, which was previously investigated using the same simulation methodologies, the KCl-UCl 3 density and mixing energy deviations from ideal solution behavior are larger by almost a factor of two. No deviation from ideal solution behavior for heat capacity was observed. The AIMD predictions for mixing energies and densities agree qualitatively with experimental data, though the spread in data obtained from the various dispersion force models utilized, measurements, and empirical estimates makes strong conclusions difficult. The dependence of thermodynamic and thermophysical properties on composition is correlated with the local chemistry of the solution phase, in particular, the tendency of UCl 3 to form network structures.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Direct observation of key aluminum hydroxide prenucleation oligomers for gibbsite nucleation and crystallization in sodium aluminate solution by liquid ToF-SIMS

The mechanism of gibbsite (aluminum hydroxide) crystallization from highly alkaline solutions such as Bayer liquors remains poorly understood, where aluminum (Al) transforms from largely tetrahedrally coordinated aluminate monomers in sodium aluminate solutions into a network of octahedra in gibbsite crystals. A variety of traditional analytical approaches applied to this system do not readily reveal the presence of higher-order oligomeric intermediates. To overcome this limitation, we employed in-situ liquid Time-of-Flight Secondary Ion Mass Spectrometry (ToF-SIMS) to examine the Al species present in concentrated sodium aluminate solutions favoring crystallization of gibbsite or sodium aluminate. A complex mixture of Al oligomers was found. By comparing the change in the relative concentration of Al oligomers with +1 and -1 charge, we were able to identify three major Al oligomer candidates, iso-tetramers, iso-pentamers, and cyclic-hexamers, for the nucleation and crystallization of gibbsite. The concentrations of iso-tetramers and iso-pentamers significantly surpass those of cyclic-hexamers. Time-dependent in-situ Raman spectroscopy analysis indicated that the appearance of gibbsite coincided with the peak concentration of these oligomers. The Density-functional theory (DFT) calculation suggests that the formation of iso-oligomers is more favorable than that of cyclic-hexamers. The combined results suggest that iso-tetramers and iso-pentamers play the most substantial role in the nucleation and growth of gibbsite in the sodium aluminate solutions. Our findings also suggest that the oligomers that promote gibbsite crystallization are more stable in dilute sodium aluminate solutions, making these solutions particularly suitable for efficient gibbsite crystallization. In conclusion, our study fills a major knowledge gap in understanding Al speciation that leads to the nucleation and crystallization of gibbsite in concentrated sodium aluminate solutions.

Bayer liquor↗
DOE OSTI2025

Structure of molten ytterbium aluminum garnet

Rare earth aluminum garnets are important materials in optical, dielectric, and thermal barrier applications. To advance the understanding of their melt processing and glass forming ability, we report the atomic structure of molten Yb 3 Al 5 O 12 over 1770–2630 K, which spans the equilibrium and supercooled liquid regimes. The melt density at T m = 2283 K is 5.50 g cm –3 , measured via silhouette imaging of electrostatically levitated drops over 1010–2420 K. Four separate structure measurements were made with aerodynamically levitated melts using x-ray and neutron diffraction with isotope substitution of Yb ( 172 Yb, 174 Yb, or nat Yb). Empirical potential structure refinement models were developed, which are in excellent agreement with the experiments. Coordination environments for Al–O are predominantly 4- and 5-coordinate, with a mean coordination of n AlO = 4.43(8), while Yb–O environments mostly range from 5- to 8-coordinate, with n YbO = 6.26(8). The cation–oxygen polyhedra are connected primarily by corner-sharing, with edge-sharing constituting up to ~1/3 of the connectivity among polyhedra with Yb or higher-coordinated Al–O. Structurally, the –Al–O– network in molten Yb 3 Al 5 O 12 appears conducive to glass formation: n OAl = 1.85(3), there are 1.86 AlO x –AlO x connections per Al atom (e.g., a mixture of Q 3 and Q 4 units), and the modal ring size is six cations. These characterize a network that is somewhat less constrained compared to SiO 2 glass, yet Yb 3 Al 5 O 12 cannot be quenched into crystal-free glass. Here, aluminum garnet compositions with larger rare earth cations do form glass, so these characterizations help reveal the structural characteristics corresponding to the limit of glass forming ability in rare earth aluminates.

36 MATERIALS SCIENCE↗