SEARCH · Engineering Papers
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.
Computer program subroutines provide the thermodynamic and transport properties of hydrogen in tabular form. The programs provide 18 combinations of input and output variables. This program is written in FORTRAN 4 for use on the IBM 7044 or CDC 3600 computers.
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.
Free energy minimization technique for computing thermodynamic properties of nonideal, complex lithium and hydrogen plasmas
We report that Q-POP-Thermo is a program designed to compute thermodynamic monodomain equilibrium states and their properties for ferroelectric single crystals and thin films based on the Landau-Ginzburg-Devonshire (LGD) Theory. Utilizing symbolic manipulation with the SymPy Library, the governing equations along with appropriate boundary conditions are solved for speedy minimization of the free energy of a crystal. Utilizing the popular Differential Evolution algorithm, with appropriate hybridization, multiple phase diagrams, such as the pressure-temperature phase diagram for bulk single crystals and the common strain-temperature phase diagram for monodomain thin-film systems can be readily generated. Furthermore, a variety of material properties of stable ferroelectric phases, including dielectric, piezoelectric, and electrocaloric properties, can simultaneously be calculated. Validation studies are presented for both thin-film and single crystal systems to test the effectiveness and capability of the open-source program.
Near horizons, quantum fields of low spin exhibit densities of states that behave asymptotically like 1 + 1 dimensional conformal field theories. In effective field theory, imposing some short-distance cutoff, one can compute thermodynamic quantities associated with the horizon, and the leading cutoff sensitivity of the heat capacity is found to equal to the leading cutoff sensitivity of the entropy. One can also compute contributions to the thermodynamic quantities from the gravitational path integral. For the cosmological horizon of the static patch of de Sitter space, a natural conjecture for the relevant heat capacity is shown to equal the Bekenstein-Hawking entropy. These observations allow us to extend the well-known notion of the generalized entropy to a generalized heat capacity for the static patch of de Sitter (dS). The finiteness of the entropy and the nonvanishing of the generalized heat capacity suggest it is useful to think about dS as a state in a finite dimensional quantum gravity model that is not maximally uncertain. Published by the American Physical Society 2024
Here, the accurate mapping of potential energy surfaces (PESs) is crucial to our understanding of the numerous physical and chemical processes mediated by atomic rearrangements, such as conformational changes and chemical reactions, and the thermodynamic and kinetic feasibility of these processes. Stochastic electronic structure theories, e.g., Quantum Monte Carlo (QMC) methods, enable highly accurate total energy calculations that in principle can be used to construct the PES. However, their stochastic nature poses a challenge to the computation and use of forces and Hessians, which are typically required in algorithms for minimum-energy pathway (MEP) and transition state (TS) identification, such as the nudged elastic band (NEB) algorithm and its climbing image formulation. Here, we present strategies that utilize the surrogate Hessian line-search method, previously developed for QMC structural optimization, to efficiently identify MEP and TS structures without requiring force calculations at the level of the stochastic electronic structure theory. By modifying the surrogate Hessian algorithm to operate in path-orthogonal subspaces and at saddle points, we show that it is possible to identify MEPs and TSs by using a force-free QMC approach. We demonstrate these strategies via two examples, the inversion of the ammonia (NH 3 ) molecule and the nucleophilic substitution (S N 2) reaction F – + CH 3 F → FCH 3 + F – . We validate our results using Density Functional Theory (DFT)- and Coupled Cluster (CCSD, CCSD(T))-based NEB calculations. We then introduce a hybrid DFT-QMC approach to compute thermodynamic and kinetic quantities, free energy differences, rate constants, and equilibrium constants that incorporates stochastically optimized structures and their energies, and show that this scheme improves upon DFT accuracy. Our methods generalize straightforwardly to other systems and other high-accuracy theories that similarly face challenges computing energy gradients, paving the way for highly accurate PES mapping, transition state determination, and thermodynamic and kinetic calculations at significantly reduced computational expense.
Algorithmic complexity is discussed as a computational counterpart to the second law of thermodynamics. It is shown that algorithmic complexity, which is a measure of randomness, sets limits on the thermodynamic cost of computations and casts a new light on the limitations of Maxwell's demon. Algorithmic complexity can also be used to define distance between binary strings.
Abstract Phase diagrams offer substantial predictive power for materials synthesis by identifying the stability regions of target phases. However, thermodynamic phase diagrams do not offer explicit information regarding the kinetic competitiveness of undesired by-product phases. Here we propose a quantitative and computable thermodynamic metric to identify synthesis conditions under which the propensity to form kinetically competing by-products is minimized. We hypothesize that thermodynamic competition is minimized when the difference in free energy between a target phase and the minimal energy of all other competing phases is maximized. We validate this hypothesis for aqueous materials synthesis through two empirical approaches: first, by analysing 331 aqueous synthesis recipes text-mined from the literature; and second, by systematic experimental synthesis of LiIn(IO 3 ) 4 and LiFePO 4 across a wide range of aqueous electrochemical conditions. Our results show that even for synthesis conditions that are within the stability region of a thermodynamic Pourbaix diagram, phase-pure synthesis occurs only when thermodynamic competition with undesired phases is minimized.
Computer program for calculating thermodynamic performance of Brayton cycle space power systems
SPECIES program computes thermodynamic and transport properties and equilibrium constants. Uses least-squares polynomial curve fits of most accurate data available to provide requested values more quickly than table-lookup methods. Subroutines computing transport coefficients and collision cross sections include additional code to correct for any electron pressure when working with ionic species. Species represented: diatomic nitrogen, diatomic oxygen, atomic nitrogen, atomic oxygen, nitric oxide, ionized nitric oxide, free electron, ionized atomic nitrogen, ionized atomic oxygen, ionized diatomic nitrogen, and ionized diatomic oxygen. Written in FORTRAN 77.
Thermodynamic modeling of the molten salt systems is critical in the design of advanced nuclear power plants, materials recycling, and many other important clean energy applications. As one of the most capable computational thermodynamic approaches, the CALPHAD (Calculation of Phase Diagrams) method has been widely used for multicomponent thermodynamic database development and further support more physics-based modeling. Although considerable CALPHAD modeling efforts for molten salt systems have been made, we found that no evaluation has been performed regarding the unary thermodynamic models for the molten salt systems despite the fact that more than one pure substance database is available for the CALPHAD modeling. Therefore, the thermodynamic descriptions of the pure salts should be critically evaluated to ensure a high-fidelity molten salt thermodynamic database. In this work, we comprehensively analyzed some selected molten salts using different available thermodynamic descriptions from two databases as an example. One is the SSUB database released by the SGTE (Scientific Group Thermodata Europe), and the other is the FactPS database from the FactSage. Finally, the significant discrepancies observed in the thermodynamic description between the existing CALPHAD pure substance databases and experiments call for an urgent effort to improve.
Batch-to-batch variation in powder compositions for binder jet additive manufacturing (BJAM) can significantly deter defining an “ideal” sintering window for a given alloy. One way to overcome the problem is by running sintering experiments at various temperatures for each batch of the powder. However, such an approach increases the time required to achieve large-scale production of parts. The predictive capabilities of computational thermodynamic tools like CALPHAD can be leveraged to overcome the challenge, especially for binder jet additive manufacturing, since the process occurs under near-equilibrium conditions. However, calculating the sintering window using CALPHAD can be computationally expensive, considering many possible feedstock compositions within “specification”. Here, we generate high throughput CALPHAD data for nickel-based superalloys to develop machine learning models to predict the sintering window rapidly. The predictive capability of the models has been validated using published results on BJAM of Inconel 718 and 625. Further, validated models are lightweight and can be deployed in an industrial setting to get sintering window in an accelerated manner.
This is the final technical report for the DOE MURI project on developing high-operating temperature heat transfer fluids for solar thermal power generation. This project aimed to evaluate liquid metals as a heat-transport fluid and potential energy storage medium for concentrated solar powered (CSP) systems. The project consisted of five phases, starting with fluid selection and screen and ending with potential CSP integration. The initial fluid selection arrived at heavy liquid metals as candidate fluids such as lead bismuth eutectic (LBE) as the prime liquid-alloy candidate to be studied, given the conditions above. Fifty-seven liquid compositions were initially identified with liquidus temperatures lower than 200°C, composed of elements that provide no driving force to the reduction of Cr 2 O 3 or Al 2 O 3 passivation layers in pipe material which is required for corrosion prevention. LBE was identified as the primary candidate fluid after these compositions were compared against computational thermodynamic modeling to screen all of the selected systems that would be compatible with our design objectives. Twenty-eight of these compositions were also investigated with the rapid characterization methods; the accuracies of the liquidus temperature and the specific heats predicted by the computations were better than 3 percent and 9 percent, respectively. We conducted numerus corrosion tests of various Fe-Cr-Al based steels in LBE and its variants at temperate as high as 800°C. It was found that minor additions to the liquid did not provide sufficient corrosion protection to warrant their use in real applications. Furthermore, we investigated the detailed mechanism of corrosion between heavy liquid metal and the pipe material at these very high temperatures. A liquid metal loop was constructed to operate at these high temperatures. The loop construction and operation lead to valuable lessons learned that will be needed for industrial scale loops. Furthermore, we learned that electromagnetic pumps at these extreme temperatures will have to be reengineered and redesigned to be used for CSP applications.
This paper presents data on a preliminary analysis of the thermal dynamic characteristics of the Airborne Information Management System (AIMS), which is a continuing design project at NASA Dryden. The analysis established the methods which will be applied to the actual AIMS boards as they become available. The paper also describes the AIMS liquid cooling system design and presents a thermodynamic computer model of the AIMS cooling system, together with an experimental validation of this model.
Synthesis remains a challenge for advancing materials science. A key focus of this challenge is how to enable selective synthesis, particularly as it pertains to metastable materials. This perspective addresses the question: how can “spectator” elements, such as those found in double ion exchange (metathesis) reactions, enable selective materials synthesis? By observing reaction pathways as they happen (in situ) and calculating their energetics using modern computational thermodynamics, we observe transient, crystalline intermediates that suggest that many reactions attain a local thermodynamic equilibrium dictated by local chemical potentials far before achieving a global equilibrium set by the average composition. Furthermore, using this knowledge, one can thermodynamically “shortcut” unfavorable intermediates by including additional elements beyond those of the desired target, providing access to a greater number of intermediates with advantageous energetics and selective phase nucleation. Ultimately, data-driven modeling that unites first-principles approaches with experimental insights will refine the accuracy of emerging predictive retrosynthetic models for complex materials synthesis.
The current project was based on a combined experimental and computational approach, which can help recommend better cathode materials under multiple impurities conditions. The PI will mainly take in charge of experimental and computational thermodynamics of the selected cathode materials for the SOFC applications under multiple impurities. While the co-PI will run the electrochemical tests of the cathodes recommended and eventually the long-term degradation tests. At the end of the project, a multiple tolerant cathode material based on the combined experimental and computational approach will be recommended and the reliability of the commonly used accelerated testing will be evaluated. It will address multiple impurities poisoning effect of SO2, CO2, Cr and H2O on the LSM, LSCF and LNO cathodes by identifying the formation of the detrimental secondary phases by XRD, SEM and TEM techniques. And further recommended cathode material will be subjected to electrochemical testing and the most promising ones will be applied to long-term tests. The hybrid approach the PI proposed will not only be applied to the design of multiple impurities tolerant cathodes but will also be considered in the future oxygen electrode applications in SOECs or reversible SOCs. This hybrid computational and experimental approach includes four sections: 1) Investigation of single impurity poisoning on LSM, LSCF and LNO cathodes in the presence of SO2, CO2, Cr. In this section, LSM, LSCF cathodes from FuelCellMaterials and LNO cathodes from Sol-Gel synthesis will be heat-treated in the above single impurity. And the formation of the secondary phases as well as the corresponding simulations will be cross compared, which shows good agreement between each other. Meanwhile, the accelerated testing approach will be evaluated in these systems compared with the previous published work to understand the reliability of the approach. 2) Investigation of multiple impurities poisoning on LSM, LSCF and LNO cathodes in the presence of Cr+H2O, SO2+Cr and SO2+Cr+H2O conditions. We have also applied these 3 candidate cathodes under these multiple impurities’ conditions and the long-term degradation mechanism of the multiple impurities will be understood with the help of the combined experimental and computational approach. Meanwhile, the synergistic effect of those impurities will be compared with the individual ones in the same cathode system to further reveal the actual operating conditions. 3) Electrochemical testing and polarization of the recommended cathode. We have demonstrated very low polarization resistance in LSCF (core)-LSM (shell) electrodes using MSD process and analyzed impedance spectra using DRT analysis and confirming that the low polarization resistance in LSM infiltrated MSD cells is due to reduction in polarization resistance associated with O2-adsorption process. 4). Long-term degradation testing. We have developed a versatile MSD-based process to deposit various continuous coatings onto porous scaffolds and established baseline for longer term Cr-impurity testing in future projects.
This computational tool provides the standard Gibbs energy of formation for aqueous species, the dielectric constant and ionization constant of water, over a wide range of temperatures, pressures, and densities both above and below the critical point of water.
The TChem open-source software is a toolkit for computing thermodynamic properties, source term, and source term’s Jacobian matrix for chemical kinetic models that involve gas and surface reactions.