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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Solvent Phase Optimizations Improve Correlations with Experimental Stability Constants for Aqueous Lanthanide Complexes

Stability constants provide insight into ion complexation in water. While computational studies have been shown to model the energy of the complexation successfully using a thermodynamic cycle approach, it does not extend to calculating the stability constants for 1:1 lanthanide to ligand complexes in solution. Using B3LYP and 6-31+G* Pople basis with small core effective core potential (ECP) on the lanthanum ion, and a solvent model based on the full solute electron density (SMD) solvation model we computed and compared with previously published stability constants of the ligands: acetate, acetohydroximate, acetylacetonate, methanoate, tropolonate, hydroxide, catecholate, malonate, oxalate, phthalate, and sulfate. The best R 2 values for the thermodynamic cycle can only be determined by separating the mono and divalent ions to achieve an R 2 value of 0.86 and 0.74 for mono and divalent ions, respectively. We show that by optimizing the lanthanide-ligand structures in implicit solvent, we achieve an improved correlation between experimental and computed stability constants of R 2 value of 0.89 for the combined mono and divalent ions.

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

Computational methods in solution-based plastics purification

Plastic waste can be recycled into resins with near-virgin properties by solution-based purification processes that selectively dissolve polymers, remove contaminants, or detach printing residues. Here, in this review, we examine computational methods for predicting the behavior governing solution-based plastic purification, motivated by the vast polymer–solvent–contaminant compositional space. We discuss thermodynamic and machine learning methods for predicting polymer–solvent and polymer–contaminant interaction and review physics-based molecular dynamics simulations that resolve molecular-scale phenomena within polymer matrices inaccessible to screening methods. We highlight how these methods have informed experimental design for dissolution-based recycling and solvent-based contaminant removal. Finally, we discuss the prospective role of agentic AI in integrating these computational tools with real-time sorting data to adapt purification conditions to the compositional variability of real post-consumer feedstocks. This review charts a path toward computationally guided solution-based purification workflows that can respond to the complexity inherent in plastic waste streams.

Altamimi, Ali [Univ. of Wisconsin, Madison, WI (Un↗
DOE OSTI2020

Exploring competitive metal binding and crystallization of UO 2 2+ and Cu 2+ tetrahydrofuran-2,3,4,5-tetracarboxylic acid complexes

Solvent extractions are used to separate actinide elements from fission products in nuclear waste streams and the successful isolation of specific species utilize subtle differences in metal ligand binding. Metal ligand binding is also important in crystallization of metal organic materials and herein we explore the importance of competitive binding in the tetrahydrofuran-2,3,4,5-tetracarboxylic acid (THFTCA) system. This ligand has been previously evaluated for the selective extraction of uranium from other lanthanides and actinides and in the current research, we evaluate the crystallization of uranyl-THFTCA complexes in the presence of Cu 2+ , Sr 2+ , Th 4+ , and Ce 3+ . Three major phases were formed in the crystallization experiments and characterized with single crystal X-ray diffraction, powder X-ray diffraction, thermogravimetric analysis, and Raman spectroscopy. Two of the resulting phases were novel (UTHF1 ((C 4 H 10 N 2 )[UO 2 (C 8 H 8 O 9 ) 2 ]∙2H 2 O) and UTHF2 (Na[(UO 2 )(C 8 H 5 O 9 )(H 2 O)]∙ 3.5 H 2 O)), whereas the third (CuTHF1) was previously reported in the literature. Raman spectroscopy was utilized to evaluate spectral changes in the mother liquor of UTHF1, UTHF2, and CuTHF1 over time to assess the crystallization process. Further, isothermal titration calorimetry was used to determine the binding constants for UO 2 2+ and Cu 2+ to the THFTCA ligand in solution and evaluate the role of competitive metal binding in this system. The thermodynamic parameters and the crystallographic data were used to justify the formation of a weaker UO 2 2+ -THFTCA complex. Formation of a weaker UO 2 2+ -THFTCA complex supports our findings that UTHF1 only forms in homomeric systems, but suggests that the crystallization of UTHF2 and CuTHF1 is reliant on the presence of additional counter ions and ligands to change the amount of available THFTCA ligand.

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

Carbon dioxide capture with aqueous amino acids: Mechanistic study of amino acid regeneration by guanidine crystallization and process intensification

CO 2 capture from powerplant-generated flue gas via a phase-changing process involving absorption with aqueous amino acids (e.g., glycine or sarcosine) and bicarbonate crystallization with bis-iminoguanidines (e.g., glyoxal-bis-iminoguanidine or GBIG) is investigated in this paper. This process is of high interest due to its potential to decrease the energy penalty for CO 2 capture by significantly reducing the solvent regeneration energy typically associated with aqueous amine solvents. A critical step in the proposed CO 2 capture mechanism is the regeneration of the amino acid by removal of protons and bicarbonate ions from solution through crystallization of GBIGH 2 2+ bicarbonate salt. Here, we investigated the thermodynamics and kinetics of glycine regeneration by crystallization of GBIGH 2 2+ (HCO 3 – ) 2 (H 2 O) 2 . A theoretical model was developed and compared to experimental data to simulate and predict the glycine regeneration and determine its reaction mechanism. This combined experimental and theoretical study led to the conclusion that, while the GBIGH 2 2+ bicarbonate crystallization step provides most of the thermodynamic driving force for the glycine regeneration, the rate-limiting step is the protonation of GBIG prior to crystallization. The CO 2 loading and amino acid regeneration steps were combined into a single, intensified process using a bubble column reactor. The CO 2 loading capacity of GBIG was experimentally determined to be roughly 1.36 mol CO 2 per mol GBIG. These results provide the fundamental basis for developing an effective carbon capture technology with phase-changing amino acid/guanidine absorbents.

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

Molecular Programming of Diorganyl Dichalcogenides for Rational Nanocrystal Design

Soft-chemistry nanocrystal synthesis leverages low-temperature, solution-phase reactions to access materials that can be kinetically stabilized rather than thermodynamically favored. Under such mild conditions, reaction pathways are governed not only by precursor composition but also by the molecular details that dictate how reactive atomic species are generated and delivered. Furthermore, harnessing this kinetic sensitivity offers a powerful opportunity: by deliberately programming precursor reactivity, nanocrystal composition, structure, and morphology can be rationally designed rather than empirically discovered.

Chalcogenides↗
DOE OSTI2020

Understanding Cation Selectivity in Carbon Nanopores with Hybrid First-Principles/Continuum Simulations: Implications for Water Desalination and Separation Technologies

Understanding ion adsorption in porous carbons is crucial for a range of technologies, including water desalination and energy storage. In this work, we combined density functional theory with a continuum solvation model to investigate thermodynamics and kinetics of the adsorption process of alkali metal ions from aqueous solutions into carbon nanopores with different sizes and geometries. We found that cations with a larger ionic radius are more favorable to enter the nanopores because of a lower energy penalty of dehydration. In addition, the pore size and geometry were found to have a significant impact on the ion–pore interaction under confinement and cation selectivity. Our study highlights a complex interplay among nanopore geometry, ion size, and hydration on the cation adsorption selectivity, suggesting that tuning the porosity could represent a general strategy for improving ion separations.

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

Impact of Mg Substitution on the Structure, Stability, and Properties of the Na 2 Fe 2 F 7 Weberite Cathode

Of the few weberite-type Na-ion cathodes explored to date, Na 2 Fe 2 F 7 exhibits the best performance, with capacities up to 184 mAh/g and energy densities up to 550 Wh/kg reported for this material. However, the development of robust structure–property relationships for this material is complicated by its tendency to form as a mixture of metastable polymorphs, and transform to a lower-energy Na y FeF 3 perovskite compound during electrochemical cycling. Our first-principles-guided exploration of Fe-based weberite solid solutions with redox-inactive Mg 2+ and Al 3+ predicts an enhanced thermodynamic stability of Na 2 Mg x Fe 2–x F 7 as the Mg content is increased, and the x = 0.125 composition is selected for further exploration. We demonstrate that the monoclinic polymorph (space group C2/c) of Na 2 Fe 2 F 7 (Mg0) and of a new Mg-substituted weberite composition, Na 2 Mg 0.125 Fe 1.875 F 7 (Mg0.125), can be isolated using an optimized synthesis protocol. The impact of Mg substitution on the stability of the weberite phase during electrochemical cycling, and on the extent and rate of Na (de)intercalation, is examined. Irrespective of the Mg content, we find that the weberite phase is retained when cycling over a narrow voltage window (2.8–4.0 V vs Na/Na+). Over a wider voltage range (1.9–4.0 V), Mg0 shows steady capacity fade due to its transformation to the NayFeF3 perovskite phase, while Mg0.125 displays more reversible cycling and a reduced phase transformation. Yet, Mg incorporation also leads to kinetically limited Na extraction and a reduced overall capacity. These findings highlight the need for the continued compositional optimization of weberite cathodes to improve their structural stability while maximizing their energy density.

25 ENERGY STORAGE↗
DOE OSTI2020

Tutorial review on structure – dendrite growth relations in metal battery anode supports

This tutorial review explains the emerging understanding of the surface and bulk chemistry – electrochemical performance relations in anode supports (aka secondary current collectors, substrates, templates, hosts) for lithium, sodium and potassium metal batteries (LMBs, SMBs or NMBs, and KMBs or PMBs). In relation to each section, the possible future research directions that may yield both new insight and improved cycling behavior are explored. Representative case studies from Li, Na and K metal anode literature are discussed. The tutorial starts with an overview of the solid electrolyte interphase (SEI), covering both the “classic” understanding of the SEI structure and the “modern” insights obtained by site-specific cryogenic stage TEM analysis. Next, the multiple roles of supports in promoting cycling stability are detailed. Without an optimized support architecture, the metal–electrolyte interface becomes geometrically unstable at a lower current density and cycle number. Taking into consideration the available literature on LMBs, SMBs and KMBs, it is concluded that effective architectures are geometrically complex and electrochemically lithiophilic, sodiophilic or potassiophilic, so as to promote conformal electrochemical wetting of the metal during plating/stripping. One way that philicity is achieved is through support oxygen surface chemistry, which yields a reversibly reactive metal–support interface. Examples of this include the well-known oxygen–carbon moieties in reduced graphene oxide (rGO), as well as classic ion battery reversible conversion reaction oxides such as SnO 2 . Unreactive surfaces lead to dewetted island growth of the metal, which is a precursor to dendrites, and possibly to non-uniform dissolution. Surveying the literature on various Li, Na and K metal supports, it is concluded that the key bulk thermodynamic property that will predict electrochemical wetting behavior is the enthalpy of infinite solution (Δ sol H ∞ ) of the metal (solute) into the support (solvent). Large and negative Δ sol H ∞ promotes uniform metal wetting on the support surface, corresponding to relatively low plating overpotential. Positive Δ sol H ∞ promotes dewetted islands and a relatively high overpotential. This simple rule explains a broad range of studies on Li, Na and K metal – support interactions, including the previously reported correlation between mutual solubility and wetting.

25 ENERGY STORAGE↗
DOE OSTI2021

Thermodynamically Stable, Plasmonic Transition Metal Oxide Nanoparticle Solar Selective Absorbers towards 95% Optical-to-Thermal Conversion Efficiency at 750 °C

Generation 3 (Gen) concentrating solar power (CSP) systems requires a high operation temperature ≥750°C to increase the power-cycle efficiency towards ≥50%. However, such a high operation temperature poses a notable challenge to high temperature materials. On the receiver side, a critical challenge is the lack of solar selective absorbers that demonstrate both high optical-to-thermal energy conversion efficiency η therm approaching 95% and long-term thermal stability at >750°C in air. Existing solar absorbers either have limited η therm ≤89% or deteriorate significantly within 500 h at 750°C; some of these also require costly vacuum deposition for stringent thickness control. Therefore, it is highly desirable to simultaneously achieve η therm ~95% AND high thermal stability at 750°C for future generations of CSP receivers. This project has investigated and developed low-cost, highly scalable spray-coated transition metal oxide nanoparticle (NP) pigmented solar selective coatings on various types of Inconel alloy tube sections that are thermodynamically stable at 750-800°C in air, maintaining η therm >94.3%(93.2%) under a solar concentration ratio of C=1000 after 60 simulated day-night cycles at 750ºC (800ºC) (1 cycle=12 h at 750 or 800°C and 12h ramping down to 25°C). We have also coated up to 48 inches long receiver tubes for preliminary solar testing under Norwich Technology’s parabolic trough systems, demonstrating notably improved performance in solar heating compared to benchmark Pyromark 2500 coatings. Two key innovations have been developed in this project: (1) We achieved an unprecedented high thermal efficiency >94% by optimizing the d-band optical absorption spectra of transition metal ions, engineering their valences and stoichiometry; (2) We were able to maintain or even slightly increase the efficiency when operating at 750°C in air by engineering the interdiffusion of transition metal ions between the coating and the Inconel substrate to our advantage. Featuring high-temperature solar spectral selectivity and thermodynamic stability in air, as well as low-cost, highly scalable solution-chemical synthesis and spray coating techniques, this innovation in solar selective absorber technology alone accounts for nearly 40% of the targeted reduction in the levelized cost of electricity (LCOE) from the receiver, thermal energy storage, and operation & maintenance combined by 2030, as proposed by the U.S. Department of Energy. For a 110 MWe CSP power plant, this LCOE reduction from the solar selective coating alone transfers to ~$\$ $1.9 M increase in annual profit based on a sales price of 10¢/kWh (or an annual sale of $50 M/year). The high solar absorptance (~98%) NP pigment materials developed in this project are equally applicable to volumetric receivers either as a coating or as bigger ceramic microspheres after sintering, potentially benefiting Gen3 falling particle CSP technologies. All these contributions facilitate the larger scale deployment of CSP systems with energy storage capability to address the intermittency issue of solar energy towards dispatchable solar electricity, bridging the temporal gap between peak solar electricity production and peak electricity consumption. The collaboration with Norwich Technologies and Brayton Energy in this project will also facilitate the future commercialization of this new solar selective coating technology developed in this SIPS project.

14 SOLAR ENERGY↗
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↗
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↗