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,585 records · Page 88
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Modeling Equilibrium Solid–Liquid Interfaces under Effective Constant Chemical Potential Using Machine Learning Interatomic Potentials

The chemical potential (μ) of species in solution is essential for understanding various chemical processes at interfaces. Molecular dynamics (MD) simulations, constrained by fixed compositions, cannot maintain constant chemical potential with reference to a targeted concentration or chemical potential under nonequilibrium or dynamic conditions, as solute species can migrate to the interface and deplete (or enrich) the bulk due to solute-interface interactions. In this study, we introduce a simple and computationally efficient approach named iterative quasi-constant chemical potential molecular dynamics (iqCμMD) simulation, which helps simulate targeted molar concentrations of species in solution. iqCμMD overcomes the limitations of conventional MD by adjusting the number of species in the solution to reach a target bulk concentration (chemical potential), which allows simulation of the interface under the bulk conditions comparable to experiment. We demonstrate our approach using machine learning interatomic potential (MLIP)-based MD simulations of the Na 2 SO 4,aq –graphene interface, and to show the transferability of our approach, we also perform classical force field-based MD simulations of NaCl aq –air and NaCl aq –graphite interfaces, which produce comparable results to previous CμMD simulations. Our results also show that the iqCμMD approach efficiently achieves the desired bulk ion concentration within two iterations, and by utilizing MLIPs, we can achieve converged results using relatively small-scale simulations compared to previous CμMD simulations. By combining iqCμMD with MLIP-driven simulations, solid–liquid interfaces can be modeled under an effective constant chemical potential with DFT-level accuracy. Here, we show that iqCμMD offers a robust and simple computational framework for constant chemical potential simulations, as its only requirement is to be able to converge interfacial simulations with a measurable bulk region.

Chemical structure↗
DOE OSTI2024

Water–Hydrocarbon Interactions in Anionic Pyrene Monohydrate

Interactions between water and polycyclic aromatic hydrocarbons are essential in many aspects of chemistry, from interstellar and atmospheric processes to interfacial hydrophobicity and wetting phenomena. Despite their growing importance, the intermolecular potentials of the water-hydrocarbon interactions are underdeveloped compared to water-water potentials, and there are similarly few experimental probes that are sensitive to the details of the water-hydrocarbon potential. We present a combined experimental and computational study of anionic pyrene monohydrate, one of the simplest water/hydrocarbon clusters. The action spectrum in the OH region of the mass-selected cluster ion provides a rigorous benchmark for intermolecular potentials and computational methodologies. We identify missing intermolecular interactions and shortcomings in conventional dynamics calculations by comparing experimental data to density functional theory and classical molecular dynamics calculations. Kinetic trapping is prevalent, even for one water molecule and one pyrene molecule, leading to slow equilibration in conventional molecular dynamics calculations, even on nanosecond timescales and at low temperatures (50 K). At constant energy, temperature fluctuations for the pair of molecules are substantial. Immersing the system in a bath of soft spheres and employing parallel tempering alleviates kinetic trapping and dampens temperature fluctuations, bringing the system closer to the thermodynamic limit. With such augmented sampling, a simple, flexible water model reproduces the linewidth and the asymmetric broadening of the symmetric OH stretching mode, which we assign to spectral diffusion. In the OH stretching region, dynamics calculations predict a more intense antisymmetric peak than experiments observe but do not predict the bimodal split symmetric peak that the experiments show. Furthermore, our work suggests that electronic polarization, missing in the empirical force field, is responsible for the first discrepancy and that quantum nuclear effects, captured neither in density functional theory nor in classical dynamics, may be responsible for the second.

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

Benchmarking nitrous oxide adsorption and activation in metal–organic frameworks bearing coordinatively unsaturated metal centers

Anthropogenic emissions of N 2 O, the third most abundant greenhouse gas after CO 2 and CH 4 , are contributing to global climate change. Although metal–organic frameworks (MOFs) have been widely studied as adsorbents for CO 2 and CH 4 , less effort has focused on the use of MOFs to remove N 2 O from emission streams or from air. Further, N 2 O activation would enable its use as an inexpensive oxidant for fine chemical synthesis. Here, in this paper, we identify features that contribute to strong binding and high uptake of N 2 O at coordinatively unsaturated metal sites in the M 2 Cl 2 (btdd) (M = Mn, Co, Ni, Cu; btdd 2– = bis(1,2,3-triazolo[4,5-b],[4',5'-i])dibenzo[1,4]dioxin) and M 2 (dobdc) (M = Mg, Mn, Fe, Co, Ni, Cu, Zn; dobdc 4– = 2,5-dioxido-1,4-benzenedicarboxylate) series of MOFs. Combined experimental and computational studies suggest that N 2 O adsorption at open-metal-sites is primarily based on electrostatic interactions, rather than π-backbonding, causing MOFs with more Lewis acidic metal centers to be superior N 2 O adsorbents. As a result, Mg 2 (dobdc) demonstrates strong binding and record-setting N 2 O uptake (8.75 mmol g –1 at 1 bar and 298 K). Using density functional theory (DFT) to characterize reactive intermediates and transition states, we demonstrate that N 2 O activation to form a M(IV)–oxo species and N 2 is thermodynamically favorable in Mn 2 (dobdc) and Fe 2 (dobdc) but appears to be kinetically limited in Mn 2 (dobdc). Our work lays a foundation for understanding N 2 O adsorption and activation in MOFs, paving the way for the design of promising next-generation materials for N 2 O capture and utilization.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Synergistic Effect of Diatomic Mo-B Site Confined in Graphene-Like C 2 N Enables Electrocatalytic Nitrogen Reduction via Novel Mechanism

Structural modulation of active site with atomic-level precision is of great importance to meet the activity and selectivity challenges that electrocatalysts are commonly facing. In this work, we have designed a metal (M) - nonmetal (NM) diatomic site embedded in graphene-like C 2 N (denoted as Mo-B@C 2 N), where electrocatalytic N 2 reduction reaction (eNRR) was thoroughly explored using density functional theory combined with computational hydrogen electrode method. Compared to M-M diatomic sites, Mo-B site can generate pronounced synergistic effect that led to eNRR proceeding via a novel quasi-dissociative reaction mechanism that has not been reported relative to the conventional enzymatic, consecutive, distal, and alternating associative mechanism. This newly uncovered mechanism in which N-N bond scission takes place immediately after the first proton-coupled electron transfer (PCET) step (i.e., *NH-*N + H + + e – ® *NH 2 *N) has demonstrated much advantage in PCET process over the four conventional mechanism in terms of thermodynamic barrier, except that the adsorption of side-on *N 2 seemed thermodynamically unfavorable (DG ads = 0.61eV). Our results have revealed that the activation of inert NºN triple bond is dominated by π*-backdonation mechanism as a consequence of charge transfers from both B and Mo site, and unexpectedly, from the substrate C 2 N itself as well. Moreover, the hybrid Mo-B diatomic site demonstrated superior performance over either Mo-Mo or B-B site for driving eNRR. Furthermore, our study could provide insight into the delicate relationships among atomic site, substrate and electrocatalytic performance.

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

Consistent thermodynamic properties for alicyclic components of jet fuels: Experimental data, estimation methods, and homologous series trends

Alkylcycloalkanes represent a significant fraction of jet fuel components. An evaluation of their thermodynamic properties, enthalpies of formation in liquid and gas phases and enthalpies of vaporization, was conducted. A combination of available experimental data, up-to-date group-contribution methods, high-level quantum-chemical calculations, and homologous series trends was used to identify outliers and to recommend the most reliable values. The group-contribution approach was found to work well for the enthalpies of vaporization. Its performance for the enthalpies of formation in the liquid and gas phases was found to be substantially less effective, especially considering notable differences in this property among stereoisomers. Computationally affordable high-level ab initio results and homologous series trend analysis appeared more reliable. In conclusion, the recommended property values for 212 individual compounds and their isomeric mixtures were provided.

09 BIOMASS FUELS↗
DOE OSTI2024

Bipartite mutual information in classical many-body dynamics

Information theoretic measures have helped to sharpen our understanding of many-body quantum states. As perhaps the most well-known example, the entanglement entropy (or more generally, the bipartite mutual information) has become a powerful tool for characterizing the dynamical growth of quantum correlations. By contrast, although computable, the bipartite mutual information (MI) is almost never explored in classical many particle systems; this owes in part to the fact that computing the MI requires keeping track of the evolution of the full probability distribution, a feat which is rarely done (or thought to be needed) in classical many-body simulations. Here, we utilize the MI to analyze the spreading of information in 1D elementary cellular automata (CA). Broadly speaking, we find that the behavior of the MI in these dynamical systems exhibits a few different types of scaling that roughly correspond to known CA universality classes. Of particular note is that we observe a set of automata for which the MI converges parametrically slowly to its thermodynamic value. We develop a microscopic understanding of this behavior by analyzing a two-species model of annihilating particles moving in opposite directions. Furthermore, our work suggests the possibility that information theoretic tools such as the MI might enable a more fine-grained characterization of classical many-body states and dynamics.

Cellular automata↗
NASA NTRS1992

Proposed radiometric measurement of the wake of a blunt aerobrake

This paper describes the aerothermal environment in the afterbody region of a blunt entry body. Recent ground-based experiments and computational predictions of the afterbody flow structure and radiation are presented. The similarity between the flowfield structures observed in the ground-based experiments and that obtained by calculation is encouraging. Approximate calculations of the radiative heating rate to the base are presented. Many of the phenomena associated with the expanding flow at the corner and the formation of the wake neck, however, are not well understood and require further study. A flight experiment is described that would use spectral and total measurements of the wake radiation as a nonintrusive diagnostic method to provide insight into the thermodynamic state of the wake gas.

Strawa, A. W.↗
NASA NTRS2017

Computational and Experimental Characterization of the Mach 6 Facility Nozzle Flow for the Enhanced Injection and Mixing Project at NASA Langley Research 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.

Drozda, Tomasz G.↗
NASA NTRS1987

A thermodynamic model of the Hg(0.8)Cd(0.2)Te-iodine transport system. I - Te-saturated source material. II - Source material composition within the homogeneity range

A thermochemical analysis of the Hg(0.8)Cd(0.2)Te-iodine vapor transport system is presented, and theoretical calculations of diffusion-controlled mass transport rates are made. The predicted mass fluxes are compared with experimental data obtained from transport experiments under vertical, stabilizing conditions reported earlier and with results of additional transport experiments conducted during the present study. Experimental mass transport rate studies of the transport system for fixed amount of excess Hg as a function of transport agent pressure are presented. The mass fluxes are determined for the vertical, stabilizing orientation of the density gradient relative to the gravitational vector. In order to compare experimental mass transport rates with computed values, the thermochemical analysis is extended to take the formation of Hg vacancies in the above compound into account along with their effect on the partial pressure of the system.

Chandra, D.↗
DOE OSTI2025

Lanczos Algorithm, the Transfer Matrix, and the Signal-to-Noise Problem

This Letter introduces a method for determining the energy spectrum of lattice quantum chromodynamics by applying the Lanczos algorithm to the transfer matrix and using a bootstrap generalization of the Cullum-Willoughby method to filter out spurious eigenvalues. Proof-of-principle analyses of the simple harmonic oscillator and the lattice quantum chromodynamics proton mass demonstrate that this method provides faster ground-state convergence than the “effective mass,” which is related to the power-iteration algorithm. Lanczos provides more accurate energy estimates than multistate fits to correlation functions with small imaginary times while achieving comparable statistical precision. Two-sided error bounds are computed for Lanczos results and guarantee that excited-state effects cannot shift Lanczos results far outside their statistical uncertainties.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Influence of Real Gas Effects on Chemical Kinetics in Oxycombustion in Supercritical Carbon Dioxide

Oxycombustion in supercritical carbon dioxide is an integral part of the Allam Cycle, a technology that enables carbon-neutral use of fossil-fuels and carbon-negative use of biofuels. We simulate oxycombustion in a realistic combustor geometry for two sets of fuel conditions: pure methane, and a 40 percent methane/60 percent carbon dioxide blend, both at 343.15 K. The fuel jet mixes with a preheated swirler of 20 percent oxygen and 80 percent carbon dioxide at 1005.35 K, with a 100 percent carbon dioxide coflow at 783.15 K. The entire system operates at a pressure of 300 bar, putting the entire system above the critical temperature and pressure of carbon dioxide. Simulations are performed using PeleC, a compressible block-structured adaptive-mesh refinement (AMR) reacting flow code. Direct comparison of results using thermodynamically self-consistent implementations of the ideal gas equation of state and the Soave-Redlich-Kwong (SRK) equation of state allow quantification of the real gas impacts on combustion kinetics, flow temperatures, and pollutant formation.

carbon dioxide↗
DOE OSTI2022

Structure, Stability, and Electronic Properties of Boron Suboxide: A Density Functional Theory Study

Boron suboxide (B 6 O) is a boron-rich compound derived from the α-rhombohedral boron lattice with extreme hardness and unusual semiconducting properties. In this work, density functional theory (DFT) was used to show that unit cell volume, mechanical strength, band gaps, and thermodynamic stabilities of B 6 O were influenced by the interstitial elements and point defects at the icosahedral sites. While the hexagonal unit cell volume (HUCV) varies with interstitial occupancy, it is the icosahedral defect that weakens the intrinsic bulk modulus of B 6 O. Using the hybrid HSE functional, we confirmed that the perfect B 6 O bulk is a p-type semiconductor with a direct band gap of 2.8 eV. Furthermore, by screening α-boron compounds systematically, we found that a simple octet rule may offer a consistent explanation for the variations in the computed electronic structures. The formation free energies calculated over a wide range of temperatures (0–2500 K) and pressures (0–80 GPa) predict that formations of interstitial defects become favorable only at higher temperatures (ca. 1800 K) in bulk B 6 O lattices. The nudged elastic band (NEB) method was employed to identify the minimum energy pathways for the diffusions of dislocated B and O atoms. The diffusion of icosahedral B atoms has an energy barrier of 0.16 eV. More complex B diffusion paths involving the reorganization of icosahedral boron atoms incur higher barriers (>1 eV). In contrast, the diffusion of interstitial O atoms is facile with a barrier of 0.4 eV. Lastly, successive O insertions into the α-B lattice were performed using DFT to generate a basic understanding of the oxidation process. Furthermore, these calculations provide fundamental atomistic insights into the growth of B 6 O crystals and control of their point defects.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Transport coefficients of the quark-gluon plasma at the critical point and across the first-order line

A bottom-up Einstein-Maxwell-dilaton holographic model is used to compute, for the first time, the behavior of several transport coefficients of the hot and baryon-rich strongly coupled quark-gluon plasma at the critical point and also across the first-order phase transition line in the phase diagram. The observables under study are the shear and bulk viscosities, the baryon and thermal conductivities, the baryon diffusion, the jet quenching parameter $\hat{q}$, as well as the heavy-quark drag force and the Langevin diffusion coefficients. These calculations provide a phenomenologically promising estimate for these coefficients, given that our model quantitatively reproduces lattice QCD thermodynamics results, both at zero and finite baryon density, besides naturally incorporating the nearly perfect fluidity of the quark-gluon plasma. Here we find that the diffusion of baryon charge, and also the shear and bulk viscosities, are suppressed with increasing baryon density, indicating that the medium becomes even closer to perfect fluidity at large densities. On the other hand, the jet quenching parameter and the heavy-quark momentum diffusion are enhanced with increasing density. The observables display a discontinuity gap when crossing the first-order phase transition line, while developing an infinite slope at the critical point. The transition temperatures associated with different transport coefficients differ in the crossover region but are found to converge at the critical point.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

BattPhase—A Convergent, Non-Oscillatory, Efficient Algorithm and Code for Predicting Shape Changes in Lithium Metal Batteries Using Phase-Field Models: Part I. Secondary Current Distribution

Electrochemical models at different scales and varying levels of complexity have been used in the literature to study the evolution of the anode surface in lithium metal batteries. This includes continuum, mesoscale (phase-field approaches), and multiscale models. Thermodynamics-based equations have been used to study phase changes in lithium batteries using phase-field approaches. However, grid convergence studies and the effect of additional parameters needed to simulate these models are not well-documented in the literature. In this paper, using a motivating example of a moving boundary model in one- and two-dimensions, we show how one can formulate phase-field models, implement algorithms for the same and analyze the results. An open-access code with no restrictions is provided as well. This article concludes with some thoughts on the computational efficiency of phase-field models for simulating dendritic growth.

25 ENERGY STORAGE↗
DOE OSTI2022

Simulating Effective QED on Quantum Computers

In recent years simulations of chemistry and condensed materials has emerged as one of the preeminent applications of quantum computing, offering an exponential speedup for the solution of the electronic structure for certain strongly correlated electronic systems. To date, most treatments have ignored the question of whether relativistic effects, which are described most generally by quantum electrodynamics (QED), can also be simulated on a quantum computer in polynomial time. Here we show that effective QED, which is equivalent to QED to second order in perturbation theory, can be simulated in polynomial time under reasonable assumptions while properly treating all four components of the wavefunction of the fermionic field. In particular, we provide a detailed analysis of such simulations in position and momentum basis using Trotter-Suzuki formulas. We find that the number of T -gates needed to perform such simulations on a 3 D lattice of n s sites scales at worst as O ( n s 3 / ϵ ) 1 + o ( 1 ) in the thermodynamic limit for position basis simulations and O ( n s 4 + 2 / 3 / ϵ ) 1 + o ( 1 ) in momentum basis. We also find that qubitization scales slightly better with a worst case scaling of O ~ ( n s 2 + 2 / 3 / ϵ ) for lattice eQED and complications in the prepare circuit leads to a slightly worse scaling in momentum basis of O ~ ( n s 5 + 2 / 3 / ϵ ) . We further provide concrete gate counts for simulating a relativistic version of the uniform electron gas that show challenging problems can be simulated using fewer than 10 13 non-Clifford operations and also provide a detailed discussion of how to prepare multi-reference configuration interaction states in effective QED which can provide a reasonable initial guess for the ground state. Finally, we estimate the planewave cutoffs needed to accurately simulate heavy elements such as gold.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1972

Some Concepts in Reaction Dynamics

In 1929 London 1 published a very approximate solution of the Schroedinger equation for a system of chemical interest: H3. To the extent that chemistry can be regarded as existing separately from physics, this was a landmark in the history of chemistry, comparable in importance to the landmark in the history of physics marked by the appearance of the Heitler-London equation for H2. The expression for H3, was, of necessity, even less accurate than that for H2, but chemists, like the habitual poor, were accustomed to this sort of misfortune. Together with the physicists they enjoyed the sensation of living in a renaissance. The physicists still could not calculate a great deal that was of interest to them, and the chemists could calculate less, but both could now dream. It would be too easy to say that their dreams were dreams of unlimited computer time. Their dreams were a lot more productive than that. Two years after London published his equation, H. Eyring and M. Polanyi obtained the first numerical energy surface for H3. They infused the London equation with a measure of empiricism to produce an energy surface which, whether or not it was correct in its details, provided a basis for further speculations of an important sort. The existence of a tangible energy surface in 1931 stimulated speculation along two different lines. The following year Pelzer and Wigner used this London-Eyring-Polanyi (LEP) energy surface for a thermodynamic treatment of the reaction rate in H + H2. This important development reached its full flowering a few years later. In these remarks I shall be concerned with another line of development. A second more-or-less distinct category of speculation that began with (and, indeed, in) the 1931 paper has to do with the dynamics of individual reactive encounters under the influence of specified interaction potentials.

Polannyi, John C.↗
DOE OSTI2024

Engineered Nanoporous Frameworks for Adsorption Cooling Applications

The energy demand for traditional vapor-compressed technology for space cooling continues to soar year after year due to global warming and the increasing human population’s need to improve living and working conditions. Thus, there is a growing demand for eco-friendly technologies that use sustainable or waste energy resources. Here this review discusses the properties of various refrigerants used for adsorption cooling applications followed by a brief discussion on the thermodynamic cycle. Next, sorbents traditionally used for cooling are reviewed to emphasize the need for advanced capture materials with superior properties to improve refrigerant sorption. The remainder of the review focus on studies using engineered nanoporous frameworks (ENFs) with various refrigerants for adsorption cooling applications. The effects of the various factors that play a role in ENF–refrigerant pair selection, including pore structure/dimension/shape, morphology, open-metal sites, pore chemistry and possible presence of defects, are reviewed. Next, in-depth insights into the sorbent–refrigerant interaction, and pore filling mechanism gained through a combination of characterization techniques and computational modeling are discussed. Finally, we outline the challenges and opportunities related to using ENFs for adsorption cooling applications and provide our views on the future of this technology.

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

Liquid Structure and Transport Properties of the Deep Eutectic Solvent Ethaline

A range of techniques including physical property measurements, neutron scattering experiments, ab initio molecular dynamics, and classical molecular dynamics simulations are used to probe the structural, thermodynamic, and transport properties of a deep eutectic solvent comprised of a 1:2 molar ratio of choline chloride and ethylene glycol. This mixture, known as Ethaline, has many desirable properties for use in a range of applications, and therefore, understanding its liquid structure and transport properties is of interest. Simulation results are able to capture experimental densities, diffusivities, viscosities, and structure factors extremely well. Here, the solvation environment is dynamic and dominated by different hydrogen bonding interactions. Dynamic heterogeneities resulting from hydrogen bonding interactions are quantified. Rotational dynamics of molecular dipole moments of choline and ethylene glycol are computed and found to exhibit a fast and slow mode.

25 ENERGY STORAGE↗