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 1,621 records · Page 90
All sourcesNASA NTRSDOE OSTI
NASA NTRS2006

Thermodynamic Analysis of Dual-Mode Scramjet Engine Operation and Performance

Recent analytical advances in understanding the performance continuum (the thermodynamic spectrum) for air-breathing engines based on fundamental second-law considerations have clarified scramjet and ramjet operation, performance, and characteristics. Second-law based analysis is extended specifically in this work to clarify and describe the performance characteristics for dual-mode scramjet operation in the mid-speed range of flight Mach 4 to 7. This is done by a fundamental investigation of the complex but predictable interplay between heat release and irreversibilities in such an engine; results demonstrate the flow and performance character of the dual mode regime and of dual mode transition behavior. Both analytical and computational (multi-dimensional CFD) studies of sample dual-mode flow-fields are performed in order to demonstrate the second-law capability and performance and operability issues. The impact of the dual-mode regime is found to be characterized by decreasing overall irreversibility with increasing heat release, within the operability limits of the system.

Riggins, David↗
DOE OSTI2021

First principles calculations in support of Pu aging: calculating the effects of lattice imperfections on thermodynamics

Plutonium (Pu) has, in theory, well defined crystal structures: its atoms are arranged in regular spatial patterns. But Pu is radioactive, and as its nuclei decay those regular spatial patterns are interrupted. The interruptions are lattice imperfections, which are known to affect how materials respond to their environment. To adequately model Pu, we need to know which lattice imperfections are present, how they interact with each other, and how they affect the material’s response to its environment. The work presented here aims to use density functional theory (DFT) calculations to begin to answer the latter, in particular, how individual lattice imperfections affect measurable effects including thermal expansion (the change in volume in response to a change in temperature), heat capacity (the amount of thermal energy needed to change a material’s temperature), and elastic moduli (a material’s resistance to applied stresses). Pu poses many computational challenges. The $\textit{f}$ electrons require special attention. Of all the elements, Pu has the largest number of electrons that must be included in the calculations. Thermal effects require calculating the phonons (the lattice vibrations), which for systems containing lattice imperfections demand large, complex computational cells - but computational resources limit the size and complexity. With careful restructuring of how the calculations are performed, all these challenges have been met to enable calculations that provide insight into how lattice imperfections affect Pu’s response to its environment. Reported here are the computational challenges and the advances developed to meet them, along with the first results showing the strong effect that one prototype lattice imperfection (an interstitial Pu atom in a delta-phase Pu lattice) has on Pu’s response to its environment.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Neural Network for Principle of Least Action

The principle of least action is the cornerstone of classical mechanics, theory of relativity, quantum mechanics, and thermodynamics. Here, we describe how a neural network (NN) learns to find the trajectory for a Lennard-Jones (LJ) system that maintains balance in minimizing the Onsager–Machlup (OM) action and maintaining the energy conservation. The phase-space trajectory thus calculated is in excellent agreement with the corresponding results from the “ground-truth” molecular dynamics (MD) simulation. Furthermore, we show that the NN can easily find structural transformation pathways for LJ clusters, for example, the basin-hopping transformation of an LJ 38 from an incomplete Mackay icosahedron to a truncated face-centered cubic octahedron. Unlike MD, the NN computes atomic trajectories over the entire temporal domain in one fell swoop, and the NN time step is a factor of 20 larger than the MD time step. The NN approach to OM action is quite general and can be adapted to model morphometrics in a variety of applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2004

Aqueous Chemical Modeling of Sedimentation on Early Mars with Application to Surface-Atmosphere Evolution

This project was to investigate models for aqueous sedimentation on early Mars from fluid evaporation. Results focused on three specific areas: (1) First, a fluid evaporation model incorporating iron minerals was developed to compute the evaporation of a likely solution on early Mars derived from the weathering of mafic rock. (2) Second, the fluid evaporation model was applied to salts within Martian meteorites, specifically salts in the nakhlites and ALH84001. Evaporation models were found to be consistent with the mineralogy of salt assemblages-anhydrite, gypsum, Fe-Mg-Ca carbonates, halite, clays-- and the concentric chemical fractionation of Ca-to Mg-rich carbonate rosettes in ALH84001. We made progress in further developing our models of fluid concentration by contributing to updating the FREZCHEM model. (3) Third, theoretical investigation was done to determine the thermodynamics and kinetics involved in the formation of gray, crystalline hematite. This mineral, of probable ancient aqueous origin, has been observed in several areas on the surface of Mars by the Thermal Emission Spectrometer on Mars Global Surveyor. The "Opportunity" Mars Exploration Rover has also detected gray hematite at its landing site in Meridiani Planum. We investigated how gray hematite can be formed via atmospheric oxidation, aqueous precipitation and subsequent diagenesis, or hydrothermal processes. We also studied the geomorphology of the Aram Chaos hematite region using Mars Orbiter Camera (MOC) images.

Catling, David C.↗
NASA NTRS2009

Uncertain Representations of Sub-Grid Pollutant Transport in Chemistry-Transport Models and Impacts on Long-Range Transport and Global Composition

Sub-grid transport, by convection and turbulence, is known to play an important role in lofting pollutants from their source regions. Consequently, the long-range transport and climatology of simulated atmospheric composition are impacted. This study uses the Goddard Earth Observing System, Version 5 (GEOS-5) atmospheric model to study pollutant transport. The baseline model uses a Relaxed Arakawa-Schubert (RAS) scheme that represents convection through a sequence of linearly entraining cloud plumes characterized by unique detrainment levels. Thermodynamics, moisture and trace gases are transported in the same manner. Various approximate forms of trace-gas transport are implemented, in which the box-averaged cloud mass fluxes from RAS are used with different numerical approaches. Substantial impacts on forward-model simulations of CO (using a linearized chemistry) are evident. In particular, some aspects of simulations using a diffusive form of sub-grid transport bear more resemblance to space-biased CO observations than do the baseline simulations with RAS transport. Implications for transport in the real atmosphere will be discussed. Another issue of importance is that many adjoint/inversion computations use simplified representations of sub-grid transport that may be inconsistent with the forward models: implications will be discussed. Finally, simulations using a complex chemistry model in GEOS-5 (in place of the linearized CO model) are underway: noteworthy results from this simulation will be mentioned.

Pawson, Steven↗
NASA NTRS2015

Global Ocean Evaporation: How Well Can We Estimate Interannual to Decadal Variability?

Evaporation from the world's oceans constitutes the largest component of the global water balance. It is important not only as the ultimate source of moisture that is tied to the radiative processes determining Earth's energy balance but also to freshwater availability over land, governing habitability of the planet. Here we focus on variability of ocean evaporation on scales from interannual to decadal by appealing to three sources of data: the new MERRA-2 (Modern-Era Retrospective analysis for Research and Applications -2); climate models run with historical sea-surface temperatures, ice and atmospheric constituents (so-called AMIP experiments); and state-of-the-art satellite retrievals from the Seaflux and HOAPS (Hamburg Ocean-Atmosphere Parameters and Fluxes from Satellite) projects. Each of these sources has distinct advantages as well as drawbacks. MERRA-2, like other reanalyses, synthesizes evaporation estimates consistent with observationally constrained physical and dynamical models-but data stream discontinuities are a major problem for interpreting multi-decadal records. The climate models used in data assimilation can also be run with lesser constraints such as with SSTs and sea-ice (i.e. AMIPs) or with additional, minimal observations of surface pressure and marine observations that have longer and less fragmentary observational records. We use the new ERA-20C reanalysis produced by ECMWF embodying the latter methodology. Still, the model physics biases in climate models and the lack of a predicted surface energy balance are of concern. Satellite retrievals and comparisons to ship-based measurements offer the most observationally-based estimates, but sensor inter-calibration, algorithm retrieval assumptions, and short records are dominant issues. Our strategy depends on maximizing the advantages of these combined records. The primary diagnostic tool used here is an analysis of bulk aerodynamic computations produced by these sources and uses a first-order Taylor series analysis of wind speed, SST, near-surface stability and relative humidity variations around climatology to gauge the importance of these components. We find that the MERRA-2 evaporation record is strongly influenced by the availability of wind speed and humidity from passive microwave imagers beginning in the late 1980s as well as by the SST record. The trend over the period 1980 to present is nearly 10%. AMIP or the ERA-20C trends are much smaller. We find that ENSO-related signals involving both wind speed and thermodynamic variability remain the primary signal in the latter and are confirmed by satellite retrievals. We present uncertainty estimates based on the various data sources and discuss the implications for GEWEX water and energy budget science challenges.

Robertson, Franklin R.↗
DOE OSTI2020

A perspective on challenges and opportunities in developing high entropy-ultra high temperature ceramics

The emergence of a new class of high entropy ultra-high temperature ceramics (HE-UHTCs) is expected to possess exceptionally superior mechanical, oxidation, and erosion resistance as compared to conventional UHTCs. The improvement in the properties of HE-UHTCs is attributed to the high configurational entropy of the multi-component system stabilized into a single-phase, but with little discussion on entropy contributions. This perspective article presents a viewpoint on the critical scientific question regarding the actual role entropy plays in stabilizing and enhancing the thermo-mechanical properties of a multi-component system. The high throughput first-principle calculations, which can provide the entropy formation ability (EFA) of thermodynamically stable HE-UHTCs, can prove to be valuable toolset to this puzzle. HE-UHTCs show enhancement in the mechanical properties much higher than the predicted rule of the mixture, which cannot traditionally be explained by a mere solid solution formation effect. Additionally, isothermal oxidation studies of HE-UHTCs limited to furnace testing state remarkable oxidation resistance which is attributed to the passivity of complex, layered structure, and composition of the oxide scale. However, whether all the constituents in the HE-UHTC system react simultaneously or preferential oxidation takes place needs to be addressed. To address all these issues, an integrated computational and experimental approach has been postulated towards designing a thermal protection system (TPS) for re-entry applications. Succeeding the approach will provide new opportunities in the compositional space of HE-UHTCs with tailored properties.

36 MATERIALS SCIENCE↗
NASA NTRS1971

A computer program for two-dimensional and axisymmetric nonreacting perfect gas and equilibrium chemically reacting laminar, transitional and-or turbulent boundary layer flows

A computer program is described in detail for laminar, transitional, and/or turbulent boundary-layer flows of non-reacting (perfect gas) and reacting gas mixtures in chemical equilibrium. An implicit finite difference scheme was developed for both two dimensional and axisymmetric flows over bodies, and in rocket nozzles and hypervelocity wind tunnel nozzles. The program, program subroutines, variables, and input and output data are described. Also included is the output from a sample calculation of fully developed turbulent, perfect gas flow over a flat plate. Input data coding forms and a FORTRAN source listing of the program are included. A method is discussed for obtaining thermodynamic and transport property data which are required to perform boundary-layer calculations for reacting gases in chemical equilibrium.

Miner, E. W.↗
NASA NTRS1996

Computer Program for Calculation of Complex Chemical Equilibrium Compositions and Applications II. Users Manual and Program Description: Users Manual and Program Description - 2

This users manual is the second part of a two-part report describing the NASA Lewis CEA (Chemical Equilibrium with Applications) program. The program obtains chemical equilibrium compositions of complex mixtures with applications to several types of problems. The topics presented in this manual are: (1) details for preparing input data sets; (2) a description of output tables for various types of problems; (3) the overall modular organization of the program with information on how to make modifications; (4) a description of the function of each subroutine; (5) error messages and their significance; and (6) a number of examples that illustrate various types of problems handled by CEA and that cover many of the options available in both input and output. Seven appendixes give information on the thermodynamic and thermal transport data used in CEA; some information on common variables used in or generated by the equilibrium module; and output tables for 14 example problems. The CEA program was written in ANSI standard FORTRAN 77. CEA should work on any system with sufficient storage. There are about 6300 lines in the source code, which uses about 225 kilobytes of memory. The compiled program takes about 975 kilobytes.

McBride, Bonnie J.↗
DOE OSTI2022

Low-Pressure Electrochemical Synthesis of Complex High-Pressure Superconducting Superhydrides

There is great current interest in multicomponent superhydrides due to their unique quantum properties under pressure. A remarkable example is the ternary superhydride Li 2 MgH 16 computationally identified to have an unprecedented high superconducting critical temperature T c of ~470 K at 250 GPa. However, the very high synthesis pressures required remains a significant hurdle for detailed study and potential applications. In this Letter, we evaluate the feasibility of synthesizing ternary Li-Mg superhydrides by the recently proposed pressure-potential (P 2 ) method that uniquely combines electrochemistry and applied pressure to control synthesis and stability. Furthermore, the results indicate that it is possible to synthesize Li-Mg superhydrides at modest pressures by applying suitable electrode potentials. Using pressure alone, no Li-Mg ternary hydrides are predicted to be thermodynamically stable, but in the presence of electrode potentials, both Li 2 MgH 16 and Li 4 MgH 24 can be stabilized at modest pressures. Three polymorphs are predicted as ground states of Li 2 MgH 16 below 300 GPa, with transitions at 33 and 160 GPa. The highest pressure phase is superconducting, while the two at lower pressures are not. Our findings point out the potentially important role of the P 2 method in controlling phase stability of complex multicomponent superhydrides.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1985

Statistical mechanics of light elements at high pressure. VII - A perturbative free energy for arbitrary mixtures of H and He

A model free energy is presented which accurately represents results from 45 high-precision Monte Carlo calculations of the thermodynamics of hydrogen-helium mixtures at pressures of astrophysical and planetophysical interest. The free energy is calculated using free-electron perturbation theory (dielectric function theory), and is an extension of the expression given in an earlier paper in this series. However, it fits the Monte Carlo results more accurately, and is valid for the full range of compositions from pure hydrogen to pure helium. Using the new free energy, the phase diagram of mixtures of liquid metallic hydrogen and helium is calculated and compared with earlier results. Sample results for mixing volumes are also presented, and the new free energy expression is used to compute a theoretical Jovian adiabat and compare the adiabat with results from three-dimensional Thomas-Fermi-Dirac theory. The present theory gives slightly higher densities at pressures of about 10 megabars.

Hubbard, W. B.↗
DOE OSTI2022

Flexocaloric effect in superelastic materials

We present a combined theoretical–experimental study of flexocaloric effects in superelastic materials exhibiting structural transitions. We study a Ginzburg–Landau model combined with a vibrational model for a beam near a ferroelastic transition loaded with a three-point bending setup. We also perform experiments on a Cu–Al–Ni single crystal undergoing a martensitic transition. We measure bent beam profiles, vertical force vs vertical deflection during a slow isothermal process, time evolution of the bending and unbending amplitudes, and the evolution of temperature profiles. We also compute the evolution of heat source and heat sink profiles. Finally, we study the location of acoustic emission events during the bending/unbending experiment. Our observations are consistent with the model predictions and allow us to identify the main physical parameters relevant for flexocaloric applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Melting curves of ice polymorphs in the vicinity of the liquid–liquid critical point

The possible existence of a liquid–liquid critical point in deeply supercooled water has been a subject of debate due to the challenges associated with providing definitive experimental evidence. The pioneering work by Mishima and Stanley [Nature 392, 164–168 (1998)] sought to shed light on this problem by studying the melting curves of different ice polymorphs and their metastable continuation in the vicinity of the expected liquid–liquid transition and its associated critical point. Based on the continuous or discontinuous changes in the slope of the melting curves, Mishima [Phys. Rev. Lett. 85, 334 (2000)] suggested that the liquid–liquid critical point lies between the melting curves of ice III and ice V. Herein we explore this conjecture using molecular dynamics simulations with a machine learning model based on ab initio quantum-mechanical calculations. We study the melting curves of ices III, IV, V, VI, and XIII and find that all of them are supercritical and do not intersect the liquid–liquid transition locus. We also find a pronounced, yet continuous, change in the slope of the melting lines upon crossing of the liquid locus of maximum compressibility. Finally, we analyze the literature in light of our findings and conclude that the scenario in which the melting curves are supercritical is favored by the most recent computational and experimental evidence. Although the preponderance of evidence is consistent with the existence of a second critical point in water, the behavior of ice polymorph melting lines does not provide strong evidence in support of this viewpoint, according to our calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

A combined experimental and machine learning exploration of Ti 2-x Zr x MnCrFeNi high entropy Laves hydrides

A series of high entropy AB 2 -type Ti 2-x Zr x MnCrFeNi alloys (x = 0.6, 0.7, 0.8, 0.9, 1.0, 1.1 and 1.2) were synthesized to investigate their potential for hydrogen storage and chemical compression. The influence of the Ti/Zr ratio was explored in terms of structural, microstructural and thermodynamic properties. The storage capacity together with the reaction enthalpy and entropy changes of the synthesized high entropy alloys were compared to predictions from Machine Learning (ML) to investigate changes in these properties across the explored composition space. The results revealed that a decreasing Zr content consistently lowered the hydride formation enthalpy and increased the plateau pressure from 8 to >90 bar H 2 at 25 °C, in good agreement with ML predictions. Selected compositions (x = 1.0 and 1.2) demonstrated reversible hydrogen storage capability over 150 cycles, with capacities of 1.34–1.40 wt % H 2 and remarkable reaction kinetics (<4 min) at ambient temperature. These experimental and computational findings highlight the potential of this Laves-HEA system as tuneable, stable, and cost-effective materials suitable for long-term operations in stationary hydrogen storage and compression applications.

36 MATERIALS SCIENCE↗
NASA NTRS2001

Simplified Analysis of Pulse Detonation Rocket Engine B1owdown Gasdynamics and Performance

Pulsed detonation rocket engines (PDREs) have generated considerable research interest in recent years as a chemical propulsion system potentially offering improved performance and reduced complexity compared to conventional rocket engines. The detonative mode of combustion employed by these devices offers a thermodynamic advantage over the constant-pressure deflagrative combustion mode used in conventional rocket engines and gas turbines. However, while this theoretical advantage has spurred a great deal of interest in building PDRE devices, the unsteady blowdown process intrinsic to the PDRE has made realistic estimates of the actual propulsive performance problematic. The recent review article by Kailasanath highlights some of the difficulties in comparing the available experimental measurements with numerical models. The goal of this paper is to improve understanding of PDRE blowdown gasdynamics and performance issues through use of a simplified model that captures the essential features of the unsteady blowdown process, and yet remains computationally inexpensive. The PDRE system studied here is highly idealized, consisting of a constant-area detonation tube with one end closed and the other end open to the environment. The tube is prefilled with a gaseous propellant mixture with no initial velocity or outflow to the environment. The detonation is initiated instantaneously at the closed end of the device. Chapman-Jouguet (C-J) post-detonation gas conditions are calculated using the CET89 version of the NASA thermochemical code. The I-D, unsteady method of characteristics is used to calculate the flowfield following the detonation front. See the compressible flow texts by Thompson and Zucrow and Hoffman for details of this method. Parametric studies of the effect of mixture stoichiometry, fill temperature, and blowdown pressure ratio on performance are reported. A comparison of the performance of an idealized straight-tube PDRE with a conventional steady-state rocket engine is provided. The effect of constant-gamma and equilibrium chemistry assumptions is also examined. Additionally, in order to form an assessment of the accuracy of the model, the flowfield time history is compared to experimental data from Stanford University.

Morris, Christopher I.↗
DOE OSTI2023

REYNO: A Reactive Hydrodynamics Modeling Suite

REYNO is a Python library that is primarily designed as a platform for the agile development and testing of novel reactive flow models. Equations of state and multi-step rate laws can be implemented using the provided class hierarchy with different thermodynamic closure conditions. The library includes a one-dimensional Lagrangian hydrodynamic solver to simulate multi-material reactive or inert problems in planar, cylindrical, or spherical coordinates. Custom time-dependent boundary conditions are also available for problems such as ramp compression. This document describes the different numerical methods and algorithms that are included in REYNO.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2003

Microelectromechanical System (MEMS) Device Being Developed for Active Cooling and Temperature Control

High-capacity cooling options remain limited for many small-scale applications such as microelectronic components, miniature sensors, and microsystems. A microelectromechanical system (MEMS) using a Stirling thermodynamic cycle to provide cooling or heating directly to a thermally loaded surface is being developed at the NASA Glenn Research Center to meet this need. The device can be used strictly in the cooling mode or can be switched between cooling and heating modes in milliseconds for precise temperature control. Fabrication and assembly employ techniques routinely used in the semiconductor processing industry. Benefits of the MEMS cooler include scalability to fractions of a millimeter, modularity for increased capacity and staging to low temperatures, simple interfaces, limited failure modes, and minimal induced vibration. The MEMS cooler has potential applications across a broad range of industries such as the biomedical, computer, automotive, and aerospace industries. The basic capabilities it provides can be categorized into four key areas: 1) Extended environmental temperature range in harsh environments; 2) Lower operating temperatures for electronics and other components; 3) Precision spatial and temporal thermal control for temperature-sensitive devices; and 4) The enabling of microsystem devices that require active cooling and/or temperature control. The rapidly expanding capabilities of semiconductor processing in general, and microsystems packaging in particular, present a new opportunity to extend Stirling-cycle cooling to the MEMS domain. The comparatively high capacity and efficiency possible with a MEMS Stirling cooler provides a level of active cooling that is impossible at the microscale with current state-of-the-art techniques. The MEMS cooler technology builds on decades of research at Glenn on Stirling-cycle machines, and capitalizes on Glenn s emerging microsystems capabilities.

Beach, Duane E.↗
DOE OSTI2025

Quantum Ornstein-Zernike theory for two-temperature two-component plasmas

Laboratory plasma production almost always preferentially heats either the ions or electrons, leading to a two-temperature state. In this state, density functional theory molecular dynamic simulation is the state of the art for modeling bulk material properties. We construct a statistical mechanics model for the two temperature limit that is theoretically consistent with the molecular dynamics method. We proceed to derive the electron-ion multi-temperature quantum Ornstein-Zernike equations for the first time. This allows the construction of a two-temperature two-component plasma model using the average atom from which we can compute bulk material properties at a fraction of the computation time of the two-temperature density functional theory simulation. The accuracy of the model is benchmarked against ion pair correlation and self-diffusion results from ab initio simulation. Here, we proceed to compute the viscosity and ion thermal conductivity as a function of both ion and electron temperature.

Ab initio molecular dynamics↗