Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermodynamics of computing”

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 559 records · Page 31

Coefficients for calculating thermodynamic and transport properties of individual species

Libraries of thermodynamic data and transport properties are given for individual species in the form of least-squares coefficients. Values of C(sup 0)(sub p)(T), H(sup 0)(T), and S(sup 0)(T) are available for 1130 solid, liquid, and gaseous species. Viscosity and thermal conductivity data are given for 155 gases. The original C(sup 0)(sub p)(T) values were fit to a fourth-order polynomial with integration constants for H(sup 0)(T) and S(sup 0)(T). For each species the integration constant for H(sup 0)(T) includes the heat of formation. Transport properties have a different functional form. The temperature range for most of the data is 300 to 5000 K, although some of the newer thermodynamic data have a range of 200 to 6000 K. Because the species are mainly possible products of reaction, the data are useful for chemical equilibrium and kinetics computer codes. Much of the data has been distributed for several years with the NASA Lewis equilibrium program CET89. The thermodynamic properties of the reference elements were updated along with about 175 species that involve the elements carbon, hydrogen, oxygen, and nitrogen. These sets of data will be distributed with the NASA Lewis personal computer program for calculating chemical equilibria, CETPC.

Mcbride, Bonnie J.↗

Real gas flow parameters for NASA Langley 22-inch Mach 20 helium tunnel

A computational procedure was developed which can be used to determine the flow properties in hypersonic helium wind tunnels in which real gas behavior is significant. In this procedure, a three-coefficient virial equation of state and the assumption of isentropic nozzle flow are employed to determine the tunnel reservoir, nozzle, throat, freestream, and post-normal shock conditions. This method was applied to a range of conditions which encompasses the operational capabilities of the LaRC 22-Inch Mach 20 Helium Tunnel. Results are presented graphically in the form of real gas correction factors which can be applied to perfect gas calculations. Important thermodynamic properties of helium are also plotted versus pressure and temperature. The computational scheme used to determine the real-helium flow parameters was incorporated into a FORTRAN code which is discussed.

Hollis, Brian R.↗

A thermodynamic approach to obtain materials properties for engineering applications

With the ever increases in the capabilities of computers for numerical computations, we are on the verge of using these tools to model manufacturing processes for improving the efficiency of these processes as well as the quality of the products. One such process is casting for the production of metals. However, in order to model metal casting processes in a meaningful way it is essential to have the basic properties of these materials in their molten state, solid state as well as in the mixed state of solid and liquid. Some of the properties needed may be considered as intrinsic such as the density, heat capacity or enthalpy of freezing of a pure metal, while others are not. For instance, the enthalpy of solidification of an alloy is not a defined thermodynamic quantity. Its value depends on the micro-segregation of the phases during the course of solidification. The objective of the present study is to present a thermodynamic approach to obtain some of the intrinsic properties and combining thermodynamics with kinetic models to estimate such quantities as the enthalpy of solidification of an alloy.

Chang, Y. Austin↗

The Navy/NASA Engine Program (NNEP89): A user's manual

An engine simulation computer code called NNEP89 was written to perform 1-D steady state thermodynamic analysis of turbine engine cycles. By using a very flexible method of input, a set of standard components are connected at execution time to simulate almost any turbine engine configuration that the user could imagine. The code was used to simulate a wide range of engine cycles from turboshafts and turboprops to air turborockets and supersonic cruise variable cycle engines. Off design performance is calculated through the use of component performance maps. A chemical equilibrium model is incorporated to adequately predict chemical dissociation as well as model virtually any fuel. NNEP89 is written in standard FORTRAN77 with clear structured programming and extensive internal documentation. The standard FORTRAN77 programming allows it to be installed onto most mainframe computers and workstations without modification. The NNEP89 code was derived from the Navy/NASA Engine program (NNEP). NNEP89 provides many improvements and enhancements to the original NNEP code and incorporates features which make it easier to use for the novice user. This is a comprehensive user's guide for the NNEP89 code.

Plencner, Robert M.↗

Numerical solution of Space Shuttle Orbiter flow field including real gas effects

The hypersonic, laminar flow around the Space Shuttle Orbiter has been computed for both an ideal gas (gamma = 1.2) and equilibrium air using a real-gas, parabolized Navier-Stokes code. This code employs a generalized coordinate transformation; hence, it places no restrictions on the orientation of the solution surfaces. The initial solution in the nose region was computed using a 3-D, real-gas, time-dependent Navier-Stokes code. The thermodynamic and transport properties of equilibrium air were obtained from either approximate curve fits or a table look-up procedure. Numerical results are presented for flight conditions corresponding to the STS-3 trajectory. The computed surface pressures and convective heating rates are compared with data from the STS-3 flight.

Prabhu, D. K.↗

A new theory of simple classical fluids

The paper presents a unified structural and thermodynamic theory of simple classical fluids in which the interactions between the particles can be represented by spherically symmetric pairwise potentials. Both the excess entropy and the gross features of the radial distribution function are determined mainly by excluded volume effects, which are in turn governed by a hard-core property intrinsic to any pair potential. The potential beyond this effective hard core is considered relatively weak and can be treated as a perturbation. It is also considered essential to sum all subclasses of diagrams to infinite order. The basic form of a diagrammatic scheme which allows both summation and the determination of the structure and excess entropy is presented, and a statement of universality is derived. The statement has been confirmed within the accuracy of present-day computer simulations, and a possible procedure for calculating both the structure and thermodynamics of every physically conceivable pair potential is presented.

Rosenfeld, Y.↗

A systematic analysis of phase stability in refractory high entropy alloys utilizing linear and non-linear cluster expansion models

We report that obtaining surrogate models for various types of materials has long been a goal of computational material scientists and physicists, as they allow the calculation of thermodynamic quantities and properties of interest far more efficiently than DFT. Many surrogate models, including interatomic potentials generated by machine learning algorthims and cluster expansions, rely on incorporating many-body interactions such as three-body or four-body interactions. When utilizing machine learning algorithms for interatomic potentials, it is common to start with two-body interactions and systematically include three- and four-body clusters to improve the fit, stopping when including the additional clusters does not measurably change the predictive power of the model. In cluster expansions for FCC and related lattices, it is often enough to include mostly two-body interactions and a handful of three-body interactions. However, it remains to be seen whether those will be enough in the case of lattices with a lower packing fraction, such as BCC.

36 MATERIALS SCIENCE↗

Adaptive Conformer Sampling for Property Prediction Using the Conductor-like Screening Model for Real Solvents

The valorization of lignocellulose-derived bioproducts requires effective separation from excessive water. Liquid–liquid extraction is a promising low-energy separation technology, but effective extraction requires solvent selection based on the thermodynamic properties of the bioproduct and solvent components. We propose a computational framework for predicting such properties by developing an adaptive conformer selection approach for use with COSMO-RS (conductor-like screening model for real solvents) calculations. In this framework, molecular dynamics simulations are used to generate many molecular structures (conformers) at representative temperatures in varying solvent environments. Conformers are then clustered based on structural metrics in a low-dimensional space and selected using a mixed-integer quadratic programming problem to iteratively insert a sampled conformer. At each iteration, we determine bioproduct properties using COSMO-RS. Here, we demonstrate the capability of the proposed framework on representative bioproducts to show convergence of the adaptive sampling toward experimentally measured properties with fewer calculations than required by random conformer sampling, enabling the improved screening of solvent systems for liquid-phase separation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

X-ray diffraction reveals two structural transitions in szomolnokite

Hydrated sulfates have been identified and studied in a wide variety of environments on Earth, Mars, and the icy satellites of the solar system. The subsurface presence of hydrous sulfur-bearing phases to any extent necessitates a better understanding of their thermodynamic and elastic properties at pressure. End-member experimental and computational data are lacking and are needed to accurately model hydrous, sulfur-bearing planetary interiors. In this work, high-pressure X-ray diffraction (XRD) and synchrotron Fourier-transform infrared (FTIR) measurements were conducted on szomolnokite (FeSO 4 ·H 2 O) up to ~83 and 24 GPa, respectively. This study finds a monoclinic-triclinic (C2/c to P1¯) structural phase transition occurring in szomolnokite between 5.0(1) and 6.6(1) GPa and a previously unknown triclinic-monoclinic (P1¯ to P2 1 ) structural transition occurring between 12.7(3) and 16.8(3) GPa. The high-pressure transition was identified by the appearance of distinct reflections in the XRD patterns that cannot be attributed to a second phase related to the dissociation of the P1¯ phase, and it is further characterized by increased H 2 O bonding within the structure. We fit third-order Birch-Murnaghan equations of state for each of the three phases identified in our data and refit published data to compare the elastic parameters of szomolnokite, kieserite (MgSO 4 ·H 2 O), and blödite (Na 2 Mg(SO 4 ) 2 ·4H 2 O). At ambient pressure, szomolnokite is less compressible than blödite and more than kieserite, but by 7 GPa both szomolnokite and kieserite have approximately the same bulk modulus, while blödite’s remains lower than both phases up to 20 GPa. Furthermore, these results indicate the stability of szomolnokite’s high-pressure monoclinic phase and the retention of water within the structure up to pressures found in planetary deep interiors.

58 GEOSCIENCES↗

Computer program for high pressure real gas effects

Computer program obtains the real-gas isentropic flow functions and thermodynamic properties of gases for which the equation of state is known. The program uses FORTRAN 4 subroutines which were designed for calculations of nitrogen and helium. These subroutines are easily modified for calculations of other gases.

Johnson, R. C.↗

Radiative transfer effects on reflected shock waves. I - Transparent gas.

Analytical and numerical calculation of the effects of radiative cooling on the flow field and heat transfer behind a reflected shock wave. The analytical solutions are developed by an expansion procedure about the Newtonian limit - i.e., the flow variables are expanded in the small parameter epsilon representing the initial density ratio across the shock front. Solutions carried out through the zeroth- and first-orders in epsilon show an accuracy to the order of epsilon squared - namely, for the specific conditions considered the analytically calculated enthalpy changes agree to within 2% with the numerically computed changes. The numerical procedure consists of a finite-difference calculation with pressure gradients neglected. The analytical calculations show that pressure-gradient effects may be neglected, except at very long times when the shocked gas has cooled by a large amount. The present calculations show much greater shock-wave attenuation than does a previous numerical computation. This difference is attributed to a better representation of the thermodynamic properties in the present calculations.

Su, F. Y.↗

The ratio of specific heats for postshock plasmas of a detached bow shock - An MHD model

The empirical relationship between the standoff distance of a detached bow shock (generated by the flow of a supersonic gas past an impenetrable obstacle), the size of the obstacle, the Mach number of the gas, and the ratio of specific heats has been generalized to include the magnetic field. The value of the ratio of specific heats (gamma-prime) in the postshock plasma has been calculated in terms of the preshock Alfvenic and sonic Mach numbers and orientation of the magnetic field. The empirical relationship is further generalized by taking into consideration the normal momentum and energy flux due to waves and/or turbulence and/or heat flow in association with high Mach number shocks. The computed value of gamma prime is substantially modified in comparison with that given by the MHD or the gas dynamic model. For this generalized model the computed gamma prime can be considered to be a more precise thermodynamic quantity, since the macroscopic parameters of the plasma have been separated out. Application of this empirical relationship to the earth's bow shock has been given.

Chao, J. K.↗

Solution influence on biomolecular equilibria - Nucleic acid base associations

Various attempts to construct an understanding of the influence of solution environment on biomolecular equilibria at the molecular level using computer simulation are discussed. First, the application of the formal statistical thermodynamic program for investigating biomolecular equilibria in solution is presented, addressing modeling and conceptual simplications such as perturbative methods, long-range interaction approximations, surface thermodynamics, and hydration shell. Then, Monte Carlo calculations on the associations of nucleic acid bases in both polar and nonpolar solvents such as water and carbon tetrachloride are carried out. The solvent contribution to the enthalpy of base association is positive (destabilizing) in both polar and nonpolar solvents while negative enthalpies for stacked complexes are obtained only when the solute-solute in vacuo energy is added to the total energy. The release upon association of solvent molecules from the first hydration layer around a solute to the bulk is accompanied by an increase in solute-solvent energy and decrease in solvent-solvent energy. The techniques presented are expectd to displace less molecular and more heuristic modeling of biomolecular equilibria in solution.

Pohorille, A.↗

A survey of spacecraft thermal design solutions

A review of activities at the NASA/Marshall Space Flight Center in the heat transfer and thermodynamics disciplines as well as attendant fluid mechanics, transport phenomena, and computer science applications is presented. Attention is focused on recent activities including the Hubble Space Telescope, and large space instruments, particularly telescope thermal control systems such as those flown aboard Spacelab 2 and the Astro missions. Emphasis is placed on defining the thermal control features, unique design schemes, and performance of selected programs. Results obtained both by ground testing and analytical means, as well as flight and postflight data are presented.

Humphries, R.↗

Characterizing Defect Dynamics in Silicon Carbide Using Symmetry-Adapted Collective Variables and Machine Learning Interatomic Potentials

Silicon carbide (SiC) divacancies are attractive candidates for spin-defect qubits possessing long coherence times and optical addressability. The high activation barriers associated with SiC defect formation and motion pose challenges for their study by first-principles molecular dynamics. In this work, we develop and deploy machine learning interatomic potentials (MLIPs) to accelerate defect dynamics simulations while retaining ab initio accuracy. We employ an active learning strategy comprising symmetry-adapted collective variable discovery and enhanced sampling to compile configurationally diverse training data, calculation of energies and forces using density functional theory (DFT), and training of an E(3)-equivariant MLIP based on the Allegro model. Here, the trained MLIP reproduces DFT-level accuracy in defect transition activation free energy barriers, enables the efficient and stable simulation of multidefect 216-atom supercells, and permits an analysis of the temperature dependence of defect thermodynamic stability and formation/annihilation kinetics to propose an optimal annealing temperature to maximally stabilize VV divacancies.

Computer simulations↗

FY21 Status Report on the Molten Salt Thermal Properties Database (MSTDB) Development

This report describes the thermodynamic assessments for the ZrF 4 -BeF 2 , AlCl 3 -NaCl, AlCl3-KCl and the PuCl 3 -CsCl pseudo-binary subsystems. Calculation results are compared to the inputs used to optimize the adjustable model parameters. The data came from experimental studies reported in the open literature, from the MSR research community, and/or computational results generated within the NEAMS program and from collaborators. The models from the thermodynamic assessment of these salt systems are integrated into the MSTDB-TC (Molten Salt Thermal Properties Database-Thermochemical). A brief description of the companion MSTDB-TP (Molten Salt Thermal Properties Database-Thermophysical) is given along with an overview of the approach for predicting both thermodynamic and thermophysical property behavior of multicomponent systems.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

TROPHY: A Topologically Robust Physics-Informed Tracking Framework for Tropical Cyclones

Tropical cyclones (TCs) are among the most destructive weather systems. Realistically and efficiently detecting and tracking TCs are critical for assessing their impacts and risks. In particular, the eye is a signature feature of a mature TC. Therefore, knowing the eyes’ locations and movements is crucial for both operational weather forecasts and climate risk assessments. Recently, a multilevel robustness framework has been introduced to study the critical points of time-varying vector fields. The framework quantifies the robustness (i.e., structural stability) of critical points across varying neighborhoods. By relating the multilevel robustness with critical point tracking, the framework has demonstrated its potential in cyclone tracking. An advantage is that it identifies cyclonic features using only 2D wind vector fields, which is encouraging as most tracking algorithms require multiple dynamic and thermodynamic variables at different altitudes. A disadvantage is that the framework does not scale well computationally for datasets containing a large number of cyclones. Herein this paper introduces a topologically robust physics-informed tracking framework (TROPHY) for TC tracking. The main idea is to integrate physical knowledge of TC to drastically improve the computational efficiency of multilevel robustness framework for large-scale climate datasets. First, during preprocessing, we propose a physics-informed feature selection strategy to filter 90% of critical points that are short-lived and have low stability, thus preserving good candidates for TC tracking. Second, during in-processing, we impose constraints during the multilevel robustness computation to focus only on physics-informed neighborhoods of TCs. We apply TROPHY to 30 years of 2D wind fields from reanalysis data in ERA5 and generate a number of TC tracks. In comparison with the observed tracks, we demonstrate that TROPHY can capture TC characteristics (e.g., frequency, intensity, duration, latitudes with maximum intensity, and genesis) that are comparable to and sometimes even better than a well-validated TC tracking algorithm that requires multiple dynamic and thermodynamic scalar fields.

97 MATHEMATICS AND COMPUTING↗

DFTTK: Density Functional Theory ToolKit for high-throughput lattice dynamics calculations

In this work, we present a software package in Python for high-throughput first-principles calculations of thermodynamic properties at finite temperatures, which we refer to as DFTTK (Density Functional Theory ToolKit). DFTTK is based on the atomate package and integrates our experiences in the last decades on the development of theoretical methods and computational softwares. It includes task submissions on all major operating systems and task executions on high-performance computing environments. Furthermore, the distribution of the DFTTK package comes with examples of calculations of phonon density of states, heat capacity, entropy, enthalpy, and free energy under the quasi-harmonic phonon scheme for the stoichiometric phases of Al, Ni, Al 3 Ni, AlNi, AlNi 3 , Al 3 Ni 4 , and Al 3 Ni 5 , and the fcc solution phases treated using the special quasirandom structures at the compositions of Al 3 Ni, AlNi, and AlNi 3 .

97 MATHEMATICS AND COMPUTING↗