Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Thermal state”

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

State preparation of lattice field theories using quantum optimal control

Here, we explore the application of quantum optimal control (QOC) techniques to state preparation of lattice field theories on quantum computers. As a first example, we focus on the Schwinger model, quantum electrodynamics in 1+1 dimensions. We demonstrate that QOC can significantly speed up the ground state preparation compared to gate-based methods, even for models with long-range interactions. Using classical simulations, we explore the dependence on the interqubit coupling strength and the device connectivity, and we study the optimization in the presence of noise. While our simulations indicate potential speedups, the results strongly depend on the device specifications. In addition, we perform exploratory studies on the preparation of thermal states. Our results motivate further studies of QOC techniques in the context of quantum simulations for fundamental physics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Quantum Joule expansion of one-dimensional systems

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

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Synchrotron Firehose Instability

We demonstrate using linear theory and particle-in-cell (PIC) simulations that a synchrotron-cooling collisionless plasma acquires pressure anisotropy and, if the plasma beta is sufficiently high, becomes unstable to the firehose instability, in a process that we dub the synchrotron firehose instability (SFHI). The SFHI channels free energy from the pressure anisotropy of the radiating, relativistic electrons (and/or positrons) into small-amplitude, kinetic-scale, magnetic-field fluctuations, which pitch-angle scatter the particles and bring the plasma to a near-thermal state of marginal instability. The PIC simulations reveal a nonlinear cyclic evolution of firehose bursts interspersed by periods of stable cooling. We compare the SFHI for electron–positron and electron–ion plasmas. As a byproduct of the growing electron-firehose magnetic-field fluctuations, magnetized ions gain a pressure anisotropy opposite to that of the electrons. If these ions are relativistically hot, we find that they also experience cooling due to collisionless thermal coupling with the electrons, which we argue is mediated by a secondary ion-cyclotron instability. We suggest that the SFHI may be activated in a number of astrophysical scenarios, such as within ejecta from black hole accretion flows and relativistic jets, where the redistribution of energetic electrons from low to high pitch angles may cause transient bursts of radiation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

HFIR LEU High Density Silicide Dispersion Optimized Design Neutronics Analyses with PHAME

A high-fidelity neutronics model of the Oak Ridge National Laboratory High Flux Isotope Reactor (HFIR) with the low-enriched uranium (LEU) high-density silicide dispersion Optimized fuel design was updated and analyzed to generate reactor physics-based metrics to support follow-on thermal hydraulic and transient analyses of this design. The Python HFIR Analysis and Measurement Engine (PHAME) was also updated to enhance the automation capabilities of the framework developed and maintained to perform these reactor physics modeling and simulation efforts. The automated framework significantly increases the efficiency and reproducibility to design and thoroughly analyzes HFIR LEU core designs, changes, and uncertainties. Reactor physics metrics evaluated include but are not limited to fuel depletion, cycle length, fission rate density distributions, axial power peaking factors, kinetics data, reactivity coefficients, control element worths, heat deposition rates, and decay heat. These neutronics results provide essential input to follow-on steady state thermal, thermal hydraulic and reactor transient analyses, which are subject of other reports. The Optimized design operates at 95 MW to maintain HFIR’s current highly enriched uranium core performance level at 85 MW.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Eigenstate entanglement scaling for critical interacting spin chains

With increasing subsystem size and energy, bipartite entanglement entropies of energy eigenstates cross over from the groundstate scaling to a volume law. In previous work, we pointed out that, when strong or weak eigenstate thermalization (ETH) applies, the entanglement entropies of all or, respectively, almost all eigenstates follow a single crossover function. The crossover functions are determined by the subsystem entropy of thermal states and assume universal scaling forms in quantum-critical regimes. This was demonstrated by field-theoretical arguments and the analysis of large systems of non-interacting fermions and bosons. Here, we substantiate such scaling properties for integrable and non-integrable interacting spin-1/2 chains at criticality using exact diagonalization. In particular, we analyze XXZ and transverse-field Ising models with and without next-nearest-neighbor interactions. Indeed, the crossover of thermal subsystem entropies can be described by a universal scaling function following from conformal field theory. Furthermore, we analyze the validity of ETH for entanglement in these models. Even for the relatively small system sizes that can be simulated, the distributions of eigenstate entanglement entropies are sharply peaked around the subsystem entropies of the corresponding thermal ensembles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Coupled Multiphysics Modeling of Lithium-Ion Batteries for Automotive Crashworthiness Applications

Considerable advances have been made in battery safety models, but achieving predictive accuracy across a wide range of conditions continues to be challenging. Interactions between dynamically evolving mechanical, electrical, and thermal state variables make model prediction difficult during mechanical abuse scenarios. In this study, we develop a physics-based modeling approach that allows for choosing between different mechanical and electrochemical models depending on the required level of analysis. We demonstrate the use of this approach to connect cell-level abuse response to electrode-level and particle-level transport phenomena. A pseudo-two-dimensional model and simplified single-particle models are calibrated to electrical-thermal cycling data and applied to mechanically induced short-circuit scenarios to understand how the choice of electrochemical model affects the model prediction under abuse scenarios. These models are implemented using user-defined subroutines on ls-dyna finite element software and can be coupled with existing automotive crash safety models.

analysis and design of components↗
DOE OSTI2021

Adiabatic quantum decoherence in many non-interacting subsystems induced by the coupling with a common boson bath

Highlights: • System–environment quantum correlation: a main solid state NMR decoherence channel. • Non-separable system–environment model yields realistic spin decoherence rates. • New open quantum system approach explains irreversible decay of refocused NMR echoes. • Adiabatic quantum decoherence is inherently irreversible and eigen-selective. This work addresses adiabatic quantum decoherence of many-body spin systems coupled with a boson field in the framework of open quantum systems theory. We generalize the traditional spin-boson model by considering a system–environment interaction Hamiltonian that represents a partition of non-interacting subsystems and highlights the collective correlation that appears exclusively due to the coupling with a common environment. Remarkably, this simple, exactly solvable model encompasses relevant aspects of a many-body open quantum system and features the subtle quantum effects that arise when the size scales up to a macroscopic level. We derive an analytical expression for the time dependence of the density matrix elements (in the preferred basis) without assuming coarse-graining. The resulting decoherence function is eigen-selective and is a complex exponential whose exponent has a real part that introduces a decay similar to that in the spin-boson model. On the contrary, the imaginary part depends on the quantum numbers and geometry of the whole partition and does not reflect the system temperature. Motivated by decoherence in solid-state NMR, and in search of realistic numerical estimations, we apply the theoretical results to a partition of dipole-coupled spin pairs in contact with a common phonon bath, using typical parameters of hydrated salts. The proposal allows estimating the decoherence time scale in terms of the system physical constants: sound velocity and eigenvalue distribution width. As a significant novelty, the decoherence function phase depends on the eigenvalue distribution throughout the sample. It plays the leading role, overshadowing the mechanism associated with the bath thermal state. Finally, we apply the formalism to describe decoherence in the “magic echo” NMR reversal experiment. We find that the system–environment correlation explains the origin of irreversibility, and both the decoherence rate value and its dependence on the dipolar frequency, are remarkably similar to the experiment.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Behavior of an air-cooled thermosiphon during a non-air gas ingress event

Here, this study investigates the transient response of an air-cooled thermosiphon subjected to a sudden non-air gas ingress event. Results presented are from experimental testing conducted on a large-scale thermal hydraulic facility that examines heat removal performance of the Reactor Cavity Cooling System (RCCS) concept, a passive safety system for advanced nuclear reactors. Relying on buoyancy driven natural circulation flow, these systems feature a series of riser standpipes and coupled chimney network that provide a pathway for heat rejection to the atmosphere. In this work, the facility was modified to allow a transition of the draft intake from normal ambient air, to draft intake from an open-top volume containing high purity argon. Testing began by establishing normal operation and allowing the facility to reach steady-state thermal hydraulic flow conditions with natural air draft. The transitioned sequence was then initiated which simultaneously changed the inlet boundary from open air to a pathway from the quiescent argon volume. This event caused 1,200 cu-feet of the heavy gas, twice the internal volume of total facility flow path, to be drawn into the inlet plenum and ingress into the heated riser standpipes. With the upper chimney network still containing residual air, there was insufficient density difference to maintain buoyancy driven natural circulation. After a period of only 90 s past the ingress event, the facility experienced complete flow stagnation before entering an extended period of severely degraded system flow. Due to the cessation of bulk fluid movement and subsequent failure of its heat removal function, fluid and structural temperatures began to rise sharply. Re-circulation patterns developed within the multiple parallel riser standpipes, where hot gas near the outlets was observed to travel downward and re-enter the inlet of adjacent channels. After approximately 18-minutes, fluid temperatures and their associated density difference rose to a level sufficient to allow reestablishment of buoyancy driven system flow, and ultimately, recovery of facility operation to normal behavior.

42 ENGINEERING↗
DOE OSTI2022

Phase transitions of correlations in black hole geometries

We study the holographic realization of optimized correlation measures—measures of quantum correlation that generalize elementary entropic formulas—in two-dimensional thermal states dual to spacetimes with a black hole horizon. We consider the symmetric bipartite optimized correlation measures: the entanglement of purification, Q-correlation, R-correlation, and squashed entanglement, as well as the mutual information, a nonoptimized correlation measure, and identify the bulk surface configurations realizing their geometric duals over the parameter space of boundary region sizes and the black hole radius. This parameter space is divided into phases associated with given topologies for these bulk surface configurations, and first-order phase transitions occur as a new topology of bulk surfaces becomes preferred. The distinct phases can be associated with different degrees of correlation between the boundary regions and the thermal environment. The Q-correlation has the richest behavior, with a structure of nested optimizations leading to two topologically distinct bulk surface configurations being equally valid as geometric duals at generic points in the phase diagram.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Phonon-assisted formation of an itinerant electronic density wave

Abstract Electronic instabilities drive ordering transitions in condensed matter. Despite many advances in the microscopic understanding of the ordered states, a more nuanced and profound question often remains unanswered: how do the collective excitations influence the electronic order formation? Here, we experimentally show that a phonon affects the spin density wave (SDW) formation after an SDW-quench by femtosecond laser pulses. In a thin film, the temperature-dependent SDW period is quantized, allowing us to track the out-of-equilibrium formation path of the SDW precisely. By exploiting its persistent coupling to the lattice, we probe the SDW through the transient lattice distortion, measured by femtosecond X-ray diffraction. We find that within 500 femtoseconds after a complete quench, the SDW forms with the low-temperature period, directly bypassing a thermal state with the high-temperature period. We argue that a wavevector-matched phonon launched by the quench changes the formation path of the SDW through the dynamic pinning of the order parameter.

42 ENGINEERING↗
DOE OSTI2024

SISGR: The regulation of carbon fixation in plant and green algae: Rubisco activase and the origin of heat inactivation of CO 2 assimilation

Rubisco activase (Rca) is a critical AAA+ ATPase protein complex that remodels and promotes the Rubisco enzyme, a key player in photosynthetic performance and carbon fixation. The assembly and function of the Rca protein complex are regulated by a range of factors, including subunit concentration, nucleotide-binding states, thermal conditions, metal-ion coordination, and post-translational modifications, such as phosphorylation. Despite its importance in photosynthesis, the detailed molecular mechanisms underlying the regulation of plant Rca and how it activates Rubisco remain elusive. This project aims to bridge this knowledge gap by integrating sophisticated enzymology tools with single-molecule methods and high-resolution electron microscopy to elucidate the structure and function of plant Rca. Through these multiple approaches, we have systematically investigated how the activity of plant Rca is impacted by various factors, such as phosphorylation and metal-ion coordination. The Rca complex assembly/disassembly dynamics were captured using anti-Brownian electrokinetic (ABEL) trap-based measurements, providing unprecedented insight into its structural flexibility and diverse assembly states. Furthermore, the structural analysis of Rca through electron crystallography and single-particle cryogenic electron microscopy (cryo-EM) reveals novel assembly states of the spinach Rca, providing insight into the mechanistic action for Rubisco remodeling. By combining cutting-edge tools and approaches, this work uncovers critical aspects of Rca’s regulation and assembly, paving the way for a deeper understanding of its role in photosynthetic efficiency and the potential for enhancing carbon fixation in crops.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Measurements of the thermal and ionization state of the intergalactic medium during the cosmic afternoon

We perform the first measurement of the thermal and ionization state of the intergalactic medium (IGM) across 0.9 < z < 1.5 using 301 Ly α absorption lines fitted from 12 archival Hubble Space Telescope Space Telescope Imaging Spectrograph quasar spectra. We employ the machine-learning-based inference method that uses joint Doppler parameter–column density (⁠b-N HI ⁠) distributions obtained from Ly α forest decomposition. Our results show that the Γ HI photoionization rates, ⁠, agree with recent ultraviolet background synthesis models, with log(Γ HI /s -1 ) = $-11.79^{+0.18}_{-0.15}$, $-11.98^{+0.09}_{-0.09}$⁠, and $-12.32^{+0.10}_{-0.12}$⁠, at z = 1.4, 1.2, and 1, respectively. We obtain the IGM temperature at the mean density, T 0 ⁠, and the adiabatic index, γ⁠, as [log(T 0 /K), γ] = $[4.13^{+0.12}_{-0.10}, 1.34^{+0.10}_{-0.15}]$, $[3.79^{+0.11}_{-0.11}, 1.70^{+0.09}_{-0.09}]$, and $[4.12^{+0.15}_{-0.25}, 1.34^{+0.21}_{-0.26}]$ at z = 1.4⁠, 1.2, and 1. Our measurements of T 0 at z = 1.4 and 1.2 are consistent with the trend predicted from previous z < 3 temperature measurements and theoretical expectations, where the IGM cools down after $He\tiny{II}$ reionization in the absence of any non-standard heating. However, our T 0 measurement at z = 1 unexpectedly high IGM temperature. Given the relatively large uncertainty in these measurements, where σ T$_0$ ~ 5000 K, mostly emanating from the limited size of our data set, we cannot conclude whether the IGM cools down as expected. Lastly, we generate mock data sets to test the constraining power of future measurement with larger data sets. The results demonstrate that, with redshift path-length Δz ~ 2 for each redshift bin, three times the current data set, we can constrain the T 0 of IGM within 1500 K, which would be sufficient to constrain the IGM thermal history at z < 1.5 conclusively.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Gravitational index of the heterotic string

The fundamental heterotic string has a tower of BPS states whose supersymmetric index has an exponential growth in the charges. We construct the saddle-point of the gravitational path integral corresponding to this index. The saddle-point configuration is a supersymmetric rotating non-extremal Euclidean black hole. This configuration is singular in the two-derivative theory. We show that the addition of higher-derivative terms in four-dimensional N = 2 supergravity resolves the singularity. In doing so, we extend the recently-developed “new attractor mechanism” to include the effect of higher-derivative terms. Remarkably, the one-loop, four-derivative F-term contribution to the prepotential leads to a precise match of the gravitational and microscopic index. We also comment, using the effective theory near the horizon, on the possibility of a string-size near-extremal black hole. Our results clarify the meaning of different descriptions of this system in the literature. The thermal state transitions to a winding condensate and a gas of strings without ever reaching a small black hole, while the index is captured by the rotating Euclidean black hole solution and is constant and thus smoothly connected to the microscopic ensemble.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

STAT7 v2.0 Verification and Validation Report

This report documents the software testing which has been performed for STAT7 Version 2.0 (hereafter STAT7 v2.0), a code for the statistical propagation of uncertainties in steady-state thermal hydraulics (TH) analysis of plate-fueled research and test reactors. The key capabilities addressed in this testing were selected based on the use of the software for analysis work supporting the conversion to low-enriched uranium fuel of U.S. High-Performance Research Reactors (USHPRR) such as the Massachusetts Institute of Technology Research Reactor (MITR-II). This testing scope includes comparison of the code’s statistical processing capabilities using statistical tests, comparison of coolant property evaluations to NIST data tables, and systematic comparison of key TH solver capabilities to the results of hand calculations and to predictions from PLTEMP/ANL Version 4.4. The

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Physics-Infused AI/ML Based Digital-Twin Framework for Flow-Induced-Vibration Damage Prediction in a Nuclear Reactor Heat Exchanger

This report summarizes some of the ongoing work related to the development of an expert-elicitation-digital-twin framework for real time damage state prediction in heat exchanger components of a nuclear reactor. The framework is targeted towards predicting damage associated with coupled low cycle fatigue (associated with regular heat-up, cool-down and power operation transients) and high cycle fatigue (associated with flow induced vibration transients). The overall framework will be based on a NoSQL based database, physics-infused-geometry-dependent virtual-sensor data, different AI/ML techniques-based data-driven-predictive-model applications (Apps) and real-time plant sensor measurements available through few existing sensors. Towards this overall goal, this report updates some of the ongoing work, such as on implementation of a NoSQL Database (such as MongoDB), FE based heat transfer analysis of a heat exchanger (e.g. of a PWR steam generator) for generating geometry-dependent virtual sensor data and evaluation of various AI/ML models such as based on multivariate linear regression, ensembled decision-tree based Random-Forest and Gradient-Boosting regression and high-dimensional-kernel-function-transformation based Support-Vector-Machine regression models. The AI/ML models were evaluated for predicting multi-time-series thermal states at thousands of 3D point-clouds

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Enhancement of PyARC for Westinghouse Electric Company’s Lead Fast Reactor Design and Modeling (Final TCF Report)

Westinghouse Electric Company is a nuclear reactor vendor headquartered in the U.S. that is developing advanced reactor technology for the U.S. and global markets. Westinghouse has been relying on the neutronics Argonne Reactor Codes (ARC) executed through the NEAMS Workbench and its PyARC module that are developed under the DOE-NE Nuclear Energy Advanced Modeling and Simulation (NEAMS) and Advanced Reactor Technology (ART) – Fast Reactor programs. Through this user experience, Westinghouse identified several enhancements that would benefit the ARC codes’ usability by the US industry and therefore its commercialization potential. The enhancements were proposed to deliver both improvements in workflow and analysis capabilities to better support effective fast reactor core design and analysis to the nuclear industry. The PyARC workflow was extended in this project by integrating non-neutronic ARC codes DASSH and NUBOW-3D. The Ducted Assembly Steady-State Heat equation (DASSH) code is developed at ANL to perform steady-state thermal hydraulic sub-channel analysis in liquid metal fast reactor assemblies to determine optimized coolant flow and temperature distributions, which in this project was updated and validated for lead fast reactor (LFR) applications. The interface between REBUS and NUBOW-3D were improved in this project to assess the impact of the core restraint design and thermal induced expansion effects on the reactivity of the core, and to model the deformations of the fuel assemblies induced by temperature and irradiation. Finally, the ARC models that were extensively verified and validated through various SFR-based modeling benchmarks are extended in this project through code-to-code comparison on relevant LFR-specific neutronics benchmarks against Monte-Carlo neutronic solutions. Overall, this work enables verification of the capability of the ARC codes for a wide range of Generation-IV reactor designs. The outcome of this project is the release of a comprehensive modeling toolkit of validated, robust and efficient codes, as well as their user interface, that enables industry to perform a wide range of fast reactor analyses for design and licensing of their concepts.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Towards Quantum Computing Phase Diagrams of Gauge Theories with Thermal Pure Quantum States

The phase diagram of strong interactions in nature at finite temperature and chemical potential remains largely theoretically unexplored due to inadequacy of Monte-Carlo–based computational techniques in overcoming a sign problem. Quantum computing offers a sign-problem-free approach, but evaluating thermal expectation values is generally resource intensive on quantum computers. To facilitate thermodynamic studies of gauge theories, we propose a generalization of the thermal-pure-quantum-state formulation of statistical mechanics applied to constrained gauge-theory dynamics, and numerically demonstrate that the phase diagram of a simple low-dimensional gauge theory is robustly determined using this approach, including mapping a chiral phase transition in the model at finite temperature and chemical potential. Quantum algorithms, resource requirements, and algorithmic and hardware error analysis are further discussed to motivate future implementations. Thermal pure quantum states, therefore, may present a suitable candidate for efficient thermal simulations of gauge theories in the era of quantum computing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Concurrent Inter-Model Spread of Boreal Winter Westerly Jet Meridional Positions Between the Northern and Southern Hemispheres in CMIP6 Models

Here, this study investigates the inter-model spread of climatological extratropical westerly jets in boreal winter, using the historical simulation of 52 Coupled Model Intercomparison Project phase 6 (CMIP6) models from 1851 to 2014. The results show that there is a substantial spread in the latitude of the upper-tropospheric westerly jet across models, characterised by large inter-model standard deviations to both the poleward and equatorward sides of the jet axis, although the multi-model ensemble mean (MME) performs well in simulating meridional position of westerly jets. Furthermore, we detect the consistency of inter-model jet position spread between the Northern and Southern Hemispheres, based on the inter-model empirical orthogonal function (EOF) decomposition and correlation of regional-averaged zonal winds. Specifically, the models that simulate the westerly jets poleward/equatorward relative to the MME position in one hemisphere also tend to simulate the jets poleward/equatorward in the other hemisphere. Accordingly, we define a global jet spread index to depict the concurrence of jet shift in the two hemispheres. The results of inter-model regression analyses based on this index indicate that the models positioning the jets poleward than the MME tend to simulate a wider Hadley Cell, a poleward-shifted Ferrel Cell in the Southern Hemisphere, enhanced precipitation in the subtropics and suppressed precipitation in the tropics, and warmer sea surface temperatures in the subtropics and mid-latitudes. The present results suggest that improving the simulation of jet positions in climate models requires a comprehensive consideration of thermal states in the tropics and subtropics/mid latitudes.

54 ENVIRONMENTAL SCIENCES↗