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

Holographic Quantum Simulation of Strongly Correlated Electron Systems

The project aimed to demonstrate a new holographic quantum simulation approach and co‐ designed quantum hardware to tackle three specific problems that fall within the broad umbrella of unraveling the physics of strongly correlated electron systems (SCES). These tasks were: (1) holographic preparation of ground‐ and thermal‐ states of correlated magnetic and electronic systems including quasi‐2d frustrated‐spin, Fermi‐Hubbard, and fractional quantum Hall (FQH) systems, (2) holographic‐simulation of long‐time out‐of‐equilibrium dynamics and (3) holographic analogs of embedding methods such as dynamical mean‐ field theory (DMFT) and density‐matrix embedding theory (DMET) to solve systems with complex structure or long‐range interactions. These tasks are prototypes for the kinds of material simulation problems of interest to BES, such as the simulation of multiferroic materials, perovskite photovoltaics and high‐temperature superconductors, that tax the capabilities of the most powerful classical supercomputers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Gas-Cooled High-Temperature Pebble-Bed Reactor Reference Plant Model Updates

This work presents the latest improvements to, and investigations performed with, the pebble-bed high-temperature gas-cooled reactor (PB-HTGR) reference plant models for the United States Nuclear Regulatory Commission. These models serve as the foundation for the future development of detailed design evaluation models based on license applications. The reference plant models have been developed with the Comprehensive Reactor Analysis Bundle, or BlueCRAB, which is the code suite proposed for non-light-water reactor systems safety analysis. It incorporates various simulation tools developed by the Nuclear Energy Advanced Modeling and Simulation program, including the Griffin code for reactor physics, the Pronghorn and SAM codes for core thermal fluids, the BISON code for solid conduction and fuel performance, and the SAM code for system analysis. The primary objective of the work that was performed was to assess BlueCRAB’s level of readiness for modeling a PB-HTGR. To do so, we first developed numerical models in BlueCRAB that include the key physics for this technology to ensure an adequate level of fidelity for modeling PB-HTGR core performance and for performing multiphysics simulations for equilibrium core conditions and different accident scenarios. Then we simulated transient scenarios, including depressurized and pressurized loss of forced cooling accidents, over-cooling, and control rod withdrawal events with delayed and prompt supercritical reactivity insertions. The analysis in this report includes comparisons of the 2D thermal fluid porous media models in Pronghorn and SAM, and comparisons of coupled SAM/Griffin/SAM and coupled Pronghorn/Griffin for depressurized and pressurized loss of forced cooling, over-cooling, and control rod withdrawal events. In addition, we compare 3D, 2D, and 0D/PKE neutronic models for the two control rod withdrawal scenarios with coupled Pronghorn/Griffin. The comparisons show that the BlueCRAB models lead to physically intuitive solutions for the scenarios examined. The changes in the various scalar and vector fields, such as neutron flux, power, temperature, density, pressure, and velocity, are within the expected ranges, and their distributions can be explained by the system response of the transients and the geometric and material variations. Several comparisons suggest that the porous media models in Pronghorn and SAM can lead to similar solutions, even though they are based on different methodologies. This work further highlights the need for flexible tools with various levels of fidelity to cover the breadth and depth of needs that may arise in future technical evaluations of the PB-HTGR. We believe that the BlueCRAB capabilities will be a significant asset for confirmatory analyses in order to resolve important safety questions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Examining the Reactions of Ethanolamine’s Thermal Degradation Compounds in Carbon Capture through 1 H NMR and 13 C NMR

In amine scrubbing carbon capture, concerns about amine solvent degradation include whether it can affect the ability of the solvent to capture CO 2 . This study examines the interactions between the three most reported monoethanolamine (MEA) thermal degradation compounds namely, oxazolidine-2-one (OZD), N-(2-hydroxyethyl)-ethylenediamine (HEEDA), and N-(2- hydroxyethyl)-imidazoline-2-one (HEIA) with CO 2 in the presence and absence of MEA. We compared 1 H NMR, 13 C NMR, and heteronuclear single-quantum coherence (HSQC) NMR for neat OZD, HEEDA, and HEIA samples with CO 2 -loaded samples to observe changes in protonation and unique carbamate species formation. We found that OZD and HEIA did not directly react with CO 2 or undergo proton shifting based on our comparison of the neat OZD and HEIA samples with CO 2 -loaded spectra. However, we observed that the OZD can protonate when sparged with CO 2 in the presence of MEA, which suggests that the OZD acts as an intermediate. The NMR spectra for HEEDA indicated that HEEDA directly reacts with CO 2 at both amino groups and can protonate. In the presence of MEA, HEEDA and MEA can act like a solvent blend, resulting in multiple carbamates forming within the solvent. The neat HEIA spectrum, compared with CO 2 -loaded HEIA spectra, revealed similar results as OZD, where it does not react with CO 2 or protonate directly. However, HEIA does not protonate in the presence of MEA. These degradation compound reactions can increase the number of general equilibrium reactions during carbon capture and impact the model MEA solvent. This work helps provide a more complete picture of the reactions as the solvent degrades. Although this study examines CO 2 effects on thermal degradation products for MEA, other amines such as piperazine or 1- amino-3-propanol will degrade, and the degradation may speciate similarly with CO 2 . Furthermore, this study can impact and improve process models by assessing the degradation compounds’ reactions with CO 2 and potentially incorporating them into the model based on a neat solvent.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2022

Current challenges in the physics of white dwarf stars

White dwarfs are a class of stars with unique physical properties. They present many challenging problems whose solution requires advanced theories of dense matter, state-of-the-art experimental techniques, and extensive computing efforts. New ground- and space-based observatories will soon provide an increasingly detailed view of white dwarf stars and reveal new phenomena that will challenge our models. This review is an introduction to the nature of white dwarfs, the physical processes that determine their structure and evolution, and the physical conditions they span. Here we discuss a wide variety of currently unsolved or partially resolved problems in their constitutive physics that are broadly related to equations of state, transport processes and opacities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

New horizons in the holographic conformal phase transition

Abstract We describe 5D dynamical cosmological solutions of the stabilized holographic dilaton and their role in completion of the conformal phase transition. This analysis corresponds, via the AdS/CFT dictionary, to a study of out-of-equilibrium dynamics where trajectories of the dilaton do not depend solely on thermodynamic quantities in the early universe, but have sensitivity also to initial conditions. Unlike the well-studied thermal transition, which requires quantum tunneling of an infrared brane through the surface of an AdS-Schwarzschild horizon, our approach instead invokes an early epoch in which the cosmology is fully 5-dimensional, with highly relativistic brane motion and with Rindler horizons obscuring the infrared brane at early times. In this context, we demonstrate the existence of a large class of natural initial conditions that seed trajectories where the brane simply passes through the Rindler horizon and into the basin of attraction of the stabilized dilaton potential. This corresponds to successful completion of the phase transition without sacrificing perturbativity of the 5D theory.

Physics↗
DOE OSTI2023

Staggered scheme for the compressible fluctuating hydrodynamics of multispecies fluid mixtures

Here, we present a numerical formulation for the solution of nonisothermal, compressible Navier-Stokes equations with thermal fluctuations to describe mesoscale transport phenomena in multispecies fluid mixtures. The novelty of our numerical method is the use of staggered grid momenta along with a finite volume discretization of the thermodynamic variables to solve the resulting stochastic partial differential equations. The key advantages of the numerical scheme are that it significantly simplifies the discretization of diffusive and stochastic momentum fluxes into a more compact form, and it provides an unambiguous prescription of boundary conditions involving pressure. The staggered grid scheme more accurately reproduces the equilibrium static structure factor of hydrodynamic fluctuations in gas mixtures compared to a collocated scheme described previously by Balakrishnan et al. [Phys. Rev. E 89, 013017 (2014)1539-375510.1103/PhysRevE.89.013017]. The numerical method is tested for ideal noble gases mixtures under various nonequilibrium conditions, such as applied thermal and concentration gradients, to assess the role of cross-diffusion effects, such as Soret and Dufour, on the long-ranged correlations of hydrodynamic fluctuations, which are also more accurately reproduced compared to the collocated scheme. We numerically study giant nonequilibrium fluctuations driven by concentration gradients and fluctuation-driven Rayleigh-Taylor instability in gas mixtures. Wherever applicable, excellent agreement is observed with theory and measurements from the direct simulation Monte Carlo method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Coherent Reaction between Molecular and Atomic Bose-Einstein Condensates: Integrable Model

Here, we solve a model that describes a stimulated conversion between ultracold bosonic atoms and molecules. The reaction is triggered by a linearly time-dependent transition throughout the Feshbach resonance. Our solution predicts a dependence, with a dynamic phase transition, of the reaction efficiency on the transition rate for both atoms-to-molecule pairing and molecular dissociation processes. We find that for the latter process with a linear energy dispersion of atomic modes, the emerging phase can have a thermalized energy distribution of noninteracting bosons with the temperature defined by the rate of the transition. This provides a simple interpretation of the phase transition in terms of the creation of equilibrium Bose-Einstein condensate.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Fundamentals of Spark-Plasma Sintering (Final Report)

This project was focused on the development of the general theory of spark-plasma sintering (SPS) taking into account the role of both thermal and non-thermal factors in the acceleration of SPS mass transport. SPS is a particular kind of field-assisted sintering, which provides potentially revolutionary capabilities to the processing of materials into configurations previously unattainable. This approach significantly improves the processing time- and quality-wise. It carries the potential of maintaining the nano and sub-micron structure in nano-powder-based materials after consolidation. SPS gains particular prominence in connection with its exceptional potential of rapid and ultra-rapid processing of very hard-to-deform materials, which would typically require lengthy consolidation times at significantly elevated temperatures under conditions of conventional powder pressing or sintering. The achievement of the ultimate goal of the formulation of the general theory of spark plasma sintering required the identification of the contributions of all the thermal and non-thermal factors in the enhancement of mass transport under SPS processing conditions. As our research indicated, this enhancement is mostly reduced to the creation of the conditions of controlled non-equilibrium, which is especially important in ultra-rapid field-assisted sintering techniques. Thus, the main project objective was the analysis of the SPS physical basis at multiple scales specifically exploring the role of electric current in the acceleration of mass transport and with an emphasis on the conditions of controlled non-equilibrium.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Tropical Cyclone Precipitation Response to Surface Warming in Aquaplanet Simulations With Uniform Thermal Forcing

While many modeling studies have attempted to estimate how tropical cyclone (TC) precipitation is impacted by climate change, the multitude of analysis techniques and methodologies have resulted in varying conclusions. Simplified models may be able to help overcome this problem. Radiative-convective equilibrium (RCE) model simulations have been used in various configurations to study fundamental aspects of Earth's climate. While many RCE modeling studies have focused on TC genesis, intensification, and size, limited work has been done using RCE to study TC precipitation. Here, in this study, the response of TC precipitation to sea surface temperature (SST) change is analyzed in global Community Atmosphere Model (CAM) aquaplanet simulations run with Radiative-Convective Equilibrium Model Intercomparison Project protocols, with the addition of planetary rotation. We expect that the insight gained about how TC precipitation responds to SST warming will help predict how TCs in the real world respond to climate change. In the CAM RCE simulations, the warmer SST simulations have less TCs on average, but the TCs tend to be larger in outer size and more intense. As simulation SST increases, more extreme precipitation rates occur within TCs, and more of the TC precipitation comes from these extreme rates. For extreme (99th percentile) TC precipitation, SST, and TC intensity increases dominate the 8.6% per K increase, while TC outer size changes have little impact. For accumulated TC precipitation, SST, and TC intensity contributions are still the majority, but TC outer size changes also contribute to the 6.6% per K increase.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Electronic and optical characterization of bulk single crystals of cubic boron nitride (cBN)

Cubic boron nitride (cBN) is a relatively less studied wide bandgap semiconductor despite its many promising mechanical, thermal, and electronic properties. We report on the electronic, structural, and optical characterization of commercial cBN crystal platelets. Temperature dependent transport measurements revealed the charge limited diode behavior of the cBN crystals. The equilibrium Fermi level was determined to be 0.47 eV below the conduction band, and the electron conduction was identified as n-type. Unirradiated dark and amber colored cBN crystals displayed broad photoluminescence emission peaks centered around different wavelengths. RC series zero phonon line defect emission peaks were observed at room temperature from the electron beam irradiated and oxygen ion implanted cBN crystals, making this material a promising candidate for high power microwave devices, next generation power electronics, and future quantum sensing applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Fast ion stabilization of tilt in large radius FRCs

The field reversed configuration (FRC) has been a curious case in plasma physics research in that early MHD analysis suggested FRCs should be grossly unstable, while experimental results contradicted that prediction. Later, this theory was able to resolve this contradiction by understanding that finite Larmor radius effects largely negated the MHD predictions. Similarly, previous theoretical studies of beam driven FRCs predicted that such system would be unstable to beam driven modes while, again, experimental results indicated the contradiction. In this paper, we reconcile the theoretical understanding of beam driven modes with experimental observations of stability in these systems. By self-consistently capturing fast ion generation from neutral beam injection and its impact on the plasma equilibrium, we show that low amplitude perturbations in the magnetic field, driven by betatron particles, modify the precession frequencies of the betatron particles such that the drive for compressional Alfvén waves in the thermal plasma is reduced. Finally, we are able to demonstrate, for the first time, stable beam driven FRC evolution at high S*/E in 3D kinetic simulations.

Betatrons↗
DOE OSTI2025

Machine learned potential for high-throughput phonon calculations of metal—organic frameworks

Metal–organic frameworks (MOFs) are highly porous and versatile materials studied extensively for applications such as carbon capture and water harvesting. However, computing phonon-mediated properties in MOFs, like thermal expansion and mechanical stability, remains challenging due to the large number of atoms per unit cell, making traditional Density Functional Theory (DFT) methods impractical for high-throughput screening. Recent advances in machine learning potentials have led to foundation atomistic models, such as MACE-MP-0, that accurately predict equilibrium structures but struggle with phonon properties of MOFs. In this work, we developed a workflow for computing phonons in MOFs within the quasi-harmonic approximation with a fine-tuned MACE model, MACE-MP-MOF0. The model was trained on a curated dataset of 127 representative and diverse MOFs. The fine-tuned MACE-MP-MOF0 improves the accuracy of phonon density of states and corrects the imaginary phonon modes of MACE-MP-0, enabling high-throughput phonon calculations with state-of-the-art precision. The model successfully predicts thermal expansion and bulk moduli in agreement with DFT and experimental data for several well-known MOFs. These results highlight the potential of MACE-MP-MOF0 in guiding MOF design for applications in energy storage and thermoelectrics.

Elena, Alin Marin↗
DOE OSTI2024

Thermochemical Nonequilibrium Modeling in a Continuous-Galerkin, Finite-Element Framework

The presented work discusses the implementation, verification, and validation of Park's two-temperature model in a scalable, computational fluid dynamics (CFD) code developed at the US Department of Energy's Oak Ridge National Laboratory (ORNL). The implementation of Park's two-temperature model was verified through 0D test cases involving an adiabatic reactor and a nitrogen thermal bath. The implementation was then validated through comparisons with other validated CFD codes and experimental data on a hypersonic cylinder and double cones. These are standard benchmark test cases for thermochemical non-equilibrium (TCNE) modeling, and all data are shared publicly. The verification and validation results showed that ORNL's in-house CFD code could model complex, high-speed flow problems with and without TCNE modeling. This work is essential for future research involving 3D shock wave/boundary layer interactions (SBLIs).

Nutter, Nicole↗
DOE OSTI2024

UNCERTAINTIES FOR IN-REACTOR TESTING OF TRITIUM LOADED GETTERS WITH TEMPERATURE CONTROL

The fourth experiment in TPBAR Materials Irradiation Separate-Effects Test (TMIST) series (TMIST-4) is planned to experimentally determine the effect of neutron irradiation on the equilibrium vapor pressure of tritium over a tritium loaded Nickel-Plated Zirconium (NPZ) getter in a simulated TPBAR service environment. This pressure will be used to estimate the tritium permeation rate. In the thermal design in the TMIST-4, two evaluations are required. 1) The cladding temperature is calculated to estimate the feasibility of detecting potential irradiation enhancement of the partial pressure for an NPZ getter. This work is required to obtain a precise cladding temperature for measurement of a low tritium partial pressure. This work has been completed and indicated that the required/target cladding temperatures based on the TMIST-4 test plan can be satisfied under the test condition range by modifying the capsule component shapes and inserting a thermally conductive metal sleeve. 2) The evaluation of uncertainties associated with the test design and operation is performed. This includes fabrication tolerances, and power and nuclear operating conditions. These factors will impact the component temperatures inside capsules and resulting the tritium partial pressure. The tritium permeation rate is a function of temperature and tritium partial pressure. The prediction of tritium permeation rates affects the feasibility of obtaining measurement resolutions required to reach reasonable conclusion for potential irradiation enhanced performance. This work will investigate the uncertainties.

TTP, TPBAR, TMIST-4, TPBAR Materials Irradiation S↗
DOE OSTI2025

Structure of molten ytterbium aluminum garnet

Rare earth aluminum garnets are important materials in optical, dielectric, and thermal barrier applications. To advance the understanding of their melt processing and glass forming ability, we report the atomic structure of molten Yb 3 Al 5 O 12 over 1770–2630 K, which spans the equilibrium and supercooled liquid regimes. The melt density at T m = 2283 K is 5.50 g cm –3 , measured via silhouette imaging of electrostatically levitated drops over 1010–2420 K. Four separate structure measurements were made with aerodynamically levitated melts using x-ray and neutron diffraction with isotope substitution of Yb ( 172 Yb, 174 Yb, or nat Yb). Empirical potential structure refinement models were developed, which are in excellent agreement with the experiments. Coordination environments for Al–O are predominantly 4- and 5-coordinate, with a mean coordination of n AlO = 4.43(8), while Yb–O environments mostly range from 5- to 8-coordinate, with n YbO = 6.26(8). The cation–oxygen polyhedra are connected primarily by corner-sharing, with edge-sharing constituting up to ~1/3 of the connectivity among polyhedra with Yb or higher-coordinated Al–O. Structurally, the –Al–O– network in molten Yb 3 Al 5 O 12 appears conducive to glass formation: n OAl = 1.85(3), there are 1.86 AlO x –AlO x connections per Al atom (e.g., a mixture of Q 3 and Q 4 units), and the modal ring size is six cations. These characterize a network that is somewhat less constrained compared to SiO 2 glass, yet Yb 3 Al 5 O 12 cannot be quenched into crystal-free glass. Here, aluminum garnet compositions with larger rare earth cations do form glass, so these characterizations help reveal the structural characteristics corresponding to the limit of glass forming ability in rare earth aluminates.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Implicit and coupled fluid plasma solver with adaptive Cartesian mesh and its applications to non-equilibrium gas discharges

In this work, we present a new fluid plasma solver with adaptive Cartesian mesh (ACM) based on a full-Newton (nonlinear, implicit) scheme for non-equilibrium gas discharge plasma. The electrons and ions are described using drift-diffusion approximation coupled to Poisson equation for the electric field. The electron-energy transport equation is solved to account for electron thermal conductivity, Joule heating, and energy loss of electrons in collisions with neutral species. The rate of electron-induced ionization is a function of electron temperature and could also depend on electron density (important for plasma stratification). The ion and gas temperature are kept constant. The transport equations are discretized using a non-isothermal Scharfetter-Gummel scheme to resolve possible large temperature gradients in the sheaths. We demonstrate the new solver for simulations of direct current (DC) and radiofrequency (RF) discharges. The implicit treatment of the coupled equations allows using large time steps. The full-Newton method (FNM) enables fast nonlinear convergence at each time step, offering significantly improved simulation efficiency. We discuss the selection of time steps for solving different plasma problems. The new solver enables solving several problems we could not solve before with existing software: two- and three-dimensional structures of the entire DC discharges including cathode and anode regions, electric field reversals and double-layer formation, the normal cathode spot and an anode ring, moving striations in diffuse and constricted DC discharges, and standing striations in RF discharges. The developed FNM-ACM technique offers many benefits for tackling the disparity of gas discharge plasma systems' time scales and nonlinearity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Thermal fluctuations in the Landau-Lifshitz-Bloch model

A formulation for thermal noise in the stochastic form of the Landau-Lifshitz-Bloch equation used for modeling the magnetization dynamics at elevated temperatures is presented here. The diffusion coefficients for thermal fluctuations are obtained via the Fokker-Plank equation using the mean-field approximation of the field in defining the free energy. The presented model leads to a mean magnetization consistent with the equilibrium magnetization for small and large particles. The distribution of the magnetization magnitude is of the Maxwell-Boltzmann type. The presented model was tested by studying the equilibrium magnetization in macrospin particles at high temperatures. The model is appealing for multiscale modeling, such as modeling heat assisted magnetic recording systems and all-optical magnetization reversal.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Transient Modeling and Simulation of a Generic Stable Salt Reactor

A SAM system-level model of a generic stable salt reactor has been developed to investigate thermal-hydraulic behavior and safety performance under steady and transient conditions. The model integrates information generated from a reactor physics analysis using PROTEUS and PERSENT, and a computation fluid dynamics (CFD) analysis using STAR-CCM+. A loose, iterative coupling scheme between PROTEUS and SAM is implemented to calculate the equilibrium power and temperature distributions in the steady-state critical core condition. The converged steady-state model is then used in PERSENT to calculate the four reactivity feedback temperature coefficients (Doppler, fuel density, coolant density, and core radial expansion) and kinetic parameters that are needed in SAM to model the temperature feedback effects in transient simulations. Within the fully enclosed liquid fuel pins, natural convection is the dominant heat transfer mechanism. The STAR-CCM+ model of the fuel pin considers conjugate heat transfer from the liquid fuel salt to the pin cladding and external reactor coolant. The CFD results of the axial and radial temperature profiles are used to empirically determine an effective fuel salt thermal conductivity in the SAM fuel pin model so that the temperatures predicted by the SAM model match as closely as possible the CFD results. In the central region of the fuel pin, the effective thermal conductivity is as high as similar to 60 times the physical fuel salt thermal conductivity. The whole-plant SAM model is then used to simulate an unprotected station blackout transient. The results of this simulation showed that the large negative fuel axial expansion reactivity feedback reduces fission power to similar to 2.4% nominal power. The core is cooled by natural circulation, which removes heat in the core to the emergency heat removal system, and ultimately, to the ambient. However, peak fuel salt and cladding temperatures can potentially reach as high as 1500 K, albeit briefly, if the shutdown mechanism fails to operate.

stable salt reactor; transient simulations; system↗