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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

Real-Gas Effects on Binary Mixing Layers

This paper presents a computational study of real-gas effects on the mean flow and temporal stability of heptane/nitrogen and oxygen/hydrogen mixing layers at supercritical pressures. These layers consist of two counterflowing free streams of different composition, temperature, and density. As in related prior studies reported in NASA Tech Briefs, the governing conservation equations were the Navier-Stokes equations of compressible flow plus equations for the conservation of total energy and of chemical- species masses. In these equations, the expressions for heat fluxes and chemical-species mass fluxes were derived from fluctuation-dissipation theory and incorporate Soret and Dufour effects. Similarity equations for the streamwise velocity, temperature, and mass fractions were derived as approximations to the governing equations. Similarity profiles showed important real-gas, non-ideal-mixture effects, particularly for temperature, in departing from the error-function profile, which is the similarity solution for incompressible flow. The temperature behavior was attributed to real-gas thermodynamics and variations in Schmidt and Prandtl numbers. Temporal linear inviscid stability analyses were performed using the similarity and error-function profiles as the mean flow. For the similarity profiles, the growth rates were found to be larger and the wavelengths of highest instability shorter, relative to those of the errorfunction profiles and to those obtained from incompressible-flow stability analysis. The range of unstable wavelengths was found to be larger for the similarity profiles than for the error-function profiles

Okong'o, Nora↗
DOE OSTI2021

Molecular Properties and Chemical Transformations Near Interfaces

The properties of bulk water and aqueous solutions are known to change in the vicinity of an interface and/or in a confined environment, including the thermodynamics of ion selectivity at interfaces, transition states and pathways of chemical reactions, and nucleation events and phase growth. Here we describe joint progress in identifying unifying concepts about how air, liquid, and solid interfaces can alter molecular properties and chemical reactivity compared to bulk water and multicomponent solutions. Furthermore, we also discuss progress made in interfacial chemistry through advancements in new theory, molecular simulation, and experiments.

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

Thermodynamic and Kinetic Modulation of Methylammonium Lead Bromide Crystallization Revealed by In Situ Monitoring

Hybrid organic–inorganic perovskite (HOIP) crystals are promising optoelectronic materials, but little is known about either the thermodynamic and kinetic controls on crystal growth or the underlying growth mechanism(s). Herein, we use fluid cell atomic force microscopy (AFM) and solution nuclear magnetic resonance (NMR) spectroscopy to investigate growth of the model HOIP crystal CH 3 NH 3 PbBr 3 (MAPbBr 3 ) and to determine how formic acid (HCOOH) modulates the thermodynamics and kinetics of growth. The results show that growth of MAPbBr 3 in dimethylformamide (DMF) proceeds through the classical pathway by the spreading of atomic crystal steps generated at screw dislocations on the {100} surface. Temperature dependent step velocity measurements demonstrate that with increasing concentration, HCOOH decreases both the solubility of MAPbBr 3 and the kinetic coefficient (b) of step movement. 1 H-NMR measurements indicate that HCOOH increases the lifetime of the methylammonium (MA + ) ions and promotes the association of MAPbBr 3 , thus tuning the solubility of the perovskite. HCOOH also alters the molecular tumbling motion and bulk diffusion of the MA + ions, possibly via H-bonding. Further, our findings establish a direct correlation between the mesoscale crystal growth kinetics and the molecular-scale interactions between organic additives and constituent ions, providing unprecedented insights for developing predictive syntheses of HOIP crystals with defined size, crystal habit and shape, and defect distribution.

36 MATERIALS SCIENCE↗
NASA NTRS2013

Prediction of Ablation Rates from Solid Surfaces Exposed to High Temperature Gas Flow

A mathematical model and a solution algorithm is developed to study the physics of high temperature heat transfer and material ablation and identify the problems associated with the flow of hydrogen gas at very high temperatures and velocities through pipes and various components of Nuclear Thermal Rocket (NTR) motors. Ablation and melting can be experienced when the inner solid surface of the cooling channels and the diverging-converging nozzle of a Nuclear Thermal Rocket (NTR) motor is exposed to hydrogen gas flow at temperatures around 2500 degrees Kelvin and pressures around 3.4 MPa. In the experiments conducted on typical NTR motors developed in 1960s, degradation of the cooling channel material (cracking in the nuclear fuel element cladding) and in some instances melting of the core was observed. This paper presents the results of a preliminary study based on two types of physics based mathematical models that were developed to simulate the thermal-hydrodynamic conditions that lead to ablation of the solid surface of a stainless steel pipe exposed to high temperature hydrogen gas near sonic velocities. One of the proposed models is one-dimensional and assumes the gas flow to be unsteady, compressible and viscous. An in-house computer code was developed to solve the conservations equations of this model using a second-order accurate finite-difference technique. The second model assumes the flow to be three-dimensional, unsteady, compressible and viscous. A commercial CFD code (Fluent) was used to solve the later model equations. Both models assume the thermodynamic and transport properties of the hydrogen gas to be temperature dependent. In the solution algorithm developed for this study, the unsteady temperature of the pipe is determined from the heat equation for the solid. The solid-gas interface temperature is determined from an energy balance at the interface which includes heat transfer from or to the interface by conduction, convection, radiation, and ablation. Two different ablation models are proposed to determine the heat loss from the solid surface due to the ablation of the solid material. Both of them are physics based. Various numerical simulations were carried out using both models to predict the temperature distribution in the solid and in the gas flow, and then predict the ablation rates at a typical NTR motor hydrogen gas temperature and pressure. Solid mass loss rate per foot of a pipe was also calculated from these predictions. The results are presented for fully developed turbulent flow conditions in a sample SS pipe with a 6 inch diameter.

Akyuzlu, Kazim M.↗
DOE OSTI2021

Thermodynamic stabilization of nanocrystalline aluminum

Nanocrystalline metals are generally unstable due to a large volume fraction of high-energy grain boundaries associated with a small grain size. Preferential dopant segregation to the high-energy grain boundaries is observed to enhance the stability of the material’s microstructure by minimizing its energy. Nanocrystalline aluminum-dopant systems were evaluated for thermodynamic stability against grain growth and phase precipitation via the mechanism of grain boundary segregation according to a modified regular nanocrystalline solution model. Fifty-one potential dopant elements have been evaluated for their efficacy in stabilizing nanostructures with three potential candidates, magnesium, lanthanum, and silicon, identified possessing the characteristics to promote grain boundary segregation and a state of thermodynamic stability in aluminum’s nanocrystalline regime. Here, the minimum dopant content required to achieve nanocrystalline microstructure stability is identified for each of the three candidate elements. Beyond this minimum content, further addition of the dopant elements decreased the final microstructure’s stability with no effects on the existence of a stable nanocrystalline state.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A carbon-efficient bicarbonate electrolyzer

Carbon efficiency is one of the most pressing problems of carbon dioxide electroreduction today. While there have been studies on anion exchange membrane electrolyzers with carbon dioxide (gas) and bipolar membrane electrolyzers with bicarbonate (aqueous) feedstocks, both suffer from low carbon efficiency. In anion exchange membrane electrolyzers, this is due to carbonate anion crossover, whereas in bipolar membrane electrolyzers, the exsolution of carbon dioxide (gas) from the bicarbonate solution is the culprit. Here, we first elucidate the root cause of the low carbon efficiency of liquid bicarbonate electrolyzers with thermodynamic calculations and then achieve carbon-efficient carbon dioxide electroreduction by adopting a near-neutral-pH cation exchange membrane, a glass fiber intermediate layer, and carbon dioxide (gas) partial pressure management. We convert highly concentrated bicarbonate solution to solid formate fuel with a yield (carbon efficiency) of greater than 96%. A device test is demonstrated at 100 mA cm -2 with a full-cell voltage of 3.1 V for over 200 h.

42 ENGINEERING↗
DOE OSTI2025

Toward Establishing Uniqueness of Experimentally Determined Transference Numbers

The passage of current through a battery results in the development of concentration gradients in the electrolytic phase. For a fully characterized binary electrolyte, where the conductivity, salt diffusion coefficient, cation transference number, and the thermodynamic factor are known, concentration and potential gradients in the electrolytic phase can be modeled using Newman’s concentrated solution theory. We report two methods for measuring the transference number: the standard method based on electrochemical measurements ( t + , echem 0 ) and electrophoretic NMR ( t + , eNMR 0 ). The electrochemical approach requires combining measurements from multiple experiments; the equations used to determine the cation transference number and the thermodynamic factor are coupled, nonlinear algebraic equations. In the electrophoretic-NMR-based approach, however, the equations used to determine the cation transference number and the thermodynamic factor are decoupled. We find for a liquid electrolyte comprised of a lithium salt dissolved in tetraglyme, the values of the transference numbers obtained by these two methods are distinct. For example, at 30 °C, t + , echem 0 = −1.02 ± 1.11 and t + , eNMR 0 = 0.25 ± 0.04. The corresponding thermodynamic factors are also different. While the magnitude of the predicted concentration gradients based on the two sets of parameters are different, the predicted current-voltage relationships are similar.

Hickson, Darby T. (ORCID:0000000251339755)↗
DOE OSTI2021

Solution and Interface Structure and Dynamics in Geochemistry: Gateway to Link Elementary Processes to Mineral Nucleation and Growth

Our understanding of the process of mineral nucleation is starting to evolve from the solely thermodynamical aspects in classical nucleation theory toward a mechanistic study of reactions influenced by solvation structures of ions and their dynamics. A knowledge of how these atomic and molecular interactions respond to changing solution composition and/or the nature of an interface will form a fundamental foundation to understand rate-limiting reactions and processes for nucleation and growth. Here, we review recent advances from multidisciplinary experimental and computational approaches of the structure and dynamics of aqueous solutions and mineral–liquid interfaces that will influence rates of nucleation (and growth behavior) both homogeneously in aqueous solutions and heterogeneously at mineral–liquid interfaces. Furthermore, these processes include structure (and extent) of solute clustering, solvent degrees of freedom, and the presence of noninnocent charge-balancing ions, all of which can drive the selectivity and exert control over nucleation mechanisms. Similar controlling reactions can also occur at interfaces with the added complexities of how solute ions interact with and assemble on surface sites and use different mechanisms from bulk solutes, such as lateral solvent exchanges between distinct surface sites. To gain a predictive understanding of nucleation chemistry, we emphasize several key knowledge gaps that will need to be addressed in this field and highlight the importance of combining computations and experiments.

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

Unseeded, spontaneous nucleation of spherulitic magnesium calcite

Most of the sedimentary carbonates deposited in the marine environments are composed of calcium carbonate minerals with varying amounts of incorporated Mg 2+ . However, understanding how interactions of impurities with carbonate and their incorporation affect sediments behavior remains a challenge. Here, a new insight is obtained by monitoring solution composition, morphology, and electrokinetic potential of carbonate particles formed in a spontaneous unseeded batch precipitation experiment using electrochemical and scanning electron microscopy methods. The solid composition and growth rate are extracted from changes in the bulk composition and fitted to chemical affinity rate law, revealing that the precipitation pathway consists of second-order dissolution and first-order precipitation. The molecular dynamics simulations show that the lattice strain induced by randomly substituting Ca 2+ by Mg 2+ stabilizes spherical nanoparticles and reduces their surface area and volume. Combining kinetics and thermodynamics insight, we conclude that variation in the carbonate bulk and interfacial energies, along with the solution supersaturation, lead to the dissolution-precipitation transformation pathway from Mg-rich to Mg-poor carbonate phase that preserves spherulitic morphology. Furthermore, our findings are relevant for long-standing questions of how impurities influence diagenesis of carbonate sediments and spherulitic carbonate particles' origin.

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

Structure of Complex Liquid–Liquid Extraction Organic Phases for Rare Earth Separations

Complex, multicomponent liquids with hierarchical structure and phase transitions are encountered in many natural and industrial processes, including in chemical separations. One notable example is aggregation and organic phase splitting in liquid–liquid extraction (LLE) of metal ions. While these two phenomena that have long been closely associated, a mechanistic link between mesoscale structure and the capacity-limiting organic phase splitting remains elusive due to complexity of these systems. Here, in this study, we combine small-angle X-ray scattering (SAXS), X-ray photon correlation spectroscopy (XPCS), and molecular dynamics simulation to reveal a comprehensive picture of structure at the nano- and mesoscale in these complex solutions. For the representative case of rare earth extraction from an acidic aqueous phase by a malonamide extractant in dodecane, we investigate a wide range of process-relevant extractant and acid concentrations to provide a complete picture of how aggregation depends on composition. We decompose organic phase structure from SAXS into two contributions, which together can capture the scattering at all compositions: composition fluctuations described by the Ornstein–Zernike equation at low wavenumber Q, and nanostructure modeled by a “pre-peak” at intermediate Q. The former contains information about the thermodynamics of demixing, while the latter reflects nanoscopic self-assembly of the extractant and extracted solutes. While fluctuations have typically not been considered in the literature, we find they in fact dominate the total structure for nearly all practical conditions. As only the fluctuations have a strong temperature response, we confirm this attribution with temperature-dependent SAXS measurements, including for extracted europium nitrate complexes. SAXS and XPCS measurements near the critical point find static and dynamic scaling consistent with theory. Overall, this new paradigm for understanding LLE organic phases connects composition, nanoscale, and mesoscale structuring to phase behavior, providing both a comprehensive picture of solution structure and a quantitative link between aggregation and third phase formation.

Peroutka, Allison A. [Argonne National Laboratory ↗
DOE OSTI2022

A fast and accurate modeling approach for water and steam thermodynamics with practical applications in district heating system simulation

In U.S. district heating (DH) systems, steam is the most common heat transport medium. Industry demand for new advanced modeling capabilities of complete steam DH systems is increasing; however, the existing models for water/steam thermodynamics are too slow for large system simulations because of computationally expensive algebraic loops that require the solution to nonlinear systems of equations. For practical applications, this work presents a novel split-medium approach that implements numerically efficient liquid water models alongside various water/steam models, breaking costly algebraic loops by decoupling mass and energy balance equations. New component models for steam DH systems are also presented. We implemented the models in the equation based Modelica language and evaluated accuracy and computing speed across multiple scales: from fundamental thermodynamic properties to complete districts featuring 10 to 200 buildings. Compared to district models with the IF97 water/steam model and equipment models from the Modelica Standard Library, the new implementation improves the scaling rate for large districts from cubic to quadratic with negligible compromise to accuracy. Additionally, for an annual simulation with 180 buildings, this translates to a computing time reduction from 33 to 1-1.5 h. These results are critically important for industry practitioners to simulate steam DH systems at large scales.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Mechanism of ion pairing–induced co-ion penetration into the stern layer

In classical electric double layer (EDL) theory, co-ions—ions carrying the same charge as the electrode—are assumed to be excluded from the Stern layer because of electrostatic repulsion. However, ion pairing with adsorbed counterions may enable co-ions to access the EDL, a phenomenon that remains underexplored. Key questions therefore remain unresolved, including how deeply co-ions penetrate, how ion-specific properties govern them, and how ion pairing within the EDL differs from bulk solution. Here, in this study, we employ molecular dynamics simulations combined with well-tempered metadynamics to quantify co-ion penetration and ion-pairing thermodynamics at positively charged Mg-Al layered double hydroxide interfaces (σ = +20.6 μC/cm 2 ) in aqueous chloride electrolytes across 0.005–3 M range. We reveal a concentration-driven crossover: at dilute conditions, counterion–co-ion contact ion pairing (CIP) is thermodynamically unfavorable or only weakly favorable (ΔG CIP = +0.12 kcal/mol for Li + to −0.42 kcal/mol for K + ), leading to co-ion exclusion. As concentration increases (≥0.5 M), enhanced screening reduces surface repulsion and stabilizes CIP thermodynamics (ΔG CIP = −0.89 kcal/mol for Li + to −1.05 kcal/mol for K + at 3 M), enabling monovalent co-ions to penetrate the Stern layer following the hydration hierarchy K + > Na + > Li + , while divalent co-ions (Mg 2+ , Ca 2+ ) remain in the diffuse layer, forming only solvent-separated pairs. Comparison with bulk solution shows interfacial ion pairing is suppressed at low concentrations but converges to bulk-like thermodynamics at high concentrations, indicating surface effects vanish under strong screening. These results establish a quantitative framework linking ion-pairing thermodynamics to co-ion penetration and providing strategies to design electrochemical interfaces.

Co-ion penetration↗
DOE OSTI2025

Mechanistic Transformation of CuI Nanoparticles Into Oxidation‐Resistant 2D Copper Nanoplates

Unconventional phase transformations reveal new crystallization mechanisms, yet direct observation of such pathways during nanoscale solution-phase synthesis remains challenging. This study uncovers an atypical growth process in which thermodynamically stable CuI nanoparticles (NPs) transform into high-energy 2D Cu plates. Using a combination of in situ transmission electron microscopy, ex situ structural analysis, and density functional theory calculations shows that the formation of structural defects induced by hexadecylamine and chloride ions facilitates the transformation by promoting surface iodine vacancies. The resulting Cu{111} nanoplates, with ultrathin thicknesses (≈4 nm) and exceptionally high aspect ratios (≈450), display enhanced oxidation resistance and long-term stability under ambient conditions. This resistance is attributed to the close-packed {111} facets, which suppress chemical oxidation even after extended exposure to air over 100 days. These findings provide new insights into non-classical crystallization pathways in metal nanomaterials and suggest a versatile approach for preparing oxidation-resistant, structurally defined Cu nanostructures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Predicting zeolites’ stability during the corrosion of nuclear waste immobilization glasses: Comparison with glass corrosion experiments

During the long-term corrosion of nuclear waste glasses under nuclear waste disposal conditions, the precipitation of zeolitic phases has been linked to a delayed acceleration in glass corrosion (known as “Stage III”). Hence, predicting the thermodynamic propensity for zeolites to form upon the dissolution of nuclear waste glasses is key to ensure their long-term performance. Here, we compile a unified, internally-consistent thermodynamic database “clay20” to estimate the stability of clay and feldspar phases relevant to nuclear waste immobilization glasses, including beidellite(Mg, Ca, Na, K), kaolinite, montmorillonite(Mg, Ca, Na, K), nontronite(Mg, Ca, Na, K), saponite(Ca, Na, K), and albite. Based on this, we report a geochemical modeling method allowing us to predict the stability of secondary phases (including zeolites, calcium–silicate–hydrate gels, and clays) upon the dissolution of nuclear waste immobilization glasses. We show that this approach offers a realistic description of the stability of the secondary phases forming during the dissolution of two archetypical model nuclear glasses (namely, the International Simple Glass, ISG, and WVUTh-203) under conditions relevant to nuclear waste disposal (T = 90°C, p = 1 bar) as a function of pH. We find that the formation of silica and clay secondary phases is thermodynamically favored at low pH (pH < 10), whereas, in contrast, zeolite (analcime) and calcium–silicate–hydrate phases are favored at high pH (pH > 10.5). This suggests that thermodynamics (i.e., not solely kinetics) plays a key role in determining the range of solution pH wherein stage III corrosion may occur, i.e., when zeolite formation is favored.

Zhen-Wu, Bi Yun↗
DOE OSTI2024

Role of Metal–Organic Framework Topology on Thermodynamics of Polyoxometalate Encapsulation

Polyoxometalates (POMs) are discrete anionic clusters whose rich redox properties, strong Brønsted acidity, and high availability of active sites poise them as potent catalysts for oxidation reactions. Here, metal–organic frameworks (MOFs) have emerged as tunable, porous platforms to immobilize POMs, thus increasing their solution stability and catalytic activity. While POM@MOF composite materials have been widely used for a variety of applications, little is known about the thermodynamics of the encapsulation process. Here, we utilize an up-and-coming technique in the field of heterogeneous materials, isothermal titration calorimetry (ITC), to obtain full thermodynamic profiles (ΔH, ΔS, ΔG, K a ) of POM binding. Six different 8-connected hexanuclear Zr-MOFs were investigated to determine the impact of MOF topology (csq, scu, the) on POM encapsulation thermodynamics.

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

Synthesis of core/shell nanocrystals with ordered intermetallic single-atom alloy layers for nitrate electroreduction to ammonia

Structurally ordered intermetallic nanocrystals (NCs) and single-atom catalysts (SACs) are two emerging catalytic motifs for sustainable chemical production and energy conversion. However, both have synthetic limitations which can lead to the aggregation of NCs or metal atoms. Single-atom alloys (SAAs), which contain isolated metal atoms in a host metal, can overcome the aggregation concern because of the thermodynamic stabilization of single atoms on host metal surfaces. We report a direct solution-phase synthesis of Cu/CuAu core/shell NCs with tunable SAA layers. This synthesis can be extended to other Cu/CuM (M = Pt, Pd) systems, in which M atoms are isolated in the copper host. Using this method, the density of SAAs on a copper surface can be controlled, resulting in both low and high densities of single atoms. Alloying gold into the copper matrix introduced ligand effects that optimized the chemisorption of *NO 3 and *N. As a result, the densely packed Cu/CuAu material demonstrated a high selectivity toward NH 3 from the electrocatalytic nitrate reduction reaction with an 85.5% Faradaic efficiency while maintaining a high yield rate of 8.47 mol h -1 g -1 . This work advances the design of atomically precise catalytic sites by creating core/shell NCs with SAA atomic layers, opening an avenue for broad catalytic applications.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Numerical fluid dynamics for FRG flow equations: Zero-dimensional QFTs as numerical test cases. II. Entropy production and irreversibility of RG flows

Herein we demonstrate that the reformulation of renormalization group (RG) flow equations as nonlinear heat equations has severe implications on the understanding of RG flows in general. We demonstrate by explicitly constructing an entropy function for a zero-dimensional Z 2 -symmetric model that the dissipative character of generic nonlinear diffusion equations is also hard-coded in the functional RG equation. This renders RG flows manifestly irreversible, revealing the semigroup property of RG transformations on the level of the flow equation itself. Additionally, we argue that the dissipative character of RG flows, its irreversibility and the entropy production during the RG flow may be linked to the existence of a so-called C– / A-function. In total, this introduces an asymmetry in the so-called RG time—in complete analogy to the thermodynamic arrow of time—and allows for an interpretation of infrared actions as equilibrium solutions of dissipative RG flows equations. The impossibility of resolving microphysics from macrophysics is evident in this framework. Furthermore, we directly link the irreversibility and the entropy production in RG flows to an explicit numerical entropy production, which is manifest in diffusive and non-linear partial differential equations (PDEs) and a standard mathematical tool for the analysis of PDEs. Using exactly solvable zero-dimensional Z 2 -symmetric models, we explicitly compute the (numerical) entropy production related to the total variation nonincreasing property of the PDE during RG flows toward the infrared limit. Finally, we discuss generalizations of our findings and relations to the C– / A-theorem as well as how our work may help to construct truncations of RG flow equations in the future, including numerically stable schemes for solving the corresponding PDEs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Singularity-EOS XCAP Report

Solving the Euler equations is a fundamental component of simulating many physical phenomena ranging from high explosives to astrophysics. An equation of state (EOS) is a required piece that relates any two thermodynamic quantities to all other thermodynamic values. The presence of multiple materials within a control volume further complicates the solution requiring additional equations to describe the interaction of materials at a sub-grid level. Equations of state themselves can also come in many forms ranging from simple algebraic relations to more complicated differential equation models that describe material interactions over a broad range of physical conditions. In the latter case, the EOS is often pre-computed at a given set of grid points and provided in a tabular form where additional properties can be derived from the interpolation functions.

97 MATHEMATICS AND COMPUTING↗