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 91 records · Page 5

Review of computational approaches to predict the thermodynamic stability of inorganic solids

Improvements in the efficiency and availability of quantum chemistry codes, supercomputing centers, and open materials databases have transformed the accessibility of computational materials design approaches. Thermodynamic stability predictions play a central role in the efficacy of these approaches and should be considered carefully. This review covers the fundamentals of calculating thermodynamic stability using first-principles methods. Stability is delineated into two main topics — stability with respect to decomposition into competing phases and stability with respect to phase transition into alternative structures at fixed composition. For each topic, a summary of the state-of-the-art is provided along with a tutorial overview of practical considerations. The application of machine learning to both kinds of stability predictions is also covered. Finally, the limitations of thermodynamic stability predictions are discussed within the context of predicting the synthesizability of materials.

36 MATERIALS SCIENCE↗

Sampling lattices in semi-grand canonical ensemble with autoregressive machine learning

Calculating thermodynamic potentials and observables efficiently and accurately is key for the application of statistical mechanics simulations to materials science. However, naive Monte Carlo approaches, on which such calculations are often dependent, struggle to scale to complex materials in many state-of-the-art disciplines such as the design of high entropy alloys or multi-component catalysts. To address this issue, we adapt sampling tools built upon machine learning-based generative modeling to the materials space by transforming them into the semi-grand canonical ensemble. Furthermore, we show that the resulting models are transferable across wide ranges of thermodynamic conditions and can be implemented with any internal energy model U, allowing integration into many existing materials workflows. We demonstrate the applicability of this approach to the simulation of benchmark systems (AgPd, CuAu) that exhibit diverse thermodynamic behavior in their phase diagrams. Finally, we discuss remaining challenges in model development and promising research directions for future improvements.

36 MATERIALS SCIENCE↗

Generating Tables of Thermodynamic Properties

Thermodynamic Property Generation program produces tables for any fluid for which state equation, vapor pressure equation, specific heat equation, and critical constants are known. Program calculates thermodynamic properties at any state providing way of obtaining properties at extreme environmental conditions.

Nguyen, S. N.↗

Selenite and Selenate Sequestration during Coprecipitation with Barite: Insights from Mineralization Processes of Adsorption, Nucleation, and Growth

Coprecipitation of selenium oxyanions with barite is a facile way to sequester Se in the environments. However, the chemical composition of Se-barite coprecipitates usually deviates from that predicted from thermodynamic calculations. This discrepancy was resolved by considering variations in nucleation and growth rates controlled by ion–mineral interactions, solubility, and interfacial energy. For homogeneous precipitation, ~10% of sulfate, higher than thermodynamic predictions (<0.3%), was substituted by Se(IV) or Se(VI) oxyanion, which was attributed to adsorption-induced entrapment during crystal growth. For heterogeneous precipitation, thiol- and carboxylic-based organic films, utilized as model interfaces to mimic the natural organic-abundant environments, further enhanced the sequestration of Se(VI) oxyanions (up to 41–92%) with barite. Such enhancement was kinetically driven by increased nucleation rates of selenate-rich barite having a lower interfacial energy than pure barite. In contrast, only small amounts of Se(IV) oxyanions (~1%) were detected in heterogeneous coprecipitates mainly due to a lower saturation index of BaSeO 3 and deprotonation degree of Se(IV) oxyanion at pH 5.6. Finally, these roles of nanoscale mineralization mechanisms observed during composition selection of Se-barite could mark important steps toward the remediation of contaminants through coprecipitation.

54 ENVIRONMENTAL SCIENCES↗

Thermochemistry of Rare Earth Silicates for Environmental Barrier Coatings

Rare earth silicates are promising candidates as environmental protective coatings (EBCs) for silica-forming ceramics and composites in combustion environments since they are predicted to have lower reactivity with the water vapor combustion products. The reactivity of rare earth silicates is assessed by the thermodynamic activity of the silica component which is best measured by Knudsen effusion mass spectrometry (KEMS). Here, we discuss a novel method based on a reducing agent to increase the partial pressure of SiO(g) which is then used to calculate thermodynamic activity of silica in Y2O3-SiO2 and Yb2O3-SiO2 systems. After the KEMS measurements, samples were probed by X-ray diffraction and their phase content was calculated from Rietveld refinement.

Rare Earth Silicates↗

Cycling phosphorus on the Archean Earth: Part I. Continental weathering and riverine transport of phosphorus

Phosphorus (P) is the key nutrient thought to limit primary productivity on geological timescales. Phosphate levels in Archean marine sediments are low, but quantification of the P cycle and how it changed through a billion years of recorded Archean history remain a challenge, hindering our understanding of the role played by P in biosphere/geosphere co-evolution on the early Earth. Here, we design kinetic and thermodynamic models to quantitatively assess one key component of the early P cycle – continental weathering – by considering the emergence and elevation of continents, as well as the evolution of climate, the atmosphere, and the absence of macroscopic vegetation during the Archean Eon. Our results suggest that the weathering rate of apatite, the major P-hosting mineral in the rocks, was at least five times higher in the Archean Eon than today, attributable to high levels of pCO2,g. Despite this, the weathering flux of P to the oceans was negligible in the early Archean Eon, increasing to a level comparable to or greater than the modern by the end of eon, a consequence of accelerating continental emergence. Furthermore, our thermodynamic calculations indicate high solubilities of primary and secondary P-hosting minerals in the acidic weathering fluids on land, linking to high Archean pCO2,g. Thus, weathering of P was both kinetically and thermodynamically favorable on the Archean Earth, and river water could transport high levels of dissolved P to the oceans, as also supported by the observed P-depletion in our new compilation of Archean paleosols. Lastly, we evaluated the relative rates of physical erosion and chemical weathering of silicates during the Archean Eon. The results suggest that continental weathering on the early and middle Archean Earth might have been transport-limited due to low erosion rates associated with limited subaerial emergence and low plateau elevations; by the late Archean, however, continental weathering would have transited to kinetically-limited state because of continental emergence, increased plateau elevation, and weakening weathering rates. Overall, our weathering calculations together with paleosol evidence indicate an increasing flux of bioavailable P to the oceans through time, associated with late Archean continental emergence, reaching levels comparable to or higher than modern values by the end of the eon. Increased P fluxes could have fueled increasing rates of primary production, including oxygenic photosynthesis, through time, facilitating the irreversible oxidation of the Earth’s atmosphere early in the Proterozoic Eon.

Phosphorus↗

Density functional theory calculations of the thermodynamic and kinetic properties of point defects in β -U

Density functional theory (DFT) calculations of the thermodynamic and kinetic properties of point defects in the β phase of uranium are reported. Defect energies and entropies were calculated using 2 x 2 x 2 supercells and the Generalized Gradient Approximation (GGA) of Perdew-Burke-Ernzerhof (PBE) for the exchange-correlation potential. Due to computational cost, calculations of the vibrational properties governing entropies were performed by only displacing atoms within (roughly) the 3rd or 4th nearest neighbor shell of the defect, which implicitly assumes that atoms beyond this distance are unaffected by the defect. Migration barriers were estimated by nudged elastic band (NEB) calculations. The low symmetry of the β-U phase (the unit cell is tetragonal and contains 30 atoms) results in many point defect configurations and even more migration pathways. A connectivity map, starting from the most stable point defects, was developed in order to identify the rate-limiting step controlling the net diffusion rate in each crystallographic direction. The uranium self-diffusivity tensor was calculated by combining the defect formation energies, entropies, migration barriers and attempt frequencies. The fastest diffusion rate was determined to be a vacancy mechanism in the z crystallographic direction. The predicted uranium self-diffusivity for this mechanism agrees well with available experimental data. Finally, the diffusion mechanisms and rates identified in this study will inform fuel performance models of swelling and gas evolution.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Surface Phase Stability of Fe 2 O 3 (001) in Hydrogen Reducing Environments: A DFT and XPS Analysis

Here, this study combines density functional theory (DFT) and ab initio thermodynamics calculations with X-ray photoelectron spectroscopy (XPS) investigations to identify the reduction properties of the Fe 2 O 3 (001) surface with implications for corrosion resistance, hydrogen transport, and energy safety. Ab initio thermodynamics modeling predicts fully hydroxylated surface stability across a broad range of pressures (1 × 10 –23 to 1 × 10 5 mbar) and temperatures below 700 K, consistent with previous experimental studies. Above 800 K, exposures to 1 × 10 –4 mbar H 2 , 1 × 10 –4 mbar O 2 , or 1 × 10 –4 mbar H 2 + 1 × 10 –4 mbar O 2 each yield unique XPS signals indicating a loss of −OH coverage, aligning with DFT predictions. Insight into the mechanism of reduction as a function of H 2 exposure is provided, as well as conditions that promote further reduction toward Fe 3 O 4 . Theoretical and experimental investigations indicate the ability to maintain the Fe 2 O 3 protective layer of iron oxides that have been exposed to H 2 environments by including trace amounts of aqueous O2.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Evidence from massive siderite beds for a CO2-rich atmosphere before approximately 1.8 billion years ago

It is generally thought that, in order to compensate for lower solar flux and maintain liquid oceans on the early Earth, methane must have been an important greenhouse gas before approximately 2.2 billion years (Gyr) ago. This is based upon a simple thermodynamic calculation that relates the absence of siderite (FeCO3) in some pre-2.2-Gyr palaeosols to atmospheric CO2 concentrations that would have been too low to have provided the necessary greenhouse effect. Using multi-dimensional thermodynamic analyses and geological evidence, we show here that the absence of siderite in palaeosols does not constrain atmospheric CO2 concentrations. Siderite is absent in many palaeosols (both pre- and post-2.2-Gyr in age) because the O2 concentrations and pH conditions in well-aerated soils have favoured the formation of ferric (Fe3+)-rich minerals, such as goethite, rather than siderite. Siderite, however, has formed throughout geological history in subsurface environments, such as euxinic seas, where anaerobic organisms created H2-rich conditions. The abundance of large, massive siderite-rich beds in pre-1.8-Gyr sedimentary sequences and their carbon isotope ratios indicate that the atmospheric CO2 concentration was more than 100 times greater than today, causing the rain and ocean waters to be more acidic than today. We therefore conclude that CO2 alone (without a significant contribution from methane) could have provided the necessary greenhouse effect to maintain liquid oceans on the early Earth.

Greenhouse Effect↗

A statistical mechanical calculation of the thermodynamic properties of interstitial solid solutions involving second nearest neighbor interactions.

A previous calculation of the thermodynamic properties of interstitial solid solutions based on the technique of Kirkwood expansions has been extended to include the effects of second nearest neighbor solute atom mutual interactions. The error inherent in the first order (or quasi-chemical) counting of the degeneracy of the solution crystal is avoided. It is shown that, at high temperatures, even strong second nearest neighbor solute mutual interactions have a negligible effect on the entropy of the solution and a small, temperature-dependent effect on the solute partial enthalpy.

Alex, K.↗

The Formation of Wassonite: A New Titanium Monosulfide Mineral in the Yamato 691 Enstatite Chondrite

Wassonite, ideally stoichiometric TiS, is a titanium monosulfide not previously observed in nature, that was discovered within the Yamato 691 EH3 enstatite chondrite [1]. Because of the submicrometer size of the wassonite grains, it was not possible to determine conventional macroscopic properties. However, the chemical composition and crystal structure were well constrained by extensive quantitative energy dispersive x-ray analysis and electron diffraction using transmission electron microscopy (TEM). The crystal system for wassonite is rhombohedral (a = 3.42 plus or minus 0.07, c = 26.50 plus or minus 0.53 Angstroms) with space group: R(sup 3 raised bar) m (R9 type), cell volume: 268.4 plus or minus 0.53 Angstroms(sup 3), Z=9, density (calculated): 4.452 grams per cubic centimeter, empirical formula: (Ti(sub 0.93), Fe(sub 0.06), Cr(sub 0.01))S. In this study, we discuss possible formation mechanisms of wassonite and its associated minerals based on the petrology, mineralogy, crystallography, thermodynamic calculations, Al/Mg isotopic systematics and the O-isotopic composition of the wassonite-bearing BO chondrule.

Nakamura-Messenger, K.↗

Thermodynamics of rock forming crystalline solutions

Analysis of phase diagrams and cation distributions within crystalline solutions as means of obtaining thermodynamic data on rock forming crystalline solutions is discussed along with some aspects of partitioning of elements in coexisting phases. Crystalline solutions, components in a silicate mineral, and chemical potentials of these components were defined. Examples were given for calculating thermodynamic mixing functions in the CaW04-SrW04, olivine-chloride solution, and orthopyroxene systems.

Saxena, S. K.↗

Can enstatite meteorites form from a nebula of solar composition

Doubts have been expressed as to whether solar matter is sufficiently reducing to explain the minerals occurring in the enstatite chondrites and in the enstatite achondrites. Thermodynamic calculations on the stabilities of TiN, Si2N2O, CaS and silicon-bearing iron metal show that these substances can form under equilibrium conditions from a nebula of solar composition provided that the total pressure exceeds about 1 atm and that thermodynamic equilibria are frozen in at near-formation temperatures.

Herndon, J. M.↗

MaterialsMap: A CALPHAD-based tool to design composition pathways through feasibility map for desired dissimilar materials, demonstrated with resistance spot welding joining of Ag-Al-Cu

Assembly of dissimilar metals can be achieved by different methods, for example, casting, welding, and additive manufacturing (AM). However, undesired phases formed in liquid-phase assembling processes due to solute segregation during solidification diminish mechanical and other properties of the processed parts. In the present work, an open-source software named MaterialsMap, has been developed based on the CALculation of Phase Diagrams (CALPHAD) approach. The primary objective of MaterialsMap is to facilitate the design of an optimal composition pathway for assembling dissimilar alloys with liquid-phases based on the formation of desired and undesired phases along the pathway. In MaterialsMap, equilibrium thermodynamic calculations are used to predict equilibrium phases formed at slow cooling rate, while Scheil-Gulliver simulations are employed to predict non-equilibrium phases formed during rapid cooling. Additionally, by combining these two simulations, MaterialsMap offers a thorough guide for understanding phase formation in various manufacturing processes, assisting users in making informed decisions during material selection and production. As a demonstration of this approach, a compositional pathway was designed from pure Al to pure Cu through Ag using MaterialsMap. The design was experimentally verified using resistance spot welding (RSW).

36 MATERIALS SCIENCE↗

The bearing of spinel cataclasites on the crust-mantle structure of the moon

Subsolidus thermodynamic calculations have been made to define the temperature and pressure conditions required to equilibrate lunar spinel cataclasites (olivine + high alumina orthopyroxene + pleonaste spinel + plagioclase + or - cordierite) that occur as clasts in 15445, 73263, and 72435. The results, which are subject to modification by improved thermodynamic data and experiment, indicate that those samples that are cordierite-free and of high Mg/(Mg + Fe) were derived from the lower crust and possibly from a high-velocity zone of the uppermost mantle. However, the cordierite-bearing type in 72435,8 /low Mg/(Mg + Fe)/ resided in the upper levels of the crust prior to excavation by impact. Consideration of the relevant supersolidus phase equilibria indicates that the whole-rock chemistry of all spinel cataclasites can only be explained by pleonaste spinel accumulation. These materials are interpreted to be primordial cumulate rocks formed during the differentiation of the lunar magma ocean.

Herzberg, C. T.↗

Assessment of thermodynamic stability of sapphire in eutectic molten chloride environment

The continued development of molten salt reactors requires highly efficient and stable instruments that can efficiently track monitor conditions of the molten salt during long-term application operation in both the fuel and coolant/heat transfer fluid loops. Sapphire (Al 2 O 3 ) fibers have shown tremendous potential by enabling radiation resistance and a broader operational range of temperature . Thermodynamic modeling using the ThermoCalc and the SGTE Molten Salts (SALT1) and Pure Substances (Pure5) database is done performed to understand the compatibility of Al 2 O 3 fibers with NaCl-MgCl 2 eutectic molten salt in the temperature range of 600 to 1000°C. The thermodynamic calculations show that sapphire fibers are not expected to be compatible with long term exposure to molten chloride salts within this temperature range. Furthermore, various pure metallic elements were evaluated as potential cladding materials for part of the Al 2 O 3 fibers. Based on the thermodynamic analysis, molybdenum (Mo) and nickel (Ni) could be effective cladding materials to enhance the stability of Al 2 O 3 in NaCl and MgCl 2 chloride salt molten bath in the desired temperature range. The presence of a metallic cladding leads to the suppression of aluminum dissolution from Al 2 O 3 , thus improving its stability in chloride salts. Additionally, it is also shown that Al 2 O 3 remains stable up to 2300 K in the presence of pre-existing Al 2 MgO 4 and Al 2 NiO 4 in the eutectic molten chloride bath environment.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Understanding the microstructural stability in a y’-strengthened Ni-Fe-Cr-Al-Ti alloy

Ni-Fe-Cr-Al-Ti alloys, with Ni levels near 50 wt.%, have the potential to develop a microstructure consisting of a face-centered cubic ? matrix with the homogeneous precipitation of fine ordered ?’ precipitates similar to traditional Ni-based superalloys with significantly greater Ni content. Scanning electron microscopy (SEM), transmission electron microscopy (TEM), atom probe tomography (APT), and CALPHAD -based thermodynamic modeling were employed to understand the phase stabilities and microstructural evolution in an age-hardenable Ni-27Fe-18Cr-1Co-1.6Al-3.75Ti-1.2Mo-0.03C (wt%) alloy. The primary heat-treatment of solution annealing at 1121°C for 4h and age-hardening treatment at 760 °C for 16h resulted in a microstructure consisting of fine ?’ precipitates in an austenitic matrix along with grain boundary carbides, consistent with thermodynamic calculations. Long-term aging at 900 °C for 250h resulted in the coarsening of ?’ along with a change in the morphology of the precipitates from spherical to a more cuboidal shape. In addition, ? phase formation was observed concomitant with the partial dissolution of the ?’ phase. The ability of computational thermodynamic models to predict microstructural characteristics is discussed.

Gwalani, Bharat↗

Gas-Phase Combustion Synthesis of Aluminum Nitride Powder

Due to its combined properties of high electrical resistivity and high thermal conductivity aluminum nitride (AlN) is a highly desirable material for electronics applications. Methods are being sought for synthesis of unagglomerated, nanometer-sized powders of this material, prepared in such a way that they can be consolidated into solid compacts having minimal oxygen content. A procedure for synthesizing these powders through gas-phase combustion is described. This novel approach involves reacting AlCl3, NH3, and Na vapors. Equilibrium thermodynamic calculations show that 100% yields can be obtained for these reactants with the products being AlN, NaCl, and H2. The NaCl by-product is used to coat the AlN particles in situ. The coating allows for control of AlN agglomeration and protects the powders from hydrolysis during post-flame handling. On the basis of thermodynamic and kinetic considerations, two different approaches were employed to produce the powder, in co-flow diffusion flame configurations. In the first approach, the three reactants were supplied in separate streams. In the second, the AlCl3 and NH3 were premixed with HCl and then reacted with Na vapor. X-ray diffraction (XRD) spectra of as-produced powders show only NaCl for the first case and NaCl and AlN for the second. After annealing at 775 C tinder dynamic vacuum, the salt was removed and XRD spectra of powders from both approaches show only AlN. Aluminum metal was also produced in the co-flow flame by reacting AlCl3 with Na. XRD spectra of as-produced powders show the products to be only NaCl and elemental aluminum.

Axelbaum, R. L.↗