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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
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↗
NASA NTRSDate not supplied

Liquefaction of Cryogenic Fluids for Production and Storage of Commodities on Extra-Terrestrial Surfaces

The Cryogenic Fluid In-situ Liquefaction for Landers (CryoFILL) activity has been investigating concepts for the liquefaction of cryogenic fluids produced on the surface of the Moon and Mars. CryoFILL consists of four activities: liquid nitrogen liquefaction testing that was completed in 2019, liquid oxygen liquefaction testing planned to begin in early 2022, fiber-optic sensor temperature measurement system development, and lightweight vacuum jacketed systems designed for the atmosphere of Mars. Additional activities include the development and validation of numerical modeling systems at three stages: thermodynamic models of the liquefaction process for basic energy/mass estimates, nodal models to predict liquefaction tank system level performance, and computational fluid dynamics to assess fluid phenomena occurring within the tank, specifically condensation and stratification within the ullage. These developments form a combined liquefaction and storage solution to support Lunar and Martian exploration. Key testing results from liquid nitrogen testing as well as plans for liquid oxygen testing will be discussed including key sensitivities from analytical evaluation of completed test results. Progress on the development and improvements on fiber optic sensor testing and developments will also be provided. Finally, work on the progress of the lightweight vacuum jacketed systems being co-developed by NASA and industry will be summarized.

Liquefaction↗
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↗
NASA NTRS2005

Development of a Stirling System Dynamic Model With Enhanced Thermodynamics

The Stirling Convertor System Dynamic Model developed at NASA Glenn Research Center is a software model developed from first principles that includes the mechanical and mounting dynamics, the thermodynamics, the linear alternator, and the controller of a free-piston Stirling power convertor, along with the end user load. As such it represents the first detailed modeling tool for fully integrated Stirling convertor-based power systems. The thermodynamics of the model were originally a form of the isothermal Stirling cycle. In some situations it may be desirable to improve the accuracy of the Stirling cycle portion of the model. An option under consideration is to enhance the SDM thermodynamics by coupling the model with Gedeon Associates Sage simulation code. The result will be a model that gives a more accurate prediction of the performance and dynamics of the free-piston Stirling convertor. A method of integrating the Sage simulation code with the System Dynamic Model is described. Results of SDM and Sage simulation are compared to test data. Model parameter estimation and model validation are discussed.

Regan, Timothy F.↗
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↗
NASA NTRS1989

Climatic equilibrium of the atmospheric convective boundary layer over a tropical ocean

A radiative-convective boundary layer model was developed by coupling a thermodynamic model of a partially mixed convective boundary layer (CBL) with a radiation model, and energy balance constraints were used to study coupled boundary layer (CBL) equilibrium over three timescales (about 1 day, about 10 days, and more than 100 days). It is shown that the variation in cloud top decreases with greater coupling to the atmosphere and the ocean. The slope of the latent heat flux with increasing SST decreases with more tropospheric coupling, and reverses sign with a coupled ocean.

Ridgway, W.↗
DOE OSTI2022

A modified thermodynamic sea ice model and its application

This paper presents a modified thermodynamic sea ice model based on the Winton’s three-layer model framework. Several improvements focused on the vertical thermodynamics have been made in the model. Results from a series of one-dimensional experiments show that equilibrium ice thickness in our modified model is increased by about 45 cm when compared with the original Winton’s model. Sensitivity tests indicate that all modifications mentioned above contribute to this change of ice thickness, among which the increase of ice layer has the most significant effect followed by the sea ice conductivity parameterization. The modified thermodynamic sea ice model is then coupled into the sea ice component SIS in GFDL MOM4 to exert long-term integration experiments forced by CORE2.0 data. Results indicate that the modified model shows significant improvement in simulating the Arctic sea ice, including an increase in both the sea ice thickness and volume over the whole Arctic region. This confirms the above founding from 1-D simulations.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1992

Modeling of heavy-gas effects on airfoil flows

Thermodynamic models were constructed for a calorically imperfect gas and for a non-ideal gas. These were incorporated into a quasi one dimensional flow solver to develop an understanding of the differences in flow behavior between the new models and the perfect gas model. The models were also incorporated into a two dimensional flow solver to investigate their effects on transonic airfoil flows. Specifically, the calculations simulated airfoil testing in a proposed high Reynolds number heavy gas test facility. The results indicate that the non-idealities caused significant differences in the flow field, but that matching of an appropriate non-dimensional parameter led to flows similar to those in air.

Drela, Mark↗
DOE OSTI2023

Thermodynamic Band Gap Model for Photoinduced Phase Segregation in Mixed-Halide Perovskites

Provided is a comprehensive description of a band gap thermodynamic model, which predicts and explains key features of photosegregation in lead-based, mixed-halide perovskites. The model gives a prescription for illustrating halide migration driven by photocarrier energies. Where possible, model predictions are compared with experimental results. Free energy derivations are provided for three assumptions: (1) halide mixing in the dark, (2) a fixed number of photogenerated carriers funneling to and localizing in low band gap inclusions of the alloy, and (3) the statistical occupancy of said inclusions from a bath of thermalized photocarriers in the parent material. Model predictions include excitation intensity (I exc )-dependent terminal halide stoichiometries (x terminal ), temperature-independent x terminal , excitation intensity thresholds (I exc,threshold ) below which photosegregation is suppressed, and reduced segregation in nanocrystals as compared to thin films. Here, the model also offers insight into kinetically manipulating photosegregation rates via control of underlying mediators and rationalizes asymmetries in forward and reverse photosegregation rate constants/activation energies. What emerges is a cohesive framework for understanding ubiquitous photosegregation in mixed-halide perovskites and a rational basis by which to manage the phenomenon.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗