Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Classical thermodynamics”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Estimation of isosteric heat of adsorption from generalized Langmuir isotherm

Simulation and design of adsorptive separation units demand accurate estimation of thermodynamic properties. Isosteric heat of adsorption as calculated from generalized Langmuir (gL) isotherm coupled with Clausius-Clapeyron expression for pure component and mixed-gas adsorption equilibria is presented in this work. The estimated isosteric heat of adsorption as functions of surface loading and composition is validated against the experimental data for various adsorption systems. Furthermore, the gL results are compared against classical Langmuir (cL) and Toth isotherm for pure components and with Ideal Adsorbed Solution Theory (IAST) for mixed-gas adsorption equilibria. Furthermore, the comparison highlights that gL outperforms cL and Toth for pure component adsorption and IAST for mixed-gas adsorption, and gL reliably captures the loading dependence and the composition dependence for isosteric heat of adsorption.

42 ENGINEERING↗
NASA NTRS1977

Concepts and methods for describing critical phenomena in fluids

The predictions of theoretical models for a critical-point phase transistion in fluids, namely the classical equation with third-degree critical isotherm, that with fifth-degree critical isotherm, and the lattice gas, are reviewed. The renormalization group theory of critical phenomena and the hypothesis of universality of critical behavior supported by this theory are discussed as well as the nature of gravity effects and how they affect cricital-region experimentation in fluids. The behavior of the thermodynamic properties and the correlation function is formulated in terms of scaling laws. The predictions of these scaling laws and of the hypothesis of universality of critical behavior are compared with experimental data for one-component fluids and it is indicated how the methods can be extended to describe critical phenomena in fluid mixtures.

Sengers, J. V.↗
DOE OSTI2021

High-pressure phase diagram of beryllium from ab initio free-energy calculations

In this report we use first-principles molecular dynamics simulations coupled with the thermodynamic integration method to study the hexagonal close-packed (hcp) to body-centered cubic (bcc) transition and melting of beryllium up to a pressure of 1600 GPa. We derive the melting line by equating solid and liquid Gibbs free energies and represent it by a Simon-Glatzel fit T m = 1564 K [1 + P/(15.6032 GPa)] 0.383 , which is in good agreement with previous two-phase simulations <6000 K. We also derive the hcp-bcc solid-solid phase boundary and show that the quasiharmonic approximation underestimates the stability of the hcp structure, predicting lower transition pressures between hcp and bcc phases. Our results are consistent with the stability regime predicted by the phonon quasiparticle method. We also predict that the hcp-bcc-liquid triple point is located at 164.7 GPa and 4314 K. In addition, we compute the shock Hugoniot curve and show that it is in good agreement with experiments, intersecting our derived melting curve at ~235 GPa and 4900 K. Finally, we make predictions for future ramp compression experiments. Starting with an isentropic compression of the liquid, we predict the path to intersect the melting line at low pressure and temperature, then to continue along the melting line over a large temperature interval of 7000 K as the sample remains in the mixed solid-liquid state before it enters the solid phase.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Time-reversal of an unknown quantum state

For decades, researchers have sought to understand how the irreversibility of the surrounding world emerges from the seemingly time-symmetric, fundamental laws of physics. Quantum mechanics conjectured a clue that final irreversibility is set by the measurement procedure and that the time-reversal requires complex conjugation of the wave function, which is overly complex to spontaneously appear in nature. Building on this Landau-Wigner conjecture, it became possible to demonstrate that time-reversal is exponentially improbable in a virgin nature and to design an algorithm artificially reversing a time arrow for a given quantum state on the IBM quantum computer. However, the implemented arrow-of-time reversal embraced only the known states initially disentangled from the thermodynamic reservoir. Here we develop a procedure for reversing the temporal evolution of an arbitrary unknown quantum state. This opens the route for general universal algorithms sending temporal evolution of an arbitrary system backward in time.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Single T gate in a Clifford circuit drives transition to universal entanglement spectrum statistics

Clifford circuits are insufficient for universal quantum computation or creating t t -designs with t\ge 4 t ≥ 4 . While the entanglement entropy is not a telltale of this insufficiency, the entanglement spectrum of a time evolved random product state is: the entanglement levels are Poisson-distributed for circuits restricted to the Clifford gate-set, while the levels follow Wigner-Dyson statistics when universal gates are used. In this paper we show, using finite-size scaling analysis of different measures of level spacing statistics, that in the thermodynamic limit, inserting a single T (\pi/8) ( π / 8 ) gate in the middle of a random Clifford circuit is sufficient to alter the entanglement spectrum from a Poisson to a Wigner-Dyson distribution.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Temperature and density dependent pair potential for deuterium under shock

Large-scale classical molecular dynamics (CMD) simulations naturally include the microscopic physics necessary for atomistic modeling of shock release at the ablator-fuel interface in an inertial confinement fusion (ICF) capsule. Here, the multi-megabar shocks utilized in ICF experiments can drive the deuterium fuel from ambient to electron volt temperatures (T) and multi-fold compression. Modeling interatomic interactions over such an extreme range of conditions is challenging for empirical bond order potentials. We generate a pair potential for deuterium with explicit temperature and mass density dependence from ab initio density functional theory molecular dynamics using the iterative Boltzmann inversion method. This potential accurately reproduces the radial distribution functions and pressures from DFT in CMD equilibrium simulations across a wide range of thermodynamic conditions, yet fails to return the expected Hugoniot relations when used in direct CMD shock simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Rate constants and third-body collision efficiencies for recombination of Na with OH and O2: Implications for flame inhibition by alkali salts

Flame inhibition by alkali metals has implications for solid fuel combustion, fire safety, and a number of industrial processes. While the mechanism of inhibition is fairly well understood, details related to thermodynamic properties and rate constants are still in question. In the present work, the recombination of Na with OH (R5) and O2 (R7), respectively, was characterized theoretically, and the implications for modeling laminar premixed hydrogen flames doped with sodium species were examined. Third-body collision efficiencies and low-pressure-limit rate constants were obtained using newly fitted ab initio-based potential energy surfaces, classical trajectories, and one-dimensional master equation calculations. The results are consistent with available experimental results and aid in the present modeling study by providing rate information for bath gases and conditions that remain unexplored experimentally. Chemical kinetic modeling of relative Na and absolute H and OH profiles in H2-fueled laminar, premixed flames doped with a sodium salt shows that the most important radical removal cycle is the sequence Na + OH (+M) → NaOH (+M) (R5), NaOH + H → Na + H2O (R12), even under oxidizing conditions. A secondary cycle, Na + O2 (+M) ⇄ NaO2 (+M) (R7), NaO2 + OH → NaOH + O2 (R14), is less important due to the low thermal stability of NaO2. In most flames, reaction R7 is partially equilibrated, and reaction R14 becomes rate-limiting for the second cycle. The flame analysis supports a lower value of k14 than indicated by recent work on KO2 + OH, but more work is required to confirm this.

Jasper, Ahren W.↗
DOE OSTI2020

Quantum Joule expansion of one-dimensional systems

Here, we investigate the Joule expansion of nonintegrable quantum systems that contain bosons or spinless fermions in one-dimensional lattices. A barrier initially confines the particles to be in half of the system in a thermal state described by the canonical ensemble and is removed at time t = 0. We investigate the properties of the time-evolved density matrix, the diagonal ensemble density matrix, and the corresponding canonical ensemble density matrix with an effective temperature determined by the total energy conservation using exact diagonalization. The weights for the diagonal ensemble and the canonical ensemble match well for high initial temperatures that correspond to negative effective final temperatures after the expansion. At long times after the barrier is removed, the time-evolved Rényi entropy of subsystems bigger than half can equilibrate to the thermal entropy with exponentially small fluctuations. The time-evolved reduced density matrix at long times can be approximated by a thermal density matrix for small subsystems. Few-body observables, like the momentum distribution function, can be approximated by a thermal expectation of the canonical ensemble with strongly suppressed fluctuations. The negative effective temperatures for finite systems go to non-negative temperatures in the thermodynamic limit for bosons, but there is a true thermodynamic effect for fermions, which is confirmed by finite-temperature density matrix renormalization group calculations. We propose the Joule expansion as a way to dynamically create negative temperature states for fermion systems with repulsive interactions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Experimental and Computational Studies of Crystal Nucleation in Composition Gradients (Final Report)

The major goals of this project were to develop and validate a predictive nucleation model that incorporates composition gradients and is applicable to a large range of metallic systems which form both stochiometric and non-stochiometric compounds. Computationally we planned to expand the extent of thermodynamics-based theories and lay the groundwork to improve the general understanding and predictive capabilities of the role that gradients play in phase formation, glass formability and stability, and nucleation and growth events. To address these goals, we used Molecular Dynamics (MD) simulations and both isothermal and isochronal nanocalorimetric experiments on amorphous phases with a controlled composition gradient with the hopes of validating and improving the classical nucleation model and its use in solid solutions. While the computations were successful and identified ways to improve the classical nucleation model, the in situ nanocalorimetry studies proved very challenging due to unexpected difficulties in fabricating effective calorimeters and samples. Thus, we could not experimentally validate our predicted influence of composition gradients on nucleation. Nonetheless, important insights were gained and modifications to the classical nucleation theory were suggested.

36 MATERIALS SCIENCE↗
NASA NTRS2020

Ultrafast Melting at Metal-Semiconductor Interfaces: A Molecular Dynamics Study

Quantifying the thermodynamics of interfacial melting is important for several applications such as laser machining and designing radiation hardened electronic devices. Due to the short time and length scales that are involved in the process, molecular simulations provide an ideal path for understating the interfacial melting phenomena. In this work we investigate tungsten (W) – silicon carbide (SiC) interfaces which are relevant in modern high-power electronic devices. Classical molecular dynamics simulations along with the two-temperature model (TTM) were used to study the radiation induced melting of the interface while a neural network-based machine learning framework was used to track the transient re-crystallization process. The damage endured by the semiconductor at the interface following the radiation event is discussed in detail while the re-crystallization of the metal is described in terms of the classical nucleation theory.

Ashwin Ravichandran↗
DOE OSTI2022

Trajectory Ensemble Methods Provide Single-Molecule Statistics for Quantum Dynamical Systems

The emergence of experiments capable of probing quantum dynamics at the single-molecule level requires the development of new theoretical tools capable of simulating and analyzing these dynamics beyond an ensemble-averaged description. In this article, we present an efficient method for sampling and simulating the dynamics of the individual quantum systems that make up an ensemble and apply it to study the nonequilibrium dynamics of the ubiquitous spin-boson model. We generate an ensemble of single-system trajectories, and we analyze this trajectory ensemble using tools from classical statistical mechanics. Our results demonstrate that the dynamics of quantum coherence is highly heterogeneous at the single-system level due to variations in the initial bath configuration, which significantly affects the transient exchange of coherence between the system and its bath. Here, we observe that single systems tend to retain coherence over time scales longer than that of the ensemble. We also compute a novel thermodynamic entanglement entropy that quantifies a thermodynamic driving force favoring system–bath entanglement.

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

Liquid state properties of SEI components in dimethoxyethane

The solid-electrolyte interphase (SEI) layer is a critical constituent of battery technology, which incorporates the use of lithium metals. Since the formation of the SEI is difficult to avoid, the engineering and harnessing of the SEI are absolutely critical to advancing energy storage. One problem is that much fundamental information about SEI properties is lacking due to the difficulty in probing a chemically complex interfacial system. One such property that is currently unknown is the dissolution of the SEI. This process can have significant effects on the stability of the SEI, which is critical to battery performance but is difficult to probe experimentally. Here, we report the use of ab initio computational chemistry simulations to probe the solution state properties of SEI components LiF, Li 2 O, LiOH, and Li 2 CO 3 in order to study their dissolution and other solution-based characteristics. Ab initio molecular dynamics was used to study the solvation structures of the SEI with a combination of radial distribution functions, discrete solvation structure maps, and vibrational density of states, which allows for the determination of free energies. From the change in free energy of dissolution, we determined that LiOH is the most likely component to dissolve in the electrolyte followed by LiF, Li 2 CO 3 , and Li 2 O although none were favored thermodynamically. Furthermore, this indicates that dissolution is not probable, but Li 2 O would make the most stable SEI with regard to dissolution in the electrolyte.

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

Potentials of mean force fail to describe chemical bond-breaking in solution

Many liquid phase studies assume that the potential energy surfaces of reacting molecules are the same as in the gas phase, neglecting complex solvent dynamics that can completely alter the nature of chemical reactivity. Even studies that include solvent effects typically only consider them in an average, equilibrium way as part of a potential of mean force (PMF). In this work, we use mixed quantum/classical simulations to compare how equilibrium and non-equilibrium solvent motions affect the photodissociation of a simple diatomic molecule, NaK + , in liquid tetrahydrofuran. A PMF analysis shows that as the excited-state molecule dissociates with the solvent at equilibrium, the bonding electron remains associated with K + at short bond distances but eventually localizes on Na + at the end of dissociation. When we examine non-equilibrium dynamical photodissociation trajectories, however, we find that they fall into three distinct categories: about a quarter of them have the bonding electron mainly associated with Na + , another quarter stay mainly associated with K + , and about half have the bonding electron shared roughly equally between the two ions. The results show that equilibrium PMFs cannot accurately describe the dynamics of bond-breaking chemical reactions in solution because there is insufficient time for the solvent to reach equilibrium on the time scale over which bond dissociation occurs. Furthermore, our analysis shows that the solvent coupling between the electronic energy surfaces is similar at and away from equilibrium, suggesting that other factors, such as solute velocity-driven solvent memory effects, play a more important role in explaining the failure of the equilibrium PMF to predict the non-equilibrium dynamics.

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

Computing Free Energies with Fluctuation Relations on Quantum Computers

As a central thermodynamic property, free energy enables the calculation of virtually any equilibrium property of a physical system, allowing for the construction of phase diagrams and predictions about transport, chemical reactions, and biological processes. Thus, methods for efficiently computing free energies, which in general is a difficult problem, are of great interest to broad areas of physics and the natural sciences. The majority of techniques for computing free energies target classical systems, leaving the computation of free energies in quantum systems less explored. Recently developed fluctuation relations enable the computation of free energy differences in quantum systems from an ensemble of dynamic simulations. While performing such simulations is exponentially hard on classical computers, quantum computers can efficiently simulate the dynamics of quantum systems. Here, we present an algorithm utilizing a fluctuation relation known as the Jarzynski equality to approximate free energy differences of quantum systems on a quantum computer. In this work we discuss under which conditions our approximation becomes exact, and under which conditions it serves as a strict upper bound. Furthermore, we successfully demonstrate a proof of concept of our algorithm using the transverse field Ising model on a real quantum processor. As quantum hardware continues to improve, we anticipate that our algorithm will enable computation of free energy differences for a wide range of quantum systems useful across the natural sciences.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2003

Small Signal Point-to-Point Tracking of a Propellant Mixer

This paper addresses some theoretical modelling and control issues for a mixing chamber used in rocket engine testing at NASA Stennis Space Center. The mixer is responsible for combining high pressure LH2 and GH2 to produce a hydrogen flow that meets certain thermodynamic properties before it is fed into a test article. The desired properties are maintained by precise control of the LH2 and GH2 flows. The mixer is modelled as a general multi-flow lumped volume for single constituent fluids using density and internal energy as states. The set of nonlinear differential equations is modelled in the SIMULINK environment including a table look-up feature of the fluid thermodynamic properties. a small-signal (linear) model is developed based on the nonlinear model and simulated as well. Pulse disturbances are introduced to the valve positions and the quality of the linear model is ascertained by comparing its behavior against the nonlinear model simulations. Valve control strategies that simulate an operator-in-the-loop scenario are then explored demonstrating the need for automatic feedback control. Finally, classical optimal single-output and multi-output Proportional/Integral controllers are designed based on the linear model and applied to the nonlinear model with excellent results to track simultaneous, constant setpoint changes in desired exit flow, exit temperature, and mixer pressure, as well as to reject unmeasurable but bounded additive step perturbations in the valve positions.

Barbieri, Enrique↗
NASA NTRS2003

Comparisons and Contrasts Between Mediterranean Flood-Producing Orographic Storms and Tropical Cyclones Modified by Complex Orography

Often, late-year, flood-producing storms in the Mediterranean Alps region arise from baroclinic disturbances developing on westerly longwaves, with the incipient storm circulations appearing well west of the flood zone. Regardless of the genesis locations and before the storms begin their flow interactions with the orography of the Mediterranean Alps, a process which can produce intense long-lived rainfall over the mountainous terrain of Italy and France, most of their moisture is drawn from warm Mediterranean waters while important elements of their thermodynamic structure are derived from airmass properties originating in north Africa and the Mediterranean basin itself. In fact, because of strong thermodynamic contrasts between the baroclinic zone producing the incipient disturbance and high potential temperature air from the south advecting into and wrapping around the storm center, a Mediterranean storm whose origins are from an evolving shortwave baroclinic perturbation, can become a quasi-warm core hybrid in term of its energetics composition, a process somewhat akin to the development of polar lows. This process begs the question of whether there are similarities between the autumn-winter season flood producing storms of the Mediterranean Alps and warm core tropical cyclones which become altered dynamically and thermodynamically by landfall over steep orographic terrain. This study investigates, from a comparison and contrast perspective, how the similarities and differences between the two types of storms manifest themselves over their respective life cycles. Before the orographically-induced precipitation-release stage in which both types of storm have evolved to the same hydrometeorlogical conclusion, the two storm categories are transforming from entirely distinctive origins, i.e., the classic eastward traveling, baroclinic mid-latitude cyclone forming in western Europe, and the classic eastward-traveling, warm core tropical cyclone developing in the western Pacific. The analyses use high resolution, nonhydrostatic model simulations, satellite data, and various types of ground data to conduct the study. The flood events of Piemonte- 2000 in northern Italy and Typhoon Rusa-2002 in South Korea are selected to focus the analysis on well documented meteorological events.

Smith, Eric A.↗
DOE OSTI2022

Mechanisms of in-plane magnetic anisotropy in superconducting NbSe 2

Here we present a unifying picture of the magnetic in-plane anisotropies of two-dimensional superconductors based on transition metal dichalcogenides. The symmetry considerations are first applied to constrain the form of the conductivity tensor. We hence conclude that the twofold periodicity of transport distinct from the planar Hall-related contributions requires a tensor perturbation. At the same time, the sixfold periodic variation of the critical field results from the Rashba spin-orbit coupling on a hexagonal lattice. We have considered the effect of a weak tensor perturbation on the critical field, gap function, and magnetoconductivity. The latter is studied using the time-dependent Ginzburg-Landau phenomenology. The common origin of the π periodicity in transport and thermodynamics properties is identified. The scheme constructed here is applied to describe the existing theoretical scenarios from a unified point of view. This allows us to single out the differences and similarities between the suggested approaches.

36 MATERIALS SCIENCE↗
NASA NTRS1973

A general relaxation theory of simple liquids

A relatively simple relaxation theory to account for the behavior of liquids under dynamic conditions was proposed. The general dynamical equations are similar in form to the phenomenological relaxation equations used in theories of viscoelasticity, however, they differ in that all the coefficients of the present equations are expressed in terms of thermodynamic and molecular quantities. The theory is based on the concept that flow in a liquid distorts both the radial and the velocity distribution functions, and that relaxation equations describing the return of these functions to their isotropic distributions, characterizing a stationary liquid, can be written. The theory was applied to the problems of steady and oscillatory shear flows and to the propagation of longitudinal waves. In all cases classical results are predicted for strain rates, and an expression for the viscosity of a liquid, simular to the Macedo-Litovitz equation, is obtained.

Merilo, M.↗