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

The saturation mechanism of thermal instability

The literature on thermal instability (TI) reveals that even for a simple homogeneous plasma, the nonlinear outcome ranges from a gentle reconfiguration of the initial state to an explosive one, depending on whether the condensations that form evolve in an isobaric or nonisobaric manner. After summarizing the recent developments on the linear and nonlinear theory of TI, here we derive several general identities from the evolution equation for entropy that reveal the mechanism by which TI saturates; whenever the boundary of the instability region (the Balbus contour) is crossed, a dynamical change is triggered that causes the comoving time derivative of the pressure to change the sign. This event implies that the gas pressure force reverses direction, slowing the continued growth of condensation. For isobaric evolution, this “pressure reversal” occurs nearly simultaneously for every fluid element in condensation and a steady state is quickly reached. For nonisobaric evolution, the condensation is no longer in mechanical equilibrium and the contracting gas rebounds with greater force during the expansion phase that accompanies the gas reaching the equilibrium curve. The cloud then pulsates because the return to mechanical equilibrium becomes wave mediated. We show that both the contraction rebound event and subsequent pulsation behavior follow analytically from an analysis of the new identities. Our analysis also leads to the identification of an isochoric TI zone and makes it clear that unless this zone intersects the equilibrium curve, isochoric modes can only become unstable if the plasma is in a state of thermal non-equilibrium.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Assessment of SCALE and MELCOR for a generic pebble bed fluoride high-temperature reactor

This paper presents the development and application of SCALE and MELCOR models for a fluoride salt-cooled high-temperature reactor (FHR) based on publicly available specifications. SCALE version 6.3beta15 was used to generate power distributions and decay heat curves, and MELCOR version 2.2.18019 was used to calculate the thermal hydraulic response of an assumed FHR primary system. Here, an approach was developed for determining the equilibrium state of the core using depletion of a core slice model and blending of fuel compositions at different burnups to provide three-dimensional fuel composition in the core. Results compared between a core with entirely fresh fuel and one with an equilibrium fuel composition revealed that the equilibrium core led to lower steady-state temperatures but slower cooldown during a loss of flow accident (LOFA). A sensitivity study was then conducted to explore the transient response of our FHR system to variations in thermal hydraulic parameters of the system using the equilibrium core. The inlet temperature, graphite thermal conductivity, and SCRAM time all provided significant control over peak fuel and coolant temperatures. Uncertainties in radionuclide decay data in SCALE were used to perform a decay heat sensitivity study, and we found that uncertainties in decay heat led to negligible impact on peak temperatures during the course of the transient. In all cases evaluated, the observed peak fuel temperatures remained approximately 700 K below anticipated failure limits for the LOFA.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Thermal conductivity tensor of β -1,3,5,7-tetranitro-1,3,5,7-tetrazoctane ( β -HMX) as a function of pressure and temperature

Here, we have used reverse non-equilibrium molecular dynamics (RNEMD) simulations to determine the full thermal conductivity tensor for the monoclinic high explosive crystal β-1,3,5,7-tetranitro-1,3,5,7-tetrazoctane (β-HMX). In order to do so for the monoclinic crystal, four directions for heat propagation are used. Effects of the temperature and pressure are investigated between 200 and 500 K and 0 and 5 GPa, respectively, which approximately covers the range where the β polymorph is stable. Simulations are carried out with the Smith–Bharadwaj non-reactive empirical potential [Smith and Bharadwaj, J. Phys. Chem. B 103, 3570 (1999)], which is known to reproduce well the thermo-elastic properties of HMX. Our results indicate that the thermal conductivity, κ, is highly anisotropic, with 36% difference between the two extreme values at 300 K and 0 GPa. A simple function is used to interpolate κ in the pressure-temperature regime considered in this study, which can be used in continuum models. The results from RNEMD simulations compare well with available experimental results from the literature and allow the determination of κ for any direction and temperature and pressure within or around the fitting interval.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Dynamics of heavy quarks in strongly coupled $\mathcal{N}$ = 4 SYM plasma

We calculate the probability distribution P(k) for a heavy quark with velocity v propagating through strongly coupled N = 4 SYM plasma in the ’t Hooft limit (N c → ∞, λ = g 2 N c → ∞) at a temperature T to acquire a momentum k due to interactions with the plasma. This distribution encodes the well-known drag coefficient η D and the transverse and longitudinal momentum diffusion coefficients κ T and κ L . The jet quenching parameter $\hat{q}$ can be extracted from P(k) for v = 1. Going beyond these known Gaussian characteristics of P(k), our calculation determines all of the higher order and mixed moments to leading order in 1/$\sqrt{λ}$ for the first time. These non-Gaussian features of P(k) include qualitatively novel correlations between longitudinal energy loss and transverse momentum broadening at nonzero v. We show that all higher moments scale characteristically with an effective temperature of the boosted plasma in the heavy quark rest frame, and we demonstrate that these non-Gaussian characteristics can be sizable in magnitude and even dominant in physically relevant situations. We use these results to derive a Kolmogorov equation for the evolution of the probability distribution for the total momentum of a heavy quark that propagates through strongly coupled plasma. This evolution equation accounts for all higher order correlations between transverse momentum broadening and longitudinal energy loss, which we have calculated from first principles. It reduces to a Fokker-Planck equation when truncated to only include the effects of η D , κ T and κ L . Remarkably, while heavy quarks do not reach kinetic equilibrium with the plasma if evolved with this Fokker-Planck equation, by showing that the Boltzmann distribution is a static solution of the all-order Kolmogorov equation that we have derived we demonstrate that heavy quarks do reach kinetic equilibrium if evolved with this equation. Our results thus provide a dynamically complete framework for understanding the thermalization of a heavy quark that may be initially far from equilibrium in the strongly coupled N = 4 SYM plasma — as well as new insight into heavy quark transport and equilibration in quark-gluon plasma.

Holography and Hydrodynamics↗
DOE OSTI2025

On the Stability Analysis of Astrophysical Cooling Functions

To model the temperature evolution of optically thin astrophysical environments at MHD scales, radiative and collisional cooling rates are typically either pretabulated or fit into a functional form and then input into MHD codes as a radiative loss function. Thermal balance requires estimates of the analogous heating rates, which are harder to calculate, and due to uncertainties in the underlying dissipative heating processes these rates are often simply parameterized. The resulting net cooling function defines an equilibrium curve that varies with density and temperature. Such cooling functions can make the gas prone to thermal instability (TI), which will cause departures from equilibrium. There has been no systematic study of thermally unstable parameter space for nonequilibrium states. Motivated by our recent finding that there is a related linear instability, catastrophic cooling instability, that can dominate over TI, here we carry out such a study. We show that Balbus instability criteria for TI can be used to define a critical cooling rate, Λ c , that permits a nonequilibrium analysis of cooling functions through the mapping of TI zones. We furthermore illustrate how thermal conduction modifies the shape of TI zones. Upon applying a Λ c -based stability analysis to coronal loop simulations, we find that loops undergoing periodic episodes of coronal rain formation are linearly unstable to catastrophic cooling instability, while TI is stabilized by thermal conduction.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Electron thermalization and relaxation in laser-heated nickel by few-femtosecond core-level transient absorption spectroscopy

Direct measurements of photoexcited carrier dynamics in nickel are made using few-femtosecond extreme ultraviolet (XUV) transient absorption spectroscopy at the nickel M 2,3 edge. It is observed that the core-level absorption line shape of photoexcited nickel can be described by a Gaussian broadening (σ) and a red shift (ω s ) of the ground-state absorption spectrum. Theory predicts and the experimental results verify that after initial rapid carrier thermalization, the electron temperature increase (ΔT ) is linearly proportional to the Gaussian broadening factor σ, providing quantitative real-time tracking of the relaxation of the electron temperature. Measurements reveal an electron cooling time for 50 nm thick polycrystalline nickel films of 640 ± 80 fs. With hot thermalized carriers, the spectral red shift exhibits a power-law relationship with the change in electron temperature of ω s ∝ ΔT 1.5 . Rapid electron thermalization via carrier-carrier scattering accompanies and follows the nominal 4-fs photoexcitation pulse until the carriers reach a quasithermal equilibrium. Entwined with a <6 fs instrument response function, carrier thermalization times ranging from 34 fs to 13 fs are estimated from experimental data acquired at different pump fluences and it is observed that the electron thermalization time decreases with increasing pump fluence. Here, the study provides an initial example of measuring electron temperature and thermalization in metals in real time with XUV light, and it lays a foundation for further investigation of photoinduced phase transitions and carrier transport in metals with core-level absorption spectroscopy.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Structural dynamics of a thermally silent triiron( II ) spin crossover defect grid complex

The structural evolution of spin crossover (SCO) complexes during their spin transition at equilibrium and out-of-equilibrium conditions needs to be understood to enable their successful utilisation in displays, actuators and memory components. In this study, diffraction techniques were employed to study the structural changes accompanying the temperature increase and the light irradiation of a defect [2 × 2] triiron(II) metallogrid of the form [Fe II 3 L H 2 (HL H ) 2 ](BF 4 ) 4 ·4MeCN (FE3), L H = 3,5-bis{6-(2,2'-bipyridyl)}pyrazole. Although a multi-temperature crystallographic investigation on single crystals evidenced that the compound does not exhibit a thermal spin transition, the structural analysis of the defect grid suggests that the flexibility of the grid, provided by a metal-devoid vertex, leads to interesting characteristics that can be used for intermolecular cooperativity in related thermally responsive systems. Time-resolved photocrystallography results reveal that upon excitation with a ps laser pulse, the defect grid shows the first two steps of the out-of-equilibrium process, namely the photoinduced and elastic steps, occurring at the ps and ns time scales, respectively. Similar to a previously reported [2 × 2] tetrairon(II) metallogrid, FE3 exhibits a local distortion of the entire grid during the photoinduced step and a long-range distortion of the lattice during the elastic step. Although the lifetime of the pure photoinduced high spin (HS) state is longer in the tetranuclear grid than in the defect grid, suggesting that the global nuclearity plays a crucial role for the lifetime of the photoinduced species, the influence of the co-crystalising solvent on the lifetime of the photoinduced HS state remains unknown. This study sheds light on the out-of-equilibrium dynamics of a thermally silent defect triiron SCO metallogrid.

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

Priming self-assembly pathways by stacking block copolymers

Abstract Block copolymers spontaneously self-assemble into well-defined nanoscale morphologies. Yet equilibrium assembly gives rise to a limited set of structures. Non-equilibrium strategies can, in principle, expand diversity by exploiting self-assembly’s responsive nature. In this vein, we developed a pathway priming strategy combining control of thin film initial configurations and ordering history. We sequentially coat distinct materials to form prescribed initial states, and use thermal annealing to evolve these manifestly non-equilibrium states through the assembly landscape, traversing normally inaccessible transient structures. We explore the enormous associated hyperspace, spanning processing (annealing temperature and time), material (composition and molecular weight), and layering (thickness and order) dimensions. We demonstrate a library of exotic non-native morphologies, including vertically-oriented perforated lamellae, aqueduct structures (vertical lamellar walls with substrate-pinned perforations), parapets (crenellated lamellae), and networks of crisscrossing lamellae. This enhanced structural control can be used to modify functional properties, including accessing regimes that surpass their equilibrium analogs.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

In Situ Insights into Cathode Calcination for Predictive Synthesis: Kinetic Crystallization of LiNiO 2 from Hydroxides

Abstract Calcination is a solid‐state synthesis process widely deployed in battery cathode manufacturing. However, its inherent complexity associated with elusive intermediates hinders the predictive synthesis of high‐performance cathode materials. Here, correlative in situ X‐ray absorption/scattering spectroscopy is used to investigate the calcination of nickel‐based cathodes, focusing specifically on the archetypal LiNiO 2 from Ni(OH) 2 . Combining in situ observation with data‐driven analysis reveals concurrent lithiation and dehydration of Ni(OH) 2 and consequently, the low‐temperature crystallization of layered LiNiO 2 alongside lithiated rocksalts. Following early nucleation, LiNiO 2 undergoes sluggish crystallization and structural ordering while depleting rocksalts; ultimately, it turns into a structurally‐ordered layered phase upon full lithiation but remains small in size. Subsequent high‐temperature sintering induces rapid crystal growth, accompanied by undesired delithiation and structural degradation. These observations are further corroborated by mesoscale modeling, emphasizing that, even though calcination is thermally driven and favors transformation towards thermodynamically equilibrium phases, the actual phase propagation and crystallization can be kinetically tuned via lithiation, providing freedom for structural and morphological control during cathode calcination.

36 MATERIALS SCIENCE↗
DOE OSTI2023

An adjoint-based method for optimising MHD equilibria against the infinite- n , ideal ballooning mode

We demonstrate a fast adjoint-based method to optimise tokamak and stellarator equilibria against a pressure-driven instability known as the infinite-n ideal ballooning mode. We present three finite-β (the ratio of thermal to magnetic pressure) equilibria: one tokamak equilibrium and two stellarator equilibria that are unstable against the ballooning mode. Using the self-adjoint property of ideal magnetohydrodynamics, we construct a technique to rapidly calculate the change in the eigenvalue, a measure of ideal ballooning instability. Using the SIMSOPT optimisation framework, we then implement our fast adjoint gradient-based optimiser to minimise the eigenvalue and find stable equilibria for each of the three originally unstable equilibria.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Novel Poly-Si:Ga/SiOx Passivating Contacts through Non-Equilibrium Doping

Poly-Si/SiOx passivating contacts are one of the key enablers for high-efficiency, low-cost c-Si solar cells. In recent years, record devices have reached efficiencies of 26% in the laboratory and 25% in the industry. The current cell structure utilizing a tunneling oxide passivating contact structure with a back phosphorus-doped poly- Si/SiOx passivating contact and a front boron diffused emitter suffers from large emitter recombination. Thus, replacing the front B diffused emitter with a p-type passivating contact is a route to overcome this deficiency. This gives rise to front/back poly-Si based passivating contacts. To address the low passivation performance of the B-doped poly-Si passivating contacts, we replace B with Ga as a novel p-type dopant to avoid dopant accumulation in the tunneling oxide, which is known to lead to large degradation loss in passivation quality. Here, we introduce a non-equilibrium method via pulsed laser melting to thermally melt and recrystallize the poly-Si and redistribute the dopants, achieving doping concentrations above the solid solubility limit (~4E19 cm-3). We demonstrate a good passivation quality with an iVoc of 721 mV with an active Ga doping concentration in poly-Si >1020 cm-3. Furthermore, we show a low contact resistivity of 33.2 +/- 9.3 mO cm2 using a diode model calculation. Finally, cross-section scanning spreading resistance microscopy was performed to determine the resistance profile across the non-homogeneously doped poly-Si layer.

c-Si↗
DOE OSTI2022

Thermal conductivity tensor of γ and ε -hexanitrohexaazaisowurtzitane as a function of pressure and temperature

Using reverse non-equilibrium molecular dynamics simulations, we have determined the dependences on temperature and pressure of the thermal conductivity tensors for the monoclinic γ and ε polymorphs of hexanitrohexaazaisowurtzitane (HNIW or CL20). In this work, a recently developed non-reactive force field [X. Bidault and S. Chaudhuri, RSC Adv. 9, 39649–39661 (2019)], designed to study polymorphism and phase transitions in CL20, is employed. The effects of temperature and pressure are investigated between 200 and 500 K and up to 0.5 GPa for γ-CL20 and 2 GPa for ε-CL20. In order to obtain the full thermal conductivity tensor, $κ_{ij}$, for the monoclinic crystals, four distinct heat propagation directions are used. We find that $κ_{ij}$ for both polymorphs is more isotropic than for other energetic molecular crystals, including α- and γ-RDX, β-HMX, and PETN, with a maximum difference of 9.8% between orientations observed at 300 K and 0 GPa for γ-CL20 and a maximum difference of 4.8% for ε-CL20. The average thermal conductivity, $\bar{κ}$, of ε-CL20 is 6.4% larger than that of γ-CL20 at 300 K and 0 GPa. Analytic linear functions of the inverse temperature and the pressure are provided, which fit the data well and can be used to predict the thermal conductivity of both polymorphs for any orientation, pressure, and temperature in and around the fitting range. Our predictions agree reasonably well with the limited available experimental data, for which the polymorph type is unknown.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Interpretation of electromagnetic modes in the sub-TAE frequency range in JET plasmas with elevated monotonic q-profiles

Recent JET deuterium experiments with an advanced tokamak scenario using an internal transport barrier (ITB) have been performed to clearly observe destabilized toroidicity-induced Alfvén eigenmodes (TAEs) by fast ions; interestingly, these also exhibit unstable electromagnetic (EM) perturbations in the sub-TAE frequency range. We identify such EM perturbations to be beta-induced ion temperature gradient (BTG) eigenmodes and not beta-induced Alfvén eigenmodes nor beta-induced Alfvén acoustic eigenmodes, which are often unstable in such high-beta plasmas with high-power neutral beam injection (NBI). The BTG modes are the most unstable modes due to the high thermal ion temperature gradient related to the ITB, high thermal ion temperature compared to thermal electron temperature (high T i /T⁠ e ), and a high ion beta regime. BTG mode experimental characteristics match the analytical theory, i.e., location in the vicinity of a rational magnetic surface with a low magnetic shear, mode frequency scaling with the ion diamagnetic frequency (⁠ω$^*_i$⁠), and a coupling among Alfvén and drift waves. We also perform linear gyrokinetic simulations with validated plasma profiles and equilibrium, and find a mode kinetically driven by thermal ions with similar characteristics as the experimental BTG modes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Low-Temperature Plasma Assisted Kinetics Study of Ethanol

Abstract: Non- equilibrium plasma assisted kinetics of ethanol at 0.5 atm over a temperature range of 523-1203 K were studied at stoichiometric conditions by using a plasma flow reactor (PFR). The effects of plasma chemistry on fuel and its influence on neutral chemistry is decoupled with both experimental and modelling efforts. Enhanced reactivity was observed experimentally in both plasma-assisted oxidation and pyrolysis cases with fuel consumption beginning at temperatures as low as 523 K. Absolute fuel consumption was observed at temperatures nearly 200 K earlier than pure thermal reactions highlighting the effectiveness of non-equilibrium plasma to promote reactivity. Modelling results illustrated the effect of electron impact reactions towards accelerating fuel consumption and promoting low temperature chemistry by generating excited state N2 and reactive O and OH radicals.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

WIMP-like Dark Matter Without Thermalization At Freeze-Out

In the standard thermal relic scenario, dark matter remains in chemical equilibrium with the Standard Model radiation bath until freeze-out occurs at $T \sim m_X/20$, where $m_X$ is the dark matter mass. In this familiar class of models, the observed relic density is obtained for annihilation cross sections of order $σv \sim 10^{-26}$ cm$^3$/s. We show that comparable cross sections can naturally be realized in hidden-sector models in which the dark matter and Standard Model sectors decouple at a very high temperature, $T \gg m_X$, and subsequently evolve with separate thermal histories. Despite this decoupling, the two sectors have similar temperatures during freeze-out, leading to the usual thermal relic abundance. As a consequence, the coupling between the Standard Model and hidden sectors can be extremely small, potentially placing direct detection and collider signals far below foreseeable sensitivities.

Hooper, Dan [Wisconsin U., Madison] (ORCID:0009000↗
DOE OSTI2023

Holographic open quantum systems: toy models and analytic properties of thermal correlators

We present a unified picture of open quantum systems, the theory of a system probing a noisy thermal environment, distilling lessons learnt from previous holographic analyses. Our treatment is applicable both when the system is coupled to short-lived (Markovian), and long-lived (non-Markovian) environmental degrees of freedom. The thermal environment is modeled using an asymptotically AdS black hole, and the systems of interest are simple probe field theories. The effective stochastic dynamics of the system is governed by real-time thermal correlators, which we compute using the gravitational Schwinger-Keldysh (grSK) geometry. We describe the structure of arbitrary tree-level contact and exchange Witten diagrams in the grSK geometry. In particular, we argue, that all such diagrams reduce to integrals supported on a single copy of the exterior of the black hole. The integrand is obtained as a multiple discontinuity of a function comprising ingoing boundary-bulk propagators, monodromy functions which appear as radial Boltzmann weights, and vertex factors. These results allow us to deduce the analytic structure of real-time thermal n-point functions in holographic CFTs. We illustrate the general statements by a two-dimensional toy model, dual to fields in the BTZ background, which we argue captures many of the essential features of generic open holographic QFTs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Generic Hilbert space fragmentation in Kogut-Susskind lattice gauge theories

At the heart of quantum many-body physics lies the understanding of mechanisms that avoid quantum thermalization in an isolated system quenched far from equilibrium. A prominent example is Hilbert space fragmentation, which has recently emerged as an ergodicity-breaking mechanism in constrained spin models. Here, we show that Kogut-Susskind formulations of lattice gauge theories in d+1D (d spatial and one temporal dimensions) give rise to Hilbert space fragmentation, and discuss possible implications for understanding continuum physics. Lastly, our findings not only prove that lattice gauge theories are a natural platform for Hilbert space fragmentation, they also serve as a guide to the conditions under which these models can be faithfully used to infer the thermalization properties of quantum chromodynamics.

Eigenstate thermalization↗