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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Wet to Dry Controls Lanthanide Scandate Synthesis

The choice of temperature and gas conditions used in a water pressure-controlled reactor is guided by density functional theory (DFT) to synthesize nearly phase pure lanthanide scandate nanoparticles (LnScO 3 , Ln = La, Nd, Sm, Gd). In this synthetic method, low water vapor partial pressures, well below water’s gas liquidus, inhibit particle growth, while an excess of water vapor results in undesired rare earth hydroxide and oxyhydroxide secondary phases. The optimal humidity for high-purity LnScO 3 particle synthesis is shown to vary with the lanthanide; DFT is used to calculate the thermodynamics of secondary phase formation for each lanthanide tested, such that the role of water vapor may be quantified and used to maintain phase purity (greater than 96 mol%) across the series. Furthermore, the combination of thermodynamic calculation and experimental confirmation with this pressure-controlled reactor provides an opportunity to explore analogous syntheses of other inorganic perovskite nanoparticles.

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

Understanding Oxide–Metal Interactions During Hot Isostatic Pressing to Diffusion Bond Aluminum Alloy 6061 Plates

The interaction between Mg, Si, and Al 2 O 3 during hot isostatic pressing diffusion bonding of aluminum alloy 6061 (AA6061) plates was investigated through thermodynamic calculations and experimental microstructural characterization. Thermodynamic calculations as functions of temperature, pressure, and composition revealed that the interaction among Mg, Si, and Al 2 O 3 yields Mg 2 Si and either MgO + Al or MgAl 2 O 4 + Al, facilitating the reduction of Al 2 O 3 and allowing Al/Al metallic bonds to form. Total pressure variation had a negligible influence on the oxygen partial pressure, and consequently, the reaction product formation. Oxygen partial pressure variation as a function of temperature and initial amount of Al 2 O 3 determined the formation of either MgO or MgAl 2 O 4 . Experimental Hot Isostatic Pressure (HIP) bonding at 723 K and 833 K under a constant pressure of 1017 atm documented the cooling rate-dependent formation of β–Mg 2 Si precipitates. High-resolution transmission electron microscopy imaging and selected area electron diffraction patterns verified the formation of β–Mg 2 Si and MgO at the interface but did not detect MgAl 2 O 4 . In conclusion, findings from this study clarify the role of thermochemical interactions in oxide disruption and bonding mechanisms during HIP diffusion bonding of AA6061 and provide guidance for optimizing joining processes for monolithic nuclear fuel assemblies.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1973

Computer program for calculation of thermodynamic and transport properties of complex chemical systems

Program performs calculations such as chemical equilibrium for assigned thermodynamic states, theoretical rocket performance for both equilibrium and frozen compositions during expansion, incident and reflected shock properties, and Chapman-Jouget detonation properties. Features include simplicity of input and storage of all thermodynamic and transport property data on master tape.

Svehla, R. A.↗
DOE OSTI2021

Calculation of the Thermodynamic Voltage Limit of CdSeTe Solar Cells

The first step to understand the origin of losses in any photovoltaic solar cell is to determine the fundamental thermodynamic efficiency and voltage limits of such a device. In this contribution, we detail techniques to calculate the voltage limit in the case of cadmium selenium telluride (CdSeTe) solar cells, and how approaches based on bandgap alone—i.e., the Shockley-Queisser approach with step-function absorptance—can overestimate the thermodynamic open-circuit voltage limit Voc,ideal . This is particularly true for arsenic-doped samples, which tend to exhibit below-bandgap absorptance.

absorptance↗
NASA NTRS1975

GASP: A computer code for calculating the thermodynamic and transport properties for ten fluids: Parahydrogen, helium, neon, methane, nitrogen, carbon monoxide, oxygen, fluorine, argon, and carbon dioxide

A FORTRAN IV subprogram called GASP is discussed which calculates the thermodynamic and transport properties for 10 pure fluids: parahydrogen, helium, neon, methane, nitrogen, carbon monoxide, oxygen, fluorine, argon, and carbon dioxide. The pressure range is generally from 0.1 to 400 atmospheres (to 100 atm for helium and to 1000 atm for hydrogen). The temperature ranges are from the triple point to 300 K for neon; to 500 K for carbon monoxide, oxygen, and fluorine; to 600 K for methane and nitrogen; to 1000 K for argon and carbon dioxide; to 2000 K for hydrogen; and from 6 to 500 K for helium. GASP accepts any two of pressure, temperature and density as input conditions along with pressure, and either entropy or enthalpy. The properties available in any combination as output include temperature, density, pressure, entropy, enthalpy, specific heats, sonic velocity, viscosity, thermal conductivity, and surface tension. The subprogram design is modular so that the user can choose only those subroutines necessary to the calculations.

Hendricks, R. C.↗
DOE OSTI2020

Observed volatilization behavior of silicon carbide in flowing hydrogen above 2000 K

In this work, the intrinsic compatibility of silicon carbide (SiC) and hydrogen (H 2 ) at high temperatures (2000-2473 K) and pressure near one atmosphere was evaluated through a combination of thermodynamic calculations and hot hydrogen exposure testing. Thermodynamic calculations predict the decomposition of SiC in a hydrogen environment to form free silicon (Si) and free carbon (C). Free Si is predicted to vaporize from the surface as a volatile species, while free C may interact with H 2 to form the hydrocarbons CH 4 (T < 2100 K) or C 2 H 2 (T > 2100 K). Coupons of high purity chemical vapor deposition (CVD) β-SiC were exposed to slowly flowing hydrogen at temperatures ranging between 2000 and 2473 K. SiC experienced active attack as the result of H 2 exposure, exhibiting linear weight loss kinetics and an Arrhenius dependence of weight loss on exposure temperature. The linear volatilization constant was experimentally evaluated to correspond with an activation energy of 370 ± 18 kJ/mol. Due to the dependence of observed corrosion rates on gas velocity, corrosion of SiC in flowing H 2 was determined to be governed by external mass transfer of volatile Si species through the boundary layer. Experimentally derived mass losses were in good agreement with mass losses predicted by a boundary layer limited gas diffusion model.

36 MATERIALS SCIENCE↗
NASA NTRS1992

Thermodynamics, Solubility, and Diffusivity of Oxygen in Titanium and Ti-Al Alloys

Titanium aluminides and titanium aluminide-based composites are attractive candidate materials for high-temperature structural applications. As these materials may be exposed to oxidizing environments durine their use at elevated temperatures, it is essential that they possess a good oxidation resistance. Previous studies have shown that the oxidation resistance of Al-rich alloys in the Ti-Al system is superior to that of the Ti-rich alloys. The scales formed on the surface of the Al-rich and Ti-rich alloys have been reported to be predominantly Al2O3 and TiO2, respectively. Since the relative stabilities of the oxides of Al and Ti at various temperatures and oxygen pressures can be assessed from their thermodynamic data, it is possible, With the help of thermodynamic calculations, to determine the compositions of the alloys which would form scales of Al2O3, TiO(x) or a ternary oxide such as TiAl2O5 during oxidation at a given temperature. The thermodynamic calculations require reliable activity data for the Ti-Al system. These data have not been determined for the entire composition and temperature range of interest. Using the data available in the literature, recently performed thermodynamic calculations and concluded that the stable oxide changed from TiO to Al2O3 in the existence region of the tial phase. In the case of titanium aluminide-based composites, another major concern is the mutual chemical compatibility of the matrix material with the reinforcement phase. Fibers of SiC, TiB2 and Al2O3 are currently being investigated for reinforcement of titanium aluminide matrices.

Mehrotra, Gopal M.↗
DOE OSTI2023

Thermodynamic modeling with uncertainty quantification using the modified quasichemical model in quadruplet approximation: Implementation into PyCalphad and ESPEI

The modified quasichemical model in the quadruplet approximation (MQMQA) considers the first- and the second-nearest-neighbor coordination and interactions, particularly useful in describing short-range ordering (SRO) in complex liquids such as molten salts, slag in metal processing, and electrolytic solutions. Here, the present work implements the MQMQA into the Python based open-source software PyCalphad for thermodynamic calculations. This endeavor facilitates the development of MQMQA-based thermodynamic database with uncertainty quantification (UQ) and propagation (UP) using the open-source software ESPEI. A new database structure based on Extensible Markup Language (XML) is proposed for ESPEI evaluation of MQMQA model parameters. Using the KF-NiF 2 , KCl-NaCl-MgCl 2 , and CaCl 2 -CaF 2 -LiCl-LiF salt systems as examples, we demonstrate the successful implementation of MQMQA in PyCalphad through thermodynamic calculations of Gibbs energy, equilibrium quadruplet fractions, and phase diagram, as well as database development with UQ and UP using ESPEI. Furthermore, as an application of the present implementation, both the LiF–TbF 3 and LiF-HoF 3 systems have been modeled by MQMQA for the first time, which are in good agreement with experiments. The present implementation hence offers an open-source capability for performing CALPHAD modeling for complex liquids with SRO using MQMQA plus a new XML database structure.

36 MATERIALS SCIENCE↗
NASA NTRS1993

Shock-induced devolatization of calcium sulfate and implications for K-T extinctions

Calcium sulfate devolatization during the impact at Chicxulub, Mexico and dispersal in the stratosphere of the resultant sulfuric acid aerosol have been suggested as a possible mechanism for the Cretaceous-Tertiary extinctions. In this paper, we investigated two shock-induced devolatization reactions of calcium sulfate up to 42 GPa in the laboratory: CaSO4 + SiO2 yields CaSiO3 + SO3(degassed) and CaSO4 yields CaO + SO2(degassed) + 1/2 O2(degassed). We found both to proceed to a much less extent than calculated by equilibrium thermodynamic calculations. Reaction products are found to be 10(exp -2) times those calculated for equilibrium. Consequently our estimate of the amount of sulfur oxides degassed into the atmosphere from shock devolatization of CaS04 in the Chicxulub lithographic section (6x10(exp 15)-2x10(exp 16)g in sulfur mass) is lower by a factor of 70 to 400 than previous estimates; the related environmental stress arising from the resultant global cooling of approximately 4 K and fallout of acid rain does not appear to suffice to explain the widespread K-T extinctions.

Chen, Guangqing↗
NASA NTRS1994

Shock-induced devolatilization of calcium sulfate and implications for K-T extinctions

The devolatilization of calcium sulfate, which is present in the target rock of the Chicxulub, Mexico impact structure, and dispersal in the stratosphere of the resultant sulfuric acid aerosol have been suggested as a possible mechanism for the Cretaceous-Tertiary extinctions. We measured the amount of SO2 produced from two shock-induced devolatilization reactions of calcium sulfate up to 42 GPa in the laboratory. We found both to proceed to a much lower extent than calculated by equilibrium thermodynamic calculations. Reaction products are found to be approx. 10(exp -2) times those calculated for equilibrium. Upon modeling the quantity of sulfur oxides degassed into the atmosphere from shock devolatilization of CaSO4 in the Chicxulub lithographic section, the resulting 9 x 10(exp 16) to 6 x 10(exp 17) g (in sulfur mass) is lower by a factor of 10-100 than previous upper limit estimates, the related environmental stress arising from the resultant global cooling and fallout of acid rain is insufficient to explain the widespread K-T extinctions.

Chen, Guangqing↗
DOE OSTI2024

Acceleration of Solvation Free Energy Calculation via Thermodynamic Integration Coupled with Gaussian Process Regression and Improved Gelman–Rubin Convergence Diagnostics

The determination of the solvation free energy of ions and molecules holds profound importance across a spectrum of applications spanning chemistry, biology, energy storage, and the environment. Molecular dynamics simulations are powerful tools for computing this critical parameter. Nevertheless, the accurate and efficient calculation of the solvation free energy becomes a formidable endeavor when dealing with complex systems characterized by potent Coulombic interactions and sluggish ion dynamics and, consequently, slow transition across various metastable states. Here, in the present study, we expose limitations stemming from the conventional calculation of the statistical inefficiency g in the thermodynamic integration method, a factor that can hinder the determination of convergence of the solvation free energy and its associated uncertainty. Instead, we propose a robust scheme based on Gelman–Rubin convergence diagnostics. We leverage this improved estimation of uncertainties to introduce an innovative accelerated thermodynamic integration method based on the Gaussian Process regression. This methodology is applied to the calculation of the solvation free energy of trivalent rare-earth elements immersed in ionic liquids, a scenario in which the aforementioned challenges render standard approaches ineffective. The proposed method proves to be effective in computing solvation free energy in situations where traditional thermodynamic integration methods fall short.

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

Occurence and implications of silicon nitride in enstatic chondrites

Silicon nitride Si3N4, has previously been observed to be a common constituent of acid residues of Qingzhen (EH3) and Indarch (EH4). Ion probe analysis of the Si, N and C isotopic compositions of individual Si3N4 grains from Qingzhen and Indarch acid residues suggest most, if not all, grains are Solar System in origin. A few grains have isotopically anomalous C but this is probably due to small presolar SiC grains adhering to them. In situ observations of the Si3N4 in Qingzhen show that it is only present within, and probably exsolved from, host phases which contain elemental Si in solid solution. Thermodynamic calculations suggest that the Si3N4 probably formed during metamorphism and not in the nebula. Thermodynamic calculations also show that sinoite (Si2N2O) and not Si3N4 should be the stable phase during metamorphism. It appears that kinetic factors must have inhibited the formation if sinoite in Qimgzhen and Indarch.

Alexander, C. M. O'D.↗
DOE OSTI2022

Thermodynamics of electrolyte solutions near charged surfaces: Constant surface charge vs constant surface potential

Electric double layers are ubiquitous in science and engineering and are of current interest, owing to their applications in the stabilization of colloidal suspensions and as supercapacitors. While the structure and properties of electric double layers in electrolyte solutions near a charged surface are well characterized, there are subtleties in calculating thermodynamic properties from the free energy of a system with charged surfaces. These subtleties arise from the difference in the free energy between systems with constant surface charge and constant surface potential. In this work, we present a systematic, pedagogical framework to properly account for the different specifications on charged bodies in electrolyte solutions. Our approach is fully variational—that is, all free energies, boundary conditions, relevant electrostatic equations, and thermodynamic quantities are systematically derived using variational principles of thermodynamics. We illustrate our approach by considering a simple electrolyte solution between two charged surfaces using the Poisson–Boltzmann theory. Our results highlight the importance of using the proper thermodynamic potential and provide a general framework for calculating thermodynamic properties of electrolyte solutions near charged surfaces. Specifically, we present the calculation of the pressure and the surface tension between two charged surfaces for different boundary conditions, including mixed boundary conditions.

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

Experimental evaluation of the role of redox during glauconite-CO 2 -brine interactions

Greensands formations are globally abundant sedimentary rocks rich in Fe clays (typically glauconite) that commonly contain natural hydrocarbon accumulations and may be important reservoirs for geologic storage of anthropogenic CO 2 . Diagenesis in greensands is commonly accompanied by the conversion of primary glauconite to siderite (FeCO 3 ), a process that could be exploited for the permanent trapping of CO 2 . Importantly, siderite formation after glauconite requires that the mostly oxidized Fe in the primary Fe clay minerals is reduced during diagenetic interactions. Here, we explore the effect of solution redox state on the stability of glauconite in sandstones with implications for the diagenetic and/or engineered formation of siderite. We performed two flow-through experiments on intact, glauconite-rich sandstone cores at 150 °C and 150 bar. Both experiments employed a 1 mol NaCl/kg, 0.1 mol NaHCO 3 /kg solution charged with ~0.58 mol CO 2 /kg solution, but the redox state of the injected fluid was manipulated between experiments in order to compare glauconite reactivity and siderite saturation state at oxidizing and reducing end-member conditions. After reaction with the oxidizing (O 2 (aq) ≈ 6 μmol/kg) fluid, chemical and Mӧssbauer spectroscopic analyses indicate the production of Fe(III)-oxy/hydroxide minerals from glauconite, whereas, in the reducing (H 2 (aq) ≈ 5–40 mmol/kg) experiment, thermodynamic calculations and coupled chemical, mineralogical, and Mӧssbauer analyses suggest glauconite dissolution and precipitation of an Fe(II) mineral, likely siderite, and minor magnetite formation. These experimental results, along with thermodynamic calculations, confirm that solution redox state is the master variable dictating siderite formation in greensands.

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

GeoMicro3D: A novel numerical model for simulating the reaction process and microstructure formation of alkali-activated slag

Highlights: • The GeoMicro3D model was innovatively developed for simulating the reaction process and microstructure formation of AAS. • In GeoMicro3D, the dissolution of slag and reactions of aqueous ions are fundamentally described in the local lattice cell. • A novel approach is proposed to speed up the nucleation simulations. • GeoMicro3D was implemented and verified with the relevant experimental data and thermodynamic calculation results. For the first time, this study developed a novel model, named GeoMicro3D, to simulate the reaction process and microstructure formation of alkali-activated slag. The GeoMicro3D model consists of four modules that are designed to simulate, respectively: (i) the initial spatial distribution of real-shape slag particles in alkaline activator, (ii) the dissolution of slag and diffusion of ions via the transition state theory and lattice Boltzmann method, respectively, (iii) the spatial distribution of reaction products using a nucleation probability theory, and (iv) the chemical reactions with thermodynamic modelling. Afterwards the GeoMicro3D model was implemented and verified. The simulation results were discussed and compared with the relevant experimental data and thermodynamic calculation results using GEMS. A good agreement was found in the comparisons, showing the strong simulation capability of GeoMicro3D.

36 MATERIALS SCIENCE↗