Search NASA⌕ Search

SEARCH · Search NASA

Results for “Mixing thermodynamics”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Thermodynamics of electrolyte solutions near charged surfaces: Constant surface charge vs constant surface potential

Electric double layers are ubiquitous in science and engineering and are of current interest, owing to their applications in the stabilization of colloidal suspensions and as supercapacitors. While the structure and properties of electric double layers in electrolyte solutions near a charged surface are well characterized, there are subtleties in calculating thermodynamic properties from the free energy of a system with charged surfaces. These subtleties arise from the difference in the free energy between systems with constant surface charge and constant surface potential. In this work, we present a systematic, pedagogical framework to properly account for the different specifications on charged bodies in electrolyte solutions. Our approach is fully variational—that is, all free energies, boundary conditions, relevant electrostatic equations, and thermodynamic quantities are systematically derived using variational principles of thermodynamics. We illustrate our approach by considering a simple electrolyte solution between two charged surfaces using the Poisson–Boltzmann theory. Our results highlight the importance of using the proper thermodynamic potential and provide a general framework for calculating thermodynamic properties of electrolyte solutions near charged surfaces. Specifically, we present the calculation of the pressure and the surface tension between two charged surfaces for different boundary conditions, including mixed boundary conditions.

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

Development of Low Cost, Robust and Durable Cathode Materials to Support Solid Oxide Fuel Cell Commercialization (Final Technical Report)

The project is aimed at addressing critical issues related to chromia-poisoning in the cathode of solid oxide fuel cells. A technical hurdle to the commercialization of solid oxide fuel cells is that chromia sources from the metallic interconnect and the balance of plant (BoP) components will cause chromia-poisoning to the cathode performance. This presents a significant challenge since it is very difficult to eliminate chromia sources in the cathode environment. One aspect of the current project is to evaluate whether SFM (Sr 2 Fe 1.5 Mo 0.5 O 6-δ ) that we have recently developed as redox-stable mixed conducting oxide possesses tolerance to chromia-poisoning. The second focus of this project is to evaluate coating on the cathode materials to mitigate chromia poisoning effect. Both thermodynamic evaluation and experimental approaches have been applied to enhance the knowledge base and fundamental understanding of chromia-tolerance of cathode for solid oxide fuel cells. low cost, robust and durable cathode material to support SOFC commercialization.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2025

Bridging the Gap in Carbon Free Iron Making: How Hydrogen Affects the Reduction of Iron Ore between 900 and 1590 °C

Hydrogen-based reduction of iron ore for iron and steel production has emerged as a promising alternative to coal and natural gas. Unlike other hydrogen-based iron ore reduction studies, this research focuses on a wide temperature range across 900–1590 °C, encompassing reduction in solid, mixed, and liquid (slag) phases. For a 20 min exposure to hydrogen, the reduction degree increased monotonically from ∼35% at 900 °C to >90% at 1550 °C, except between 1100 °C and 1400 °C, where it stagnated ∼60%. This experimental work challenges the widely accepted notion that higher temperatures enhance the reduction process. Instead, it reveals an overlooked kinetic bottleneck, suggesting complex thermodynamic and mass transfer limitations influenced by phase transformations, diffusion barriers, and microstructural changes. Density functional theory-based molecular dynamics simulations indicate that oxygen diffusivity in BCC iron is 3.88 × 10 –5 cm 2 /s which is ∼5–10 times higher than that in FCC iron. This study reports that this stagnant reduction degree in the mixed solid–liquid phase is due to competition of multiple mechanisms, such as surface- and bulk-diffusion, pore collapse mechanisms, and crystallographic transitions.

08 HYDROGEN↗
DOE OSTI2021

An automated multi-component gas adsorption system (MC GAS)

The knowledge gap on adsorption of complex mixtures in the literature relative to single component data represents a persistent obstacle to developing accurate process models for adsorption separations. The collection of mixed gas adsorption data is an imminent need for improved understanding of the behavior of adsorbent systems in these diverse adsorption applications. Current approaches to understanding mixture adsorption using predictive theories based on pure component adsorption experiments often fail to capture the behavior of more complex, non-ideal systems. In this work, we present an automated volumetric instrument for the measurement of mixed gas adsorption isotherms. This instrument was validated by comparison to other in-house instruments and data available in the literature, and the binary adsorption measurements were found to be thermodynamically consistent. The automation of this instrument allows for rapid collection of high-quality mixture adsorption data.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Entrainment and aerosol effects on marine boundary-layer clouds: An investigation using ACE-ENA data from HOLODEC, G1, Pico and ACTOS (Final Report on Project Activities)

Marine boundary-layer clouds cover large regions of the globe and are known to strongly influence radiative balances. The microphysical properties and persistence of these clouds are tightly coupled with cloud-top entrainment and aerosol properties within both the boundary layer and the overlying free troposphere. This work addressed the microphysical response to entrainment and aerosol properties in marine stratocumulus clouds, using data from the Aerosol and Cloud Experiments in the Eastern North Atlantic (ACE-ENA) project. Specifically, emphasis was placed on 1. extensive in situ measurements taken with the Atmospheric Radiation Measurement (ARM) G1 aircraft using the Holographic Detector for Clouds (Holodec) instrument; 2. airborne measurements obtained with the helicopter-borne Airborne Cloud Turbulence Observation System (ACTOS) operated by the Leibniz Institute for Tropospheric Research (TROPOS) during the first phase of ACE-ENA; and 3. above-boundary-layer, mountain-top measurements taken at the Pico Mountain Observatory (OMP) by the research group of the co-investigator and scientists from TROPOS during the first phase of ACE-ENA. The work is aligned with the topic “Warm Boundary-Layer Atmospheric Processes.” Key results include the following. Data from the Holographic Detector for Clouds (Holodec) were reprocessed, refined, and validated through a careful instrument intercomparison, resulting in a high-quality dataset available on the ARM archive for ACE-ENA. High-resolution ACTOS measurements of thermodynamics, microphysical, and turbulence properties were analyzed to explore the behavior of the entrainment velocity at cloud top. Analysis of the vertical variability of cloud droplet size distribution shape suggests that mixing is more inhomogeneous near cloud top, and more homogeneous deeper into the cloud. Analysis of aerosol measurements from OMP explored the mixing state and cloud condensation nucleus properties, as well as the implications for radiative forcing. Finally, a machine learning algorithm was developed for identifying characteristic cloud droplet size distributions, and was employed to explore microphysical regimes in stratocumulus clouds observed during ACE-ENA. The characteristic size distributions are narrow, and only when spatially averaged do they produce the broad “gamma” distributions typically assumed in models. This implies that precipitation development should account for variability and correlations in the distribution shape as well as the number concentration of cloud droplets.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Temperature-Dependent Configurational Entropy Calculations for Refractory High-Entropy Alloys

The cluster expansion formalism for alloys is used to construct surrogate models for three refractory high-entropy alloys (NbTiVZr, HfNbTaTiZr, and AlHfNbTaTiZr). These cluster expansion models are then used along with Monte Carlo methods and thermodynamic integration to calculate the configurational entropy of these refractory high-entropy alloys as a function of temperature. Many solid solution alloy design guidelines are based on the ideal entropy of mixing, which increases monotonically with N , the number of elements in the alloy. However, our results show that at low temperatures, the configurational entropy of these materials is largely independent of N , and the assumption described above only holds in the high-temperature limit. This suggests that alloy design guidelines based on the ideal entropy of mixing require further examination.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The effects of ACRT on melt undercooling during the growth of CZT via the traveling heater method: Ekman versus Taylor-Görtler flows

For the first time, a comprehensive model for the traveling heater method, including detailed species transport and phase-change thermodynamics, is applied to understand the effects of accelerated crucible rotation (ACRT) on the growth of cadmium zinc telluride. Here, a rotation schedule that is designed from classical considerations does not completely mix the melt, rather undercooled regions of localized melt are swept periodically across the growth interface. Even without melt homogenization, these ACRT-driven flows disrupt the persistent undercooling that occurs during THM growth without rotation. An ACRT schedule designed to promote Ekman flows results in lower levels of melt undercooling than a schedule designed to promote Taylor-Görtler flows and greater mixing.

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 OSTI2021

Evaluation of the CAM6 Climate Model Using Cloud Observations at McMurdo Station, Antarctica

A comparative analysis between observational data from McMurdo Station, Antarctica and the Community Atmosphere Model version 6 (CAM6) simulation is performed focusing on cloud characteristics and their thermodynamic conditions. Ka-band Zenith Radar (KAZR) and High Spectral Resolution Lidar (HSRL) retrievals are used as the basis of cloud fraction and cloud phase identifications. Radiosondes released at 12-hour increments provide atmospheric profiles for evaluating the simulated thermodynamic conditions. Our findings show that the CAM6 simulation consistently overestimates (underestimates) cloud fraction above (below) 3 km in four seasons of a year. Normalized by total in-cloud samples, ice and mixed phase occurrence frequencies are underestimated and liquid phase frequency is overestimated by the model at cloud fractions above 0.6, while at cloud fractions below 0.6 ice phase frequency is overestimated and liquid-containing phase frequency is underestimated by the model. The cloud fraction biases are closely associated with concurrent biases in relative humidity (RH), i.e., high (low) RH biases above (below) 2 km. Frequencies of correctly simulating ice and liquid-containing phase increase when the absolute biases of RH decrease. Cloud fraction biases also show a positive correlation with RH biases. Water vapor mixing ratio biases are the primary contributor to RH biases, and hence, likely a key factor controlling the cloud biases. Finally, this diagnosis of the evident shortfalls of representations of cloud characteristics in CAM6 simulation at McMurdo Station brings new insight in improving the governing model physics therein.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Charge state distributions in dense plasmas

Charge state distributions in hot, dense plasmas are a key ingredient in the calculation of spectral quantities like the opacity. However, they are challenging to calculate, as models like Saha–Boltzmann become unreliable for dense, quantum plasmas. Here, we present a new variational model for the charge state distribution, along with a simple model for the energy of the configurations that includes the orbital relaxation effect. Comparison with other methods reveals generally good agreement with average atom-based calculations, the breakdown of the Saha–Boltzmann method, and mixed agreement with a chemical model. We conclude that the new model gives a relatively inexpensive, but reasonably high fidelity method of calculating the charge state distribution in hot dense plasmas, in local thermodynamic equilibrium.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 OSTI2022

Atomic Mixing of Initially Separated Mixtures

We consider the atomic mixing of two materials, initially separated in space, at fixed temperature– pressure conditions [T, p], and zero mass-mean velocity. Each material is composed of a specified set of atomic species whose density corresponds to the specieswise real Equation of State (EOS) $\tilde{&#961}$ i [T, p]. Because there are no spatial gradients in the mean [T, p] state, the only mechanism for mixing is a specieswise nonequilibrium thermal velocity distribution. The recipe for nonequilibrium species mass fluxes is taken directly from Hirschfelder, Curtiss, and Bird, who express the diffusion velocities in terms of gradients in the molar concentrations and pairwise coefficients of diffusion. EOS data from the SESAME library is used for all species. Data for the specieswise self diffusion coefficients is adopted from the Orbital Free Molecular Dynamics (OFMD) simulations of Ticknor et al.; the so-called Darken approximation is used to express the mutual diffusion coefficients of a mixture in the manner used by White et al. A one-dimensional spatial domain is used; with the left half occupied by one multispecies material and the right half occupied by the other, and with zero–gradient boundary conditions at the ends of the domain. Hence the total mass in the domain is a constant. At time t = 0 the materials are separated. Accordingly, spatial gradients at the separation point are infinite and the species fluxes are also infinite. As t → $\textbf{∞}$ the spatial gradients approach zero and a uniform set of species concentrations is obtained; at this point the system is said to be at “thermodynamic equilibrium.” Analytic solution of the problem is not possible in general due to the nonlinear concentration dependence of the species fluxes. Hence we utilize the numerical solution apparatus in Mathematica for obtaining accurate time-dependent solutions for the mixing of various initial species masses for conditions of interest in laser-driven Inertial Confinement Fusion. Concentration profiles and the speed of propagation of the mixing fronts are discussed in light of the nonlinearity of the problem.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Self-healing Growth of LaNiO 3 on a Mixed-Terminated Perovskite Surface

Developing atomic-scale synthesis control is a prerequisite for understanding and engineering the exotic physics inherent to transition-metal oxide heterostructures. Thus, far, however, the number of materials systems explored has been extremely limited, particularly with regard to the crystalline substrate, which is routinely SrTiO 3 . Here, we investigate the growth of a rare-earth nickelate-LaNiO 3 -on (LaAlO 3 )(Sr 2 AlTaO 6 ) (LSAT) (001) by oxide molecular beam epitaxy (MBE). Whereas the LSAT substrates are smooth, they do not exhibit the single surface termination usually assumed necessary for control over the interface structure. Performing both nonresonant and resonant anomalous in situ synchrotron surface X-ray scattering during MBE growth, we show that reproducible heterostructures can be achieved regardless of both the mixed surface termination and the layer-by-layer deposition sequence. The rearrangement of the layers occurs dynamically during growth, resulting in the fabrication of high-quality LaNiO 3 /LSAT heterostructures with a sharp and consistent interfacial structure. This is due to the thermodynamics of the deposition window as well as the nature of the chemical species at interfaces-here, the flexible charge state of nickel at the oxide surface. This has important implications regarding the use of a wider variety of substrates for fundamental studies on complex oxide synthesis.

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

Application of Scalar Filtered Density Function to Turbulent Flows Under Supercritical Condition

Abstract The scalar filtered density function (FDF) methodology is extended and employed for large eddy simulation (LES) of turbulent flows under supercritical condition. To describe real fluid behavior, the extended methodology incorporates the generalized heat and mass diffusion models along with real fluid thermodynamic relations which are derived using the cubic Peng–Robinson equation of state. These models are implemented within the stochastic differential equations comprising the scalar FDF transport. Simulations are conducted of a temporally developing mixing layer under supercritical condition and the results are assessed by comparing with data generated by direct numerical simulation (DNS) of the same layer. The consistency of the proposed FDF methodology is assessed. The LES-FDF predictions are shown to agree favorably with the DNS data and exhibit several key features pertaining to supercritical turbulent flows.

Energy & Fuels↗
DOE OSTI2022

Effect of disorder on thermodynamic instability of binary Rare-earth – Nickel – Palladium compounds

In this work, we have investigated the thermodynamic stability of disordered rare-earth phases SmX 2 and Sm 10 X 21 (X=Ni, Pd) using machine-learning based analytical descriptor and first-principles density functional theory methods. The absence of Laves phase compounds in R-Pd binary systems is a longstanding problem of rare earth science: even though Ni and Pd belong to the same group of the periodic table and have similar electronic structure, the Pd compound crystallizes in a monoclinic (C2/m) phase with 10:21 stoichiometry, i.e., Sm 10 Pd 21 , while the Ni compound adopts a cubic Laves phase (MgCu 2 ) structure. To understand this contrasting phase stability, we performed thermodynamic convex hull analysis of Sm x Ni 1-x and Sm x Pd 1-x binary systems, which is experimentally validated using powder X-ray diffraction (PXRD) analyzes of polycrystalline Sm(Ni x Pd 1-x ) 2 samples with x=0, 0.5, and 1. A detailed electronic-structure (band-structure, charge density, and Fermi-surface) analysis of the differences between SmNi 2 /SmPd 2 and Sm 10 Ni 21 /Sm 10 Pd 21 compounds provides the quantum mechanical origin of the unfavorable mixing of Pd with Ni in cubic Laves phase. We show that the stability of Sm-Pd in 10:21 stoichiometry arises from improved intra-/inter-layer 5d-4d bonding compared to the 1:2 stoichiometry. Our work emphasizes the importance of ab-initio methods and computationally inexpensive analytical descriptors for the detailed analysis of thermodynamic and electronic properties of hard-to-prepare rare-earth compounds.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Dispatch optimization of a concentrating solar power system under uncertain solar irradiance and energy prices

The integration of thermal energy storage into a concentrating solar power system allows for mitigating some of the risk associated with uncertain solar irradiance and uncertain energy prices. We solve a 48 h dispatch optimization model with continually updated conditional point forecasts of both direct normal irradiance (DNI) and electricity prices with a rolling-horizon scheme at hourly resolution over the course of a year. Joint, conditional forecasts for DNI and prices are formed using an autoregressive moving-average time series model with exogenous weather predictors. We guide dispatch using a mixed-integer programming model, but in order to evaluate performance we use the System Advisor Model (SAM) of the National Renewable Energy Laboratory. SAM is a techno-economic simulation model that accounts for plant thermodynamics with higher fidelity. Our conditional DNI forecasts improve annual revenue by 4%–12% over using historical forecasts based on data from previous years. Conditional price forecasts improve annual revenue by 6%–19% in the real-time market over analogous historical forecasts. Updating these forecasts every six hours, rather than every 24 h, further improves annual revenue by 5%–6%. Here, we also investigate a method that values terminal inventory in our dispatch optimization model, again when used in a rolling-horizon scheme.

14 SOLAR ENERGY↗
DOE OSTI2022

Liquid and Gamma-BCC Phase Equilibria in the Uranium-Niobium System: A Simple Thermodynamic Model for Phase Diagram and Other Predictions

A simple thermodynamic model for the uranium-niobium (U-Nb) system was created using a Calphad-type approach, circa 2002. Only the liquid and gamma-BCC phases were modeled, with their subregular solution mixing terms optimized by trial and error. The custom (user-defined) database file (.tdb) and experimental data file (.exp) are listed, both of which are compatible with Thermo-Calc software. This model is used to calculate the solidus, liquidus, miscibility gap on the U-Nb phase diagram and also the enthalpy of melting in the U-6Nb alloy.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗