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

Langevin dynamics of generalized spins as SU( N ) coherent states

Classical models of spin systems traditionally retain only the dipole moments, but a quantum spin state will frequently have additional structure. Spins of magnitude S have N = 2⁢S + 1 levels. Alternatively, the spin state is fully characterized by a set of N 2 - 1 local physical observables, which we interpret as generalized spin components. For example, a spin with S = 1 has three dipole components and five quadrupole components. Further, these components evolve under a generalization of the classical Landau-Lifshitz dynamics, which can be extended with noise and damping terms. In this paper, we reformulate the dynamical equations of motion as a Langevin dynamics of SU(N) coherent states in the Schrödinger picture. This viewpoint is especially useful as the basis for an efficient numerical method to sample spin configurations in thermal equilibrium and to simulate the relaxation and driven motion of topological solitons. To illustrate the approach, we simulate a nonequilibrium relaxation process that creates CP 2 skyrmions, which are topological defects with both dipole and quadrupole character.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Efficient simulation of low-temperature physics in one-dimensional gapless systems

Here, we discuss the computational efficiency of the finite-temperature simulation with minimally entangled typical thermal states (METTS). To argue that METTS can be efficiently represented as matrix product states, we present an analytic upper bound for the average entanglement Rényi entropy of METTS for a Rényi index 0 < q ≤ 1. In particular, for one-dimensional (1D) gapless systems described by conformal field theories, the upper bound scales as O⁡(cN 0 ⁢log⁡β) where c is the central charge and N is the system size. Furthermore, we numerically find that the average Rényi entropy exhibits a universal behavior characterized by the central charge and is roughly given by half of the analytic upper bound. Based on these results, we show that METTS can provide a speedup compared to employing the purification method to analyze thermal equilibrium states at low temperatures in 1D gapless systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Generalized Kompaneets formalism for inelastic neutrino-nucleon scattering in supernova simulations

Based on the Kompaneets approximation, we develop a robust methodology to calculate spectral redistribution via inelastic neutrino-nucleon scattering in the context of core-collapse supernova simulations. The resulting equations conserve lepton number to machine precision and scale linearly, not quadratically, with number of energy groups. The formalism also provides an elegant means to derive the rate of energy transfer to matter which, as it must, automatically goes to zero when the neutrino radiation field is in thermal equilibrium. Furthermore, we derive the next higher order in $ϵ/mc^2$ correction to the neutrino Kompaneets equation. Unlike other Kompaneets schema, ours also generalizes to the case of anisotropic angular distributions, while retaining the conservative form that is a hallmark of the classical Kompaneets equation. In this work, our formalism enables immediate incorporation into supernova codes that follow the spectral angular moments of the neutrino radiation fields.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Asymmetric temperature equilibration with heat flow from cold to hot in a quantum thermodynamic system

A model computational quantum thermodynamic network is constructed with two variable temperature baths coupled by a linker system, with an asymmetry in the coupling of the linker to the two baths. It is found in computational simulations that the baths come to “thermal equilibrium” at different bath energies and temperatures. In a sense, heat is observed to flow from cold to hot. Additionally, a description is given in which a recently defined quantum entropy S univ Q for a pure state “universe” continues to increase after passing through the classical equilibrium point of equal temperatures, reaching a maximum at the asymmetric equilibrium. Thus, a second law account Δ S univ Q ≥ 0 holds for the asymmetric quantum process. In contrast, a von Neumann entropy description fails to uphold the entropy law, with a maximum near when the two temperatures are equal, then a decrease Δ S v N < 0 on the way to the asymmetric equilibrium.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Neuroevolutionary Learning of Particles and Protocols for Self-Assembly

Within simulations of molecules deposited on a surface we show that neuroevolutionary learning can design particles and time-dependent protocols to promote self-assembly, without input from physical concepts such as thermal equilibrium or mechanical stability and without prior knowledge of candidate or competing structures. Furthermore, the learning algorithm is capable of both directed and exploratory design: it can assemble a material with a user-defined property, or search for novelty in the space of specified order parameters. In the latter mode it explores the space of what can be made, rather than the space of structures that are low in energy but not necessarily kinetically accessible.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Nonthermal Bonding Origin of a Novel Photoexcited Lattice Instability in SnSe

Lattice dynamics measurements are often crucial tools for understanding how materials transform between different structures. Here, we report time-resolved x-ray scattering-based measurements of the nonequilibrium lattice dynamics in SnSe, a monochalcogenide reported to host a novel photoinduced lattice instability. By fitting interatomic force models to the fluence dependent excited-state dispersion, we determine the nonthermal origin of the lattice instability to be dominated by changes of interatomic interactions along a bilayer-connecting bond, rather than of an intralayer bonding network that is of primary importance to the lattice instability in thermal equilibrium.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Elevated Electron Temperature Coincident with Observed Fusion Reactions in a Sheared-Flow-Stabilized Z Pinch

The sheared-flow-stabilized Z pinch concept has been studied extensively and is able to produce fusion-relevant plasma parameters along with neutron production over several microseconds. In this work, we present elevated electron temperature results spatially and temporally coincident with the plasma neutron source. An optical Thomson scattering apparatus designed for the FuZE device measures temperatures in the range of 1–3 keV on the axis of the device, 20 cm downstream of the nose cone. The 17-fiber system measures the radial profiles of the electron temperature. Scanning the laser time with respect to the neutron pulse time over a series of discharges allows the reconstruction of the T e temporal response, confirming that the electron temperature peaks simultaneously with the neutron output, as well as the pinch current and inductive voltage generated within the plasma. Comparison to spectroscopic ion temperature measurements suggests a plasma in thermal equilibrium. The elevated T e confirms the presence of a plasma assembled on axis, and indicates limited radiative losses, demonstrating a basis for scaling this device toward net gain fusion conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Biophysical and Structural Features of αβT ‐Cell Receptor Mechanosensing: A Paradigmatic Shift in Understanding T‐Cell Activation

ABSTRACT αβT cells protect vertebrates against many diseases, optimizing surveillance using mechanical force to distinguish between pathophysiologic cellular alterations and normal self‐constituents. The multi‐subunit αβT‐cell receptor (TCR) operates outside of thermal equilibrium, harvesting energy via physical forces generated by T‐cell motility and actin‐myosin machinery. When a peptide‐bound major histocompatibility complex molecule (pMHC) on an antigen presenting cell is ligated, the αβTCR on the T cell leverages force to form a catch bond, prolonging bond lifetime, and enhancing antigen discrimination. Under load, the αβTCR undergoes reversible structural transitions involving partial unfolding of its clonotypic immunoglobulin‐like (Ig) domains and coupled rearrangements of associated CD3 subunits and structural elements. We postulate that transitions provide critical energy to initiate the signaling cascade via induction of αβTCR quaternary structural rearrangements, associated membrane perturbations, exposure of CD3 ITAMs to phosphorylation by non‐receptor tyrosine kinases, and phase separation of signaling molecules. Understanding force‐mediated signaling by the αβTCR clarifies long‐standing questions regarding αβTCR antigen recognition, specificity and affinity, providing a basis for continued investigation. Future directions include examining atomistic mechanisms of αβTCR signal initiation, performance quality, tissue compliance adaptability, and T‐cell memory fate. The mechanotransduction paradigm will foster improved rational design of T‐cell based vaccines, CAR‐Ts, and adoptive therapies.

Immunology↗
DOE OSTI2020

Cryo-EM and MD infer water-mediated proton transport and autoinhibition mechanisms of V o complex

Rotary vacuolar adenosine triphosphatases (V-ATPases) drive transmembrane proton transport through a Vo proton channel subcomplex. Despite recent high-resolution structures of several rotary ATPases, the dynamic mechanism of proton pumping remains elusive. Here, we determined a 2.7-Å cryo–electron microscopy (cryo-EM) structure of yeast V o proton channel in nanodisc that reveals the location of ordered water molecules along the proton path, details of specific protein-lipid interactions, and the architecture of the membrane scaffold protein. Moreover, we uncover a state of V o that shows the c-ring rotated by ~14°. Molecular dynamics simulations demonstrate that the two rotary states are in thermal equilibrium and depict how the protonation state of essential glutamic acid residues couples water-mediated proton transfer with c-ring rotation. Our cryo-EM models and simulations also rationalize a mechanism for inhibition of passive proton transport as observed for free V o that is generated as a result of V-ATPase regulation by reversible disassembly in vivo.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Observation of the magnonic Dicke superradiant phase transition

Two-level atoms ultrastrongly coupled with single-mode cavity photons are predicted to exhibit a quantum phase transition, entering a phase in which both the atomic polarization and the photonic field are finite even without external driving. However, this phenomenon, the superradiant phase transition (SRPT), is forbidden by a no-go theorem due to the existence of the diamagnetic term. Here, we present spectroscopic evidence for a magnonic SRPT in ErFeO 3 , where the role of the photonic mode (two-level atoms) in the photonic SRPT is played by an Fe 3+ magnon mode (Er 3+ spins). The absence of the diamagnetic term in the Fe 3+ -Er 3+ exchange coupling ensures that the no-go theorem does not apply. Ultrabroadband terahertz and gigahertz magnetospectroscopy experiments revealed the signatures of the SRPT in thermal equilibrium, a kink and a softening, respectively, of two spin-magnon hybridized modes at the critical point. Systems near this phase are expected to harbor large-scale squeezing, which will potentially provide a route to next-generation quantum technologies.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Path to Increasing p-B11 Reactivity via ps and ns Lasers

The Lawson criterion for proton-boron (p-11B) thermonuclear fusion is substantially higher than that for deuterium-tritium (DT) because the fusion cross section is lower and peaks at higher ion energies. The Maxwellian averaged p-11B reactivity peaks at several hundred keV, where bremsstrahlung radiation emission may dominate over fusion reactions if electrons and ions are in thermal equilibrium and the losses are unrestricted. Nonequilibrium burn has often been suggested to realize the benefits of this aneutronic reaction, but the predominance of elastic scattering over fusion reactivity makes this difficult to achieve. The development of ultrashort pulse lasers (USPL) has opened new possibilities for initiating nonequilibrium thermonuclear burns and significant numbers of p-11B alpha particles have been reported from several experiments. We present an analysis that shows that these significant alpha yields are the result of beam fusion reactions that do not scale to net energy gain. We further find that the yields can be explained by experimental parameters and recently updated cross sections such that a postulated avalanche mechanism is not required. We use this analysis to understand the underlying physics of USPL-driven nonequilibrium fusion reactions and whether they can be used to initiate fusion burns. We conclude by outlining a path to increasing the p-11B reactivity towards the goal of achieving ignition and describing the design principles that we will use to develop a computational point design.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Replica symmetry breaking in spin glasses in the replica-free Keldysh formalism

We show that the algebra of Parisi ultrametric matrices is recovered by the real-time, replica-free, Dyson-Keldysh equations of infinite-range quantum spin glasses in the late time glassy limit. This connects to earlier results on classical and quantum systems showing how ultrametricity emerges from the persistent slow aging dynamics of the glass phase. The stationary spin glass state thereby spontaneously breaks thermal symmetry, or the Kubo-Martin-Schwinger relation of a state in global thermal equilibrium. We describe the Keldysh path integral of the infinite-range Ising model in transverse and longitudinal fields, and in the context of the Landau expansion of the action functional, show how the long-time limit connects to the full replica symmetry breaking obtained in the equilibrium formalism. We also illustrate our formalism by applying it to the spherical quantum p p -spin model, which only exhibits one-step replica symmetry breaking.

Lang, Johannes (ORCID:0000000283533392)↗
DOE OSTI2020

Impact of ionizing radiation on superconducting qubit coherence

The practical viability of any qubit technology stands on long coherence times and high-fidelity operations, with the superconducting qubit modality being a leading example. However, superconducting qubit coherence is impacted by broken Cooper pairs, referred to as quasiparticles, with a density that is empirically observed to be orders of magnitude greater than the value predicted for thermal equilibrium by the Bardeen-Cooper-Schrieffer (BCS) theory of superconductivity. Previous work has shown that infrared photons significantly increase the quasiparticle density, yet even in the best isolated systems, it still remains higher than expected, suggesting that another generation mechanism exists. In this Letter, we provide evidence that ionizing radiation from environmental radioactive materials and cosmic rays contributes to this observed difference, leading to an elevated quasiparticle density that would ultimately limit superconducting qubits of the type measured here to coherence times in the millisecond regime. We further demonstrate that introducing radiation shielding reduces the flux of ionizing radiation and positively correlates with increased coherence time. Albeit a small effect for today’s qubits, reducing or otherwise mitigating the impact of ionizing radiation will be critical for realizing fault-tolerant superconducting quantum computers.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Explicit Modeling of Pebble Temperature in the Porous-medium Framework for Pebble-bed Reactors Applications

In this study, an explicit model for pebble temperature calculations in the framework of a porous-medium approach has been developed and implemented in the SAM code. The solid-phase energy conservation equation, including the explicit modeling of pebble temperature, is a multi-scale model, which can predict the macroscopic (pebble-bed) and microscopic (pebble) temperature distributions under both steady-state and transient conditions. Extensive and successful code verifications and demonstrations have been performed for this newly developed model. By explicitly modeling pebble temperatures, this new model addresses a major deficiency of the previously implemented model, which assumed a local solid-phase thermal equilibrium that is not appropriate for pebble-bed reactor design and safety analysis applications.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Development of a BlueCRAB/MELCOR Framework for Supporting Realistic Mechanistic Source Term Calculations in Microreactors

Efforts are currently underway to deploy microreactor modeling and simulation tools to better support vendors and regulatory authorities in submitting and reviewing licensing applications. In particular, the Nuclear Regulatory Commission is expected to rely on the Comprehensive Reactor Analysis Bundle (BlueCRAB) multiphysics toolset in performing design- and beyond-design-basis accident analyses. In addition, the Nuclear Regulatory Commission has been using the MELCOR code to estimate mechanistic source terms during accidents. As MELCOR relies on isotopic inventory and reactor temperature/power evolution profiles during accident conditions—all of which can theoretically be obtained from BlueCRAB—the ultimate goal of this activity is to establish a common BlueCRAB-MELCOR framework. However, prior to the present research, BlueCRAB had never been used to calculate such quantities of interest at the full-core level. While there are many Monte Carlo (MC) codes capable of computing such quantities of interest, they are unable to readily account for multiphysics feedback. BlueCRAB allows for the coupling of different physics codes together to perform multiphysics-informed calculations. Therefore, the purpose of this fiscal year 2023 work is to investigate the feasibility and challenges of performing such calculations within BlueCRAB so as to generate the data that MELCOR relies on. To demonstrate the methodology, the proposed workflow was applied to a prototypical heat pipe-cooled microreactor model. To predict isotopic concentrations (taking into account the ac- cumulation of fission products during operation), the necessary microscopic cross sections were generated via OpenMC and tabulated with respect to temperature and burnup. Next, a recently developed capability in Griffin (the reactor physics application in BlueCRAB) was used to convert the OpenMC output format into the ISOXML format used by Griffin. A multiphysics microscopic depletion calculation that involved performing a coupled full-core, heterogeneous neutron trans- port and thermal calculation at each depletion step was conducted to deplete the core to end of life (EOL) conditions so as to provide both isotopics and the initial condition for the transient calculation. Following a brief null-transient to verify that the initial condition had been properly restarted and was indeed in thermal equilibrium, a heat pipe failure transient was simulated. Thus, the entire workflow of using BlueCRAB to generate MELCOR inputs, from cross-section generation to producing isotopic inventory and power/temperature evolution profiles during transients, is demonstrated. This report also details the identified gaps in the workflow and how they were (for the most part) addressed. Future work should focus on directly including MEL- COR into the workflow by performing a MELCOR calculation using the BlueCRAB-generated input data. In addition, the heat pipe reactor design should be improved so as to reflect more prototypical burnup characteristics at EOL.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

SpinQuest Polarized Target System

The SpinQuest experiment at Fermilab uses a solid-state polarized ammonia target held at a magnetic field of 5 T, immersed in liquid helium-4, which is held at approximately 1K by the evaporation refrigerator. The refrigerator provides the required cooling power during the dynamic nuclear polarization (DNP) process and the high intensity interaction with the 120 GeV proton beam from the Fermilab main injector. The refrigerator was designed in compliance with the American Society of Mechanical Engineers (ASME) to operate safely at Fermilab. The high pumping capacity ($17,000 \:m^3/h$) roots stack provides the required pumping speed during DNP production data taking and a custom-made radiation hard flow control valve regulates the refrigerator temperature during the thermal equilibrium calibration measurements. The frequency of the microwave generator, an Extended Interaction Oscillator (EIO), is automated to keep the maximum polarization while the (Nuclear Magnetic Resonance) NMR system continuously measures the polarization of the target material. In this talk, an overview of the SpinQuest polarized target system will be presented as well as a brief report of recent commissioning activities and target performance during the early production runs in 2024.

Bandara, Vibodha [Colombo U.]↗
DOE OSTI2020

Impact of Vibrational Nonequilibrium on the Supersonic Mode in Hypersonic Boundary Layers

The supersonic mode in hypersonic boundary layers has recently been shown to achieve a larger peak growth rate than the traditional second mode in particular circumstances relevant to free flight. In this environment, high levels of aerodynamic heating are encountered such that incorporation of chemical and vibrational nonequilibrium effects can become critical. The impact of thermal nonequilibrium (TNE) on the supersonic mode in particular has not been thoroughly and systematically investigated. This work uses thermochemical nonequilibrium (TCNE) direct numerical simulation (DNS) and TCNE linear stability theory (LST) to obtain a more complete investigation of the supersonic mode using both nonequilibrium and frozen gas models. The simulation is Mach 5 flow over a 1-mm-nose-radius axisymmetric cone 1 m in length, representative of experimental conditions in ground test facilities. Both LST and DNS results indicate that TNE effects are destabilizing compared with frozen flow and can lead to an N-factor difference of approximately +1 at the end of the cone. LST results also indicate that TNE effects are stabilizing to the second and supersonic modes in comparison with thermal equilibrium. The supersonic mode appeared to be slightly more affected by TNE than the second mode. In conclusion, collectively, the findings in this work suggest that even in ground test facilities where TCNE effects are relatively small compared with free flight, incorporating nonequilibrium models in both the mean flow and stability solvers can have a practical effect on hypersonic boundary-layer transition prediction analyses.

42 ENGINEERING↗
DOE OSTI2020

The Effect of Intensity Fluctuations on Sequential X-ray Photon Correlation Spectroscopy at the X-ray Free Electron Laser Facilities

How materials evolve at thermal equilibrium and under external excitations at small length and time scales is crucial to the understanding and control of material properties. X-ray photon correlation spectroscopy (XPCS) at X-ray free electron laser (XFEL) facilities can in principle capture dynamics of materials that are substantially faster than a millisecond. However, the analysis and interpretation of XPCS data is hindered by the strongly fluctuating X-ray intensity from XFELs. Here we examine the impact of pulse-to-pulse intensity fluctuations on sequential XPCS analysis. We show that the conventional XPCS analysis can still faithfully capture the characteristic time scales, but with substantial decrease in the signal-to-noise ratio of the g2 function and increase in the uncertainties of the extracted time constants. We also demonstrate protocols for improving the signal-to-noise ratio and reducing the uncertainties.

36 MATERIALS SCIENCE↗