Engineering PapersSearch

SEARCH · Engineering Papers

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 109 records · Page 6

Vapor-liquid equilibrium thermodynamics of N2 + CH4 - Model and Titan applications

A thermodynamic model is presented for vapor-liquid equilibrium in the N2 + CH4 system, which is implicated in calculations of the Titan tropospheric clouds' vapor-liquid equilibrium thermodynamics. This model imposes constraints on the consistency of experimental equilibrium data, and embodies temperature effects by encompassing enthalpy data; it readily calculates the saturation criteria, condensate composition, and latent heat for a given pressure-temperature profile of the Titan atmosphere. The N2 content of condensate is about half of that computed from Raoult's law, and about 30 percent greater than that computed from Henry's law.

Thompson, W. R.

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

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

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

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

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.

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

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.

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

Colloidal Behavior of Plutonium Oxide in Concentrated Electrolyte Solutions

The Hanford Site in Washington State manages legacy high-level radioactive waste streams that display major chemistry and engineering challenges, including the high salt levels and pH values that correspond to conditions under which many classical concepts describing chemical reactivity cannot be applied. One particular challenge that needs to be tackled at Hanford is that Pu concentrations, [Pu], in the soluble phases of the tank wastes are higher than expected based on the solubility of crystalline PuO2, which is widely accepted to be caused by the formation of PuO2 colloid, consisting of nano- to submicron-sized particles (PuO2 NPs). Fundamental research underpinning the behavior of PuO2 NPs under conditions not only relevant to the Hanford tank waste but at high ionic strength in general is needed to reliably predict the chemical reactivity of PuO2 NPs and develop engineering solutions to safely and efficiently process high-level radioactive waste into forms suitable for long-term storage. In this work, we study the behavior of PuO2 NPs (particle size ~100 nm) under high ionic strength conditions by reacting it with highly concentrated (up to 5 M) salt solutions. We explore different electrolyte compositions to elucidate the impact of different anions (NO3-, Cl-, ClO4-, SO42-, C2O42-, CO32-) on the stability of PuO2 NP in the acidic and alkaline pH regime. PuO2 NP aggregation and precipitation as function electrolyte concentration is tracked by a combination of liquid scintillation counting, dynamic light scattering for determination of particle size distributions, and zeta potentials as a proxy for particle charge. At acidic pH, electrolytes containing non-coordinating anions, such as NaNO3, NaCl, and NaClO4 mostly stabilize PuO2 NPs over a large electrolyte concentration range, showing only subtle differences in their reactivity. Other electrolyte anions show a more pronounced effect on the PuO2 NP stability: SO42-, binds directly to the particles’ surface, reverses the particle charge, and precipitates the PuO2 NPs efficiently even at intermediate sulfate concentrations (>0.1 M). In contrast, C2O42- is found to lead to high [Pu] in solution, in the milli-molar range, even at mildly acidic pH (~4). Thermodynamic modeling of the dissolved Pu concentrations using PHREEQC is unable to predict the observed [Pu] in the acidic pH regime, supporting the influence of colloids in maintaining elevated [Pu]. It is noteworthy that the current thermodynamic databases do not include constants for colloidal Pu phases and cannot accurately predict many of the high ionic strength solutions relevant to this work. The mechanisms and models responsible for these observations will need further investigation in the future. At high pH values (~12), PuO2 NPs exhibits classical sol-gel chemistry, meaning that upon destabilization of the colloidal sol, for example by addition of concentrated NaOH, highly porous and viscid PuO2 coagulates are formed that consist of a three-dimensional network likely held together by physical interactions. The PuO2 NP coagulate shows no significant reversibility of the aggregation when contacted with concentrated brines; however, PuO2 NPs can be efficiently resuspended in solution by addition of diluted electrolytes, alkaline solutions containing high amounts of carbonate, or simple addition of water. Especially carbonate is shown to stabilize PuO2 NPs in solution at high pH, characterized by stable colloidal suspensions that are resistant against sedimentation during centrifugation. Thermodynamic modeling of the carbonate system was able to predict an increasing dissolved Pu concentration with increasing carbonate concentration. However, the model was profoundly sensitive to the fixed redox potential and does not include any thermodynamic constants for colloidal Pu species.

Neumann, Julia

Assessing Domains of Equilibration In Highly Metamorhposed Eucrites

This study assesses how petrographic textures and trace element distributions, combined with thermodynamic models, can act as indicators of metamorphic equilibrium and/or trace element mobility in highly metamorphosed eucrite, Elephant Moraine (EET) 90020.

Thermodynamic modelling

Advances in Modeling Capabilities for Critical Mineral Separation Technologies: A PrOMMiS Overview

This is an oral presentation at the TechConnect conference on the work developed by PrOMMiS. PrOMMiS builds on and extends capabilities developed within the Department of Energy’s (DOE) Institute for the Design of Advanced Energy Systems (IDAES), Integrated Platform, and Water Treatment Technoeconomic Assessment Platform (WaterTAP), which have been successfully leveraged by other Department of Energy research areas. The open-source toolkit facilitates validation, reproducibility, and accountability, allowing for easy extension of the framework to other systems. This talk presents an overview of the PrOMMiS capabilities, including unit model library, advances in thermophysical properties models, and capital cost libraries for simulation and optimization of mineral processing technologies. The PrOMMiS applications include (1) conceptual design and superstructure optimization for screening different process configurations and identifying promising technologies; (2) dynamic modeling and optimization to enable the creation of digital twins; (3) surrogate modeling tools to leverage data when predictive thermodynamic models are not currently available; (4) technical risk reduction via uncertainty quantification and robust optimization to identify process designs that are robust to process variability and uncertainties; and (5) deployment of uncertainty quantification tools to maximize knowledge gained from experimental campaigns, while reducing the number of experiments required

critical minerals and materials

Neural Networks for Prediction of Complex Chemistry in Water Treatment Process Optimization

Water chemistry plays a critical role in the design and operation of water treatment processes. Detailed chemistry modeling tools use a combination of advanced thermodynamic models and extensive databases to predict phase equilibria and reaction phenomena. The complexity and formulation of these models preclude their direct integration in equation-oriented modeling platforms, making it difficult to use their capabilities for rigorous water treatment process optimization. Neural networks (NN) can provide a pathway for integrating the predictive capability of chemistry software into equation-oriented models and enable optimization of complex water treatment processes across a broad range of conditions and process designs. Herein, we assess how NN architecture and training data impact their accuracy and use in equation-oriented water treatment models. We generate training data using PhreeqC software and determine how data generation and sample size impact the accuracy of trained NNs. The effect of NN architecture on optimization is evaluated by optimizing hypothetical black-box desalination processes using a range of feed compositions from USGS brackish water data set, tracking the number of successful optimizations, and testing the impact of initial guess on the final solution. Our results clearly demonstrate that data generation and architecture impact NN accuracy and viability for use in equation-oriented optimization problems.

Dudchenko, Alexander V

Unravelling Microstructure Selection in an Additively Manufactured Eutectic High‐Entropy Alloy

High-entropy alloys (HEAs) are promising candidates for advanced structural applications due to their excellent mechanical properties. Additive manufacturing (AM), with its rapid solidification conditions, enables the creation of unique nonequilibrium microstructures. To fully leverage the synergy between AM and HEAs, understanding how processing affects structure and properties is essential. Here, how solidification rate influences microstructure evolution and phase transformation pathway in laser additively manufactured AlCrFe2Ni2 eutectic HEAs is investigated. By increasing the laser scan speed and hence the solidification rate, distinct solidification modes evolving from coupled eutectic to anomalous eutectic and eventually to single-phase solidification are revealed. These transitions result in distinct microstructures and a wide range of mechanical properties. Thermodynamic modeling and molecular dynamics simulations reveal that low cooling rates allow for sufficient atomic diffusion and phase separation, facilitating coupled eutectic growth. In contrast, rapid cooling suppresses diffusion and destabilizes the solid–liquid interface, promoting anomalous or single-phase solidification. This integrated experimental and computational approach provides a multiscale understanding of solidification mechanisms in HEAs and underscores how kinetic effects can over-ride thermodynamic predictions under nonequilibrium conditions. Furthermore, these results demonstrate that AM can serve as a powerful tool to design HEAs with tailored microstructures and properties.

36 MATERIALS SCIENCE

First-Principles Thermodynamic Assessments of Sr-Containing Secondary Phase Formation in La1-xSrxMnO3±δ Perovskites for Solid Oxide Cell Applications

Sr-secondary phase formation is a potentially significant degradation mode threatening solid-oxide cell (SOC) commercial viability. A first-principles thermodynamic study was performed for rhombohedral perovskite (La1-xSrx) MnO3±δ (LSM) to assess its stability against Sr secondary phase formation in SOC applications. In this work, the Sr secondary phase formation reaction free energies were determined by combining ab initio lattice dynamics calculations for the solid phases and an ab initio thermodynamics approach for the gas phases. Furthermore, this approach goes beyond previous thermodynamic modeling studies by integrating first-principles based point-defect equilibria into the analyses. The modeling results indicate an increased tendency to form SrO oxide from LSM upon decreasing the oxygen partial pressure. Additionally, enhancing factors to form the Sr-related secondary phase from the associated SrO activity in LSM are further quantified by considering the equilibrium of SrO reacting with contaminant gas species as a function of temperature and gas pressure.

Defect and phase stability

Improved Upper Ocean/Sea Ice Modeling in the GISS GCM for Investigating Climate Change

This project built on our previous results in which we highlighted the importance of sea ice in overall climate sensitivity by determining that for both warming and cooling climates, when sea ice was not allowed to change, climate sensitivity was reduced by 35-40%. We also modified the Goddard Institute for Space Studies (GISS) 8 deg x lO deg atmospheric General Circulation Model (GCM) to include an upper-ocean/sea-ice model involving the Semtner three-layer ice/snow thermodynamic model, the Price et al. (1986) ocean mixed layer model and a general upper ocean vertical advection/diffusion scheme for maintaining and fluxing properties across the pycnocline. This effort, in addition to improving the sea ice representation in the AGCM, revealed a number of sensitive components of the sea ice/ocean system. For example, the ability to flux heat through the ice/snow properly is critical in order to resolve the surface temperature properly, since small errors in this lead to unrestrained climate drift. The present project, summarized in this report, had as its objectives: (1) introducing a series of sea ice and ocean improvements aimed at overcoming remaining weaknesses in the GCM sea ice/ocean representation, and (2) performing a series of sensitivity experiments designed to evaluate the climate sensitivity of the revised model to both Antarctic and Arctic sea ice, determine the sensitivity of the climate response to initial ice distribution, and investigate the transient response to doubling CO2.

Source record

Improved Upper Ocean/Sea Ice Modeling in the GISS GCM for Investigating Climate Change

This project built on our previous results in which we highlighted the importance of sea ice in overall climate sensitivity by determining that for both warming and cooling climates, when sea ice was not allowed to change, climate sensitivity was reduced by 35-40%. We also modified the GISS 8 deg x lO deg atmospheric GCM to include an upper-ocean/sea-ice model involving the Semtner three-layer ice/snow thermodynamic model, the Price et al. (1986) ocean mixed layer model and a general upper ocean vertical advection/diffusion scheme for maintaining and fluxing properties across the pycnocline. This effort, in addition to improving the sea ice representation in the AGCM, revealed a number of sensitive components of the sea ice/ocean system. For example, the ability to flux heat through the ice/snow properly is critical in order to resolve the surface temperature properly, since small errors in this lead to unrestrained climate drift. The present project, summarized in this report, had as its objectives: (1) introducing a series of sea ice and ocean improvements aimed at overcoming remaining weaknesses in the GCM sea ice/ocean representation, and (2) performing a series of sensitivity experiments designed to evaluate the climate sensitivity of the revised model to both Antarctic and Arctic sea ice, determine the sensitivity of the climate response to initial ice distribution, and investigate the transient response to doubling CO2.

Source record

The sensitivity of a one-dimensional thermodynamic sea ice model to changes in cloudiness

A thermodynamic sea ice-lead model is used to assess the importance of cloud cover changes to modeled ice thickness. For regions of either permanent multiyear ice or seasonal sea ice, the cloud amount variations have relatively little impact. However, for regions where the presence of summer ice is variable from year to year, the predicted ice thickness is strongly dependent on cloud cover. In general, with a snow covered surface, decreased cloud leads to surface cooling while increased cloud gives rise to a surface warming. For a melting bare ice surface, the reverse occurs. The ice model response time is too long for interannual variations in cloud amount to explain interannual variations in ice thickness and extent. Nevertheless, the implication of the results is that numerical modeling of sea ice distribution requires accurate cloud data or cloud prediction and that trends in cloud cover may lead to significant perturbations in sea ice extent and thickness.

Shine, K. P.