Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermal equilibrium”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

The structure of cluster merger shocks: turbulent width and the electron heating time-scale

We present a new 2MsChandra observation of the cluster merger Abell 2146, which hosts two huge M ~ 2 shock fronts each ~500 kpc across. For the first time, we resolve and measure the width of cluster merger shocks. The best-fitting width for the bow shock is 17 ± 1kpc and for the upstream shock is 10.7 ± 0.3 kpc. A narrow collisionless shock will appear broader in projection if its smooth shape is warped by local gas motions. We show that both shock widths are consistent with collisionless shocks blurred by local gas motions of 290 ± 30 km s –1 . The upstream shock forms later on in the merger than the bow shock and is therefore expected to be significantly narrower. From the electron temperature profile behind the bow shock, we measure the time-scale for the electrons and ions to come back into thermal equilibrium. We rule out rapid thermal equilibration of the electrons with the shock-heated ions at the 6σ level. The observed temperature profile instead favours collisional equilibration. For these cluster merger shocks, which have low sonic Mach numbers and propagate through a high β plasma, we find no evidence for electron heating over that produced by adiabatic compression. Our findings are expected to be valid for collisionless shocks with similar parameters in other environments and support the existing picture from the solar wind and supernova remnants. Furthermore, the upstream shock is consistent with this result but has a more complex structure, including a ~2ke V increase in temperature ~50 kpc ahead of the shock.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Dark matter isocurvature from curvature

Isocurvature fluctuations, where the relative number density of particle species spatially varies, can be generated from initially adiabatic, or curvature, fluctuations if the various species fall out of or were never in thermal equilibrium. The freezing of the thermal relic dark matter abundance is one such case, but for modes that are still outside the horizon the amplitude is highly suppressed and originates from the small change in the local expansion rate due to the local space curvature produced by the curvature fluctuation. We establish a simple separate-universe method for calculating this generation that applies to both freeze-in and freeze-out models, identify three critical epochs for this process, and give general scaling behaviors for the amplitude in each case: the freezing epoch, the kinetic decoupling epoch and matter-radiation equality. Freeze-out models are typically dominated by spatially modulated annihilation from the latter epochs and can generate much larger isocurvature fluctuations compared with typical freeze-in models, albeit still very small and observationally allowed by cosmic microwave background measurements. Here, we illustrate these results with concrete models where the dark matter interactions are vector or scalar mediated.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Spectroscopic Determination of Stellar Parameters and Oxygen Abundances for Hyades/Field G–K Dwarfs

It has been suggested that Fe abundances of K dwarfs derived from Fe i and Fe ii lines show considerable discrepancies, and oxygen abundances determined from high-excitation O i 7771–5 triplet lines are appreciably overestimated (the problem becoming more serious toward lower T {sub eff}), which, however, has not yet been widely confirmed. With the aim of clarifying this issue, we spectroscopically determined the atmospheric parameters of 148 G–K dwarfs (Hyades cluster stars and field stars) by assuming the classical Fe i/Fe ii ionization equilibrium as usual, and determined their oxygen abundances by applying the non-local thermal equilibrium spectrum fitting analysis to O i 7771–5 lines. It turned out that the resulting parameters did not show any significant inconsistency with those determined by other methods (for example, the mean differences in T {sub eff} and logg from the well-determined solutions of Hyades dwarfs are mostly ≲100 K and ≲0.1 dex). Likewise, the oxygen abundances of Hyades stars are around [O/H] ∼ +0.2 dex (consistent with the metallicity of this cluster) without exhibiting any systematic T {sub eff}-dependence. Accordingly, we conclude that parameters can be spectroscopically evaluated to a sufficient precision in the conventional manner (based on the Saha–Boltzmann equation for Fe i/Fe ii) and oxygen abundances can be reliably determined from the O i 7771–5 triplet for K dwarfs as far as stars of T {sub eff} ≳ 4500 K are concerned. We suspect that previously reported strongly T {sub eff}-dependent discrepancies may have stemmed mainly from overestimation of weak-line strengths and/or improper T {sub eff} scale.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Simulations of classical three-body thermalization in one dimension

One-dimensional systems, such as nanowires or electrons moving along strong magnetic field lines, have peculiar thermalization physics. The binary collision of pointlike particles, typically the dominant process for reaching thermal equilibrium in higher-dimensional systems, cannot thermalize a 1D system. We study how dilute classical 1D gases thermalize through three-body collisions. We consider a system of identical classical point particles with pairwise repulsive inverse power-law potential V ij ∝ 1/|x i –x j | n or the pairwise Lennard-Jones potential. Using Monte Carlo methods, we compute a collision kernel and use it in the Boltzmann equation to evolve a perturbed thermal state with temperature T toward equilibrium. We explain the shape of the kernel and its dependence on the system parameters. Additionally, we implement molecular dynamics simulations of a many-body gas and show agreement with the Boltzmann evolution in the low-density limit. For the inverse power-law potential, the rate of thermalization is proportional to ρ 2 ⁢T$\frac{1}{2}$ – $\frac{1}{n}$, where ρ is the number density. Furthermore, the corresponding proportionality constant decreases with increasing n.

1-dimensional systems↗
DOE OSTI2021

Disequilibrium Chemistry in Exoplanet Atmospheres Observed with the Hubble Space Telescope

Literature on the theory of exoplanet atmospheric disequilibrium chemistry is rich, although its observational counterpart has yet to emerge beyond the hints provided by a few targets in dedicated studies. We report results from an uniform data reduction and analysis for a catalog of 62 Hubble Space Telescope exoplanet transit spectra where we assess the atmospheric model preference for disequilibrium chemistry (i.e., water vapor is not the dominant absorption spectral signature) over thermal equilibrium chemistry in a comparative planetology context. Where model preference assessment is possible, we find that disequilibrium occurs in about half of the atmospheres, indicating that disequilibrium processes play an important role in the composition of exoplanet atmospheres. While very hot atmospheres, over 1800 K, prefer equilibrium chemistry, we find a clustering of preference for disequilibrium in the 1200–1800 K temperature range. We suggest that UV-augmented thermochemistry may play a significant role for those atmospheres.

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

Time evolution of a pumped molecular magnet—A time-resolved inelastic neutron scattering study

Introducing an experimental technique of time-resolved inelastic neutron scattering (TRINS), we explore the time-dependent effects of resonant pulsed microwaves on the molecular magnet Cr 8 F 8 Piv 16 . The octagonal rings of magnetic Cr 3+ atoms with antiferromagnetic interactions form a singlet ground state with a weakly split triplet of excitations at 0.8 meV. A 4.6 tesla field was applied to tune the splitting between two members of the triplet excited level |1$\rangle$ ↔ |2$\rangle$ to resonance with 105 GHz (0.434 meV) microwaves. The time-dependent occupations of the ground state |0$\rangle$, lower lying levels |1$\rangle$ and |2$\rangle$, and higher energy states |λ ≥ 3$\rangle$ were extracted during and after 20 s long microwave pulses incident along the (101) direction of a Cr 8 F 8 Piv 16 crystal held at 1.9 K. At significantly elevated spin temperatures, we found underpopulation relative to thermal equilibrium of |2$\rangle$ and spin-lattice thermalization time scales ranging from 1.6(2) s to 5.7(2) s depending on the power level. This contrasts with the relaxation time τ 1 (T → 0) = 27(5) 𝛍 s inferred for |2$\rangle$ from in situ Electron Spin Resonance measurements. By probing a broad range of excited states during intense microwave pumping, TRINS thus provides a first view of long lived excited states in a molecular antiferromagnet.

Science & Technology - Other Topics↗
DOE OSTI2021

Astromers in the Radioactive Decay of r-process Nuclei

Certain nuclear isomers are well known to affect nucleosynthesis with important observable consequences (e.g., 26 Al and 180 Ta). We study the impact of nuclear isomers in the context of rapid neutron capture process (r-process) nucleosynthesis. Here we demonstrate that nuclear isomers are dynamically populated in the r process and that some are populated far from thermal equilibrium; this makes them astrophysical isomers, or "astromers." We compute thermally mediated transition rates between long-lived isomers and the corresponding ground states in neutron-rich nuclei. We calculate the temperature-dependent β-decay feeding factors, which represent the fraction of material going to each of the isomer and ground state daughter species from the β-decay parent species. We simulate nucleosynthesis following the decay of a solar-like r-process composition and include as separate species nuclear excited states with measured terrestrial half-lives greater than 100 μs. We introduce a new metric to identify those astromers most likely to be influential and summarize them in a table. Notable entries include many second peak nuclei (e.g., the Te isotopic chain) and previously overlooked isomers in stable nuclei (e.g., 119 Sn, 131 Xe, and 195 Pt). Finally, we comment on the capacity of isomer production to alter radioactive heating in an r-process environment.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

An In-situ and Direct Confirmation of Super-Planckian Thermal Radiation Emitted From a Metallic Photonic-Crystal at Optical Wavelengths

Planck’s law predicts the distribution of radiation energy, color and intensity, emitted from a hot object at thermal equilibrium. The Law also sets the upper limit of radiation intensity, the blackbody limit. Recent experiments reveal that micro-structured tungsten can exhibit significant deviation from the blackbody spectrum. However, whether thermal radiation with weak non-equilibrium pumping can exceed the blackbody limit in the far field remains un-answered experimentally. Here, we compare thermal radiation from a micro-cavity/tungsten photonic crystal (W-PC) and a blackbody, which are both measured from the same sample and also in-situ. We show that thermal radiation can exceed the blackbody limit by >8 times at λ=1.7 μm resonant wavelength in the far-field. Our observation is consistent with a recent calculation by Wang and John performed for a 2D W-PC filament. This finding is attributed to non-equilibrium excitation of localized surface plasmon resonances coupled to nonlinear oscillators and the propagation of the electromagnetic waves through non-linear Bloch waves of the W-PC structure. This discovery could help create super-intense narrow band thermal light sources and even an infrared emitter with a laser-like input-output characteristic.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Microwave-driven heterogeneous catalysis for activation of dinitrogen to ammonia under atmospheric pressure

This paper presents an innovative approach to producing energy-dense, carbon–neutral liquid ammonia as a means for carrying energy. This approach synergistically integrates microwave reaction chemistry with novel heterogeneous catalysis that decouples dinitrogen activation from high-temperature and high-pressure reactions, altering reaction pathways and increasing ammonia formation rate. Results presented here demonstrate that ammonia synthesis can be conducted at 280 ? and ambient pressure to achieve ~1 mmol ammonia/g cat/h over supported ruthenium catalyst systems utilizing microwave irradiation. It is further shown that adding promoter ions such as potassium, cerium, and barium significantly improves the ammonia production rate over undoped ruthenium-based catalysts. This effect could be attributed to enhanced dielectric loss processes that lead to stronger microwave absorption by the catalyst. Measurement of the equilibrium constant under microwave conditions showed a higher ammonia yield than under thermal equilibrium conditions for both the iron- and ruthenium-based catalysts. Finally, this study also illustrates the advantages of using a variable-frequency microwave reactor for ambient-pressure ammonia synthesis. Mechanistically, investigators believed that the oscillating electric fields of the radiation can couple with adsorbed nitrogen on the surface and accelerate its dissociation. Since dinitrogen dissociation on the surface is rate limiting, this effectively accelerates the reaction.

Hu, Jianli (John)↗
DOE OSTI2021

Interdiffusion of Elements During Ultrasonic Additive Manufacturing

In this paper, we report evidence for enhanced elemental interdiffusion during ultrasonic additive manufacturing (UAM) across metal boundaries of copper-aluminum, nickel-gold, and nickel-gold-aluminum. The high solute interdiffusion measured by energy dispersive X-ray spectroscopy line scans is rationalized with calculated vacancy concentrations orders of magnitude larger than thermal equilibrium values. The above estimates are supported by existing knowledge related to defect physics and UAM thermal cycles. The observation of pronounced elemental mixing are evidence for the presence of enhanced non-equilibrium immiscible metal interdiffusion during UAM processing.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A simple and effective five-equation two-phase numerical model for liquid-vapor phase transition in cavitating flows

Numerical difficulties, notably the non-monotonic behavior of the Wood speed of sound and the volume fraction positivity, associated with the reduced five-equation two-phase flow model of Kapila et al. (2001) [A.K. Kapila, R. Menikoff, J.B. Bdzil, S.F. Son, D.S. Stewart, 2001. Two-phase modeling of deflagration-to-detonation transition in granular materials: reduced equations, Physics of Fluids 13(10), 3002–3024] have been resolved in the past through the introduction of a frozen speed of sound and an algebraic approach for mechanical relaxation afforded by a pressure non-equilibrium six-equation model proposed by [R. Saurel, F. Petitpas, R.A. Berry, 2009. Simple and efficient relaxation methods for interfaces separating compressible fluids, cavitating flows and shocks in multiphase mixture, J. Comput. Phys. 228, 1678–1712]. By contrast, it is explored and demonstrated in this work that these difficulties can in fact still be resolved within the numerical scheme for solving the reduced five-equation model by numerically replacing the Wood speed of sound for the estimates of wave speeds in the approximate Riemann solver HLLC with the monotonic mixture speed of sound for a transport five-equation model. For shock interface (artificial mixture separating pure or nearly pure fluids) interaction problems, with the apparent advantage of monotonic behavior of the speed of sound in the interface, the effect of the numerical replacement is also confined to the interface. Differences other than the behavior of the speed of sound within the interface in the solutions due to the replacement diminish with increasing resolution when reasonable solution can be obtained with Wood speed of sound. For cavitating/expansion problems in physical fluid mixture, it is pointed out and explained why acoustics in the numerical solutions still propagate at the Wood speed of sound (therefore consistent with the reduced five-equation model) even though in some cases a much higher speed of sound like the numerical replacement above for solving the reduced five-equation model or the frozen speed of sound for solving a six-equation model is used for the estimates of wave speeds in the HLLC scheme. A variant of the five-equation two-phase flow model by Saurel et al. (2008) [R. Saurel, F. Petitpas, R. Abgrall, 2008. Modelling phase transition in metastable liquids: application to cavitating and flashing flows, J. Fluid Mech. 607, 313–350] is then constructed for liquid-vapor phase transition in cavitating flows. The relaxation toward thermo-chemical equilibrium during phase transition is achieved by solving a simple system of algebraic equations for the equilibrium state variables for better efficiency, following Pelanti and Shyue (2014) [M. Pelanti, K.-M. Shyue, 2014. A mixture-energy-consistent six-equation two-phase numerical model for fluid with interfaces, cavitation and evaporation waves. J. Comput. Phys. 259, 331–357]. Therefore, the current model retains both the simplicity afforded by the five-equation model and the efficiency of the algebraic relaxation solver. An alternative algebraic approach for handling the non-conservative term (the so-called K∇ · u term) in the reduced five-equation model for mechanical equilibrium of a liquid-vapor mixture is also explored by enforcing the thermal equilibrium at the same time. Finally, numerical results of sample tests in both one and two dimensions in the literature as well as that in three dimensions demonstrate the effectiveness and ability of the proposed model to simulate cavitating flows. An interesting mechanism of shock generation by acoustics in water due to phase transition is then found by the numerical simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Transverse momentum broadening of a jet in quark-gluon plasma: an open quantum system EFT

We utilize the technology of open quantum systems in conjunction with the recently developed effective field theory for forward scattering to address the question of massless jet propagation through a weakly-coupled quark-gluon plasma in thermal equilibrium. We discuss various possible hierarchies of scales that may appear in this problem, by comparing thermal scales of the plasma with relevant scales in the effective field theory. Starting from the Lindblad equation, we derive and solve a master equation for the trans- verse momentum distribution of a massless quark jet, at leading orders both in the strong coupling and in the power counting of the effective field theory. Markovian approximation is justified in the weak coupling limit. Using the solution to the master equation, we study the transverse momentum broadening of a jet as a function of the plasma temperature and the time of propagation. We discuss the physical origin of infrared sensitivity that arises in the solution and a way to handle it in the effective field theory formulation. We suspect that the final measurement constraint can only cut-off leading infrared singularities and the solution to the Markovian master equation resums a logarithmic series. This work is a stepping stone towards understanding jet quenching and jet substructure observables on both light and heavy quark jets as probes of the quark-gluon plasma.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Equilibrium spin polarization arising from chirality

Chirality-induced spin selectivity (CISS) describes how chiral molecules and materials generate spin polarization even at thermal equilibrium. This observation has challenged established principles of microscopic reversibility and Onsager reciprocity. We resolve this paradox by formulating a pseudo-Hermitian quantum framework that separates thermodynamic equilibrium from time-reversal symmetry. Within this approach, structural chirality and electron correlations, irrespective of their microscopic origin, are sufficient to produce CISS observables. Chirality enters through a non-local metric η that couples spin and spatial motion, leading to real spectra, unitary evolution, and thermodynamic consistency. The framework predicts a chirality-induced spin magnetic ordering characterized by a spin-displacement order, which reconciles equilibrium spin polarization with detailed balance and explains the persistence of CISS in materials composed of light elements. We derive generalized Onsager-Casimir relations that respect the observed CISS symmetry, i.e., parity $(\mathscr{P})$-odd and time-reversal $(\mathscr{T})$-odd, but exhibiting $\mathscr{PT}$-even symmetry. This approach establishes a coherent foundation for equilibrium CISS and provides a route to link chemical chirality with measurable spin-to-charge conversion effects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

The Lifshitz Regime and its Experimental Signals

We discuss the possibility of a Lifshitz regime, where the dispersion relation for Goldstone bosons and related fields has a minimum at nonzero momenta. Studies with the Functional Renormalization Group suggest that this occurs over a wide region in the plane of temperature and baryon chemical potential. Conversely, the FRG finds that the region in which fluctuations from a critical endpoint are significant is rather small. Here, we suggest that this is due generically to the narrowness of the tricritical region in the chiral limit. Even if particles are produced in thermal equilibrium, a dispersion relation which is non-monotonic in momenta produces what appears to be non-thermal behavior.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Quantum algorithms from fluctuation theorems: Thermal-state preparation

Fluctuation theorems provide a correspondence between properties of quantum systems in thermal equilibrium and a work distribution arising in a non-equilibrium process that connects two quantum systems with Hamiltonians H 0 and H 1 = H 0 + V. Building upon these theorems, we present a quantum algorithm to prepare a purification of the thermal state of H 1 at inverse temperature β ≥ 0 starting from a purification of the thermal state of H 0 . The complexity of the quantum algorithm, given by the number of uses of certain unitaries, is $\mathcal{O}$ (e β(ΔA - w l )/2 ), where ΔA is the free-energy difference between H 1 and H 0 , and w l is a work cutoff that depends on the properties of the work distribution and the approximation error ϵ > 0. If the non-equilibrium process is trivial, this complexity is exponential i β∥V∥, where ∥V∥ is the spectral norm of V. This represents a significant improvement of prior quantum algorithms that have complexity exponential in β∥H 1 ∥ in the regime where ∥V∥$\ll$ ∥H 1 ∥. The dependence of the complexity in ϵ varies according to the structure of the quantum systems. It can be exponential in 1/ϵ in general, but we show it to be sublinear in 1/ϵ if H 0 and H 1 commute, or polynomial in 1/ϵ if H 0 and H 1 are local spin systems. The possibility of applying a unitary that drives the system out of equilibrium allows one to increase the value of w l and improve the complexity even further. To this end, we analyze the complexity for preparing the thermal state of the transverse field Ising model using different non-equilibrium unitary processes and see significant complexity improvements.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Solid-Density Ion Temperature from Redshifted and Double-Peaked Stark Line Shapes

Heβ spectral line shapes are important for diagnosing temperature and density in many dense plasmas. This work presents Heβ line shapes measured with high spectral resolution from solid-density plasmas with minimized gradients. The line shapes show hallmark features of Stark broadening, including quantifiable redshifts and double-peaked structure with a significant dip between the peaks; these features are compared to models through a Markov chain Monte Carlo framework. Line shape theory using the dipole approximation can fit the width and peak separation of measured line shapes, but it cannot resolve an ambiguity between electron density n e and ion temperature T i , since both parameters influence the strength of quasistatic ion microfields. Here a line shape model employing a full Coulomb interaction for the electron broadening computes self-consistent line widths and redshifts through the monopole term; redshifts have different dependence on plasma parameters and thus resolve the n e – T i ambiguity. The measured line shapes indicate densities that are 80–100% of solid, identifying a regime of highly ionized but well-tamped plasma. This analysis also provides the first strong evidence that dense ions and electrons are not in thermal equilibrium, despite equilibration times much shorter than the duration of x-ray emission; cooler ions may arise from nonclassical thermalization rates or anomalous energy transport. Furthermore, the experimental platform and diagnostic technique constitute a promising new approach for studying ion-electron equilibration in dense plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

BBN photodisintegration constraints on gravitationally produced vector bosons

Gravitational production of massive particles due to cosmic expansion can be significant during the inflationary and reheating period of the Universe. If the particle also has non-gravitational interactions that do not significantly affect its production, numerous observational probes open up, including cosmological probes. In this work, we focus on the gravitational production of light vector bosons that couple feebly to the Standard Model (SM) particles. Due to the very feeble coupling, the light vector bosons never reach thermal equilibrium, and if the Hubble scale at the end of inflation is above 10 8 GeV, the gravitational production can overwhelm the thermal production via the freeze-in mechanism by many orders of magnitude. As a result, much stronger constraints from the Big Bang Nucleosynthesis (BBN) can be placed on the lifetime and mass of the vector bosons compared to the scenario where only thermal production is considered. As an example, we study the sub-GeV scale dark photons, which couple to the SM only through kinetic mixing, and derive constraints on the mass and kinetic mixing parameter of the dark photon from the photodisintegration effects on the light element abundances relevant at the end of the BBN when the cosmic age was around 10 4 s.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Evidence for equilibrium exciton condensation in monolayer WTe 2

We present evidence that the two-dimensional bulk of monolayer WTe2 contains electrons and holes bound by Coulomb attraction—excitons—that spontaneously form in thermal equilibrium. On cooling from room temperature to 100 K, the conductivity develops a V-shaped dependence on electrostatic doping, while the chemical potential develops a step at the neutral point. These features are much sharper than is possible in an independent-electron picture, but they can be accounted for if electrons and holes interact strongly and are paired in equilibrium. Our calculations from first principles show that the exciton binding energy is larger than 100 meV and the radius as small as 4 nm, explaining their formation at high temperature and doping levels. Below 100 K, more strongly insulating behaviour is seen, suggesting that a charge-ordered state forms. The observed absence of charge density waves in this state is surprising within an excitonic insulator picture, but we show that it can be explained by the symmetries of the exciton wavefunction. Therefore, in addition to being a topological insulator, monolayer WTe2 exhibits strong correlations over a wide temperature range.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗