Engineering Papers⌕ Search

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Calcium and Magnesium Distribution Between Sulfide and Silicate Melts: Thermodynamic Modeling and Insights into Mercury's Mantle-Crust Differentiation

Mercury stands out as the terrestrial planet with the lowest degree of oxidation among all planets of our solar system. It is characterized by a high abundance of sulfur, which is expected to replace oxygen in silicate melts. This sulfur enrichment could lead to sulfide saturation in its mantle, potentially resulting in the formation of sulfides rich in calcium and magnesium. Utilizing experimental data from the literature, we constructed a thermodynamic model to explore the partitioning of calcium and magnesium between sulfide and silicate melts. This model aims to improve our comprehension of the chemical equilibria driving the process of mantle-crust differentiation on Mercury. Our models indicate that increased temperature, decreased pressure, and low concentrations of iron, oxygen, and carbon within the sulfide phase, along with low sulfur concentrations in the silicate phase, favor calcium’s partitioning into the sulfide. Additionally, magnesium partitioning into the sulfide melt is favored with decreased pressure, decreased temperature, low concentrations of iron and calcium within the sulfide, and low concentrations of sulfur within the silicate. Our results suggest that up 9.4% and 5.6% Ca and Mg, respectively, can be present in sulfides equilibrated with sulfide-saturated silicate melts during mantle-crust differentiation. Further investigation will focus on the distribution of sulfide phases within the mantle and crust to better constrain Mercury’s internal structure.

Taylor McCombs↗
NASA NTRSDate not supplied

Thermodynamic Modeling of the Vapor in Equilibrium With Apollo 17 Basalts.

Lunar mare basalts sampled during the Apollo 17 mission provide insight to the chemical evolution of mare basalt magmas and their associated vapor as these samples erupted and cooled rapidly at the lunar surface. We focus on minerals present on vesicle and vug surfaces in Apollo 17 basalts, such as 71036, which was opened as part of the Apollo Next Generation Sample Analysis (ANGSA) program. The mineral assemblage observed includes native Fe, SiO2, and a Mg-P mineral, possibly merrillite [1, 2]. The pet-rographic context of these phases suggests that they may have been deposited directly from the vapor phase after eruption. Here we use thermodynamic modeling to assess the plausibility of the vapor-phase deposition hypothesis.

D. L. Thompson↗
NASA NTRS2005

Aqueous Alteration and Hydrogen Generation on Parent Bodies of Unequilibrated Ordinary Chondrites: Thermodynamic Modeling for the Semarkona Composition

Ordinary chondrites are the most abundant class of meteorites that could represent rocky parts of solar system bodies. However, even the most primitive unequilibrated ordinary chondrites (UOC) reveal signs of mild alteration that affected the matrix and peripheral zones of chondrules. Major chemical changes include oxidation of kamacite, alteration of glass, removal of alkalis, Al, and Si from chondrules, and formation of phases enriched in halogens, alkalis, and hydrogen. Secondary mineralogical changes include formation of magnetite, ferrous olivine, fayalite, pentlandite, awaruite, smectites, phosphates, carbonates, and carbides. Aqueous alteration is consistent with the oxygen isotope data for magnetite. The presence of secondary magnetite, Ni-rich metal alloys, and ferrous silicates in UOC implies that H2O was the oxidizing agent. However, oxidation by H2O means that H2 is produced in each oxidative pathway. In turn, production of H2, and its redistribution and possible escape should have affected total pressure, as well as the oxidation state of gas, aqueous and mineral phases in the parent body. Here we use equilibrium thermodynamic modeling to explore water-rock reactions in UOC. The chemical composition of gas, aqueous, and mineral phases is considered.

Zolotov, M. Y.↗
NASA NTRS1987

A coupled dynamic-thermodynamic model of an ice-ocean system in the marginal ice zone

Thermodynamics are incorporated into a coupled ice-ocean model in order to investigate wind-driven ice-ocean processes in the marginal zone. Upswelling at the ice edge which is generated by the difference in the ice-air and air-water surface stresses is found to give rise to a strong entrainment by drawing the pycnocline closer to the surface. Entrainment is shown to be negligible outside the areas affected by the ice edge upswelling. If cooling at the top is included in the model, the heat and salt exchanges are further enhanced in the upswelling areas. It is noted that new ice formation occurs in the region not affected by ice edge upswelling, and it is suggested that the high-salinity mixed layer regions (with a scale of a few Rossby radii of deformation) will overturn due to cooling, possibly contributing to the formation of deep water.

Hakkinen, Sirpa↗
NASA NTRS2021

VESIcal Part I: An open-source thermodynamic model engine for mixed volatile (H2O-CO2) solubility in silicate melts

Thermodynamics has been fundamental to the interpretation of geologic data and modeling of geologic systems for decades. However, more recent advancements in computational capabilities and a marked increase in researchers’ accessibility to computing tools has outpaced the functionality and extensibility of currently available modeling tools. Here we present VESIcal (Volatile Equilibria and Saturation Identification calculator): the first comprehensive modeling tool for H 2 O, CO 2 , and mixed (H 2 O-CO 2 ) solubility in silicate melts that: a) allows users access to seven commonly used models, plus easy inter-comparison between models; b) provides universal functionality for all models (e.g., functions for calculating saturation pressures, degassing paths, etc.); c) can process large datasets (1,000’s of samples) automatically; d) can output computed data into an excel spreadsheet for simple post-modeling analysis; e) integrates advanced plotting capabilities directly within the tool; and f) provides all of these within the framework of a python library, making the tool extensible by the user and allowing any of the model functions to be incorporated into any other code capable of calling python. The tool is presented within this manuscript, which is a Jupyter notebook containing worked examples accessible to python users with a range of skill levels. The basic functions of VESIcal can also be access via a web app (https://vesical.anvil.app). The VESIcal python library is open-source and available for download at https://github.com/kaylai/VESIcal.

K. Iacovino↗
DOE OSTI2025

Thermodynamic modeling of countercurrent chemical looping reverse water gas shift process for redox material screening

The reverse water gas shift (RWGS) reaction is a key pathway for CO 2 utilization, particularly within Power-to-X process chains aimed at sustainable fuel and chemical production. Countercurrent chemical looping (CL-RWGS) using non-stoichiometric oxides can overcome equilibrium limitations of conventional RWGS reactors, enabling significantly higher CO 2 conversions. However, modeling the limiting performance of such systems is challenging due to their multiphase nature and coupled spatial and temporal variation in chemical composition. In this work, we present a discretized batch equilibrium model that simulates CL-RWGS reactors as a series of localized equilibrium exchanges between gas and solid elements. The model is numerically stable, computationally efficient, and free of kinetic source terms, making it well-suited for parametric studies and system-level integration. It is validated against established convection–diffusion models and shown to predict reasonable upper bounds on experimental results. Application of the model to a range of oxygen carrier materials identifies cerium–zirconium solid solutions, particularly Ce 0.80 Zr 0.20 O 2 , as a promising class offering superior oxygen storage characteristics compared to state-of-the-art La 0.6 Sr 0.4 FeO 3 . This framework provides a robust platform for materials screening, reactor sizing, and performance optimization in chemical looping systems. The model implementation is available as open-source software to support further research and development.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2003

Thermodynamic Modeling of the YO(l.5)-ZrO2 System

The YO1.5-ZrO2 system consists of five solid solutions, one liquid solution, and one intermediate compound. A thermodynamic description of this system is developed, which allows calculation of the phase diagram and thermodynamic properties. Two different solution models are used-a neutral species model with YO1.5 and ZrO2 as the components and a charged species model with Y(+3), Zr(+4), O(-2), and vacancies as components. For each model, regular and sub-regular solution parameters are derived fiom selected equilibrium phase and thermodynamic data.

Jacobson, Nathan S.↗
NASA NTRS2000

A Solution Methodology and Computer Program to Efficiently Model Thermodynamic and Transport Coefficients of Mixtures

A solution methodology has been developed to efficiently model multi-specie, chemically frozen, thermally perfect gas mixtures. The method relies on the ability to generate a single (composite) set of thermodynamic and transport coefficients prior to beginning a CFD solution. While not fundamentally a new concept, many applied CFD users are not aware of this capability nor have a mechanism to easily and confidently generate new coefficients. A database of individual specie property coefficients has been created for 48 species. The seven coefficient form of the thermodynamic functions is currently used rather then the ten coefficient form due to the similarity of the calculated properties, low temperature behavior and reduced CPU requirements. Sutherland laminar viscosity and thermal conductivity coefficients were computed in a consistent manner from available reference curves. A computer program has been written to provide CFD users with a convenient method to generate composite specie coefficients for any mixture. Mach 7 forebody/inlet calculations demonstrated nearly equivalent results and significant CPU time savings compared to a multi-specie solution approach. Results from high-speed combustor analysis also illustrate the ability to model inert test gas contaminants without additional computational expense.

Ferlemann, Paul G.↗
NASA NTRS1999

Thermodynamic Modeling and Analysis of Human Stress Response

A novel approach based on the second law of thermodynamics is developed to investigate the psychophysiology and quantify human stress level. Two types of stresses (thermal and mental) are examined. A Unified Stress Response Theory (USRT) is developed under the new proposed field of study called Engineering Psychophysiology. The USRT is used to investigate both thermal and mental stresses from a holistic (human body as a whole) and thermodynamic viewpoint. The original concepts and definitions are established as postulates which form the basis for thermodynamic approach to quantify human stress level. An Objective Thermal Stress Index (OTSI) is developed by applying the second law of thermodynamics to the human thermal system to quantify thermal stress or dis- comfort in the human body. The human thermal model based on finite element method is implemented. It is utilized as a "Computational Environmental Chamber" to conduct series of simulations to examine the human thermal stress responses under different environmental conditions. An innovative hybrid technique is developed to analyze human thermal behavior based on series of human-environment interaction simulations. Continuous monitoring of thermal stress is demonstrated with the help of OTSI. It is well established that the human thermal system obeys the second law of thermodynamics. Further, the OTSI is validated against the experimental data. Regarding mental stress, an Objective Mental Stress Index (OMSI) is developed by applying the Maxwell relations of thermodynamics to the combined thermal and cardiovascular system in the human body. The OMSI is utilized to demonstrate the technique of monitoring mental stress continuously and is validated with the help of series of experimental studies. Although the OMSI indicates the level of mental stress, it provides a strong thermodynamic and mathematical relationship between activities of thermal and cardiovascular systems of the human body.

Boregowda, S. C.↗
DOE OSTI2025

Thermodynamic Modeling of Complex Solid Solutions in the Lu-H-N System via Graph Neural Network Accelerated Monte Carlo Simulations

Metal hydrides are important across diverse applications, such as hydrogen storage, batteries, gas sensors, nuclear reactions, and high-temperature superconductivity. Previous computational studies of metal hydrides under extreme pressures, e.g., 𝑂⁡(10 2 ) ⁢GPa, usually treat them as stoichiometric compounds without considering interstitial lattice disorder. As pressures become more moderate in the 𝑂⁡(10 0 ) ⁢GPa and below range, hydrogen disorder at interstitial lattice sites becomes prominent, e.g., in the N-doped Lu hydride that was recently claimed superconducting near 1 GPa. Further adding compositional complexity from alloying and/or multielement interstitial occupation makes elucidating pressure- and temperature-dependent observables intractable by first-principles calculations alone. We therefore propose a lattice graph neural-network surrogate modeling approach to predict configuration- and pressure-dependent equation-of-state properties. Their efficiency permits Monte Carlo simulations to calculate Gibbs energies and pressure-dependent phase diagrams, thereby revealing insights into the synthesis conditions required for achieving desired phase equilibria. We demonstrate this concept for the compositionally complex cubic Lu(H,N,Va) 3 system where three constituents (hydrogen, nitrogen and vacancy) have disordered multielement interstitial occupancies and insights into pressure-dependent phase equilibria are critically needed, e.g., N-doping levels can significantly lower dehydrogenation temperatures and provide a new strategy to optimize hydrogen-storage alloys. This work can improve the thermodynamic understanding of the Lu-H-N system and help rational synthesis of N-doped Lu hydrides, but more generally demonstrates an efficient approach to model pressure-dependent thermodynamics of multicomponent solid solutions.

Monte Carlo methods↗
NASA NTRS2004

Thermodynamic Modeling of a Solid Oxide Fuel Cell to Couple with an Existing Gas Turbine Engine Model

The Controls and Dynamics Technology Branch at NASA Glenn Research Center are interested in combining a solid oxide fuel cell (SOFC) to operate in conjunction with a gas turbine engine. A detailed engine model currently exists in the Matlab/Simulink environment. The idea is to incorporate a SOFC model within the turbine engine simulation and observe the hybrid system's performance. The fuel cell will be heated to its appropriate operating condition by the engine s combustor. Once the fuel cell is operating at its steady-state temperature, the gas burner will back down slowly until the engine is fully operating on the hot gases exhausted from the SOFC. The SOFC code is based on a steady-state model developed by The U.S. Department of Energy (DOE). In its current form, the DOE SOFC model exists in Microsoft Excel and uses Visual Basics to create an I-V (current-voltage) profile. For the project's application, the main issue with this model is that the gas path flow and fuel flow temperatures are used as input parameters instead of outputs. The objective is to create a SOFC model based on the DOE model that inputs the fuel cells flow rates and outputs temperature of the flow streams; therefore, creating a temperature profile as a function of fuel flow rate. This will be done by applying the First Law of Thermodynamics for a flow system to the fuel cell. Validation of this model will be done in two procedures. First, for a given flow rate the exit stream temperature will be calculated and compared to DOE SOFC temperature as a point comparison. Next, an I-V curve and temperature curve will be generated where the I-V curve will be compared with the DOE SOFC I-V curve. Matching I-V curves will suggest validation of the temperature curve because voltage is a function of temperature. Once the temperature profile is created and validated, the model will then be placed into the turbine engine simulation for system analysis.

Brinson, Thomas E.↗
DOE OSTI2026

Solvent-mediated contaminant removal from plastic waste using thermodynamic modeling

Plastics recycling is hindered by the compositional complexity of plastic waste, which can include numerous polymer components as well as low concentrations of additives and non-intentionally added substances. These latter small-molecule species, which we collectively refer to as contaminants, can harm human health and will build up in recycled plastic causing environmental and downstream processing challenges if not removed. In this work, we present molecular modeling approaches using the COnductor-like Screening MOdel for Real Solvents (COSMO-RS) to guide the selection of solvents that are capable of removing targeted contaminants from plastic waste. By considering the thermodynamic partitioning of contaminant species between a solvent phase and polymer phase, we identify guidelines for solvent selection to promote either the low-temperature extraction of contaminants from plastic waste or the removal of contaminants as part of a dissolution-based plastics recycling process. We present four case studies to illustrate the application of the computational approach to the removal of brominated flame retardants, phthalates, and selected perfluoroalkyl substances, and compare to both literature and newly collected experimental data to illustrate model prediction accuracy. Furthermore, the case studies highlight the capability of the modeling approach to help design recycling processes that explicitly account for contaminant removal, thereby increasing product purity during dissolution-based recycling or facilitating chemical recycling of contaminant-free plastics.

Zhou, Panzheng [University of, Wisconsin, Madison,↗
NASA NTRS1988

Thermodynamic modeling of the no-vent fill methodology for transferring cryogens in low gravity

The filling of tanks with cryogens in the low-gravity environment of space poses many technical challenges. Chief among these is the inability to vent only vapor from the tank as the filling proceeds. As a potential solution to this problem, the NASA Lewis Research Center is researching a technique known as No-Vent Fill. This technology potentially has broad application. The focus is the fueling of space based Orbital Transfer Vehicles. The fundamental thermodynamics of the No-Vent Fill is described. The model is then used to conduct a parametric investigation of the key parameters: initial tank wall temperature, liquid-vapor interface heat transfer rate, liquid inflow rate, and inflowing liquid temperatures. Liquid inflowing temperature and the liquid-vapor interface heat transfer rate seem to be the most significant since they influence the entire fill process. The initial tank wall temperature must be sufficiently low to prevent a rapid pressure rise during the initial liquid flashing state, but then becomes less significant.

Chato, David J.↗
NASA NTRS1988

Thermodynamic modeling of the no-vent fill methodology for transferring cryogens in low gravity

The filling of tanks with cryogens in the low-gravity environment of space poses many technical challenges. Chief among these is the inability to vent only vapor from the tank as the filling proceeds. As a potential solution to this problem, the NASA Lewis Research Center is researching a technique known as No-Vent Fill. This technology potentially has broad application. The focus is the fueling of space based Orbital Transfer Vehicles. The fundamental thermodynamics of the No-Vent Fill process to develop an analytical model of No-Vent Fill is described. The model is then used to conduct a parametric investigation of the key parameters: initial tank wall temperature, liquid-vapor interface heat transfer rate, liquid inflow rate, and inflowing liquid temperatures. Liquid inflowing temperature and the liquid-vapor interface heat transfer rate seem to be the most significant since they influence the entire fill process. The initial tank wall temperature must be sufficiently low to prevent a rapid pressure rise during the initial liquid flashing stage, but then becomes less significant.

Chato, David J.↗
NASA NTRS2003

Latent Heating Retrieval from TRMM Observations Using a Simplified Thermodynamic Model

A procedure for the retrieval of hydrometeor latent heating from TRMM active and passive observations is presented. The procedure is based on current methods for estimating multiple-species hydrometeor profiles from TRMM observations. The species include: cloud water, cloud ice, rain, and graupel (or snow). A three-dimensional wind field is prescribed based on the retrieved hydrometeor profiles, and, assuming a steady-state, the sources and sinks in the hydrometeor conservation equations are determined. Then, the momentum and thermodynamic equations, in which the heating and cooling are derived from the hydrometeor sources and sinks, are integrated one step forward in time. The hydrometeor sources and sinks are reevaluated based on the new wind field, and the momentum and thermodynamic equations are integrated one more step. The reevalution-integration process is repeated until a steady state is reached. The procedure is tested using cloud model simulations. Cloud-model derived fields are used to synthesize TRMM observations, from which hydrometeor profiles are derived. The procedure is applied to the retrieved hydrometeor profiles, and the latent heating estimates are compared to the actual latent heating produced by the cloud model. Examples of procedure's applications to real TRMM data are also provided.

Grecu, Mircea↗
NASA NTRS1984

A thermodynamic model for subsolidus equilibria in the system CaO-MgO-Al2O3-SiO2

It is shown that the high temperature enthalpy of solution data for pure phases and solid solutions in the CMAS system are, for 'gabbroic' and 'peridotitic' compositions, consistent with available phase equilibrium data for the MAS, CAS and CMAS systems. A refined set of values of thermodynamic properties for these phases and solid solutions is tabulated. The small differences between the new data set and that of Helgeson et al. (1978) arise from new data on heat capacity and enthalpy being incorporated. The important constraints applied and the major difficulties in fitting which arose are summarized.

Wood, B. J.↗