Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermodynamic entropy”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Overcoming time and complexity limitations in molecular dynamics investigations of equilibrium melting

Abstract A hybrid Monte-Carlo molecular-dynamics method for determining solidus and liquidus compositions in multicomponent systems is presented that overcomes both the time limitations in conventional molecular dynamics that prevent the evolution of distinct solid and liquid compositions via diffusion and the complexity challenge that prevents use of thermodynamic assessment in systems of many components. This hybrid method is validated in the Cu–Ni system against an independent assessment of solidus and liquidus compositions based on the regular solution model. Strategies for efficient mapping of different phase diagrams, based on the thermodynamic parameter T 0 , the temperature at which two phases of the composition X 0 have equal free energies, are presented and then demonstrated for the copper-nickel fully miscible system and the gold–silicon eutectic system. A calculation of the solidus and liquidus sampled during the equilibrium melting of equiatomic CrMnFeCoNi is performed, indicating that this method has potential to be extended to the study of many component alloys.

Au-Si↗
DOE OSTI2022

Ab initio free energies of liquid metal alloys: Application to the phase diagrams of Li-Na and K-Na

Comparison of free energies between different phases and different compositions underlies the prediction of alloy phase diagrams. To allow direct comparison, consistent reference points for the energies or enthalpies are required, and the entropy must be placed on an absolute scale, yielding absolute free energies. Here we derive absolute free energies of liquids from ab-initio molecular dynamics (AIMD) by combining the directly simulated enthalpies with an entropy derived from simulated densities and pair correlation functions. Additionally, as an example of the power of this method we calculate the phase diagrams of two binary alkali metal alloys, Li-Na and K-Na, revealing a critical point and liquid-liquid phase separation in the former case, and a deep eutectic in the latter. Good agreement with experimental data demonstrates the power of this simple method.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Inverse design of equilibrium cluster fluids applied to a physically informed model

Inverse design strategies have proven highly useful for the discovery of interaction potentials that prompt self-assembly of a variety of interesting structures. However, often the optimized particle interactions do not have a direct relationship to experimental systems. In this work, we show that Relative Entropy minimization is able to discover physically meaningful parameter sets for a model interaction built from depletion attraction and electrostatic repulsion that yield self-assembly of size-specific clusters. Here, we then explore the sensitivity of the optimized interaction potentials with respect to deviations in the underlying physical quantities, showing that clustering behavior is largely preserved even as the optimized parameters are perturbed.

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

Geometric and thermodynamic characterization of binary neutron star accretion discs

Accretion disks formed in binary neutron star mergers play a central role in many astrophysical processes of interest, including the launching of relativistic jets or the ejection of neutron-rich matter hosting heavy element nucleosynthesis. In this work we analyze in detail the properties of accretion disks from 44 ab initio binary neutron star merger simulations for a large set of nuclear equations of state, binary mass ratios, and remnant fates, with the aim of furnishing reliable initial conditions for disk simulations and a comprehensive characterization of their properties. We find that the disks have a significant thermal support, with an aspect ratio decreasing with the mass ratio of the binary from ~0.7 to 0.3. Even if the disk sample spans a broad range in mass and angular momentum, their ratio is independent from the equation of state and from the mass ratio. This can be traced back to the rotational profile of the disk, characterized by a constant specific angular momentum (as opposed to a Keplerian one) of 3 – 5×10 16 cm 2 s –1 . The profiles of the entropy per baryon and of the electron fraction depend on the mass ratio of the binary. For more symmetric binaries, they follow a sigmoidal distribution as a function of the rest mass density, for which we provide a detailed description and a fit. Furthermore, the disk properties discussed in this work can be used as a robust set of initial conditions for future long-term simulations of accretion disks from binary neutron star mergers, posing the basis for progress in the quantitative study of the outflow properties.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2001

Thermodynamic Data for Fifty Reference Elements

This report is a compilation of thermodynamic functions of 50 elements in their reference state. The functions are: C(sup 0, sub p), {H(T)-H(sup 0)(0)}, S(sup 0)(T), and - {G(sup 0)(T) - H(sup 0)(0)} for the elements Ag, Al, Ar, B, Ba, Be, Br2, C, Ca, Cd, Cl2, Co, Cr, Cs, Cu, F2, Fe, Ge, H2, He, Hg, I2, K, Kr, Li, Mg, Mn, Mo, N2, Na, Nb, Ne, Ni, O2, P, Pb, Rb, S, Si, Sn, Sr, Th, Th, Ti, U, V, W, Xe, Zn, and Zr. Deuterium D, and electron gas e(sup -) are also included. The data are tabulated as functions of temperature as well as given in the form of least-squares coefficients for two functional forms for C(sup 0, sub p) with integration constants for enthalpy and entropy. One functional form for C(sup 0, sub p) is a fourth-order polynomial and the other has two additional terms, one with T(sup -1) and the other with T(sup -2). The gases Ar, D2, e(sup -), H2, He, Kr, N2, Ne, O2, and Xe are tabulated for temperatures from 100 to 20 000 K. The remaining gases Cl2 and F2 are tabulated from 100 to 6000 K. The polynomial functional form for C(sup 0, sub p) for all these gases is split into two temperature intervals of 200 to 1000 K and 1000 to 6000 K. The second functional form for (sup 0, sub p) has an additional interval from 6000 to 20 000 K for the gases tabulated to 20 000 K. The fits are constrained so that the properties match at the common temperature endpoints. The temperature ranges for the condensed species vary with range of the data, phase changes, and shapes of the C(sup 0, sub p) curves.

McBride, Bonnie J.↗
DOE OSTI2025

Experimental evidence of disordered crystalline premixing in sputter-deposited Ni(V)/Al multilayers

The sputter deposition of alternating layers of Ni(V) and Al forms a reactive multilayer known to undergo self-propagating formation reactions when ignited. The sequential deposition process leads to nanometer-scale premixing of reactants at each included interface, which ultimately affects multilayer exothermicity. This work performs the direct measurement of a disordered face-centered cubic (FCC) solid solution premixed phase at the interfaces of Ni(V)/Al multilayers via scanning transmission electron microscopy. The crystallinity of the observed phase differs from previously reported a priori predictions of an amorphous interlayer. The disordered FCC phase retains its symmetry after annealing for 16 h at 135 ± 5 °C, but the lattice parameter shifts consistently with an Al-rich composition. The existence of a crystalline premix in Ni(V)/Al is attributed to the electronic contribution to the entropy of crystallization. The importance of electronic entropy to the phase formation of energetic materials motivates its inclusion when constructing digital twins for atomistic kinetics and ignition sensitivity.

Crystal structure↗
DOE OSTI2025

How Thermodynamic, Electronic, and Steric Factors Influence Mesitylcopper Oligomers

Mesitylcopper (CuMes) is a highly versatile organocopper reagent used in both organic and inorganic syntheses. It has previously been shown that CuMes exists as a tetrameric or pentameric cyclic oligomer [CuMes] n (n = 4, 5), both in solution and in the solid state. The bonding arrangement between the [CuMes] units has qualitatively been described as localized three-center twoelectron (3c-2e) bonds. However, the electronic, structural, and thermodynamic forces driving this aggregation are still not well understood. For this reason, we employed density functional theory (DFT) calculations to study mesitylcopper as a monomeric [CuMes] unit and [CuMes]n oligomers with n = 2 to n = 7. We found that there is a strong electronic driving force for aggregation caused by strong mixing between the Cu’s d orbitals and Mes’s π orbitals in oligomers larger than the dimer. This mixing is only optimized in oligomers with n ≥ 3, where the mesityl group is no longer bonded to a single copper center but instead becomes a bridging ligand. Beyond the trimer, steric and entropic factors become relevant for determining the relative stabilities of the different aggregates, with midsized oligomers (n = 4−5) having the optimal balance between the electronic Cu−C bonding character, Cu···Cu attractive forces, entropy, reduced internal ring strain, and reduced steric interactions between the mesityl groups.

Copper↗
NASA NTRS1993

Thermodynamic data for fifty reference elements

This report is a compilation of thermodynamic functions of 50 elements in their standard reference state. The functions are C(sub p)(sup 0), (H(sup 0)(T) - H(sup 0)(0)), S(sup 0)(T), and -(G(sup 0)(T) - H(sup 0)(O)) for the elements Ag, Al, Ar, B, Ba, Be, Br2, C, Ca, Cd, Cl2, Co, Cr, Cs, Cu, F2, Fe, Ge, H2, He, Hg, I2, K, Kr, Li, Mg, Mn, Mo, N2, Na, Nb, Ne, Ni, O2, P, Pb, Rb, S, Si, Sn, Sr, Ta, Th, Ti, U, V, W, Xe, Zn, and Zr. Deuterium D2 and electron gas e(sup -) are also included. The data are tabulated as functions of temperature as well as given in the form of least-squares coefficients for two functional forms for C(sub p)(sup 0) with integration constants for enthalpy and entropy. One functional form for C(sub p)(sup 0) is a fourth-order polynomial and the other has two additional terms, one with T(exp -1) and the other with T(exp -2). The gases Ar, D2, e(sup -), H2, He, Kr, N2, Ne, O2, and Xe are tabulated for temperatures from 100 to 20,000 K. The remaining gases Cl2 and F2 are tabulated from 100 to 6000 K and 1000 to 6000 K. The second functional form for C(sub p)(sup 0) has an additional interval from 6000 to 20,000 K for the gases tabulated to 20,000 K. The fits are constrained so that the match at the common temperature endpoints. The temperature ranges for the condensed species vary with range of the data, phase changes, and shapes of the C(sub p)(sup 0) curves.

Mcbride, Bonnie J.↗
DOE OSTI2024

Phase-space entropy cascade and irreversibility of stochastic heating in nearly collisionless plasma turbulence

We consider a nearly collisionless plasma consisting of a species of “test particles” in one spatial and one velocity dimension, stirred by an externally imposed stochastic electric field—a kinetic analog of the Kraichnan model of passive advection. The mean effect on the particle distribution function is turbulent diffusion in velocity space—known as stochastic heating. Accompanying this heating is the generation of fine-scale structure in the distribution function, which we characterize with the collisionless (Casimir) invariant C 2 ∝ ∫ ∫ d x d v 〈 f 2 〉 —a quantity that here plays the role of (negative) entropy of the distribution function. We find that C 2 is transferred from large scales to small scales in both position and velocity space via a phase-space cascade enabled by both particle streaming and nonlinear interactions between particles and the stochastic electric field. We compute the steady-state fluxes and spectrum of C 2 in Fourier space, with k and s denoting spatial and velocity wave numbers, respectively. In our model, the nonlinearity in the evolution equation for the spectrum turns into a fractional Laplacian operator in k space, leading to anomalous diffusion. Whereas even the linear phase mixing alone would lead to a constant flux of C 2 to high s (towards the collisional dissipation range) at every k , the nonlinearity accelerates this cascade by intertwining velocity and position space so that the flux of C 2 is to both high k and high s simultaneously. Integrating over velocity (spatial) wave numbers, the k -space ( s -space) flux of C 2 is constant down to a dissipation length (velocity) scale that tends to zero as the collision frequency does, even though the rate of collisional dissipation remains finite. The resulting spectrum in the inertial range is a self-similar function in the ( k , s ) plane, with power-law asymptotics at large k and s . Our model is fully analytically solvable, but the asymptotic scalings of the spectrum can also be found via a simple phenomenological theory whose key assumption is that the cascade is governed by a “critical balance” in phase space between the linear and nonlinear timescales. We argue that stochastic heating is made irreversible by this entropy cascade and that, while collisional dissipation accessed via phase mixing occurs only at small spatial scales rather than at every scale as it would in a linear system, the cascade makes phase mixing even more effective overall in the nonlinear regime than in the linear one. Published by the American Physical Society 2024

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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↗
DOE OSTI2021

Enhancement and maximum in the isobaric specific-heat capacity measurements of deeply supercooled water using ultrafast calorimetry

Knowledge of the temperature dependence of the isobaric specific heat (C p ) upon deep supercooling can give insights regarding the anomalous properties of water. If a maximum in C p exists at a specific temperature, as in the isothermal compressibility, it would further validate the liquid–liquid critical point model that can explain the anomalous increase in thermodynamic response functions. The challenge is that the relevant temperature range falls in the region where ice crystallization becomes rapid, which has previously excluded experiments. Here, we have utilized a methodology of ultrafast calorimetry by determining the temperature jump from femtosecond X-ray pulses after heating with an infrared laser pulse and with a sufficiently long time delay between the pulses to allow measurements at constant pressure. Evaporative cooling of ~15-µm diameter droplets in vacuum enabled us to reach a temperature down to ~228 K with a small fraction of the droplets remaining unfrozen. We observed a sharp increase in C p , from 88 J/mol/K at 244 K to about 218 J/mol/K at 229 K where a maximum is seen. The C p maximum is at a similar temperature as the maxima of the isothermal compressibility and correlation length. From the C p measurement, we estimated the excess entropy and self-diffusion coefficient of water and these properties decrease rapidly below 235 K.

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

Critical phenomena of the layered ferrimagnet Mn 3 Si 2 Te 6 following proton irradiation

The critical phenomena and magnetic entropy of the quasi-2D ferrimagnetic crystal, Mn 3 Si 2 Te 6 (MST), is analyzed along the easy axis (H || ab) as a function of proton irradiance. The critical exponents β and γ do not fall into any particular universality class upon proton irradiation. However, for pristine and irradiated samples, the critical exponents lie closer to mean field-like interactions; therefore, long-range interactions are presumed to be sustained in MST. The effective spatial dimensionality reveals that MST remains at d =3 under proton irradiation, whereas spin dimensionality transitions from an initial n =1 to n =2 and n =3 for 1 × 10 15 and 5 × 10 15 H + /cm 2 , indicating XY and Heisenberg interactions, respectively. The spin correlation function reveals an increase in magnetic correlations at 5 × 10 15 H + /cm 2 . Maximum change in magnetic entropy at 3 T is the largest for 5 × 10 15 H + /cm 2 at 2.45 J/kg K, in comparison to 1.60 J/kg K for pristine MST. These results intriguingly align with previous findings on MST where magnetization increased by ~50% at 5 × 10 15 H + /cm 2 , in comparison to its pristine counterpart [Martinez et al., Appl. Phys. Lett. 116, 172404 (2020)]. Magnetic entropy derived from heat capacity analysis shows no large deviations across the proton irradiated samples suggesting that the antiferromagnetic (AFM) coupling between the Mn sites is stable even after proton irradiation. This implies that magnetization is enhanced through a strengthening of the super-exchange interaction between Mn atoms mediated through Te rather than a weakening of the AFM component.

2D materials↗
NASA NTRS2004

Thermodynamic Conditions Favorable to Superlative Thunderstorm Updraft, Mixed Phase Microphysics and Lightning Flash Rate

Satellite observations of lightning flash rate have been merged with proximal surface station thermodynamic observations toward improving the understanding of the response of the updraft and lightning activity in the tropical atmosphere to temperature. The tropical results have led in turn to an examination of thermodynamic climatology over the continental United States in summertime and its comparison with exceptional electrical conditions documented in earlier studies. The tropical and mid-latitude results taken together support an important role for cloud base height in regulating the transfer of Convective Available Potential Energy (CAPE) to updraft kinetic energy in thunderstorms. In the tropics, cloud base height is dominated by the dry bulb temperature over the wet bulb temperature as the lightning-regulating temperature in regions characterized by moist convection. In the extratropics, an elevated cloud base height may enable larger cloud water concentrations in the mixed phase region, a favorable condition for the positive charging of large ice particles that may result in thunderclouds with a reversed polarity of the main cloud dipole. The combined requirements of instability and cloud base height serve to confine the region of superlative electrification to the vicinity of the ridge in moist entropy in the western Great Plains.

Williams, E.↗
DOE OSTI2023

Borderline first-order phase transition and large cryogenic magnetocaloric effect in PrNdIn

We report a large cryogenic magnetocaloric effect stemming from an unconventional borderline first-order magnetic phase transition with negligibly small thermomagnetic hysteresis in a rare-earth-based intermetallic compound PrNdIn. The sample exhibits maximum magnetic field-induced entropy change as large as –10 J/Kg K (for 20 kOe magnetic field change) near the boiling point of oxygen. Magnetocaloric properties of PrNdIn are comparable to those of other known potential magnetocaloric materials with operating temperatures ranging between 50 and 125 K. Here, the magnetic properties of the present sample are qualitatively reminiscent of those of the binary Pr 2 In and Nd 2 In, including the emergence of a second low-temperature anomaly in the temperature dependence of magnetization.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Enhancing the formation of ionic defects to study the ice Ih/XI transition with molecular dynamics simulations

Ice Ih, the common form of ice in the biosphere, contains proton disorder. Its proton-ordered counterpart, ice XI, is thermodynamically stable below 72 K. However, the formation of ice XI is kinetically hindered, and experimentally it is obtained by doping with KOH. Doping creates ionic defects that promote the migration of protons and the associated change in proton configuration. In this article, we mimic the effect of doping with a bias potential that enhances the formation of ionic defects in molecular dynamics simulations. The recombination of the ions thus formed proceeds through fast migration of the hydroxide along hydrogen bond loops, providing a physical and expedite way to change the proton configuration. Here, a key ingredient of this approach is a machine learning potential trained with density functional theory data and capable of modelling molecular dissociation. We exemplify the usefulness of this idea by studying the order-disorder transition using an appropriate order parameter that distinguishes the proton environments in ice Ih and XI. We calculate the changes in free energy, enthalpy, and entropy associated with the transition. Our estimated entropy agrees with experiment within the error bars of the calculation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Fluctuations of subsystem entropies at late times

We study the fluctuations of subsystem entropies in closed quantum many-body systems after thermalization. Using a combination of analytics and numerics for both random quantum circuits and Hamiltonian dynamics, we find that the statistics of such entropy fluctuations is drastically different than in the classical setting. For instance, shortly after a system thermalizes, the probability of entropy fluctuations for a subregion is suppressed in the dimension of the Hilbert space of the complementary subregion. This suppression becomes increasingly stringent as a function of time, ultimately depending on the exponential of the Hilbert space dimension, until extremely late times when the amount of suppression saturates. We also use our results to estimate the total number of rare fluctuations at large timescales. We find that the “Boltzmann brain” paradox is largely ameliorated in quantum many-body systems, in contrast with the classical setting.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Magnetoelastic-coupling-assisted first-order antiferromagnetic transitions in the corrugated honeycomb-lattice compounds (Ca, Sr) ⁢Mn 2 ⁢ P 2

Understanding the origin of unconventional magnetic order is a central theme in the research field of functional quantum materials for spin-dependent applications. In general, a paramagnetic (PM) to antiferromagnetic (AFM) transition is of thermodynamically second order in nature in the absence of a magnetic field. In the magnetically frustrated corrugated-honeycomb lattice compounds (Ca, Sr) ⁢Mn 2 ⁢ P 2 , electrical resistivity 𝜌 measurements in the 𝑎⁢𝑏 plane are reported to exhibit sharp discontinuities at the respective PM to AFM ordering temperatures 𝑇 N = 69.5 and 53 K, respectively, signifying the unusual first-order nature of the AFM transition. To reveal the origin of such anomalous transitions, high-resolution linear thermal-expansion measurements are reported for both CaMn 2 ⁢P 2 and SrMn 2 ⁢ P 2 in the 𝑎⁢𝑏 plane and along the 𝑐 axis using capacitance dilatometry. The thermal expansion measurements show a sharp discontinuity in Δ⁢𝐿/𝐿 at 𝑇 N , yielding a diverging behavior in the thermal expansion coefficient 𝛼⁡(𝑇) associated with strong magnetoelastic coupling, driving the first-order AFM transition in CaMn 2 ⁢P 2 . However, this effect is weaker in SrMn 2 ⁢P 2 , consistent with the magnetic entropy changes 𝑆⁡(𝑇) obtained from the heat-capacity measurements. Temperature-dependent 𝑎⁢𝑏-plane x-ray diffraction measurements for CaMn 2 ⁢ P 2 near its 𝑇 N yield a volume thermal expansion coefficient consistent with literature data for other materials. The uniaxial pressure derivatives of 𝑇 N in the two compounds are found to be respectively opposite in sign, presumably due to the different magnetic structures.

Pakhira, Santanu [Maulana Azad National Institute ↗
DOE OSTI2025

Roadmap on thermodynamics and thermal metamaterials

Thermal metamaterials represent a transformative paradigm in modern physics, synergizing thermodynamic principles with metamaterial engineering to master heat flow at will. As next-generation technologies demand multi-scale thermal control, this field urgently requires systematic frameworks to unify its multidisciplinary advances. Curated through a global collaboration involving over 50 specialists across 25 subdisciplines, this review primarily summarizes two decades of advancements, ranging from theoretical breakthroughs to functional implementations. The review reveals groundbreaking innovations in heat manipulation through the exploration of both classical and non-classical transport regimes, topological thermal control mechanisms, and quantum-informed phonon engineering strategies. By bridging physical insights like non-Hermitian thermal dynamics and valleytronic phonon transport with cutting-edge applications, we demonstrate paradigm-shifting capabilities: environment-adaptive thermal cloaks, AI-optimized metamaterials, and nonlinear thermal circuits enabling heat-based computation. Experimental milestones include 3D thermal null media with reconfigurable invisibility and thermal designs breaking classical conductivity limits. Here, this collaborative effort establishes an indispensable roadmap for physicists, highlighting pathways to quantum thermal management, entropy-controlled energy systems, and topological devices. As thermal metamaterials transition from laboratory marvels to technological cornerstones, this work provides the foundational lexicon and design principles for the coming era of intelligent thermal matter.

heat conduction control↗