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 91 records · Page 5

A whole-core steady-state thermal-hydraulic model for annular fuel type fluoride-salt-cooled reactors

A whole-core, steady-state thermal-hydraulic model is developed for the fluoride-salt-cooled small modular advanced high-temperature reactor (SmAHTR) that employs an annular fuel configuration. This pre-conceptual reactor design by Oak Ridge National Laboratory (ORNL) has the annular fuel and moderator pins arranged in a hexagonal layout. The FLiBe coolant flows from the bottom to the top of the core, parallel to the hexagonal bundle. The fuel and moderator pins in the core are discretized into finite volumes and the 3-D heat conduction equation is solved to obtain the temperature profile. Inter-fuel assembly conduction is also addressed. For this fuel assembly configuration, the coolant flows through two distinct regions – the hexagonal pin bundle and the annulus between the fuel pin and the tie rod. The fluid flow through the hexagonal bundles is modeled using the subchannel approach, in which the coolant region is discretized into corner, edge and interior subchannels and the resulting conservation equations are systematically solved. The 1-D mass, momentum and energy equations are solved for the annulus channels between the fuel pin and the tie rod. Pertinent closure models from the literature are employed to close the system of equations. We also performed a preliminary code-to-code comparison between the present model and a CFD model.. The resulting thermal-hydraulic model can provide temperature, flow rate and pressure drop profiles for the different solid and fluid regions throughout the entire core. Whole-core thermal-hydraulic results for a representative power profile are presented and discussed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Thermal State-of-Charge in Solar Heat Receivers

A theoretical framework is developed to determine the so-called thermal state-of-charge (SOC) in solar heat receivers employing encapsulated phase change materials (PCMS) that undergo cyclic melting and freezing. The present problem is relevant to space solar dynamic power systems that would typically operate in low-Earth-orbit (LEO). The solar heat receiver is integrated into a closed-cycle Brayton engine that produces electric power during sunlight and eclipse periods of the orbit cycle. The concepts of available power and virtual source temperature, both on a finite-time basis, are used as the basis for determining the SOC. Analytic expressions for the available power crossing the aperture plane of the receiver, available power stored in the receiver, and available power delivered to the working fluid are derived, all of which are related to the SOC through measurable parameters. Lower and upper bounds on the SOC are proposed in order to delineate absolute limiting cases for a range of input parameters (orbital, geometric, etc.). SOC characterization is also performed in the subcooled, two-phase, and superheat regimes. Finally, a previously-developed physical and numerical model of the solar heat receiver component of NASA Lewis Research Center's Ground Test Demonstration (GTD) system is used in order to predict the SOC as a function of measurable parameters.

Hall, Carsie A., Jr.↗

Thermal State-of-Charge in Solar Heat Receivers

A theoretical framework is developed to determine the so-called thermal state-of-charge (SOC) in solar heat receivers employing encapsulated phase change materials (PCMS) that undergo cyclic melting and freezing. The present problem is relevant to space solar dynamic power systems that would typically operate in low-Earth-orbit (LEO). The solar heat receiver is integrated into a closed-cycle Brayton engine that produces electric power during sunlight and eclipse periods of the orbit cycle. The concepts of available power and virtual source temperature, both on a finite-time basis, are used as the basis for determining the SOC. Analytic expressions for the available power crossing the aperture plane of the receiver, available power stored in the receiver, and available power delivered to the working fluid are derived, all of which are related to the SOC through measurable parameters. Lower and upper bounds on the SOC are proposed in order to delineate absolute limiting cases for a range of input parameters (orbital, geometric, etc.). SOC characterization is also performed in the subcooled, two-phase, and superheat regimes. Finally, a previously-developed physical and numerical model of the solar heat receiver component of NASA Lewis Research Center's Ground Test Demonstration (GTD) system is used in order to predict the SOC as a function of measurable parameters.

Hall, Carsie, A., III↗

Steady-State Thermal-Hydraulic Analysis of the LEU-Fueled Dalat Nuclear Research Reactor

This paper presents results of steady-state thermal-hydraulic analysis for the designed working core of the Dalat Nuclear Research Reactor (DNRR) using the PLTEMP/ANL code. The core was designed to be loaded with 92 low-enriched uranium (LEU) VVR-M2 fuel bundles (FBs) and 12 beryllium rods surrounding a neutron trap at the core center, for replacement of the previous core with 104 high-enriched uranium (HEU) VVR-M2 FBs. Before using this code for thermohydraulic analysis of the designed LEU working core, it was validated by comparing calculation results with experimental data collected from the HEU working core of the DNRR. The discrepancy between calculated results and measured data was at the maximum about 0.8°C and 1.5°C of fuel cladding and outlet coolant temperatures, respectively. In the design calculation, thermohydraulic safety was confirmed through evaluation of the fuel cladding and coolant temperatures, as well as of other safety parameters such as Departure from Nucleate Boiling Ratio (DNBR) and Onset of Nucleate Boiling Ratio (ONBR). The calculation results showed that, in normal operation conditions at full nominal thermal power of 500 kW without uncertainty parameters, the maximum fuel cladding temperature of the hottest FB was about 90.4°C, which is lower than its limit value of 103°C, the minimum DNBR was 32.0, which is much higher than the recommended value of 1.5, and the minimum ONBR was 1.43, which is higher than the recommended value of 1.4 for VVR-M2 LEU fuel type. When the global and local hot channel factors were taken into account, the maximum temperature of fuel cladding at the hottest FB was about 98.4 °C, for global only, and 114.3°C, for global together with local hot channel factors. The calculation results confirm the safety operation of the designed LEU core loaded with 92 fresh VVR-M2 FBs.

Nguyen, Kien-Cuong↗

Dominant effect of anharmonicity on the equation of state and thermal conductivity of MgO under extreme conditions

Understanding heat transfer in the earth mantle region is of great scientific interest but highly challenging due to dissimilar conditions from our daily lives. Thermal conductivity of highly pressured and thermally activated materials can be completely different from that at ambient condition. In this study we calculate the equation of state and thermal conductivity of cubic crystalline MgO, a major component of the mantle, in such an environment. It is shown that the material properties are not accurately captured unless rigorous treatment of anharmonicity of phonon vibrational motion is explicitly applied. To predict renormalized phonon dispersion relations, interatomic force constants of high-order expansion are required, especially at elevated temperature. The anharmonicity plays an increasingly substantial role for regulating volume expansion as temperature rises, because the Gibbs free energy of MgO considerably deviates from what was predicted by the quasiharmonic approximation. We find that the lattice thermal conductivity of MgO and the heat transfer of the earth's mantle are reliably estimated only when both the temperature-dependent renormalization of phonon frequencies and four-phonon scattering effects are taken into consideration. We map the characteristics at Geotherm into contour diagrams in which equilibrium volumes and the thermal conductivities are described as functions of temperature and pressure. Furthermore, the thermal conductivity values of MgO in this study are smaller than those in previous theoretical studies due to the additional incorporation of four-phonon scatterings in the calculations. Specifically, deviation between the values in previous studies and in this study becomes larger as thermodynamic variables reach those near the core-mantle boundary. For instance, thermal conductivities of 52.15 and 23.67 Wm -1 K -1 are estimated at ambient and extreme conditions (T = 4000K and P = 135 GPa), respectively.

36 MATERIALS SCIENCE↗

Application of Satellite Frost Forecast Technology to Other Parts of the United States

Thermal infrared data taken from the GOES satellite over a period of several hours was color enhanced by computer according to temperature. The varying temperatures were then used to assist in frost forecasting. Input from Michigan and Pennsylvania to the cold climate mapping project is emphasized in the report of the second year's activities of a two year effort.

Martsolf, J. D.↗

A review on the state of thermal hydraulics research on air ingress scenarios in High-Temperature Gas-cooled Reactors following a D-LOFC

With the expectation of near-immediate carbon neutrality, widespread implementation of proven High-Temperature Gas-cooled Reactors (HTGRs) embodies a viable solution pathway given their inherent, passive safety features and high thermal efficiency. This study provides an overview of the current state of research involving the thermal hydraulics associated with air ingress from a depressurized loss of forced cooling (D-LOFC) in HTGRs. Accurately characterizing and predicting the physical phenomena underlying air ingress is of paramount concern, as the integrity of the fuel and core graphite support structures are threatened by the presence of oxygen. Broadly speaking, the air ingress scenario can be delineated into three main stages: (1) Depressurization, (2) Density-Driven Flow, and (3) Natural Convection. In tandem with the underlying fundamental theory, this review collates and synthesizes the existing body of contemporary research concerning the air ingress scenario following a D-LOFC. As evinced by this review, our current understanding and predictive abilities have benefited from extensive research, predominantly concentrated on the rate of air ingestion into the core. Here, additional research is necessary to holistically capture the phenomenology of an air ingress scenario following a D-LOFC by considering an additional variable: the oxygen content of the ingressing air. The latter variable requires investigation into the complex interactions of the fully integrated system. Additionally, while numerical tools are evolving domestically through the Nuclear Energy Advanced Modeling and Simulation program, a sufficiently validated code remains absent.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Applications of an Energy Transfer Model to Three Problems in Planetary Regoliths: The Solid-State Greenhouse, Thermal Beaming, and Emittance Spectra

Several problems of interest in planetary infrared remote sensing are investigated using a new radiative-conductive model of energy transfer in regoliths: the solid-state greenhouse effect, thermal beaming, and reststrahlen spectra. The results of the analysis are as follows: (1) The solid-state greenhouse effect is self-limiting to a rise of a few tens of degrees in bodies of the outer solar system. (2) Non-Lambertian directional emissivity can account for only about 20% of the observed thermal beaming factor. The remainder must have another cause, presumably surface roughness effects. (3) The maximum in a reststrahlen emissivity spectrum does not occur exactly at the Christiansen wavelength where, by definition, the real part of the refractive index equals one, but rather at the first transition minimum in reflectance associated with the transition from particle scattering being dominated by volume scattering to that dominated by strong surface scattering. The transparency feature is at the second transition minimum and does not require the presence of a second band at longer wavelength for its occurance. Subsurface temperature gradients have only a small effect on emissivity bands.

Hapke, Bruce↗

Observation of solid-state bidirectional thermal conductivity switching in antiferroelectric lead zirconate (PbZrO 3 )

Materials with tunable thermal properties enable on-demand control of temperature and heat flow, which is an integral component in the development of solid-state refrigeration, energy scavenging, and thermal circuits. Although gap-based and liquid-based thermal switches that work on the basis of mechanical movements have been an effective approach to control the flow of heat in the devices, their complex mechanisms impose considerable costs in latency, expense, and power consumption. As a consequence, materials that have multiple solid-state phases with distinct thermal properties are appealing for thermal management due to their simplicity, fast switching, and compactness. Thus, an ideal thermal switch should operate near or above room temperature, have a simple trigger mechanism, and offer a quick and large on/off switching ratio. In this study, we experimentally demonstrate that manipulating phonon scattering rates can switch the thermal conductivity of antiferroelectric PbZrO 3 bidirectionally by –10% and +25% upon applying electrical and thermal excitation, respectively. Our approach takes advantage of two separate phase transformations in PbZrO 3 that alter the phonon scattering rate in different manners. In this study, we demonstrate that PbZrO 3 can serve as a fast (<1 second), repeatable, simple trigger, and reliable thermal switch with a net switching ratio of nearly 38% from ~1.20 to ~1.65 W m –1 K –1 .

25 ENERGY STORAGE↗

Performance assessment of active insulation systems in residential buildings for energy savings and peak demand reduction

Active insulation systems (AISs) in buildings are envelopes that integrate thermal insulation, thermal energy storage, and controls. Although different designs for AISs have been proposed in the literature, a comprehensive analysis of feasible AISs is lacking. This paper discusses the energy performance, peak demand reduction potential, and performance characteristics of an AIS that uses a concrete wall as thermal mass sandwiched between two solid-state thermal switches (STSs). These STSs change their thermal conductivity using an on/off metal switch to create or break a thermal bridge across the STS. This paper first describes the experimental setup, used to determine the ratio of thermal resistance during R-high (low thermal conductivity) and R-low (high thermal conductivity) states of the STSs. This ratio was then used in whole-building energy simulations to evaluate the performance of AIS walls across different climate zones with/without a freeze timer of 60 min. The timer was added to reduce the number of switches of STSs from one state to another, and hence the energy needed for these switches. Analysis of the switching frequency and interval of STSs, thermal conductivity of walls, impact of wall orientation, and heat transfer through the wall from the use of AIS at different climate zones/locations were performed. The simulation results show that the AIS can achieve energy savings ranging from ~980 to 2,290 kWh in a single-family home with a floor area of ~220 m 2 compared with an IECC 2018 baseline. We found the energy savings was higher in dry climate zones which represent 17% of residential buildings in the United States, compared to humid or marine climate.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Characterization of SiC/SiC Composites Subjected to Sequential Tensile Creep and Dwell Fatigue in Air under Thermal Gradient Conditions

Uncoated SiC/SiC composites, with 2D and 3D fiber architectures, were investigated under steady state thermal gradient conditions with sequential tensile creep and dwell fatigue. Melt infiltration (MI), chemical vapor infiltration (CVI), polymer impregnation and pyrolysis (PIP), and hybrid (CVI and PIP) techniques were used to fabricate ceramic matrix composites (CMCs) with a 2D balanced fiber architecture. CMCs with 3D fiber weaves were also fabricated using MI and hybrid methods. Steady state through the thickness thermal gradients were generated for the tensile specimens using laser heating on the front (hot side) and forced air cooling on the back (cold side). Hot side and cold side temperatures of 2400 and 2200 °F, respectively, were used for testing CMCs fabricated with the MI process, whereas hot side and cold side temperatures of 2700 and 2400 °F, respectively, were used for CMCs fabricated with CVI, PIP, and hybrid processes. Each CMC specimen was initially tested under tensile creep at 69 MPa (10 ksi) for 80 hours, which was then followed by tensile dwell fatigue testing at a maximum stress of 69 MPa (10 ksi) for another 80 hours for a total of 160 hours. This type of sequential tensile creep followed by tensile dwell fatigue testing was performed to ascertain the integrity of uncoated CMCs under steady state thermal gradient conditions. Tensile creep and dwell fatigue test results, generated under steady state thermal gradient conditions on different types of SiC/SiC composites, are compared with corresponding creep behaviors observed under isothermal conditions to determine the influence of both monotonic creep (static fatigue) and dwell fatigue on accumulated strains and deformation behaviors of the CMCs.

Kalluri, Sreeramesh↗

Classical and quantum computing of shear viscosity for ( 2 + 1 ) D SU(2) gauge theory

We perform a nonperturbative calculation of the shear viscosity for ( 2 + 1 )-dimensional SU(2) gauge theory by using the lattice Hamiltonian formulation. The retarded Green’s function of the stress-energy tensor is calculated from real time evolution via exact diagonalization of the lattice Hamiltonian with a local Hilbert space truncation, and the shear viscosity is obtained via the Kubo formula. When taking the continuum limit, we account for the renormalization group flow of the coupling but no additional operator renormalization. We find the ratio of the shear viscosity and the entropy density η s is consistent with a well-known holographic result 1 4 π at several temperatures on a 4 × 4 honeycomb lattice with the local electric representation truncated at j max = 1 2 . We also find the ratio of the spectral function and frequency ρ x y ( ω ) ω exhibits a peak structure when the frequency is small. Both the exact diagonalization method and simple matrix product state classical simulation method beyond j max = 1 2 on bigger lattices require exponentially growing resources. So we develop a quantum computing method to calculate the retarded Green’s function and analyze various systematics of the calculation including j max truncation and finite size effects, Trotter errors and the thermal state preparation efficiency. Our thermal state preparation method still requires resources that grow exponentially with the lattice size, but with a very small prefactor at high temperature. We test our quantum circuit on both the Quantinuum emulator and the IBM simulator for a small lattice and obtain results consistent with the classical computing ones. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Separation of infrared and bulk in thermal QCD

A new thermal regime of QCD, featuring decoupled scale-invariant infrared glue, has been proposed to exist both in pure-glue (N f =0) and “real-world” (N f =2+1 at physical quark masses) QCD. In this IR phase, elementary degrees of freedom flood the infrared, forming a distinct component independent from the bulk. This behavior necessitates non-analyticities in the theory. In pure-glue QCD, such non-analyticities have been shown to arise via Anderson-like mobility edges in Dirac spectra (λ IR = 0, ±λ A ≠ 0), as manifested in the dimension function d IR (λ). Here, we present the first evidence, based on lattice QCD calculation at a=0.105 fm, that this mechanism is also at work in real-world QCD, thus supporting the existence of the proposed IR regime in nature. An important aspect of our results is that, while at T = 234 MeV we find a dimensional jump between zero modes and lowest near-zero modes very close to unity (d IR = 3 to d IR ≃ 2), similar to the IR phase of pure-glue QCD, at T = 187 MeV we observe a continuous λ-dependence. This suggests that thermal states just above the chiral crossover are non-analytically (in T) connected to thermal state at T = 234 MeV, supporting the key original proposition that the transition into the IR regime occurs at a temperature strictly above the chiral crossover.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Quantum Markov chain Monte Carlo with digital dissipative dynamics on quantum computers

Modeling the dynamics of a quantum system connected to the environment is critical for advancing our understanding of complex quantum processes, as most quantum processes in nature are affected by an environment. Modeling a macroscopic environment on a quantum simulator may be achieved by coupling independent ancilla qubits that facilitate energy exchange in an appropriate manner with the system and mimic an environment. This approach requires a large, and possibly exponential number of ancillary degrees of freedom which is impractical. In contrast, we develop a digital quantum algorithm that simulates interaction with an environment using a small number of ancilla qubits. By combining periodic modulation of the ancilla energies, or spectral combing, with periodic reset operations, we are able to mimic interaction with a large environment and generate thermal states of interacting many-body systems. We evaluate the algorithm by simulating preparation of thermal states of the transverse Ising model. Our algorithm can also be viewed as a quantum Markov chain Monte Carlo process that allows sampling of the Gibbs distribution of a multivariate model. To illustrate this we evaluate the accuracy of sampling Gibbs distributions of simple probabilistic graphical models using the algorithm.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Exploring finite temperature properties of materials with quantum computers

Abstract Thermal properties of nanomaterials are crucial to not only improving our fundamental understanding of condensed matter systems, but also to developing novel materials for applications spanning research and industry. Since quantum effects arise at the nano-scale, these systems are difficult to simulate on classical computers. Quantum computers can efficiently simulate quantum many-body systems, yet current quantum algorithms for calculating thermal properties of these systems incur significant computational costs in that they either prepare the full thermal state on the quantum computer, or they must sample a number of pure states from a distribution that grows with system size. Canonical thermal pure quantum (TPQ) states provide a promising path to estimating thermal properties of quantum materials as they neither require preparation of the full thermal state nor require a growing number of samples with system size. Here, we present an algorithm for preparing canonical TPQ states on quantum computers. We compare three different circuit implementations for the algorithm and demonstrate their capabilities in estimating thermal properties of quantum materials. Due to its increasing accuracy with system size and flexibility in implementation, we anticipate that this method will enable finite temperature explorations of relevant quantum materials on near-term quantum computers.

36 MATERIALS SCIENCE↗

Recent progress on in-medium properties of heavy mesons from finite-temperature EFTs

Mesons with heavy flavor content are an exceptional probe of the hot QCD medium produced in heavy-ion collisions. In the past few years, significant progress has been made toward describing the modification of the properties of heavy mesons in the hadronic phase at finite temperature. Ground-state and excited-state thermal spectral properties can be computed within a self-consistent many-body approach that employs appropriate hadron-hadron effective interactions, providing a unique opportunity to confront hadronic Effective Field Theory predictions with recent and forthcoming lattice QCD simulations and experimental data. In this article, we revisit the application of the imaginary-time formalism to extend the calculation of unitarized scattering amplitudes from the vacuum to finite temperature. These methods allow us to obtain the ground-state thermal spectral functions. The thermal properties of the excited states that are dynamically generated within the molecular picture are also directly accessible. We present here the results of this approach for the open-charm and open-bottom sectors. We also analyze how the heavy-flavor transport properties, which are strongly correlated to experimental observables in heavy-ion collisions, are modified in hot matter. In particular, transport coefficients can be computed using an off-shell kinetic theory that is fully consistent with the effective theory describing the scattering processes. The results of this procedure for both charm and bottom transport coefficients are briefly discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗