Search NASA⌕ Search

SEARCH · Search NASA

Results for “Solution 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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Thermodynamics and kinetics of solution combustion synthesis: Ni(NO 3 ) 2 + fuels systems

Solution combustion synthesis (SCS) utilizes exothermic self-propagating reactions to prepare nanoscale materials that can be used widely in energy, electronics, and biomedical technologies and other applications. SCS is a specific variety of a more general combustion synthesis (CS) method. Investigations of the thermodynamics, kinetics, and the mechanisms of SCS reactions, are not as well studied as the other CS processes. This work reports on a systematic study of the thermodynamics and kinetics of SCS reactions involving Ni(NO 3 ) 2 , an oxidizer, and either glycine (C 2 H 5 NO 2 ) or hexamethylenetetramine (HMT, C 6 H 12 N 4 ) as fuels. A thermodynamic modeling approach, based on the Gibbs free energy minimization principle, is applied to the simultaneous calculations of the adiabatic temperatures and compositions of the equilibrium products. Our calculations reveal the influence of fuel-to-oxidizer ratio, amount of water, and the oxygen in air on the combustion temperature under adiabatic conditions and the composition of the resulting products. We have, in turn, measured the combustion temperature and phase composition of products and compared them with the calculations. Variations of drying times for the solutions yield precursor gels with varying water contents. This approach enables the manipulation of combustion parameters and confirms the use of calculated activation energies for reactions using the Merzhanov-Khaikin method. Here, the results show that SCS reactions in fuel-lean solutions producing NiO have higher activation energy in contrast to reactions with fuel-rich solutions that form Ni. Reduction of activation energies due to the increase in the fuel-to-oxidizer ratio could be related to the observed change of the rate-limiting stages of the endothermic decomposition of the individual reactants to the exothermic decomposition of coordinate compounds formed between the reactants.

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

Thermodynamic analysis of hydrogel swelling in aqueous sodium chloride solutions

Here, a comprehensive thermodynamic model is presented for swelling and salt partitioning of poly (n-isopropyl acrylamide) hydrogels in aqueous sodium chloride (NaCl) solutions at 298 K. Together with a Helmholtz energy expression for network elastic energy, a modified electrolyte Nonradom Two-Liquid activity coefficient model is used to represent both the long-range electrostatic interactions and the short-range van der Waals interactions present in the gel phase. The model parameters include one network parameter per hydrogel system plus two binary parameters per interaction pair in the polymer–solvent-salt systems. With semiquantitave agreement with experimental data for both hydrogel swelling and salt partitioning up to salt saturation, the analysis suggests the ion hydration between NaCl and water is the root cause for hydrogel deswelling in aqueous NaCl solution at high NaCl concentrations, i.e., NaCl weight fraction >0.03. On the other hand, the strong repulsive interaction between NaCl and hydrogel polymer is the dominant factor for salt partitioning.

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

Molecular-scale origins of solution nanostructure and excess thermodynamic properties in a water/amphiphile mixture

The molecular and nanoscale origins of nonideality in excess thermodynamic properties are essential to understanding cosolvent mixtures, yet they remain challenging to determine. Here, we consider a binary mixture of water and an amphiphile, N,N,N',N'-tetramethylmalonamide (TMMA), which is characterized by strong hydrogen bonding between the two components and no hydrogen bonding between amphiphiles. Using molecular dynamics simulation, validated with excess volume measurements and X-ray scattering, we identify three distinct solution regimes across the composition range of the binary mixture and find that the transition between two of these regimes, marked by the water percolation threshold, is closely correlated with minima in the excess volume and excess enthalpy. Structural analysis of the simulations reveals an interplay between local interactions and solution nanostructure, determined by the relative strength of the water–water and water–amphiphile hydrogen bonding interactions. By comparison with other amphiphiles, such as linear alcohols, the relative strength of like and unlike interactions between water and amphiphile affects the relationship between thermodynamics and structural regimes. This provides insight into how molecular forces of mutual solvation interact across length scales and how they manifest in excess thermodynamic properties.

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

Equilibrium Thermodynamic Properties of Aqueous Solutions of Ionic Liquid 1-Ethyl-3-Methylimidazolium Methanesulfonate [EMIM][MeSO3]

Abstract The ionic liquid 1-ethyl-3-methylimidazolium methanesulfonate ([EMIM][MeSO 3 ]) has been considered as a promising alternative desiccant to triethylene glycol and lithium bromide commonly used in the industry. In this paper, the water activity coefficient of this binary system was measured from 303 K to 363 K with water concentration from 18% to 92%. The interaction energies between the ionic liquid molecules ( $${g}_{22}$$ g 22 ) and between the ionic liquid and water molecules ( $${g}_{12}$$ g 12 ) for the [EMIM][MeSO 3 ]/water binary system were determined from the water activity coefficient data using the Non-Random Two-Liquid (NRTL) model. The magnitude of the interaction energy between the [EMIM][MeSO 3 ] and water molecules ( $${g}_{12}$$ g 12 ) was found to be in the range of 45~49 kJ/mol, which was about 20% larger than that between the water molecules ( $${g}_{11}$$ g 11 ) in the [EMIM][MeSO 3 ]/water system. The large ( $${g}_{12}$$ g 12 ) can explain many observed macroscopic thermodynamic properties such as strong hygroscopicity in the ionic liquid [EMIM][MeSO 3 ]. These interaction energies were used to determine the heat of desorption of the [EMIM][MeSO 3 ]/water system, and the obtained heat of desorption was in good agreement with that calculated from the conventional Clausius-Clapeyron Equation.

42 ENGINEERING↗
DOE OSTI2024

A simple method for obtaining heat capacity coefficients of minerals

Abstract Heat capacity data are unavailable or incomplete for many minerals at geologically relevant temperatures. Despite the availability of entropy and enthalpy values in numerous thermodynamic tables (even sometimes at elevated temperatures), there remains need for extrapolation beyond, or interpolation between, temperatures. This approach inevitably results in estimates for entropy and enthalpy values because the heat capacity coefficients required for optimal thermodynamic treatment are less frequently available. Here we propose a simple method for obtaining heat capacity coefficients of minerals. This method requires only the empirically measured temperature-specific heat capacity for calculation via a matrix algorithm. The system of equations solver is written in the Python computing language and has been made accessible in an online repository. Thermodynamically, the solution to a system of equations represents the heat capacity coefficients that satisfy the mineral-specific polynomial. Direct coefficient calculation will result in more robust thermodynamic data, which are not subject to fitting uncertainties. Using hematite as an example, this method provides results that are comparable to conventional means and is applicable to any solid material. Coefficients vary within the traditional large 950 K temperature interval, indicating that best results should instead utilize a smaller 400 K temperature interval. Examples of large-scale implications include the refinement of geothermal gradient estimation in rapidly subsiding sedimentary basins or metamorphic and hydrothermal evolution.

Geochemistry & Geophysics↗
DOE OSTI2022

Thermodynamic modeling of aqueous lithium salt solutions with association electrolyte nonrandom two-liquid activity coefficient model

The high charge density of lithium ion and the resulting strong association phenomena make thermodynamic modeling of aqueous lithium electrolyte solutions extremely challenging. In this study, the association electrolyte nonrandom two-liquid activity coefficient model of Lin et al. (AIChE J. 2022, 68(2), e17422) is utilized to correlate and predict thermodynamic properties and solubility behavior of aqueous single electrolyte solutions of LiCl, LiBr, LiI, and LiNO 3 , and their mixed electrolyte solutions. Capturing self-association of water, cross-association of ion and water for hydration, and cross-association of cation and anion for ion-pairing, the association model accurately represents the literature experimental data up to saturation concentrations and at the temperature ranging from 263 K to 523 K. Here, this study further investigated the effect of anions of the lithium salts, and re-confirmed that the order of solution non-ideality as LiI > LiBr > LiCl > LiNO 3 because the anions with stronger association strengths are more likely to form ion pairs and thus lower the mean ionic activity coefficients.

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

Affinity of small-molecule solutes to hydrophobic, hydrophilic, and chemically patterned interfaces in aqueous solution

Significance Increasing global population, urbanization, and climate change pose tremendous constraints on water availability. While membrane processes offer attractive, energy-efficient options, transformative advances in membrane technology are needed to enable reuse of poor-quality waters from diverse sources. This work brings fundamental understanding of solute–surface interactions to membrane phenomena, including selectivity and fouling. Molecular simulations reveal limitations of traditional predictive measures of solute–surface partitioning, suggesting alternative thermodynamic signatures for solute binding that highlight the contextual nature of hydrophobicity. An inverse design algorithm provides a tool to pattern membrane nanoscale chemistry to optimize solute affinity and selectivity. This work demonstrates how simulations can inform practical materials design for sustainable water treatment and, more broadly, many technologies involving water-mediated solute–surface interactions.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Atomistic simulations and machine learning of solute grain boundary segregation in Mg alloys at finite temperatures

Understanding solute segregation thermodynamics is the first step in investigating grain boundary (GB) properties, such as strong yttrium (Y) effects on grain growth and texture evolution in micro-scale polycrystalline magnesium (Mg) alloys. To estimate the average GB segregation behavior in low-solute-concentration Mg alloys (e.g., 2 at.% Y), a state-of-the-art spectral approach is applied based on a per-site segregation energy spectrum for Y solute atoms at zero K obtained from molecular statistics (MS) simulations of ~10 4 GB sites in Mg symmetric tilt GBs (STGBs). Although selected MS simulation results are consistent with verification by density functional theory (DFT) calculations, estimates of average segregation tendency based on the zero-K energy spectrum deviate from experimental observations. To resolve this problem, thermodynamic integration (TI) methods based on molecular dynamics (MD) simulations are used to determine the per-site segregation free energies of Y at representative GB sites, which show contributions beyond harmonic approximations can be important for certain GB sites at high temperatures. A surrogate model of per-site segregation free energy is constructed from a small set of TI data points using stacking cross-validation regressors and physics-informed descriptors. This model is applied to predict the Y segregation free energy spectra for thousands of GB sites in Mg STGBs with uncertainty quantification. Finally, the average segregation tendency predicted by the spectral approach based on the free energy spectra agrees well (within the uncertainty range) with experimental observations for micro-scale polycrystalline Mg alloys at typical thermomechanical processing temperatures (500 ~ 800 K), where thermodynamic equilibrium states are likely to be achieved due to fast diffusion.

Atomistic simulations↗
DOE OSTI2020

Engineering Multi-scale B2 Precipitation in a Heterogeneous FCC Based Microstructure to Enhance the Mechanical Properties of a Al0.5Co1.5CrFeNi1.5 High Entropy Alloy

While ordered L12 or gamma prime precipitates in face centered cubic (FCC) based microstructures have been extensively used for strengthening Nickel or Cobalt base superalloys, and more recently in high entropy alloys (HEAs) or complex concentrated alloys (CCAs), the possibility of exploiting ordered B2 precipitates in FCC-based systems has been relatively less investigated. The present study shows the propensity of developing a heterogeneous microstructure, consisting of two different distributions of FCC grain sizes, and two different size scales of B2 precipitates, within an FCC-based Al0.5Co1.5CrFeNi1.5 HEA/CCA. This alloy composition has been designed using solution thermodynamics-based modeling such that it has a high phase fraction and solvus temperature of the B2 phase. The resulting heterogenous microstructure exhibited an approximately 400% increase in yield strength with respect to the single-phase FCC solid solution condition of the same alloy while maintaining very good tensile ductility ~20 %.

Dasari, Sriswaroop↗
DOE OSTI2025

Stability of aqueous neodymium complexes in carbonate-bearing solutions from 100–600 °C

Rare earth element exploration requires a quantitative understanding of factors governing their mobilization and economic concentration. However, the behavior of rare earth elements in carbonate- bearing hydrothermal fluids associated with carbonatite-hosted deposits is poorly understood, and conflicting mechanisms of rare earth transport by anionic ligands and alkali behavior have been described. Here, we report quantitative data to characterize the role of carbonate-bearing solutions in the hydrothermal mobilization of neodymium. Solubility studies of neodymium phosphate were performed at temperatures ranging from 100 to 600 °C in carbonate-bearing solutions. The thermodynamic data determined for the predominant complex were used to model the separation of neodymium from thorium in a simple flow-through system based on fluid and mineral compositions characteristic of carbonatite deposits. Our data suggest that neodymium transport is controlled by the stability of the carbonate species NdCO 3 OH o , and at temperatures of 500–600 °C, the concentrations of neodymium in solutions can reach ~1000 ppm.

58 GEOSCIENCES↗
DOE OSTI2026

KBKit: A Python Toolkit for Kirkwood–Buff Theory from Molecular Dynamics

Thermodynamic properties of liquid mixtures govern processes that range from drug delivery to energy storage, yet extracting these properties from molecular simulations remains challenging. Kirkwood–Buff (KB) theory offers a rigorous route by linking microscopic pair distribution functions to macroscopic free energies, but practical use of the theory has been hindered by two obstacles: (i) the long simulations needed to obtain well-converged Kirkwood-Buff integrals (KBIs) and (ii) the specialized corrections required to translate finite-size data to the thermodynamic limit. $\texttt{KBKit}$ is an open-source Python package that removes these barriers. It automatically computes KBIs and derived thermodynamic quantities from GROMACS input files, applies state-of-the-art finite-size corrections, and provides built-in diagnostic tools to quantify statistical uncertainty. Written with modern software-engineering practices—continuous integration, extensive unit testing, and thorough documentation—$\texttt{KBKit}$ is both reliable and easy to extend. By condensing complex KBI analysis into a few intuitive commands, $\texttt{KBKit}$ enables researchers to incorporate KB theory into routine simulation workflows and accelerate the discovery of solution-phase thermodynamics.

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

Formation kinetics and thermodynamic stability of the Aurivillius compounds in Bi 4 Ti 3 O 12 –BiFeO 3 system

Abstract The Aurivillius compounds in the Bi 2 O 3 –Fe 2 O 3 –TiO 2 system, combining ferroelectric, semiconducting, and ferromagnetic properties, have attracted particular interest. Formation kinetics and thermodynamic stability are the fundamental knowledge needed for modeling and predicting the temporal microstructure and property evolution during materials processing but have not yet been addressed by quantitative experimental measurement. This article focuses on the Bi n +1 Fe n –3 Ti 3 O 3 n +3 Aurivillius compounds on the Bi 4 Ti 3 O 12 –BiFeO 3 tie‐line to elucidate the mechanisms and thermodynamic controls responsible for phase formation of compounds with various perovskite‐like layers. Five high‐purity Aurivillius compounds Bi 4 Ti 3 O 12 , Bi 5 FeTi 3 O 15 , Bi 6 Fe 2 Ti 3 O 18 , Bi 7 Fe 3 Ti 3 O 21 , and Bi 8 Fe 4 Ti 3 O 24 with integer n = 3–7 values were synthesized and their phase transformation properties and enthalpies of formation were studied by X‐ray diffraction in situ, high temperature differential scanning calorimetry, and high temperature oxide melt solution calorimetry. Thermodynamic stability of the compounds decreases with increasing n , and formation kinetics gradually slow down, demonstrating the inherent difficulty to synthesize pure Aurivillius compounds with n larger than 8. This difficulty was confirmed by an impurity phase coexisting with Bi 9 Fe 5 Ti 3 O 27 .

Gong, Weiping↗
DOE OSTI2022

Heterogenous activation of dynamic recrystallization and twinning during friction stir processing of a Cu-4Nb alloy

An interplay between high degree of shear deformation and deformation-induced heating occurs during friction stir processing (FSP) of metals. In medium-to-low stacking fault energy Cu alloys, this can lead to a complex spatially heterogenous activation of dynamic recrystallization (DRX) and twinning mechanisms. Within the Cu-Nb system, the presence of Nb is further expected to influence the DRX mechanism of the Cu matrix. However, the microstructural changes induced by the co-deformation of Nb during FSP are still not well understood. So, this study uses a combination of multimodal microstructural characterization, solution thermodynamics-based predictions, and computational crystal plasticity simulation to reveal the various microstructural evolution mechanisms that can occur during FSP of a Cu-4at.%Nb binary model alloy. The formation of softer DRX zones, and harder shear localization regions are revealed using electron backscatter diffraction, transmission electron microscopy, atom probe tomography, and crystal plasticity modeling.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The “good,” the “bad,” and the “hidden” in neutron scattering and molecular dynamics of ionic aqueous solutions

We characterize a concentrated 7.3 m CaCl 2 solution, combining neutron diffraction with chloride isotopic substitution (Cl-NDIS) in null water and molecular dynamics (MD) simulations. We elucidate the solution structure, thermodynamic properties, and extent of ion pairing previously suggested as concentration-dependent and often not observed at lower concentrations. Our Cl-NDIS measurements designate the solvent-shared ion pairing as dominant and the contact ion pairing (CIP) as insignificant even under conditions close to the solubility limit. The MD models parameterized against neutron diffraction with calcium isotopic substitution (Ca-NDIS) overestimate CIP despite successfully reproducing most of the Cl-NDIS signal. This drawback originates from the fact that Ca 2+ –Cl - interactions were primarily “hidden” in the Ca-NDIS signal due to overlapping with Ca 2+ –O w and Ca 2+ –H w contributions to the total scattering. Contrary, MD models with moderate CIP and possessing generally good performance at high concentrations fail to reproduce the NDIS measurements accurately. Therefore, the electronic polarization, introduced in most of the recent MD models via scaling ionic charges, resolves some but not all parameterization drawbacks. We conclude that despite improving the quality of MD models “on average,” the question “which model is the best” has not been answered but replaced by the question “which model is better for a given research.” An overall “good” model can still be inappropriate or, in some instances, “bad” and, unfortunately, produce erroneous results. Additionally, the accurate interpretation of several NDIS datasets, complemented by MD simulations, can prevent such mistakes and help identify the strengths, weaknesses, and convenient applications for corresponding computational models.

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

Nuclear Thermal Energy Storage Configurations for Industrial Combined Heat and Power Supply: Conceptual Study and Engineering Designs

The industries examined in this report primarily rely on moderate-temperature heat provided by gas- or coal-fired boilers and combined heat and power (CHP) plants, delivered through standard process steam systems. High-temperature energy demands are often industry-specific and typically exceed the capabilities of high-temperature gas-cooled reactors (HTGRs). While it is technically feasible to replace process steam from fossil-based heat sources with nuclear energy, certain industries, such as methanol production and pulp and paper, face technoeconomic challenges in integrating nuclear energy without major changes or a technological shift. This is mainly due to the limited external energy demand remaining after the use of internal byproducts, waste heat recovery, and simple efficiency improvements. Achieving full decarbonization of these processes with nuclear energy would require significant technological advancements, involving experimental technology and substantial investments, making widespread adoption in existing industrial plants unlikely in the near term. This study reviews TES options in the context of enabling a flexible CHP supply while maintaining a steady nuclear heat input. Heat storage systems that interface between the reactor primary fluid and the CHP system offer superior performance and flexibility. Specifically, steam extraction downstream of the reheater with a two-tank molten-salt TES appears as the best solution regarding thermodynamic system benefits and system drawbacks. Using selected system configurations, a conceptual design of an industrial energy park was developed for industries with varying energy demands, such as steel production plants utilizing electric arc furnaces (EAFs) and chemical plants, as well as for those with constant energy demands, like petroleum refineries. This design highlights the capabilities of TES and explores its potential business cases. The study also conceptually develops the potential for integrating additional energy sources with nuclear systems through the implementation of TES. The potential of the HTGR-TES-CHP system was also evaluated considering key uncertainties such as industrial demand profiles, external grid access availability, and eligible tax credit levels, using the Holistic Energy Resource Optimization Network. Sensitivity of net present value to these uncertainties was analyzed to determine the optimal number of nuclear reactors (and CHP systems) and the suitable TES capacity. The results were interpreted from a decision-maker’s perspective, focusing on three key areas: deployment strategy (oversized units vs. undersized units with TES support), industrial process characteristics (thermal-intensive single profiles vs. electricity-intensive combined profiles), and operational goals (maximizing profits vs. minimizing natural gas (NG) consumption or external grid dependence). The optimization results indicate that the HTGR-TES-CHP system significantly reduces reliance on NG boilers for individual industrial processes by 9-60% (in NG capacity factor), with an average reduction of 38%, compared to standalone NG boiler operation case (Business As Usual [BAU]). For combined industrial processes, the reduction ranges from 37-77%, with an average of 60%. Additionally, the system greatly reduces dependence on external grids. In meeting industrial electrical demands, a 33-100% self-sufficient internal electricity supply is achieved for single industrial process, with an average of 74%, compared to the BAU scenario, where 100% of electricity is imported. For combined processes, 35-100% of internal electricity demands are met by the reactor, with an average of 73%. At last, the relative NG price levels at which the proposed HTGR-TES-CHP system can cost-effectively enter the market currently dominated by existing NG boilers were estimated. For a moderate HTGR CAPEX level ($\$$2500/kWth, $\$$6329/kWe), the analysis suggests that NG prices must be 2.5 to 7 times higher than HTGR variable operating and maintenance costs for single industrial process, and 5.5 to 9.5 times higher for a combined process scenario. Tax credit modeling shows that the Investment Tax Credit significantly reduces the price threshold needed to break even, making the system competitive with NG boilers in certain cases.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗