Search NASA⌕ Search

SEARCH · Search NASA

Results for “Thermodynamic Modeling”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

On the competitive antagonism effect in combined chloride-sulfate attack: A numerical exploration

Highlights: • A chemo-physical-mechanical method based on thermodynamic modelling is implemented. • The contradicting findings from reported experiments are numerically explored. • The competitive antagonism effect is time-dependent rather than a conclusive term. • The competitive antagonism effect is quantifiable with a novel assessment criterion. It is of tremendous difficulty to draw a conclusive verdict on the complex process of mutual interference regarding the combined chloride-sulfate attack on concrete, where many have explored the topic, but contradicting results have been reported. In this paper, by taking advantages of the chemo-physical-mechanical method, based on the robust technique of thermodynamic modelling, the competitive antagonism effect in combined chloride-sulfate attack is investigated. Firstly, a range of reported tests are numerically modelled to demonstrate the effectiveness of the employed method in interpreting those opposing and diverse experimental outcomes. More importantly, a fresh conjecture incorporating a novel assessment criterion for the competitive antagonism effect is proposed and supported through a carefully designed illustrative example. Based on the numerical exploration, unique findings and implications for future engineering practices are revealed.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Dual-module humidity pump for efficient air dehumidification: Demonstration and performance limitations

Condensation dehumidification in conventional air conditioning technologies is energy-intensive, accounting for up to 50% of building cooling energy used in some climates. Selective vacuum membrane dehumidification (VMD) is one of the leading alternative dehumidification technologies due to its potential for significant energy savings, and the “dual-module humidity pump” is one of the most promising VMD concepts. Here, this work is the first to provide experimental proof-of-concept for the dual-module humidity pump system and provides the first thermodynamic modeling framework that accounts for realistic steady-state operating limitations, both of which are lacking in the current literature. Additionally, this work is the first to provide a system design solution that overcomes practical challenges associated with air accumulation in the vacuum channels. The experimental results in this work show that the current prototype can remove up to 45% of the water vapor in the air stream, and the vapor pressure difference in the vapor rejection module needs to be approximately 2–4 times greater than that of the dehumidification module in order to maintain balanced mass transfer. The thermodynamic model applied to typical air conditioning conditions shows that the ideal dehumidification (latent) COPs can reach up to 40, but practical COPs are limited to approximately 10. Furthermore, the model shows that the overall energy efficiency increases as the membrane air selectivity increases, though this improvement gradually starts to diminish when the membrane selectivity is increased above 10,000.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Influence of ethylenediaminetetraacetic acid on the long-term oxidation state distribution of plutonium

Here, spectrophotometry was used to study the effect of EDTA on plutonium oxidation state distribution as a function of time, pH, and ligand-to-metal ratio (L/M) under anoxic conditions. Novel Pu(V)-EDTA absorption bands were identified at 571, 993, 1105, and 1150 nm with molar absorption coefficients of 15 ± 1, 6 ± 1, 10 ± 1, and 10 ± 1 cm –1 M –1 , respectively. Pu(V)-EDTA spectral changes occurred at L/M < 1, indicating only Pu V O 2 (EDTA) 3- formed with logK = 3.6 ± 0.3. Time-resolved experiments showed EDTA drastically increased the Pu(V/VI) reduction rate, which we propose is driven by amine lone-pair electron donation and the oxidative decarboxylation of EDTA. Oxidation of Pu(III)-EDTA to Pu(IV)-EDTA occurred on a slower time scale (110–237 days) than previously reported (<15 min) and is hypothesized to be radiolysis driven. Pu(V/VI)-EDTA and Pu(III)-EDTA both approached Pu(IV)-EDTA stabilization over time, yet Pu(V/VI)-EDTA solubility data was ≥ 1.0 log 10 units higher than predicted by Pu(IV)-EDTA solubility models, indicating that current thermodynamic models are incomplete. Ultimately, the data show EDTA preferentially stabilizes Pu(IV) over time regardless of initial oxidation state, but Pu(V)-EDTA can persist under environmentally-relevant conditions, emphasizing the need to continue investigating redox reactions, speciation, and behavior of these complexes to support the transuranic waste disposal and surface remediation/containment efforts.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2023

Ab initio modeling and thermodynamics of hydrated plutonium oxalates

An ab initio study on the plutonium oxalate hydrates: Pu 2 (C 2 O 4 ) 3 ·10H 2 O and Pu(C 2 O 4 ) 2 ·6H 2 O, using PBE exchange-correlation with D3 dispersion correction and Hubbard correction for the plutonium atoms, was performed and compared to experimental vibrational spectral and thermodynamic property values. Here we demonstrated that this technique can accurately predict the experimental infrared spectra Pu(III) oxalate hydrate, as well as the Raman peak of PuO 2 (used to calculate the thermodynamic properties of the oxalates). For Pu(IV) oxalate hydrate we found that our predicted structure agreed qualitatively with PXRD measurements, the only available experimental determination of the structure. Using this method at standard temperature and pressure, we predicted standard enthalpies of formation of -6,755 kJ mol -1 and -3,923 kJ mol -1 and standard Gibbs free energy of formation of -5,899 kJ mol -1 and -3,386 kJ mol -1 for Pu 2 (C 2 O 4 ) 3 ·10H 2 O and Pu(C 2 O 4 ) 2 ·6H 2 O, respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Thermal properties and Thermodynamic Equilibrium Modeling of Cementitious Waste Forms

INTRODUCTION Cementitious matrices are used in the US DOE complex and worldwide to solidify aqueous radioactive, hazardous, mixed salt solutions, and sludges to meet low-level radioactive waste (LLW) disposal requirements. Blended formulations are used for most waste form mix designs. Substitution of pozzolans for Class F fly ash has the potential to alter the processing properties and stabilization properties of the resulting waste forms because they add chemical and mineralogical complexity to the final material. The Savannah River Site saltstone waste form was selected as the test case material. CemGEMS software was chosen as the model for predicting hydrated cementitious phases assemblages and phase evolution. Isothermal calorimetry was used as the method for evaluating processing and properties of fresh, uncured waste forms. The test cases consisted of reference case saltstone, and five natural pozzolan substituted saltstone mixes.

Bustamante, Michael E. [Savannah River National La↗
DOE OSTI2023

Modeling the Photoelectrochemical Evolution of Lead-Based, Mixed-Halide Perovskites Due to Photosegregation

Lead-based, mixed-halide perovskites such as methylammonium lead iodide-bromide [MAPb(I 1-x Br x ) 3 ] undergo anion photosegregation under illumination. This is observed as low band gap photoluminescence from photogenerated iodine-rich domains due to favorable band offsets that induce carrier funneling into them. Unfortunately, theoretical rationalizations of mixed-halide photosegregation are complicated by biases inherent to photoluminescence-based observations. Recent compositionally-weighted X-ray diffraction (XRD) measurements now reveal broad distributions of photosegregated stoichiometries not captured by existing photosegregation models. To better bridge experiment and theory, we perform kinetic Monte Carlo (KMC) simulations of photosegregation within the context of a band gap-based thermodynamic model, which has previously accounted for numerous experimental observations. Our KMC simulations are modified to consider high carrier density Fermi-Dirac statistics that result from carrier funneling and accumulation within photosegregated I-rich domains. Obtained KMC results reproduce broad x terminal distributions seen experimentally and illustrate how they are characterized by a central, heavily I-enriched stoichiometry. I-rich domain “drifting” during photosegregation rationalizes the long photosegregation timescales seen experimentally with drifting simultaneously producing a wake of variable stoichiometry I-rich inclusions that form the lion’s share of stoichiometric heterogeneities seen in compositionally-weighted XRD measurements. These simulations and accompanying rationalizations further reveal a general criterion for realizing favorable free energies to induce demixing. Central to the criterion is the statistical occupation of low gap inclusions in the parent alloy by excitations. Furthermore, the resulting model thus provides a general framework for conceptualizing mixed-halide perovskite light and temperature sensitivities, mediated by photocarriers.

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

Fundamental Issues Identified for Thermodynamic Description of Molten Salt Systems

Thermodynamic modeling of the molten salt systems is critical in the design of advanced nuclear power plants, materials recycling, and many other important clean energy applications. As one of the most capable computational thermodynamic approaches, the CALPHAD (Calculation of Phase Diagrams) method has been widely used for multicomponent thermodynamic database development and further support more physics-based modeling. Although considerable CALPHAD modeling efforts for molten salt systems have been made, we found that no evaluation has been performed regarding the unary thermodynamic models for the molten salt systems despite the fact that more than one pure substance database is available for the CALPHAD modeling. Therefore, the thermodynamic descriptions of the pure salts should be critically evaluated to ensure a high-fidelity molten salt thermodynamic database. In this work, we comprehensively analyzed some selected molten salts using different available thermodynamic descriptions from two databases as an example. One is the SSUB database released by the SGTE (Scientific Group Thermodata Europe), and the other is the FactPS database from the FactSage. Finally, the significant discrepancies observed in the thermodynamic description between the existing CALPHAD pure substance databases and experiments call for an urgent effort to improve.

42 ENGINEERING↗
DOE OSTI2023

Assessing the Robustness of Arctic Sea Ice Bi‐Stability in the Presence of Atmospheric Feedbacks

Abstract Arctic sea‐ice loss is influenced by multiple positive feedbacks, sparking concerns of accelerated loss in the coming years or even a tipping point, where a sea‐ice equilibrium disappears at a given CO 2 value and sea ice rapidly evolves to a new steady state. Such a tipping point would imply a bi‐stability of the Arctic climate—where multiple steady‐state Arctic climates are possible at the same CO 2 value. Previous works have sought to establish the existence of bi‐stability using a range of models, from zero‐dimensional sea ice thermodynamic models to fully coupled global climate models, with conflicting results. Here, we present a new model of the Arctic that includes both sea‐ice thermodynamics and key atmospheric feedbacks in a simple framework. We exploit the model's simplicity to identify physical mechanisms that control the timing and extent of sea‐ice bi‐stability, and the abruptness of ice loss. We show that longwave radiation feedbacks can have a strong influence on Arctic surface climate from atmospheric temperature increases alone, even without major contributions from clear‐sky moisture or convective clouds suggested previously. While winter sea‐ice bi‐stability is robust to changes in uncertain model parameters in this study, summer sea ice is more sensitive. Finally, our model indicates that positive feedbacks may modulate the CO 2 threshold of sea‐ice loss and the width of bi‐stability much more strongly than the abruptness of loss. These results lead to a comprehensive understanding of the conditions that favor Arctic sea‐ice bi‐stability, particularly the role of atmospheric feedbacks, in both future and past climates.

Meteorology & Atmospheric Sciences↗
DOE OSTI2025

Predicting Multicomponent Aqueous Phase Adsorption Equilibria of Organic Acids Using the Generalized Brunauer–Emmett–Teller Isotherm Model

Here, to support process development of adsorptive separation of organic acids, this work presents a thermodynamic modeling methodology to predict multicomponent aqueous phase organic acid adsorption equilibria from single-component adsorption isotherms using the generalized Brunauer–Emmett–Teller isotherm model (gBET). With the organic acid fugacities rigorously accounting for the aqueous phase organic acid solution chemistry and solution nonideality, gBET precisely represents the single-component adsorption isotherms and accurately predicts the binary and ternary mixed-acid adsorption equilibria for the ranges of initial pH (∼3–7), acid concentration (100–400 mmol/L), and temperature (298.15–328.15 K) with less than 10% average absolute relative deviation. In addition, gBET with pH-independent parameters provides insights into the underlying adsorption phenomena, including the adsorbate loadings and compositions in the monolayer and subsequent layers under varying initial pH, temperature, concentration, and composition. The gBET model predictions outperform the predictions from the overloading model and the Ideal Dilute Solution Theory.

09 BIOMASS FUELS↗
DOE OSTI2020

A thermodynamics-based damage model for the non-linear mechanical behavior of SiC/SiC ceramic matrix composites in irradiation and thermal environments

A damage model is developed and validated with experimental data for the non-linear mechanical behavior of SiC/SiC composite materials in nuclear applications. Cyclic thermal and mechanical loading associated with neutron irradiation effects of these composites leads to wide-spread and progressive micro-cracking that leads to loss of thermal conductivity and further enhancement of thermo-mechanical damage. A physics-based model of wide-spread micro-cracking is developed within the thermodynamic framework of continuum damage mechanics. Evolution equations for damage parameters that describe the growth of continuum damage are developed, where the material variables are obtained from experiments. The model novelty is in coupling mechanical, thermal, and irradiation damage through a consistent thermodynamic framework, including loss of thermal conductivity due to the evolution of mechanically induced micro-cracks. A number of thermo-mechanical experiments were conducted to confirm model assumptions. The model is shown to be validated with out-of-pile experiments, and then implemented using commercial finite element code COMSOL to the fuel cladding problem with normal and high radiation dose cases.

Materials Science↗
DOE OSTI2022

Development of a Thermodynamic Database for Corrosion in Chloride MSRs

Control of the salt redox condition is the salient method for corrosion inhibition in molten salts. For chloride salt-fueled molten salt reactors of primary concern is the U(IV)/U(III) ratio, influenced directly by the Cl2 activity in the salt. Although thermodynamic assessments of MSR-relevant salts abound, none have included consideration of UCl4 which is essential to correctly describe the salt redox condition. This work reports the assessment of the NaCl-KCl-MgCl2 base salt with UCl3-UCl4 and corrosion chlorides involving Cr, Fe, and Ni as minor components and the resultant thermodynamic models. The effort utilized available thermodynamic measurements and phase equilibria complemented by differential scanning calorimetry determinations of phase equilibria. Example calculations explore the temperature and composition space to understand its influence on corrosion. The foregoing models are incorporated into an update of the Molten Salt Thermal Properties Database – Thermochemical which is publicly available.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Manufacturing of Al-Ce/Metal-Matrix-Composites (MMC) for Harsh Environments

An integrated modeling/simulation and experimental approach will be used to design and synthesize Al-Ce/MMC for use in key technologies where improved energy efficiency and performance is harsh environments is required. Modeling and simulation work including thermodynamic modeling, DFT calculations, and continuum modeling of casting processes and mechanical behavior, will be combined with lab-scale synthesis to identify potential matrix and reinforcing phase compositions for large-scale synthesis. The experimental synthesis work will rely on lab-scale casting and gas atomization to identify the optimal matrix compositions for pilot/commercial-scale synthesis is Al-Ce/MMC. Selected Al-Ce/MMC will be synthesized by commercial scale melt processing, powder metallurgy and thermomechanical processing. The corrosion resistance, mechanical properties as a function of temperature, and thermal conductivity of the composites will be characterized to validate target properties are met or exceeded. This work will focus on two selected applications; however, the design framework will have broad applicability for composites for numerous harsh service condition applications. We will have synthesized Al-Ce/MMC that exhibit better performance in harsh environments than currently used materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Dynamic Modeling of Latent Heat Thermal Battery

A dynamic thermodynamic model is constructed in Modelica for the thermal battery portion of a heat pipe integrated thermal battery and is demonstrated using a shakedown test. The model is also leveraged to inform aspects of an experiment configuration. The theory section discusses modeling assumptions of the thermal battery system, establishes how melting and fusing is calculated, and how internal convection is modeled within the thermal battery. The shakedown test uses three charge-standby-discharge-standby cycles in which the charging cycles are 10 hours in length and the discharge cycles are 8 hours. System behavior including heat loss, temperature distributions and divergences, and internal convection behavior are presented. The model is then leveraged to show how experiment design can be informed through testing the model for insulation thickness impact on heat loss, round trip efficiency calculation, and heat tracing control.

25 ENERGY STORAGE↗
DOE OSTI2022

Evaluation of the energy balance of chemical looping combustion of solid fuels using CuO-based oxygen carriers

Chemical looping combustion (CLC) can inherently capture the CO2 generated during thermal energy production. When using CuO-based oxygen carriers, the reactions in both the fuel reactor (FR) and air reactor (AR) are exothermic, facilitating autothermal operation. To prevent the oxygen carrier from becoming agglomerated while ensuring a sufficient oxygen release rate, it is necessary to predict the FR temperature and remove an appropriate amount of heat from the AR when burning high-rank fuels. Here in this work, the energy balance of a CuO-based CLC system firing various solid fuels (biomass, coal, and petroleum coke) was investigated by establishing a thermodynamic model using FactSage. The energy distribution characteristics of the FR were quantified based on a virtual five-step thermodynamic sequence. The effects of fuel properties, oxygen carrier properties, and operating conditions on the energy balance were evaluated systematically. The results indicate that the temperature difference between the FR and AR is not very sensitive to the heating value of the fuel, while the CuO loading of the oxygen carrier, the heat capacity of the support material, and the flow rate of the fluidizing gas have significant impact. The insights obtained in this work will help improve the design and heat management of CLC.

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

Effect of Mn on eutectic phase equilibria in Al-rich Al-Ce-Ni alloys

Microstructural analysis of additively manufactured (AM) Al-Ce-Ni-Mn alloys has identified phases not predicted from existing ternary liquidus projections in the Al-Ce-Ni system. Because the rapid cooling rate of AM is orders of magnitude above that of traditional casting, it is unclear if these additional phases arose from the non-equilibrium processing conditions of AM, a drastic shift in phase stability in the system due to the addition of 1 wt% Mn, or some combination of these two influences. The phases and microstructure of cast samples of Al-Ce-Ni and Al-Ce-Ni-Mn alloys were characterized for several annealing conditions which revealed the equilibrium phases at different temperatures. Phase analysis confirmed that minute levels of Mn substituted for Ni in the system drastically shifts the liquidus projection in the Al-rich corner of the ternary phase diagram such that the eutectic Al 3 Ni phase is suppressed in favor of the Al 23 Ni 6 (Ce,Mn) 4 phase. Further addition of Mn promotes the formation of Al 20 Mn 2 Ce and Al 10 Mn 2 Ce phases. The phase analysis data was then used to improve the CALPHAD modeling of the liquidus projection and isothermal sections for the Al-rich Al-Ce-Ni-Mn quaternary system. Thermodynamic modeling and experimental analysis on phases in the AM sample of Al-Ce-Ni with Mn confirmed that the phases present are consistent with Mn-containing Al-Ce-Ni cast samples. Here, this investigation demonstrates the potential for using secondary alloying elements to drastically alter phase stability and microstructure in alloy systems.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Measurements of Density of Liquid Oxides with an Aero-Acoustic Levitator

Densities of liquid oxide melts with melting temperatures above 2000 °C are required to establish mixing models in the liquid state for thermodynamic modeling and advanced additive manufacturing and laser welding of ceramics. Accurate measurements of molten rare earth oxide density were recently reported from experiments with an electrostatic levitator on board the International Space Station. In this work, we present an approach to terrestrial measurements of density and thermal expansion of liquid oxides from high-speed videography using an aero-acoustic levitator with laser heating and machine vision algorithms. The following density values for liquid oxides at melting temperature were obtained: Y 2 O 3 4.6 ± 0.15; Yb 2 O 3 8.4 ± 0.2; Zr 0.9 Y 0.1 O 1.95 4.7 ± 0.2; Zr 0.95 Y 0.05 O 1.975 4.9 ± 0.2; HfO 2 8.2 ± 0.3 g/cm 3 . The accuracy of density and thermal expansion measurements can be improved by employing backlight illumination, spectropyrometry and a multi-emitter acoustic levitator.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Scaling dynamics in low-salt-rejection reverse osmosis for high-salinity produced water desalination: Mechanistic modeling and membrane autopsy

Membrane scaling remains a critical barrier to the reliable operation of desalination systems, particularly for hypersaline produced water (PW) treatment. This study fills the knowledge gap of autopsy-based model validation for PW desalination by elucidating scaling mechanisms in a Low-Salt-Rejection Reverse Osmosis (LSRRO) system through the integration of pilot-scale experimentation and complementary modeling approaches. A semi-empirical modeling framework was developed and applied to a multistage pilot LSRRO system equipped with nanofiltration and RO membranes treating high-salinity PW from the Permian Basin. Water quality analysis showed that total dissolved solids decreased from ~130,000 mg/L to ~1900 mg/L in the permeate, then further reduced to ~300 mg/L by a second-pass RO. Two different thermodynamic modeling approaches were evaluated: the first extends the LSRRO framework by incorporating system complexity and scaling phenomena, whereas the second method explicitly captures concentration polarization in localized supersaturation. Both methods illustrate the tendency for carbonate and sulfate scaling throughout the stages. Membrane autopsies revealed a silica-dominated deposit matrix, localized CaSO 4 at Stage 2, and minor barite/celestite despite their prominence in model predictions. Quantum-chemical calculations indicated silica scaling can be rationalized by favorable adsorption of H 4 SiO 4 on Fe-oxide surfaces (ΔG ≈ −44 kJ/mol), providing a kinetic pathway for interfacial inorganic polymerization even when bulk equilibrium predictions are conservative. Overall, the thermodynamic scaling modeling and membrane autopsy revealed heterogeneous, localized deposits with limited impact on LSRRO performance, while quantum analysis rationalized the thermodynamically unfavorable precipitation formation under bulk equilibrium, reconciling model–autopsy discrepancies. These insights support targeted pretreatment and silica-specific antiscalants to extend membrane lifetime and increase recovery, providing a transferable framework for hypersaline water desalination systems. The combined experimental–computational approach provides new mechanistic insight into scaling in hypersaline membrane systems and establishes a transferable framework for predicting and mitigating scaling in next-generation desalination technologies.

Low-salt-rejection reverse osmosis↗