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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,369 records · Page 76

All sources NASA NTRS DOE OSTI

NASA NTRS1986

Acceleration effects in solid propellant rocket motors

The performance variations due to acceleration loads imposed on spinning solid propellant rocket motors are investigated. The four potentially most significant modes of acceleration-induced phenomena are identified from a study of the literature and modeled. The four modes are a mechanical mode which deals with deformations of the propellant and case: a thermodynamic mode which covers acceleration-induced combustion phenomena; a stress mode which covers the stressed propellant's effect on burn rate; and a gas dynamic mode which deals with changes in gas flow in the chamber and through the nozzle. Simplified models of each mode are developed or taken from the literature and are added to an internal ballistics evaluation computer program. The resulting analysis is the first to include all of the modes. In order to do this an original analysis of the mechanical and stress modes was necessary. However, the analysis shows that the stress mode is not important for the circular perforated grains studied. The other effects are shown to have a significant influence on solid rocket motor performance. The magnitude of the different mode effects are such that one may not be ignored over the others as has been done in the past. The results of the analysis are compared to published rocket motor data. The comparisons indicate an erosive burning effect that is a function of spin rate. A qualitative explanation of the erosive effect is presented.

Langhenry, M. T.↗

NASA NTRS1993

Analysis of the anisotropic viscoplastic-damage response of composite laminates - Continuum basis and computational algorithms

The mathematical structure underlying the rate equations of a recently-developed constitutive model for the coupled viscoplastic-damage response of anisotropic composites is critically examined. In this regard, a number of tensor projection operators have been identified, and their properties were exploited to enable the development of a general computational framework for their numerical implementation using the Euler fully-implicit integration method. In particular, this facilitated (i) the derivation of explicit expressions of the (consistent) material tangent stiffnesses that are valid for both three-dimensional as well as subspace (e.g. plane stress) formulations, (ii) the implications of the symmetry or unsymmetry properties of these tangent operators from a thermodynamic standpoint, and (iii) the development of an effective time-step control strategy to ensure accuracy and convergence of the solution. In addition, the special limiting case of inviscid elastoplasticity is treated. The results of several numerical simulations are given to demonstrate the effectiveness of the schemes developed.

Saleeb, A. F.↗

NASA NTRS1994

Nonequilibrium radiation and chemistry models for aerocapture vehicle flowfields

The primary accomplishments of the project were as follows: (1) From an overall standpoint, the primary accomplishment of this research was the development of a complete gasdynamic-radiatively coupled nonequilibrium viscous shock layer solution method for axisymmetric blunt bodies. This method can be used for rapid engineering modeling of nonequilibrium re-entry flowfields over a wide range of conditions. (2) Another significant accomplishment was the development of an air radiation model that included local thermodynamic nonequilibrium (LTNE) phenomena. (3) As part of this research, three electron-electronic energy models were developed. The first was a quasi-equilibrium electron (QEE) model which determined an effective free electron temperature and assumed that the electronic states were in equilibrium with the free electrons. The second was a quasi-equilibrium electron-electronic (QEEE) model which computed an effective electron-electronic temperature. The third model was a full electron-electronic (FEE) differential equation model which included convective, collisional, viscous, conductive, vibrational coupling, and chemical effects on electron-electronic energy. (4) Since vibration-dissociation coupling phenomena as well as vibrational thermal nonequilibrium phenomena are important in the nonequilibrium zone behind a shock front, a vibrational energy and vibration-dissociation coupling model was developed and included in the flowfield model. This model was a modified coupled vibrational dissociation vibrational (MCVDV) model and also included electron-vibrational coupling. (5) Another accomplishment of the project was the usage of the developed models to investigate radiative heating. (6) A multi-component diffusion model which properly models the multi-component nature of diffusion in complex gas mixtures such as air, was developed and incorporated into the blunt body model. (7) A model was developed to predict the magnitude and characteristics of the shock wave precursor ahead of vehicles entering the Earth's atmosphere. (8) Since considerable data exists for radiating nonequilibrium flow behind normal shock waves, a normal shock wave version of the blunt body code was developed. (9) By comparing predictions from the models and codes with available normal shock data and the flight data of Fire II, it is believed that the developed flowfield and nonequilibrium radiation models have been essentially validated for engineering applications.

Carlson, Leland A.↗

DOE OSTI2025

Mechanistic Insights into Molecular Copper Hydride Catalysis: the Kinetic Stability of CuH Monomers toward Aggregation is a Critical Parameter for Catalyst Performance

The activity of molecular copper hydride (CuH) complexes towards the selective insertion of unsaturated hydrocarbons under mild conditions has contributed significantly to versatile methodologies for upgrading these feedstocks. However, these catalysts are particularly susceptible to deleterious aggregation, leading to the depletion of active CuH species. Little is known about the mechanisms of CuH aggregation, how it influences overall catalyst performance, and how it can be controlled. We address these challenges with mechanistic studies on a model reaction of unactivated alkene hydroboration catalyzed by (IPr*CPh 3 )CuH (LCuH). Here, we report a comprehensive mechanistic investigation of this system, identifying an aggregation pathway that continuously depletes catalytically active LCuH to form inactive CuH clusters during turnover. Deactivation of LCuH is controlled primarily by the competition between the kinetics of the initial LCuH dimerization step and that of alkene insertion. We therefore propose that a more comprehensive understanding of CuH catalyst performance must account for the kinetics of the initial LCuH dimerization step, revising a previously explored thermodynamic understanding of CuH aggregation, where the concentration of active species is controlled by equilibria established between CuH dimers and monomers. With a series of (NHC)CuH congeners (NHC = N-heterocyclic carbene), we demonstrate that ostensibly minor structural modifications to the ligand peripheries can drastically affect the LCuH dimerization kinetics, while maintaining reactivity towards on–cycle alkene insertion. We employed a computational approach based on molecular dynamics simulations to provide an in-depth understanding of how specific structural ligand modifications can substantially increase the kinetic stability of monomeric CuH catalysts. Our combined experimental and computational studies suggest strategies for rational ligand design that can be broadly applied to molecular catalyst systems that are susceptible to deactivation via aggregation pathways.

Ryan, David E. [Pacific Northwest National Laborat↗

DOE OSTI2024

The Bell-Evans-Polanyi relation for hydrogen evolution reaction from first-principles

The versatile Bell-Evans-Polanyi (BEP) relation stipulates the kinetics of a reaction in terms of thermodynamics. Herein, we establish the BEP relation for the hydrogen evolution reaction (HER) from fundamental electrochemical principles leveraging the Butler-Volmer relation for a one-step, one-electron process and the transition state theory. Based on first-principles investigations of HER mechanisms on fourteen metal electrodes, we firmly justify the BEP relation solely using an easy-to compute hydrogen adsorption free energy and universal electrochemical constants.

42 ENGINEERING↗

DOE OSTI2023

Dimensionally reduced machine learning model for predicting single component octanol–water partition coefficients

Abstract MF-LOGP, a new method for determining a single component octanol–water partition coefficients ( $$LogP$$ LogP ) is presented which uses molecular formula as the only input. Octanol–water partition coefficients are useful in many applications, ranging from environmental fate and drug delivery. Currently, partition coefficients are either experimentally measured or predicted as a function of structural fragments, topological descriptors, or thermodynamic properties known or calculated from precise molecular structures. The MF-LOGP method presented here differs from classical methods as it does not require any structural information and uses molecular formula as the sole model input. MF-LOGP is therefore useful for situations in which the structure is unknown or where the use of a low dimensional, easily automatable, and computationally inexpensive calculations is required. MF-LOGP is a random forest algorithm that is trained and tested on 15,377 data points, using 10 features derived from the molecular formula to make $$LogP$$ LogP predictions. Using an independent validation set of 2713 data points, MF-LOGP was found to have an average $$RMSE$$ RMSE = 0.77 ± 0.007, $$MAE$$ MAE = 0.52 ± 0.003, and $${R}^{2}$$ R 2 = 0.83 ± 0.003. This performance fell within the spectrum of performances reported in the published literature for conventional higher dimensional models ( $$RMSE$$ RMSE = 0.42–1.54, $$MAE$$ MAE = 0.09–1.07, and $${R}^{2}$$ R 2 = 0.32–0.95). Compared with existing models, MF-LOGP requires a maximum of ten features and no structural information, thereby providing a practical and yet predictive tool. The development of MF-LOGP provides the groundwork for development of more physical prediction models leveraging big data analytical methods or complex multicomponent mixtures. Graphical Abstract

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1983

Transient flow analysis of the AEDC/HPDE MHD generator

A hybrid Lax-Wendroff/Method of Characteristics computer code has been developed for numerical simulation of flow transients associated with the operation of MHD generator facilities. The code employs the shock-fitting method, with an Eulerian formulation of the basic conservation equations and explicit tracking of shock waves. Pressure, temperature, and velocity are used as primary integration variables to simplify interfacing of the code with real-gas thermodynamic and transport property tables. Application of the code to the simulation of selected transients for the AEDC/HPDE MHD generator produced results that are in good agreement with experimental observations.

Wilson, D. R.↗

NASA NTRS2003

A Generalized Multi-Phase Framework for Modeling Cavitation in Cryogenic Fluids

A generalized multi-phase formulation for cavitation in fluids operating at temperatures elevated relative to their critical temperatures is presented. The thermal effects and the accompanying property variations due to phase change are modeled rigorously. Thermal equilibrium is assumed and fluid thermodynamic properties are specified along the saturation line using the NIST-12 databank. Fundamental changes in the physical characteristics of the cavity when thermal effects become pronounced are identified; the cavity becomes more porous, the interface less distinct, and has increased entrainment when temperature variations are present. Quantitative estimates of temperature and pressure depressions in both liquid nitrogen and liquid hydrogen were computed and compared with experimental data of Hord for hydrofoils. Excellent estimates of the leading edge temperature and pressure depression were obtained while the comparisons in the cavity closure region were reasonable. Liquid nitrogen cavities were consistently found to be in thermal equilibrium while liquid hydrogen cavities exhibited small, but distinct, non-equilibrium effects.

Dorney, Dan↗

NASA NTRS2006

Entropy Generation/Availability Energy Loss Analysis Inside MIT Gas Spring and "Two Space" Test Rigs

The results of the entropy generation and availability energy loss analysis under conditions of oscillating pressure and oscillating helium gas flow in two Massachusetts Institute of Technology (MIT) test rigs piston-cylinder and piston-cylinder-heat exchanger are presented. Two solution domains, the gas spring (single-space) in the piston-cylinder test rig and the gas spring + heat exchanger (two-space) in the piston-cylinder-heat exchanger test rig are of interest. Sage and CFD-ACE+ commercial numerical codes are used to obtain 1-D and 2-D computer models, respectively, of each of the two solution domains and to simulate the oscillating gas flow and heat transfer effects in these domains. Second law analysis is used to characterize the entropy generation and availability energy losses inside the two solution domains. Internal and external entropy generation and availability energy loss results predicted by Sage and CFD-ACE+ are compared. Thermodynamic loss analysis of simple systems such as the MIT test rigs are often useful to understand some important features of complex pattern forming processes in more complex systems like the Stirling engine. This study is aimed at improving numerical codes for the prediction of thermodynamic losses via the development of a loss post-processor. The incorporation of loss post-processors in Stirling engine numerical codes will facilitate Stirling engine performance optimization. Loss analysis using entropy-generation rates due to heat and fluid flow is a relatively new technique for assessing component performance. It offers a deep insight into the flow phenomena, allows a more exact calculation of losses than is possible with traditional means involving the application of loss correlations and provides an effective tool for improving component and overall system performance.

Ebiana, Asuquo B.↗

DOE OSTI2025

A physics-constrained neural ordinary differential equations approach for robust learning of stiff chemical kinetics

The high computational cost associated with solving for detailed chemistry poses a significant challenge for predictive computational fluid dynamics (CFD) simulations of turbulent reacting flows. While deep learning techniques have been explored to develop faster surrogate models, they often fail to integrate reliably with CFD solvers. This instability arises because traditional deep learning approaches optimize for training error without ensuring compatibility with ordinary differential equation (ODE) solvers, resulting in accumulation of errors over time. Recently, neuralODE (NODE) based approaches have been shown to be a promising technique to emulate and accelerate detailed chemistry computations. Here, in the present work, we extend this NODE framework for stiff chemical kinetics by incorporating mass conservation constraints directly into the loss function during training. This ensures that the total mass as well as the individual elemental species masses are conserved in an a-posteriori manner. Proof-of-concept studies are performed with the novel physics-constrained NODE (PC-NODE) approach for homogeneous autoignition of hydrogen-air mixture over a range of composition and thermodynamic conditions. It is demonstrated that the PC-NODE framework not only improves the physical consistency of the resulting data-driven model with respect to mass conservation criteria, but also improves training efficiency. PC-NODE is shown to achieve 2–100× speedup relative to the hydrogen-air detailed chemical mechanism depending on the type of the ODE solver (implicit or explicit) used during autoregressive inference tests. Lastly, a-posteriori studies are performed wherein the trained PC-NODE model is coupled with a CFD solver. It is shown that higher accuracy is achieved with PC-NODE relative to the purely data-driven NODE approach. Moreover, PC-NODE also exhibits robustness and generalizability to unseen initial conditions from within (interpolative capability) as well as outside (extrapolative capability) the training regime.

computational combustion↗

DOE OSTI2020

MARTINI-Compatible Coarse-Grained Model for the Mesoscale Simulation of Peptoids

Peptoids (poly-N-substituted glycines) are a class of synthetic polymers that are regioisomers of peptides (poly-C-substituted glycines), in which the point of side-chain connectivity is shifted from the backbone C to the N atom. Peptoids have found diverse applications as peptidomimetic drugs, protein mimetic polymers, surfactants, and catalysts. Computational modeling is valuable in the understanding and design of peptoid-based nanomaterials. In this work, we report the bottom-up parameterization of coarse-grained peptoid force fields based on the MARTINI peptide force field against all-atom peptoid simulation data. Our parameterization pipeline iteratively refits coarse-grained bonded interactions using iterative Boltzmann inversion and nonbonded interactions by matching the potential of mean force for chain extension. We assure good sampling of the amide bond cis/trans isomerizations in the all-atom simulation data using parallel bias metadynamics. We develop coarse-grained models for two representative peptoids—polysarcosine (poly(N-methyl glycine)) and poly(N-((4-bromophenyl)ethyl)glycine)—and show their structural and thermodynamic properties to be in excellent accord with all-atom calculations but up to 25-fold more efficient and compatible with MARTINI force fields. Here, this work establishes a new rigorously parameterized coarse-grained peptoid force field for the understanding and design of peptoid nanomaterials at length and time scales inaccessible to all-atom calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Quasi‐Binary Transition Metal Dichalcogenide Alloys: Thermodynamic Stability Prediction, Scalable Synthesis, and Application

Abstract Transition metal dichalcogenide (TMDCs) alloys could have a wide range of physical and chemical properties, ranging from charge density waves to superconductivity and electrochemical activities. While many exciting behaviors of unary TMDCs have been demonstrated, the vast compositional space of TMDC alloys has remained largely unexplored due to the lack of understanding regarding their stability when accommodating different cations or chalcogens in a single‐phase. Here, a theory‐guided synthesis approach is reported to achieve unexplored quasi‐binary TMDC alloys through computationally predicted stability maps. Equilibrium temperature–composition phase diagrams using first‐principles calculations are generated to identify the stability of 25 quasi‐binary TMDC alloys, including some involving non‐isovalent cations and are verified experimentally through the synthesis of a subset of 12 predicted alloys using a scalable chemical vapor transport method. It is demonstrated that the synthesized alloys can be exfoliated into 2D structures, and some of them exhibit: i) outstanding thermal stability tested up to 1230 K, ii) exceptionally high electrochemical activity for the CO 2 reduction reaction in a kinetically limited regime with near zero overpotential for CO formation, iii) excellent energy efficiency in a high rate Li–air battery, and iv) high break‐down current density for interconnect applications. This framework can be extended to accelerate the discovery of other TMDC alloys for various applications.

Hemmat, Zahra↗

DOE OSTI2025

Unveiling the thermite-driven lithium fire ignition in solid-state batteries

Here, this study challenges the assumption of the non-flammability of lithium metal all-solid-state batteries (LiSSBs) and other lithium metal batteries without flammable electrolytes. Through thermodynamic calculations and ex situ experiments, we reveal for the first time the risk of thermite reactions between lithium metal and LiFePO 4 in both charged and discharged states. Reactivity is worsened by excess lithium metal in the cell, reaching final maximum adiabatic temperatures of 2,500°C in the charged state, which is hot enough to boil lithium. The thermite reaction triggers spontaneously at 500°C, with poor surface contact, while increasing surface contact through mixing initiates the reaction at room temperature in an inert environment. Despite its fast kinetics, this reaction is transport limited due to lithium passivation, leading to long burn times and reignition risks. Given the risk of lithium metal contacting the cathode during failure, understanding these reactions is crucial for ensuring the safe deployment of LiSSBs.

LiFePO4↗

DOE OSTI2025

Site-specific surface reactivity on SnO 2 : Evaluating selective atomic layer deposition processes

Area selective atomic layer deposition (AS-ALD) is a bottom-up synthesis approach with potential for deposition with molecular level precision. Here, the site-specific hydration of metal oxide substrates, combined with surface H 2 O-selective ALD processes, provides a potentially powerful path to targeted synthesis. Density functional theory (DFT) calculations are used to predict the thermodynamics of ALD precursor reactivity and hydration for (001), (101), (110), and (100) rutile SnO 2 facets as a function of temperature. Trimethylaluminum (TMA) and dimethyl aluminum isopropoxide (DMAI) dimers are predicted to react with both dehydrated and hydrated SnO 2 (001), (101), and (110) facets at ALD-relevant temperatures, while the SnO 2 (100) facet is predicted to be uniquely unreactive with TMA and DMAI monomers as well as dehydrate near 177 °C making this facet more amenable to targeted ALD. In situ ellipsometric studies of Al 2 O 3 ALD on polycrystalline SnO 2 at 150 °C are consistent with the computational predictions of rapid and unselective nucleation, in stark contrast to inhibited and selective ALD on isostructural rutile TiO 2 .

Atomic Layer Deposition↗

DOE OSTI2021

Analytic and numerical vibronic spectra from quasi-classical trajectory ensembles

The truncated Wigner approximation to quantum dynamics in phase space is explored in the context of computing vibronic line shapes for monomer linear optical spectra. We consider multiple model potential forms including a shifted harmonic oscillator with both equal and unequal frequencies on the ground and excited state potentials as well as a shifted Morse potential model. For the equal-frequency shifted harmonic oscillator model, we derive an analytic expression for the exact vibronic line shape that emphasizes the importance of using a quantum mechanical distribution of phase space initial conditions. For the unequal-frequency shifted harmonic oscillator model, we are no longer able to obtain an exact expression for the vibronic line shape in terms of independent deterministic classical trajectories. We show how one can rigorously account for corrections to the truncated Wigner approximation through nonlinear responses of the line shape function to momentum fluctuations along a classical trajectory and demonstrate the qualitative improvement in the resulting spectrum when the leading-order quantum correction is included. Finally, we numerically simulate absorption spectra of a highly anharmonic shifted Morse potential model. We find that, while finite quantization and the dissociation limit are captured with reasonable accuracy, there is a qualitative breakdown of the quasi-classical trajectory ensemble’s ability to describe the vibronic line shape when the relative shift in Morse potentials becomes large. The work presented here provides clarity on the origin of unphysical negative features known to contaminate absorption spectra computed with quasi-classical trajectory ensembles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Theoretical study of CO2 capture mechanisms of SrO and Sr(OH)2.nH2O (n=0,1,8)

In this presentation, the adsorption reactions of CO2 on SrO surfaces and its hydrated counterparts are examined. The predicted energies of CO2 adsorption on follows the order of increasing H2O content. The carbonation on the monohydrate surface is barrierless, whereas the CO2 reaction on the Sr(OH)2 and the octahydrate surfaces follows pathways with finite activation barriers. The thermodynamics of bulk reaction is also addressed.

computational science and engineering↗

NASA NTRS2013

Isentropic Analysis of Convective Motions

This paper analyzes the convective mass transport by sorting air parcels in terms of their equivalent potential temperature to determine an isentropic streamfunction. By averaging the vertical mass flux at a constant value of the equivalent potential temperature, one can compute an isentropic mass transport that filters out reversible oscillatory motions such as gravity waves. This novel approach emphasizes the fact that the vertical energy and entropy transports by convection are due to the combination of ascending air parcels with high energy and entropy and subsiding air parcels with lower energy and entropy. Such conditional averaging can be extended to other dynamic and thermodynamic variables such as vertical velocity, temperature, or relative humidity to obtain a comprehensive description of convective motions. It is also shown how this approach can be used to determine the mean diabatic tendencies from the three-dimensional dynamic and thermodynamic fields. A two-stream approximation that partitions the isentropic circulation into a mean updraft and a mean downdraft is also introduced. This offers a straightforward way to identify the mean properties of rising and subsiding air parcels. The results from the two-stream approximation are compared with two other definitions of the cloud mass flux. It is argued that the isentropic analysis offers a robust definition of the convective mass transport that is not tainted by the need to arbitrarily distinguish between convection and its environment, and that separates the irreversible convective overturning fromoscillations associated with gravity waves.

Lagrangian circulation↗

DOE OSTI2021

Defying Thermodynamics: Stabilization of Alane Within Covalent Triazine Frameworks for Reversible Hydrogen Storage

The highly unfavorable thermodynamics of direct aluminum hydrogenation can be overcome by stabilizing alane within a nanoporous bipyridine-functionalized covalent triazine framework (AlH 3 @CTF-bipyridine). This material and the counterpart AlH 3 @CTF-biphenyl rapidly desorb H 2 between 95 and 154°C, with desorption complete at 250°C. Sieverts measurements, 27 Al MAS NMR and 27 Al{ 1 H} REDOR experiments, and computational spectroscopy reveal that AlH 3 @CTF-bipyridine dehydrogenation is reversible at 60°C under 700 bar hydrogen, >10 times lower pressure than that required to hydrogenate bulk aluminum. DFT calculations and EPR measurements support an unconventional mechanism whereby strong AlH 3 binding to bipyridine results in single-electron transfer to form AlH 2 (AlH 3 ) n clusters. Here the resulting size-dependent charge redistribution alters the dehydrogenation/rehydrogenation thermochemistry, suggesting a novel strategy to enable reversibility in high-capacity metal hydrides.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗