SEARCH · Search NASA
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 calculates the real fluid properties of normal or parahydrogen using a library of single function calls without initial estimates. Accurate transport and thermodynamic properties of molecular hydrogen are needed for advanced propulsion systems.
Experimental studies of the combusting mixing layers formed over a rearward facing step and between two parallel streams are reported. The steady behavior of the mean properties of the reacting flows is similar to that of nonreacting flows except for the global effects of thermodynamic property changes. Vortices form by the two-dimensional K-H instability and grow by subharmonic combination until the layer begins to interact with the walls. Unsteady behavior occurs when there are resonator interactions with the K-H instability or the instability associated with the recirculation vortex in the rearward facing step flow.
Environmental barrier coatings (EBCs) are a necessary enabling technology for the transition from superalloys to silicon carbide (SiC) ceramic matrix composites (CMCs) in gas turbine engines for increased efficiency and decreased fuel costs. SiC-based CMCs are prone to oxidation-based degradation in the engine hot section, and rare-earth (RE) silicates are promising candidates for EBCs due to their close thermal expansion match to the composite substrate and oxidation resistance. However, the design of EBCs is hindered by the large chemical space of candidate materials and the difficulty in obtaining material properties for engineering optimization. This is especially difficult as research continues into mixed-cation or “high-entropy” RE silicates. First-principles computational methods such as density functional theory (DFT) are highly effective at calculating material properties to guide coating design but are limited by their computational cost. Atomistic simulations have the potential to both accelerate property calculations and expand the properties able to be calculated due to their lower computational compared to DFT. However, they require interatomic potentials (IAPs) specific to the material system of interest, and, to our knowledge, there are no suitable IAPs for RE silicates. Machine learning (ML) is a promising technique to accelerate material property predictions indirectly by generating IAPs for atomistic simulations or via direct prediction. In this work, we present two ML approaches to accelerate the calculation of RE silicate properties relevant to EBC design: 1) a ML-derived interatomic potential (IAP) for atomistic simulations of yttrium disilicate (Y2Si2O7) from DFT training data, and 2) a neural network (NN) model to directly predict thermochemical properties of RE silicates and oxides directly from easily obtainable unit cell parameters. Classical MD simulations using the IAP yield lattice properties and bond lengths in good agreement with both DFT and experimental results from x-ray diffraction. Thermodynamic properties calculated using the finite-displacement phonon method and quasi-harmonic approximation were orders of magnitude faster than DFT with good agreement to the DFT results. The IAP was also used to calculate properties such as coefficient of thermal expansion (CTE) that require large simulation supercells and are therefore difficult with DFT. The IAP correctly predicted the anisotropic nature of the CTE in three different phases of Y2Si2O7. The NN model predicts constant pressure heat capacity, Cp, orders of magnitude faster than DFT calculations, which can enable its use as a surrogate model for multiscale simulations. The two methods presented in this work demonstrate the utility of ML for accelerating the prediction of RE silicate properties, which can in turn accelerate EBC design and optimization.
A self-consistent formalism is developed that, based upon a many-body potential, dynamically determines the thermodynamic properties of ionic clusters without an a priori designation of the equilibrium structures. Aggregates consisting of a single closed shell K(+) ion and N-1 isoelectronic argon atoms were studied. The clusters form crystallites at low temperatures, and melting transitions and spontaneous dissociations are indicated. The results confirm experimental evidence that shows that ionic clusters become less stable with increasing N. The crystallite structures formed by four different clusters are isosceles triangle, skewed form, octahedron with ion in the middle, and icosahedron with the ion in the middle.
The thermodynamic functions, transport properties, and physical properties of both liquid and gaseous hydrogen are presented. The low temperature regime is emphasized. The tabulation of the properties of normal hydrogen in both Si and engineering units is given along with the tabulation of parahydrogen.
Analyses of process system properties were continued for materials involved in the alternate processes under consideration for semiconductor silicon. Primary efforts centered on physical and thermodynamic property data for dichlorosilane. The following property data are reported for dichlorosilane which is involved in processing operations for solar cell grade silicon: critical temperature, critical pressure, critical volume, critical density, acentric factor, vapor pressure, heat of vaporization, gas heat capacity, liquid heat capacity and density. Work was initiated on the assembly of a system to prepare binary gas mixtures of known proportions and to measure the thermal conductivity of these mixtures between 30 and 350 C. The binary gas mixtures include silicon source material such as silanes and halogenated silanes which are used in the production of semiconductor silicon.
The new Web site program "ThermoBuild" allows users to easily access and use the NASA Glenn Thermodynamic Database of over 2000 solid, liquid, and gaseous species. A convenient periodic table allows users to "build" the molecules of interest and designate the temperature range over which thermodynamic functions are to be displayed. ThermoBuild also allows users to build custom databases for use with NASA's Chemical Equilibrium with Applications (CEA) program or other programs that require the NASA format for thermodynamic properties. The NASA Glenn Research Center has long been a leader in the compilation and dissemination of up-to-date thermodynamic data, primarily for use with the NASA CEA program, but increasingly for use with other computer programs.
The method of characteristics for a chemically reacting gas is used in the construction of the time-dependent, one-dimensional flow field resulting from the normal reflection of an incident shock wave at the end wall of a shock tube. Nonequilibrium chemical reactions are allowed behind both the incident and reflected shock waves. All the solutions are evaluated for oxygen, but the results are generally representative of any inviscid, nonconducting, and nonradiating diatomic gas. The solutions clearly show that: (1) both the incident- and reflected-shock chemical relaxation times are important in governing the time to attain steady state thermodynamic properties; and (2) adjacent to the end wall, an excess-entropy layer develops wherein the steady state values of all the thermodynamic variables except pressure differ significantly from their corresponding Rankine-Hugoniot equilibrium values.
When graphite is used as ablation material in heat shields, very often a prominent carbon species vaporized into the stream is C3. The thermodynamic properties of the C3 molecule are, therefore, important in calculating transport phenomena in the ablation flow field. The nature of the C3 thermodynamic functions has been in doubt because of the uncertain contribution of the bending mode vibrations to the total internal energy of the molecule. An approach for overcoming these difficulties is considered. The results of the computations are presented in the form of graphs and approximating functions.
Thermodynamic and transport combustion properties were calculated for a wide range of conditions for the reaction of hydrocarbons with air. Three hydrogen-carbon atom ratios (H/C = 1.7, 2.0, 2.1) were selected to represent the range of aircraft fuels. For each of these H/C ratios, combustion properties were calculated for the following conditions: Equivalence ratio: 0, 0.25, 0.5, 0.75, 1.0, 1.25 Water - dry air mass ratio: 0, 0.03 Pressure, kPa: 1.01325, 10.1325, 101.325, 1013.25, 5066.25 (or in atm: 0.01, 0.1, 1, 10, 50) Temperature, K: every 10 degrees from 200 to 900 K; every 50 degrees from 900 to 3000 K Temperature, R: every 20 degrees from 360 to 1600 R; very 100 degrees from 1600 to 5400 R. The properties presented are composition, density, molecular weight, enthalphy, entropy, specific heat at constant pressure, volume derivatives, isentropic exponent, velocity of sound, viscosity, thermal conductivity, and Prandtl number. Property tables are based on composites that were calculated by assuming both: (1) chemical equilibrium (for both homogeneous and heterogeneous phases) and (2) constant compositions for all temperatures. Properties in SI units are presented in this report for the Kelvin temperature schedules.
A computer program subroutine, WETAIR, was developed to calculate the thermodynamic and transport properties of air water mixtures. It determines the thermodynamic state from assigned values of temperature and density, pressure and density, temperature and pressure, pressure and entropy, or pressure and enthalpy. The WETAIR calculates the properties of dry air and water (steam) by interpolating to obtain values from property tables. Then it uses simple mixing laws to calculate the properties of air water mixtures. Properties of mixtures with water contents below 40 percent (by mass) can be calculated at temperatures from 273.2 to 1497 K and pressures to 450 MN/sq m. Dry air properties can be calculated at temperatures as low as 150 K. Water properties can be calculated at temperatures to 1747 K and pressures to 100 MN/sq m. The WETAIR is available in both SFTRAN and FORTRAN.
A FORTRAN IV computer program for the calculation of the thermodynamic and transport properties of complex mixtures is described. The program has the capability of performing calculations such as:(1) chemical equilibrium for assigned thermodynamic states, (2) theoretical rocket performance for both equilibrium and frozen compositions during expansion, (3) incident and reflected shock properties, and (4) Chapman-Jouguet detonation properties. Condensed species, as well as gaseous species, are considered in the thermodynamic calculation; but only the gaseous species are considered in the transport calculations.
The properties of the fundamental hydrogen species, proton (H+), hydrogen atom (H•), and hydride anion (H-) are critical to a vast range of chemical processes, yet their thermodynamic properties in nonaqueous solvents are not well established. A hybrid supermolecule-continuum approach is used to predict the Gibbs free energies of solvation (?G°solv) and standard redox potentials (E°) for the 2H+/H2 and H•/H- couples in acetonitrile (MeCN) and tetrahydrofuran (THF) following the approach previously used for water. Gas phase geometry optimizations were done using density functional theory (DFT) with the ?B97XD functional, Møller-Plesset second-order perturbation theory (MP2) with augmented correlation-consistent basis sets and composite-correlated molecular orbital theory (G3MP2 and G3MP2B3). Additional single point calculations in gas phase were performed for H- affinity for MeCN and THF at the coupled cluster CCSD(T) level using MP2 geometries. Solvation was included using the self-consistent reaction field SMD model. For H+, the ?G°solv values are predicted to be -252.2 kcal/mol in MeCN and -261.0 kcal/mol in THF. The redox potentials for the 2H+/H2 couple are +0.49 V in MeCN and +0.11 V in THF relative to the aqueous SHE, consistent with the available experimental data. The ?G°solv for H- is predicted to be -80 kcal/mol in MeCN and -68 in THF, and for H•, ?G°solv is predicted to be -1.3 kcal/mol in THF and -2.0 kcal/mol in MeCN. These solvation energies yield respective calculated redox potentials for the H•/H- couple of -0.14 V in MeCN and -0.68 V in THF relative to the aqueous SHE. H2 solvation is slightly exergonic in MeCN (-3.3 kcal/mol) and near thermoneutral in THF (+1.3 kcal/mol). The pKa(H2) is predicted to be 45.4 in MeCN and 47.8 in THF. The proton-coupled electron transfer (PCET) accounting term CG was evaluated giving values on the ferrocene scale within approximately 2 kcal/mol of literature in both solvents. These results provide benchmarks for hydrogen redox chemistry in MeCN and THF and improve consistency across E°, ?G°solv, pKa, and CG. These results provide benchmarks for hydrogen redox chemistry in nonaqueous media challenging previous estimates based on indirect extrapolations and assumptions.
Recent helicity-resolved magneto-Raman spectroscopy measurement demonstrates large effective phonon magnetic moments of ~2.5 μ B in monolayer MoS 2 , highlighting resonant excitation of bright excitons as a feasible route to activate Γ-point circularly polarized phonons in transition metal dichalcogenides. However, a microscopic picture of this intriguing phenomenon remains lacking. In this work, we show that an orbital transition between the split conduction bands (Δ 0 = 4 meV) of MoS 2 couples to the doubly degenerate E" phonon mode (Ω 0 = 33 meV), forming two hybridized states. Our phononic and electronic Raman scattering measurements capture these two states: (i) one with predominantly phonon contribution in the helicity-switched channels and (ii) one with primarily orbital contribution in the helicity-conserved channels. An orbital-phonon coupling model successfully reproduces the large effective magnetic moments of the circularly polarized phonons and explains their thermodynamic properties. Strikingly, the Raman mode from the orbital transition is superimposed on a strong quasi-elastic scattering background, indicating the presence of spin fluctuations. As a result, the electrons excited to the conduction bands through the exciton exhibit paramagnetic behavior although MoS 2 is generally considered as a nonmagnetic material. By depositing nanometer-thickness nickel thin films on monolayer MoS 2 , we tune the electronic structure so that the A exciton perfectly overlaps with the 633 nm laser. The optimization of resonance excitation leads to pronounced tunability of the orbital-phonon hybridized states. Our results generalize the orbital-phonon coupling model of effective phonon magnetic moments to material systems beyond the paramagnets and magnets.
We investigate the one-dimensional (1D) inviscid compressible flow equations for an ideal gas through the lens of optimization techniques. It is the case that, to our knowledge, optimization analysis applied to the so-called “linear velocity” solutions of the Euler compressible flow equations has not been previously conducted. Through both gradient-based and variational techniques, new variants of well-studied flow scenarios, i.e., self-similar, 1D, linear velocity solution class to idealized inviscid compressible flow equations, are determined, as encoded in both the kinematic and thermodynamic properties of this self-similar solution class. With the kinematics of the said solutions being driven by a self-similar “scale radius” and the thermodynamics being driven separately through the appearance of an arbitrary function, a myriad of new solution classes is possible. Acting as a guide to more realistic physical circumstances as well as discovery, it is the hope that the presented cases serve as the framework for future investigations into the intersection of self-similarity and optimization techniques. Fields of study that may find this work to be of interest include aerodynamic design, flow control, inertial confinement fusion, physics-informed neural networks, and other related areas of interest.
Electrostatic Levitation (ESL) is an advanced technique for containerless processing of metals, ceramics, and semiconductors. Because no container is required, there is no contamination from reaction with a crucible, allowing processing of high temperature, highly reactive melts. The high vacuum processing environment further reduces possible contamination of the samples. Finally, there is no container to provide heterogeneous nucleation sites, so the undercooled range is also accessible for many materials. For these reasons, ESL provides a unique environment for measuring thermophysical properties of liquid materials. The properties that can be measured in ESL include density, surface tension, viscosity, electrical and thermal conductivity, specific heat, phase diagram, TTT- and CCT- curves, and other thermodynamic properties. In this paper, we present data on surface tension and viscosity, measured by the oscillating drop technique, and density, measured by an automated photographic technique, measured in the ESL at NASA Marshall Space Flight Center.
Computational fluid dynamics analyses and experimental data are presented for the Mach 6 facility nozzle used in the Arc-Heated Scramjet Test Facility for the Enhanced Injection and Mixing Project (EIMP). This project, conducted at the NASA Langley Research Center, aims to investigate supersonic combustion ramjet (scramjet) fuel injection and mixing physics relevant to flight Mach numbers greater than 8. The EIMP experiments use a two-dimensional Mach 6 facility nozzle to provide the high-speed air simulating the combustor entrance flow of a scramjet engine. Of interest are the physical extent and the thermodynamic properties of the core flow at the nozzle exit plane. The detailed characterization of this flow is obtained from three-dimensional, viscous, Reynolds-averaged simulations. Thermodynamic nonequilibrium effects are also investigated. The simulations are compared with the available experimental data, which includes wall static pressures as well as in-stream static pressure, pitot pressure and total temperature obtained via in-stream probes positioned just downstream of the nozzle exit plane.
Pump Design (PUMPDES) is a computer program for designing a rotating pump for liquid hydrogen, liquid oxygen, liquid nitrogen, water, methane, or ethane. Using realistic properties of these fluids provided by another program called GASPAK, this code performs a station-by-station, mean-line analysis along the pump flow path, obtaining thermodynamic properties of the pumped fluid at each station and evaluating hydraulic losses along the flow path. The variables at each station are obtained under constraints that are consistent with the underlying physical principles. The code evaluates the performance of each stage and the overall pump. In addition, by judiciously choosing the givens and the unknowns, the code can perform a geometric inverse design function: that is, it can compute a pump geometry that yields a closest approximation of given design point. The code contains two major parts: one for an axial-rotor/inducer and one for a multistage centrifugal pump. The inducer and the centrifugal pump are functionally integrated. The code can be used in designing and/or evaluating the inducer/centrifugal-pump combination or the centrifugal pump alone. The code is written in standard Fortran 77.