Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Computational thermodynamics”

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 55 records · Page 3

Dissociation of gas-phase anisole induced by low-energy electron interactions: understanding patterns of aromatic bond cleavage

Abstract In view of elucidating the fragmentation patterns of aromatic systems induced by low-energy electron interactions, dissociative electron attachment (DEA) to gas-phase anisole was performed. Anionic fragments resulting from this DEA process were detected by a quadrupole mass spectrometer, and ion yields of those fragments as a function of incident electron energy were rendered. Our study showed the formation of CH 3 − , HCC − , and OCH 3 − fragments, suggesting that various dissociation channels proceed out of DEA to anisole. We employed density functional theory to compute thermodynamic threshold energies for each potential dissociation channel. Those theoretical calculations supported the prediction that the CH 3 − and OCH 3 − fragments form via mechanisms of single-bond cleavage; the HCC − fragments may form through two-, three-, or four-body dissociation channels that entail hydrogen transfers and the cleavage of multiple aromatic bonds. The experimental resonance energies that form the CH 3 − , HCC − , and OCH 3 − fragments were 6.0 eV, 5.8 and 9.7 eV, and 9.8 eV, respectively. Given the classification of anisole as a monosubstituted aromatic species, our results explain generalizable patterns of electron-mediated dissociation in aromatic systems.

Finley, Jacob (ORCID:0009000606963352)

QCD Constraints on Isospin-Dense Matter and the Nuclear Equation of State

Understanding the behavior of dense hadronic matter is a central goal in nuclear physics as it governs the nature and dynamics of astrophysical objects such as supernovae and neutron stars. Because of the nonperturbative nature of quantum chromodynamics (QCD), little is known rigorously about hadronic matter in these extreme conditions. Here, lattice QCD calculations are used to compute thermodynamic quantities and the equation of state of QCD over a wide range of isospin chemical potentials with controlled systematic uncertainties. Agreement is seen with chiral perturbation theory when the chemical potential is small. Comparison to perturbative QCD at large chemical potential allows for an estimate of the gap in the superconducting phase, and this quantity is seen to agree with perturbative determinations. Since the partition function for an isospin chemical potential μ I bounds the partition function for a baryon chemical potential μ B = 3 μ I / 2 , these calculations also provide rigorous nonperturbative QCD bounds on the symmetric nuclear matter equation of state over a wide range of baryon densities for the first time. Published by the American Physical Society 2025

73 NUCLEAR PHYSICS AND RADIATION PHYSICS

C-C Coupling Mechanism on Cu(100) A Molecular Dynamics Study at 298K

The electrochemical reduction of carbon dioxide (CO2) into valuable fuels such as C1 (syngas, methane) and C2 (ethylene, ethanol) products is a key strategy for achieving a carbon-neutral economy. Computational studies of C-C coupling, a critical step in CO2 reduction, are essential for designing more efficient catalysts. However, simulating these processes under realistic electrochemical conditions, including temperature and solvent effects, is computationally demanding. In this work, we develop a machine learning-based atomistic potential to study CO2 reduction on Cu(100) surfaces, accounting for temperature and explicit water solvent effects. We compute thermodynamic free energies of the possible C-C coupling pathways, CO*+CO*→OCCO*, CO*+CHO*→OCCHO*, CO*+COH*→OCCOH*, CHO*+COH*→OHCCOH*, COH*-COH*→HOCCOH*, and CHO*-CHO*→OHCCHO*. Our results quantify the thermodynamic tendencies of these reactions and reveal that, in addition to the well-established CO* + CO* → OCCO* pathway, CHO* is a critical intermediate in the formation of C2 products on Cu(100). Furthermore, we demonstrate that the machine learning approach offers a cost-efficient framework for studying CO2 reduction on diverse catalysts under realistic electrochemical conditions.

machine learning

Creep Property of Intermetallic Dispersive Steels for Nuclear Applications

To meet the materials performance demands of next-generation nuclear and high-temperature energy systems, a new class of intermetallic-dispersive steels (IDS) has been designed through integrated computational thermodynamics and alloy design strategies. The IDS alloys incorporate coherent L1₂ (γ′-Ni₃Al) nanoprecipitates within an Fe–Ni–Cr austenitic matrix, engineered to achieve both high temperature strength and thermodynamic stability while suppressing the formation of detrimental δ-Ni₃Nb and η-Ni₃Ti phases. Initial creep testing of Ti-IDS and TiTa-IDS alloys demonstrates rupture lives comparable to or exceeding those of Inconel 718 and ODS steels, despite being fabricated through conventional ingot metallurgy. Step-load creep tests identified a stress threshold near 300–350 MPa for the onset of tertiary creep, and in-situ neutron diffraction experiments on Ti-IDS revealed clear load partitioning between the matrix and γ′ precipitates: elastic strain is shared by both phases, while plastic strain localizes in the matrix. These results confirm that stable γ′–matrix interfaces play a dominant role in retarding dislocation motion and enhancing creep resistance. In FY26, the program will conduct repeat creep rupture testing of TiTa-IDS at 650 °C/400 MPa, initiate long-term rupture testing of TiNb-IDS, and perform TEM-based microstructural characterization to elucidate dislocation–precipitate interactions and microstructural stability. Collectively, these efforts will establish the mechanistic foundation for next-generation high-temperature structural materials based on the IDS concept.

22 GENERAL STUDIES OF NUCLEAR REACTORS

Radiative coupled viscous flow with massive blowing

An analysis of the fully-coupled viscous, radiating flow past an ablating blunt body at hyperbolic entry conditions is presented. A detailed thermodynamics computation, as well as a realistic radiation transport model, is included. A locally nonsimilar approach is employed to solve the conservation equations away from the stagnation point. The validity of the locally nonsimilar approach is demonstrated for some nonablating cases. Sample calculations are made for the typical flight condition of a Jovian entry probe. The effects of the downstream injection of the ablation products of a carbon heat shield on the flux distribution around the body are discussed in detail. It is found that most of the radiative energy is absorbed by the injected carbon gas and dumped into the wake.

Chou, Y. S.

Systems Improved Numerical Fluids Analysis Code

Systems Improved Numerical Fluids Analysis Code, SINFAC, consists of additional routines added to April, 1983, version of SINDA. Additional routines provide for mathematical modeling of active heat-transfer loops. Simulates steady-state and pseudo-transient operations of 16 different components of heat-transfer loops, including radiators, evaporators, condensers, mechanical pumps, reservoirs, and many types of valves and fittings. Program contains property-analysis routine used to compute thermodynamic properties of 20 different refrigerants. Source code written in FORTRAN 77.

Costello, F. A.

Stability of Titanium Nitride and Titanium Carbide When Exposed to Hydrogen Atoms from 298 to 1950 K

Titanium nitride and titanium carbide deposited on tungsten wires were exposed to hydrogen atoms (10(exp -4) atm pressure) produced by the action of microwave radiation on molecular hydrogen. The results of these experiments in the temperature range 298 to 1950 K indicate that no appreciable reaction takes place between atomic hydrogen and TiN or TiC. The formation of reaction products (NH3, CH4, C2H2) should be favored at lower temperatures. However, because of the high catalytic activity of Ti for H atom recombination, the rate of such reactions with H atoms is controlled by the rate of evaporation of Ti from the surface, this rate being low at temperatures below 1200 K. In order to interpret the stability of TiN and TiC in H atoms more fully, the stability of TiN and TiC in vacuum and H2 gas was also studied. The thermodynamic computations conform in order of magnitude to the experimentally found rates of decomposition of TiN and TiC in vacuum and are also consistent with the fact that no appreciable reaction is found with these compounds in molecular H2 at a pressure of 10(exp -3) atmosphere in the temperature range 2980 to 2060 K. When TiN or TiC was heated in atomic H or molecular H2, no reaction products other than those obtained from the simple decomposition of the nitride and carbide were observed. The gaseous products were analyzed in a mass spectrometer.

Philipp, Warren H.

Cloud Scaling Properties and Cloud Parameterization

Cloud liquid and cloud traction variability is studied as a function of horizontal scale in the ECMWF forecast model during several 10-day runs at the highest available model resolution, recently refined from approximately 60 km (T213) down to approximately 20 km (T639). At higher resolutions, model plane-parallel albedo biases are reduced, so that models may be tuned to have larger, more realistic, cloud liquid water amounts, However, the distribution of cloud liquid assumed -within- each gridbox, for radiative and thermodynamic computations, depends on ad hoc assumptions that are not necessarily consistent with observed scaling properties, or with scaling properties produced by the model at larger scales. To study the larger-scale cloud properties, ten locations on the Earth are chosen to coincide with locations having considerable surface data available for validation, and representing a variety of climatic regimes, scaling exponents are determined from a range or scales down to model resolution, and are re-computed every three hours, separately for low, medium and high clouds, as well as column-integrated cloudiness. Cloud variability fluctuates in time, due to diurnal, synoptic and other' processes, but scaling exponents are found to be relatively stable. various approaches are considered for applying computed cloud scaling to subgrid cloud distributions used for radiation, beyond simple random or maximal overlap now in common use. Considerably more work is needed to compare model cloud scaling with observations. This will be aided by increased availability of high-resolution surface, aircraft and satellite data, and by the increasing resolution of global models,

Cahalan, R. F.

Investigation of Conditions of Titanium Carbonization - IV

In a previous paper, results are presented of accurate investigations of the processes of titanium carbonization and the succeeding titanium carbide decarbonization as related to the phenomenon of the graphitization of soot by heating at a constant temperature in atmospheres of pure hydrogen and carbon monoxide. These tests showed that the processes of titanium carbonization-decarbonization in an atmosphere of pure gases without nitrogen proceed in the same direction as the analogous processes under the conditions of the production furnace. In this case, however, the presence of admixtures of nitrogen changes the quantitative results of the decarbonization process. Thermodynamic computations confirming the results of previous tests conducted at atmospheric pressure and additional tests of titanium carbonization at lowered pressures are presented herein.

Meerson, G. A.

Vacuum Melt Infiltration of 2D Tyranno SA3® Ceramic Matrix Composites with Cr-25(wt.%)Si Intermetallic Alloy

The potential and limitations of vacuum melt infiltrating woven preforms of 2D Tyranno SA3 fibers with Cr-25(wt.%)Si alloy were investigated. Microstructural observations revealed that the volume fractions of open voids infiltrated by the molten metal varied between 3.7 and 10.0% irrespective of the magnitudes of the melt infiltration temperatures and hold times, which was attributed in part to substantial volatilization of the silicide charge based on computational thermodynamics. In the absence of experimental data, poor wettability of the melt could not be ruled out unambiguously. Other possible reasons are carefully examined and discarded. Ring microstructures were observed around the fibers mainly at the outer surfaces of the infiltrated preform, where the BN and chemical vapor infiltrated SiC coatings had vaporized in the high vacuum environment. In contrast, they were absent in the interior of the preform, where the coatings were still present with no evidence of molten metal attack of the coatings. Details of thermodynamic calculations are presented to confirm the volatilization of the BN and CVI SiC coatings thereby resulting in the ring microstructures due to the filling of voids by the molten silicide. It is concluded that successful melt infiltration of Tyranno SA3® with Cr-25%Si can be only achieved under sufficient positive pressure to prevent volatilization of the coatings and the molten silicide.

Cr-25%Si

CEA2022: A Modernization of NASA Glenn’s Software CEA (Chemical Equilibrium with Applications)

The software program “Chemical Equilibrium with Applications” (CEA) is used to solve chemical equilibrium, and compute thermodynamic and transport properties of the resulting mixture, and also has special solvers dedicated to rocket, shock, and detonation problems. We have recently completed a full re-write of CEA with modernization and improvements, called “CEA2022”. In this paper, we will give an overview of CEA2022’s features, and discuss some of the fundamental equations used by CEA2022, as well as the fundamental assumptions, in order to provide users with a complete understanding of the software’s methodology. Several enhancements have been made to the software, including modern software development practices, interface improvements, and additional features. The feature enhancements include: running cases in parallel with thread safe solves, thermodynamic and transport database updates, and allowing for negative and inert reactants. In terms of interface improvements, we have made CEA a reusable library by adding APIs for multiple languages, including Python, Matlab, Excel, Fortran, and C. The subroutine interface allows for integration with other applications, including flow-solver integration (i.e. with CFD). We also compare results between CEA2022 and the previous version (CEA2) as a validation of the new software.

chemical equilibrium

NASA Chemical Equilibrium with Applications (CEA) Tutorial

The software program Chemical Equilibrium with Applications (CEA) is used to solve chemical equilibrium, and compute thermodynamic and transport properties of the resulting mixture, and also has special solvers dedicated to rocket, shock, and detonation problems. We will give an overview of CEA’s features and discuss some of the fundamental equations used by CEA, as well as the fundamental assumptions, in order to provide users with a better understanding of what exactly the software is doing. We have undertaken a full re-write of the software in the last two years, so we will also provide details about the new version of CEA as well.

chemical equilibrium

Computer program for thermodynamic analysis of open cycle multishaft power system with multiple reheat and intercool

A computer program to analyze power systems having any number of shafts up to a maximum of five is presented. On each shaft there can be as many as five compressors and five turbines, along with any specified number of intervening intercoolers and reheaters. A recuperator can be included. Turbine coolant flow can be accounted for. Any fuel consisting entirely of hydrogen and/or carbon can be used. The program is valid for maximum temperatures up to about 2000 K (3600 R). The system description, the analysis method, a detailed explanation of program input and output including an illustrative example, a dictionary of program variables, and the program listing are explained.

Glassman, A. J.

Computer program for calculating and fitting thermodynamic functions

A computer program is described which (1) calculates thermodynamic functions (heat capacity, enthalpy, entropy, and free energy) for several optional forms of the partition function, (2) fits these functions to empirical equations by means of a least-squares fit, and (3) calculates, as a function of temperture, heats of formation and equilibrium constants. The program provides several methods for calculating ideal gas properties. For monatomic gases, three methods are given which differ in the technique used for truncating the partition function. For diatomic and polyatomic molecules, five methods are given which differ in the corrections to the rigid-rotator harmonic-oscillator approximation. A method for estimating thermodynamic functions for some species is also given.

Mcbride, Bonnie J.

Computer program for calculating thermodynamic and transport properties of fluids

Computer code has been developed to provide thermodynamic and transport properties of liquid argon, carbon dioxide, carbon monoxide, fluorine, helium, methane, neon, nitrogen, oxygen, and parahydrogen. Equation of state and transport coefficients are updated and other fluids added as new material becomes available.

Hendricks, R. C.

Computer programs for thermodynamic and transport properties of hydrogen (tabcode-II)

The thermodynamic and transport properties of para and equilibrium hydrogen have been programmed into a series of computer routines. Input variables are the pair's pressure-temperature and pressure-enthalpy. The programs cover the range from 1 to 5000 psia with temperatures from the triple point to 6000 R or enthalpies from minus 130 BTU/lb to 25,000 BTU/lb. Output variables are enthalpy or temperature, density, entropy, thermal conductivity, viscosity, at constant volume, the heat capacity ratio, and a heat transfer parameter. Property values on the liquid and vapor boundaries are conveniently obtained through two small routines. The programs achieve high speed by using linear interpolation in a grid of precomputed points which define the surface of the property returned.

Roder, H. M.

General method and thermodynamic tables for computation of equilibrium composition and temperature of chemical reactions

A rapidly convergent successive approximation process is described that simultaneously determines both composition and temperature resulting from a chemical reaction. This method is suitable for use with any set of reactants over the complete range of mixture ratios as long as the products of reaction are ideal gases. An approximate treatment of limited amounts of liquids and solids is also included. This method is particularly suited to problems having a large number of products of reaction and to problems that require determination of such properties as specific heat or velocity of sound of a dissociating mixture. The method presented is applicable to a wide variety of problems that include (1) combustion at constant pressure or volume; and (2) isentropic expansion to an assigned pressure, temperature, or Mach number. Tables of thermodynamic functions needed with this method are included for 42 substances for convenience in numerical computations.

Huff, Vearl N