Search NASA⌕ Search

SEARCH · Search NASA

Results for “Thermodynamic calculations”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS2005

Phlogopite Decomposition, Water, and Venus

Venus is a hot and dry planet with a surface temperature of 660 to 740 K and 30 parts per million by volume (ppmv) water vapor in its lower atmosphere. In contrast Earth has an average surface temperature of 288 K and 1-4% water vapor in its troposphere. The hot and dry conditions on Venus led many to speculate that hydrous minerals on the surface of Venus would not be there today even though they might have formed in a potentially wetter past. Thermodynamic calculations predict that many hydrous minerals are unstable under current Venusian conditions. Thermodynamics predicts whether a particular mineral is stable or not, but we need experimental data on the decomposition rate of hydrous minerals to determine if they survive on Venus today. Previously, we determined the decomposition rate of the amphibole tremolite, and found that it could exist for billions of years at current surface conditions. Here, we present our initial results on the decomposition of phlogopite mica, another common hydrous mineral on Earth.

Johnson, N. M.↗
DOE OSTI2022

Machine learning molecular dynamics simulations toward exploration of high-temperature properties of nuclear fuel materials: case study of thorium dioxide

Predicting materials properties of nuclear fuel compounds is a challenging task in materials science. Their thermodynamical behaviors around and above the operational temperature are essential for the design of nuclear reactors. However, they are not easy to measure, because the target temperature range is too high to perform various standard experiments safely and accurately. Moreover, theoretical methods such as first-principles calculations also suffer from the computational limitations in calculating thermodynamical properties due to their high calculation-costs and complicated electronic structures stemming from f-orbital occupations of valence electrons in actinide elements. Here, we demonstrate, for the first time, machine-learning molecular-dynamics to theoretically explore high-temperature thermodynamical properties of a nuclear fuel material, thorium dioxide. The target compound satisfies first-principles calculation accuracy because f-electron occupation coincidentally diminishes and the scheme meets sampling sufficiency because it works at the computational cost of classical molecular-dynamics levels. We prepare a set of training data using first-principles molecular dynamics with small number of atoms, which cannot directly evaluate thermodynamical properties but captures essential atomistic dynamics at the high temperature range. Then, we construct a machine-learning molecular-dynamics potential and carry out large-scale molecular-dynamics calculations. Consequently, we successfully access two kinds of thermodynamic phase transitions, namely the melting and the anomalous λ transition induced by large diffusions of oxygen atoms. Furthermore, we quantitatively reproduce various experimental data in the best agreement manner by selecting a density functional scheme known as SCAN. Our results suggest that the present scale-up simulation-scheme using machine-learning techniques opens up a new pathway on theoretical studies of not only nuclear fuel compounds, but also a variety of similar materials that contain both heavy and light elements, like thorium dioxide.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS2024

Theoretical Prediction of Thermal Expansion Anisotropy for Y 2 Si 2 O 7 Environmental Barrier Coatings Using a Deep Neural Network Potential and Comparison to Experiment

Environmental barrier coatings (EBCs) are an enabling technology for silicon carbide (SiC)-based ceramic matrix composites (CMCs) in extreme environments such as gas turbine engines. However, development of new coating systems is hindered by the large design space and difficulty in predicting properties for these materials. Density Functional Theory (DFT) has successfully been used to model and predict some thermodynamic and thermo-mechanical properties of high-temperature ceramics for EBCs, although these calculations are challenging due to their high computational costs. In this work, we use machine learning to train a deep neural network potential (DNP) for Y 2 Si 2 O 7 , which is then applied to calculate thermodynamic and thermo-mechanical properties at near-DFT accuracy much faster and using less computational resources than DFT. We use this DNP to predict phonon-based thermodynamic properties of Y 2 Si 2 O 7 with good agreement to DFT and experiments. We also utilize the DNP to calculate the anisotropic, lattice direction-dependent coefficients of thermal expansion (CTEs) for Y 2 Si 2 O 7 . Molecular dynamics trajectories using the DNP correctly demonstrate accurate prediction of the anisotropy of the CTE in good agreement with diffraction experiments. In the future, this DNP could be applied to accelerate additional property calculations for Y 2 Si 2 O 7 compared to DFT or experiments.

rare earth silicates↗
NASA NTRS2010

Methane on Mars: Thermodynamic Equilibrium and Photochemical Calculations

The detection of methane (CH4) in the atmosphere of Mars by Mars Express and Earth-based spectroscopy is very surprising, very puzzling, and very intriguing. On Earth, about 90% of atmospheric ozone is produced by living systems. A major question concerning methane on Mars is its origin - biological or geological. Thermodynamic equilibrium calculations indicated that methane cannot be produced by atmospheric chemical/photochemical reactions. Thermodynamic equilibrium calculations for three gases, methane, ammonia (NH3) and nitrous oxide (N2O) in the Earth s atmosphere are summarized in Table 1. The calculations indicate that these three gases should not exist in the Earth s atmosphere. Yet they do, with methane, ammonia and nitrous oxide enhanced 139, 50 and 12 orders of magnitude above their calculated thermodynamic equilibrium concentration due to the impact of life! Thermodynamic equilibrium calculations have been performed for the same three gases in the atmosphere of Mars based on the assumed composition of the Mars atmosphere shown in Table 2. The calculated thermodynamic equilibrium concentrations of the same three gases in the atmosphere of Mars is shown in Table 3. Clearly, based on thermodynamic equilibrium calculations, methane should not be present in the atmosphere of Mars, but it is in concentrations approaching 30 ppbv from three distinct regions on Mars.

Levine, J. S.↗
DOE OSTI2025

First-principles thermodynamics of Al 10 ⁢V: An analytical treatment of localized anharmonic modes

Many complex intermetallic structures possess cagelike environments that can host additional guest atoms. In Al 10 ⁢V, these atoms give rise to low-frequency, localized vibrations (Einstein modes) that dominate the thermodynamic response at low temperature. They become imaginary under volume expansion as temperature rises, invalidating the harmonic approximation. Here, we develop a framework to incorporate these strongly anharmonic vibrational modes into first-principles thermodynamic calculations. By explicitly modeling the cage potential and solving the associated Schrödinger equation numerically, we compute the full anharmonic free energy contribution and demonstrate its impact on the thermodynamic behavior of Al 10⁢ V. This allows us to examine structures with different cage fillings and construct the Al-V phase diagram in the relevant composition range. Our results reproduce key experimental signatures, including the anomalous rise in the thermal expansion coefficient and heat capacity at low temperatures, and reveal that the presence and the extent of cage filling by guest atoms is essential to stabilizing the Al 10 ⁢V phase at elevated temperatures.

anharmonic lattice dynamics↗
DOE OSTI2023

An asymptotic approach for the statistical thermodynamics of certain model systems

In classical statistical thermodynamics, calculating the configuration integral is both vital and elusive. Analytic relations for configuration integrals are desirable for modeling purposes, but it is typically impossible to obtain them. Certain systems become analytically tractable after replacing steep potential energies with harmonic potentials or athermal rigid constraints, but these approximations are often inadequate, especially when modeling the stretching of molecules. It is therefore necessary to develop a systematic approach to improve upon the approximations provided by these reference systems. Here, a general asymptotic approach is introduced, where the configuration integral for the full system is obtained in terms of that of the reference system and several corrections. This asymptotic approach is first demonstrated using the simple example of a classical three-dimensional oscillator. Next, the approach is applied to modeling the stretching of single polymer chains and to modeling thermally assisted crack growth, where results are verified with respect to numerical calculations. Overall, this asymptotic approach is a valid and effective tool for statistical thermodynamics in general.

Buche, Michael Robert↗
DOE OSTI2020

Thermodynamic Modeling of Aqueous LiCl, LiBr, LiI, and LiNO 3 Solutions

Thermodynamic models are essential to facilitate the advancement of process design, optimization, and operation of electrolyte systems. In this work, a comprehensive thermodynamic framework based on the Electrolyte Nonrandom Two-Liquid model is developed to calculate phase equilibria behavior and salt solubility of aqueous LiCl, LiBr, LiI, and LiNO 3 solutions. The model describes the non-ideality of the electrolyte solutions by using two binary interaction parameters for each electrolyte-molecule pair in the system. Each binary interaction parameter is further expressed with up to three temperature coefficients which are regressed from experimental data. To take into account the hydration of lithium ion, two separate chemistries for the dissociation of lithium salts are investigated. In the first case, the lithium ion is considered as a bare ion, Li + , while in the second case hydration of the lithium ion from Li + to Li(H 2 O) + is considered. Here, the calculated thermodynamic properties compare adequately with the experimental data for both sets of chemistries for concentrations up to saturation and temperatures from 273.15 K up to 623.15 K. Moderate to significant improvements are observed with the incorporation of the hydration chemistry for aqueous LiCl, LiBr, and LiI solutions when compared to the non-hydrated lithium ion model results.

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

Phase stability in cast and additively manufactured Al-rich Al-Cu-Ce alloys

Additively manufactured (AM) eutectic Al alloy systems have been studied extensively for advantageous thermal stability and mechanical properties due to their refined microstructures. Al-Cu-Ce alloys are one subset of these AM eutectic alloys. Here we studied phase stability in AM Al-Cu-Ce alloys and compared it to that of conventionally cast ones. A new phase, Al 8 Cu 3 Ce, was identified in the microstructures of both AM and cast Al-Cu-Ce alloys. This Al 8 Cu 3 Ce phase was not previously included on experimental or thermodynamically calculated phase diagrams of the Al-Cu-Ce system. Therefore, we performed additional thermodynamic modeling of the system. These models were experimentally validated with cast and subsequently heat-treated Al-Cu-Ce alloys. We found that despite the refined microstructure of the AM alloys, the phases formed were consistent with the cast alloys, suggesting that AM processing did not significantly alter the formation and stability of phases from that in the conventional alloys. In conclusion, this work has resolved inconsistent previous descriptions of the Al-Cu-Ce ternary phase diagram in the Al-rich region and resulted in the addition of the Al 8 Cu 3 Ce as an equilibrium phase above 500 °C.

36 MATERIALS SCIENCE↗
NASA NTRS1987

Calculating Theromodynamic And Transport Properties Of Fluids

Computer program incorporates van der Waals equation and correction tables. FLUID program developed to calculate thermodynamic and transport properties of pure fluids in both liquid and gas phases. Properties calculated by use of simple gas model, empirical corrections, and efficient numerical interpolation scheme. Produces results that agree very well with measured values. Much faster than older, more complex programs developed for same purpose.

Proctor, Margaret P.↗
NASA NTRS1968

Computer program for calculation of ideal gas thermodynamic data

Computer program calculates ideal gas thermodynamic properties for any species for which molecular constant data is available. Partial functions and derivatives from formulas based on statistical mechanics are provided by the program which is written in FORTRAN 4 and MAP.

Gordon, S.↗
DOE OSTI2024

Linear Free Energy Relationships Associated with Hydride Transfer From [(6,6'-R 2 -bpy)Re(CO) 3 H]: A Cautionary Tale in Identifying Hydrogen Bonding Effects in the Secondary Coordination Sphere

Here, six rhenium hydride complexes, [(6,6'-R 2 -bpy)Re(CO) 3 H] (bpy = 2,2'-bipyridine, R = OEt, OMe, NHMe, Me, F, Br), were synthesized. These complexes insert CO 2 to form rhenium formate complexes of the type [(6,6'-R 2 -bpy)Re(CO) 3 {OC(O)H}]. All the rhenium formate species were characterized using X-ray crystallography, which revealed that the bpy ligand is not coplanar with the metal coordination plane containing the two nitrogen donors of the bpy ligand but tilted. A solid-state structure of [(6,6'-Me 2 -bpy)Re(CO) 3 H] determined using MicroED also featured a tilted bpy ligand. The kinetics of CO 2 insertion into complexes of the type [(6,6'-R 2 -bpy)Re(CO) 3 H] were measured experimentally and the thermodynamic hydricities of [(6,6'-R 2 -bpy)Re(CO) 3 H] species were determined using theoretical calculations. A Brønsted plot constructed using the experimentally determined rate constants for CO 2 insertion and the calculated thermodynamic hydricities for [(6,6'-R 2 -bpy)Re(CO) 3 H] revealed a linear free energy relationship (LFER) between thermodynamic and kinetic hydricity. This LFER is different to the previously determined relationship for CO 2 insertion into complexes of the type [(4,4'-R 2 -bpy)Re(CO) 3 H]. At a given thermodynamic hydricity, CO 2 insertion is faster for complexes containing a 6,6'-substituted bpy ligand. This is likely in part due to the tilting observed for systems with 6,6'-substituted bpy ligands. Notably, the 6,6'-(NHMe) 2 -bpy ligand could in principle stabilize the transition state for CO 2 insertion via hydrogen bonding. This work shows that if only the rate of CO 2 insertion into [(6,6'-(NHMe) 2 -bpy)Re(CO) 3 H] is compared to [(4,4'-R 2 -bpy)Re(CO) 3 H] systems, the increase in rate could be easily attributed to hydrogen bonding, but in fact all 6,6'-substituted systems lead to faster than expected rates.

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

Identifying crystal structures and chemical reactions at the interface of stanene on Bi 2 Te 3

Synthesizing monolayers and heterostructures is an enabling approach to extract new physical phenomena from bulk materials. Among the structures amenable to this approach is stanene, which is a monolayer of tin, similar to graphene, and has been predicted to host one-dimensional topological states at its edges. Stanene can be tuned by decorating with different adatoms, which makes it a promising platform on which to engineer topological devices. Here, we deposit Sn on Bi 2 Te 3 and characterize the growth using anomalous synchrotron x-ray scattering and x-ray photoelectron spectroscopy (XPS). X-ray diffraction data reveal the formation of epitaxial Sn-based structures, along with penetration of Sn into the Bi 2 Te 3 , with Sn intercalating between the upper 10 Bi 2 Te 3 quintuple layers. Additionally, XPS data show deposited Sn reacting to form SnTe and Bi at the Bi 2 Te 3 surface. The calculated heat of reaction for Sn and Bi2Te3 is consistent with an exothermic reaction to SnTe and Bi. Using thermodynamic calculations as a guide, we identify several candidate substrates that can stabilize the stanene phase.

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

Deposit geometry and oxygen concentration spatial variations due to composition change in printed functionally graded components

Functionally graded materials (FGMs) with site specific chemical composition are commonly manufactured by directed energy deposition (DED). Although previous work fabricated an FGM with a compositional variation between a ferritic and austenitic alloy, difficulties arose due to variations in deposit shape with composition change. This problem also occurs for FGMs in literature; however, unlike other cases, the thermophysical properties of these two alloys were similar throughout the build. Here, we investigate the role of chemical composition and surface-active elements on deposit geometry during the manufacture of FGMs by laser DED. Single-track experiments for the relevant FGM compositions are analyzed with results from a well-tested, three-dimensional, transient numerical heat transfer and fluid flow model and thermodynamic calculations. Experiments showed deposit shape varied as a function of composition for constant laser power and scan speed. Thermodynamic analysis indicated that the oxygen solubility in the fusion zone varied significantly for each composition used for the FGM. Numerical modeling revealed that the change in fluid flow caused by Marangoni convection due to dissolved oxygen in the fusion zone were mainly responsible for the changes in deposit shape observed in experiments. Finally, because oxygen can be introduced into the fusion zone through the feedstock as well as the surrounding atmosphere, these findings elucidate a previously unconsidered aspect of process control during DED fabrication of FGMs.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Computational thermodynamic study of SiC chemical vapor deposition from MTS-H 2

This study focuses on the computational thermodynamic analysis of the chemical vapor deposition (CVD) of SiC from the methyltrichlorosilane-hydrogen (MTS-H 2 ) using up-to-date thermodynamic databases. High-resolution computation has been performed with the fine intervals of temperature and pressure at the various H 2 /MTS ratios of interest to systematically investigate the deposition condition range (800 to 1600°C, 0 to 26 664 Pa, and H 2 /MTS ratios of 0.1 to 100) to guide experimental exploration. The influence of deposition parameters on the compositions and phase stabilities of the deposit and gas phase pertinent to vapor processing is elucidated. Low pressure and medium temperatures (1000 to 1400°C) are beneficial to reaching a higher SiC deposition efficiency and provide an optimal window for preparing a high-purity (>99 wt.% SiC) deposit. This optimal processing window expands significantly with an increasing H 2 /MTS ratio (<20). These results are supported by a number of previous theoretical and experimental observations. The mass fraction of SiC in deposit is proposed as an additional perspective to understand the discrepancy between thermodynamic calculation and experimental observation of pure CVD SiC at low H 2 /MTS ratios.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Theoretical assessments of Pd–PdO phase transformation and its impacts on H 2 O 2 synthesis and decomposition pathways

The direct synthesis of H 2 O 2 from O 2 and H 2 provides a green pathway to produce H 2 O 2 , a popular industrial oxidant. Here, in this study, we theoretically investigate the effects of Pd oxidation states, coordination environments, and particle sizes on primary H 2 O 2 selectivities, assessed by calculating the ratio of rate constants for the formation of H 2 O 2 (via OOH* reduction; k O–H ) and the decomposition of OOH* (via O–O cleavage; k O–O ). For Pd metals, the k O–H /k O–O ratio decreased from 10 -4 for Pd(111) to 10 -10 for the Pd 13 cluster at 300 K, indicating poorer H 2 O 2 selectivity as Pd particle size decreases and low primary selectivities for H 2 O 2 overall. As the oxygen chemical potential increases and metals form surface and bulk oxides, the perturbation of Pd–Pd ensemble sites by lattice O atoms results in selectivities that become dramatically higher than unity. For instance, at 300 K, the k O–H /k O–O ratio increases significantly from 10 -4 to 10 9 to 10 16 as Pd(111) oxidizes to Pd 5 O 4 /Pd(111) and to PdO(100), respectively. In contrast, such selectivity enhancements are not observed for surface and bulk oxides that persistently contain rows of more metallic, undercoordinated Pd–Pd ensemble sites, such as PdO(101)/Pd(100) and PdO(101). These Pd–Pd ensembles are also absent when smaller Pd nanoparticles fully oxidize, indicating that smaller PdO clusters can be more selective for H 2 O 2 synthesis. These trends for primary H 2 O 2 selectivities were found to inversely correlate with trends for H 2 O 2 decomposition rates via O–O bond cleavage, demonstrating that catalysts with high primary H 2 O 2 selectivity can also hinder H 2 O 2 decomposition. Ab initio thermodynamic calculations are used to estimate the thermodynamically favored phase among Pd, PdO/Pd and PdO in O 2 , H 2 O 2 /H 2 O, and O 2 /H 2 environments. These results are combined to show that smaller Pd nanoparticles are more prone to be oxidized at lower oxygen chemical potentials, upon which they become more selective than larger Pd particles for H 2 O 2 synthesis.

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

Verification Testing of OLI Systems Mixed Solvent Electrolyte Model for the Na-K-Mg-Ca-H-Cl-SO 4 -OH-HCO 3 -CO 3 -CO 2 -H 2 ) System to High Ionic Strength at 25°C

This technical report summarizes model verification results and summary statistics for 41 evaporite mineral solubility cases evaluated by Savannah River National Laboratory using OLI Systems’ aqueous electrolyte thermodynamic modeling software. The 41 verification cases containing a total of 60 solubility curves comprise mineral solubility data from low to high ionic strength at 25°C for the eight-component system Na-K-Mg-Ca-H-Cl-SO 4 -OH-HCO 3 -CO 3 -CO 2 -H 2 O as reported by Harvie et al. (1984). Thermodynamic calculations were executed using OLI Systems’ Stream Analyzer computation module within the OLI Studio software platform (Ver. 11.0, Rev. 11.0.1.9). The Mixed Solvent Electrolyte (MSE) thermodynamic framework was chosen for this investigation because of its superiority in modeling high ionic-strength inorganic salt solutions and actinide redox chemistry and solubility, both of which are relevant to the geological repository conditions at the Waste Isolation Pilot Plant in Carlsbad, New Mexico. Mineral solubility data in various inorganic salt solutions were digitized and extracted from figures generated by Harvie et al. (1984). For each of the 60 solubility curves, a case-specific chemistry model and input file were generated in OLI Studio using OLI Stream Analyzer and the MSE (H 3 O + ion) public databank provided by OLI Systems. Model simulation results were exported to Microsoft Excel to calculate summary statistics and to generate graphs comparing the OLI model predictions to the solubility data. Summary statistics include residuals (model – data) and concordance (accuracy × precision, where precision is indicated by the Pearson correlation coefficient and accuracy accounts for bias and scale differential). Private databanks were not developed, and activity coefficient model regressions were not performed to improve OLI model fits to the data. Of the 41 model verification plots, 83% have a mean of the percent residuals less than or equal to 25%. Similarly, 75% display a concordance greater than or equal to 0.75. Only seven of the 41 verification plots fail to show good agreement between the model and data. Of these seven, three are relevant to the WIPP repository because they involve the Mg-OH-Cl-SO 4 -CO 3 aqueous system. The remaining four address salt solubilities at the pH extremes (strong acid and strong base). It should be noted that in two of the three Mg-OH-Cl-SO 4 -CO 3 system cases, the regressed Harvie et al. (1984) solubility curve also deviated from the data. Lack of agreement between the OLI model-predicted solubility curves and the data is attributable to one or more of the following: specific solid species are not included in the OLI MSE databank; there is significant variation among the different solubility datasets chosen by Harvie et al. (1984); the OLI MSE model’s thermodynamic parameters were determined using different solubility datasets; and the activity coefficient parameters for certain relevant ion-ion and ion-molecule pairs have not been optimized via data regression. Two recommendations for future work are to (1) evaluate solubility data for the Mg-OH-Cl-SO 4 -CO 3 system at high ionic strength and, if necessary, develop a private OLI MSE database that includes missing species and, where necessary, regressed standard state properties and interaction parameters; (2) perform similar verification testing of the OLI model for actinide solubility data.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1992

Steady and transient performance calculation method for prediction, analysis, and identification

The detailed design and development of turbofans Involves the prediction and identification, by means of test analysis, of the performance of the engine and its components. The thermodynamic simulation and analysis codes integrate existing knowledge and Interpretations of the detailed operating procedure of the components of the engine being developed. The relevance of the predicted performance depends on the quality of the representation of the various physical phenomena affecting the characteristics of the components and, consequently, on the incorporation of experimental correlation in the modelling. In this context, the representation of compressor and turbine characteristics is particularly important. Firstly, we will analyze the ability of corrected parameters to represent MACH similitude at the component inlet under various conditions. The various measurements achieved on the engine during development are used to enhance the modelling. The methods of identifying the thermodynamic calculation code with the various measurements, considered here with their uncertainties, are then presented and described. The analysis of the tests performed on the powerplant, designed to Identify the real characteristics of the engine components, can be undertaken, considering one or more engine test points and incorporating knowledge acquired through experimentation or on the component test bench. There are many possible fields of use for these identification methods ranging from the rematching of components to the optimization of the control system.

Jean Pierre Duponchel↗