Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermal response functions”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Modeling cell venting and gas-phase reactions in 18650 lithium ion batteries during thermal runaway

A numerical model is developed to study cell venting, internal pressure, and gas-phase dynamics behavior of 18650 Li-ion cells undergoing thermal runaway. A k-e Reynolds-Averaged Navier-Stokes (RANS) model is adopted to describe the turbulent flow out of the cells, while the fluid dynamics inside the cells is described by Darcy-Forchheimer's equation. Thermal abuse reactions and gas generation kinetics are described by a single-step lumped reaction model. Then, a series of computational fluid dynamics (CFD) simulations are conducted on a single 18650 cell at various states-of-charge (100%, 50%, 25%) to study detailed flow and thermal behavior as a function of quantity of gas generated during cell venting. Venting events are categorized into two stages: i) breaching of the cell container and ii) thermal runaway reactions. It is found that the cell response is dominated by the second stage since most of the gases are generated during thermal runaway. Also, the propensity for propagation is highly affected by state-of-charge (SOC). Cells at higher SOCs produce more heat and gas during the venting event, owing to higher mass and concentrations of reacting gases, and consequently reach higher internal cell pressures which increase the risk of side-wall breaching.

25 ENERGY STORAGE↗
DOE OSTI2020

Gas Generation Models for Thermal Runaway in Lithium Ion Batteries

A numerical model is implemented to study cell venting, internal pressure, and gas-phase dynamics behavior of Li-ion cells undergoing thermal runaway. A k-ɛ Reynolds-Averaged Navier-Stokes (RANS) model is adopted to describe the turbulent flow out of the cells, while the fluid dynamics inside the cells is described by Darcy-Forchheimer’s equation. Thermal abuse reactions and gas generation kinetics are described by a single-step lumped reaction model. A series of computational fluid dynamics (CFD) simulations are conducted on a single cell at various states-of-charge (100%, 50%, 25%) to demonstrate detailed flow and thermal behavior as a function of quantity of gas generated during cell venting. Venting events are categorized into two stages: i) breaching of the cell container and ii) thermal runaway reactions. It is shown that the cell response is dominated by the second stage since most of the gases are generated during thermal runaway. Also, the propensity for propagation is highly affected by state-of-charge (SOC). Cells at higher SOCs produce more heat and gas during the venting event, owing to higher mass and concentrations of reacting gases, and consequently reach higher internal cell pressures which increase the risk of side-wall breaching.

Kim, Jinyong↗
DOE OSTI2020

Optical properties of cubic boron arsenide

The ultrahigh thermal conductivity of cubic boron arsenide (BAs) makes it a promising material for next-generation electronics and optoelectronics. Here, we report measured optical properties of BAs crystals, including the complex dielectric function, refractive index, and absorption coefficient in the ultraviolet, visible, and near-infrared wavelength range. In this work, the data were collected at room temperature using spectroscopic ellipsometry and transmission and reflection spectroscopy. We further calculated the optical response using density functional theory and many-body perturbation theory, considering quasiparticle and excitonic corrections. The computed values for the direct and indirect bandgaps (4.25 eV and 2.07 eV) agree well with the measured results (4.12 eV and 2.02 eV). Our findings pave the way for using BAs in future electronic and optoelectronic applications that take advantage of its demonstrated ultrahigh thermal conductivity and predicted high ambipolar carrier mobility.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Thermomechanical properties and microstructures of yttrium hydride

Yttrium hydride is an optimal choice for a high-temperature moderator material in advanced thermal neutron spectrum reactors that require small core volumes. However, a complete database of the thermomechanical properties of yttrium hydride is not available yet, although it is much needed to understand and predict the moderator performance during service in reactors. In this paper, we report the properties of unirradiated bulk yttrium hydride as a function of hydrogen concentration—including density, crystal structure, specific heat capacity, thermal diffusivity, thermal conductivity, hardness, elastic/shear moduli, Poisson’s ratio, fracture strength, microstructure, and thermal stability—providing a baseline measurement for the subsequent neutron irradiation response study of yttrium hydride. Furthermore, the recommended empirical treatment of the data is suggested. In addition, other properties (i.e., hydrogen retention, thermal hydrogen migration, and irradiation response) that needs to be investigated are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2020

CISME report on experiments with HMX powder, PBX 9501 prills, and heated PBX 9501

The ability to predict post-ignition explosive response is relevant to explosives safety, particularly weapons assembly/disassembly operations that must operate under a sufficiently conservative safety basis to meet risk mitigation criteria for inadvertent nuclear detonation (IND). In particular, hazard determinations hinge on whether the PBX 9501 main charge can ultimately attain detonation from an initial, non-shock, ignition stimulus. The CISME experiment uses a spherical explosive charge, which is thermally ignited in the center using a fiber-delivered infrared laser pulse. This ignition mechanism - i.e. a localized thermal hotspot - is representative of the ignition-causing insults typical in accident scenarios, including low-velocity impact, drop, skid, and punch events. In fiscal year 2018 CISME experiments were conducted with pristine PBX 9501 spheres, confined and unconfined. Sub-detonative, low-violence reactions were observed. Here we report experiments conducted since 2018. The explosive charge was varied: experiments were conducted with room temperature pour-density HMX powder and prills of PBX 9501, and heated, thermally damaged PBX 9501. Experiments were conducted in a confined configuration, with two 4 mm-thick stainless-steel hemispherical shells bolted around the sphere. Detonation, via DDT (deflagration-to-detonation transition) was observed in the 3-inch diameter powdered HMX and the 6-inch diameter heated and thermally-damaged PBX 9501 tests. No transition to DDT was observed in smaller diameters for those configurations, and no DDT was observed for pour-density prills of PBX 9501 of any diameter. The original test series demonstrated that the response of pristine, undamaged PBX 9501 to a central ignition is relatively benign. This new test series shows that increased porosity, lack of polymer binders and thermal damage can greatly increase the violence of the post-ignition response, up to and including detonation. These results bridge the explosive response of PBX 9501 from sub-detonative through detonation as a function of diameter and provide valuable test cases for modeling efforts that strive to capture the DDT threshold. By exploring explosive response in a more representative 3D geometry, we can also benchmark these new results to more plentiful DDT datasets generated from testing in 1D tube confinement (the historical test configuration)–the configuration upon which almost all DDT models were developed.

42 ENGINEERING↗
DOE OSTI2025

Phonon Olympics: Phonon property and lattice thermal conductivity benchmarking from open-source packages

Three widely used open-source packages for determining phonon properties and lattice thermal conductivities (ALAMODE, phono3py, and ShengBTE) are benchmarked by teams of expert users and the package developers. The phonons for Ge, RbBr, monolayer MoSe 2 , and AlN are modeled at zero temperature, and they scatter through three-phonon and phonon-isotope processes, with thermal conductivities obtained from the linearized Peierls–Boltzmann transport equation with input from density functional theory calculations. Over a wide range of temperatures, the thermal conductivities calculated by the teams fall within at most ±15% of their mean values for each of the four materials. The phonon frequencies, obtained from the harmonic force constants, do not show large differences between the calculations, indicating that the modal heat capacities and group velocities are not responsible for the thermal conductivity variations. It is the lifetimes associated with three-phonon scattering, obtained from the cubic force constants, that drive the variations. The many decisions required to calculate the cubic force constants (e.g., supercell size, atomic displacement, neighbor cutoff, and application of symmetries) make identification of the precise origin of the thermal conductivity variations challenging. The calculated thermal conductivities do not generally show agreement with experimental measurements, which is attributed to the limitations of the density functional theory calculations. Guidance for the development of best practices is provided, which will help to standardize protocols needed for building thermal conductivity databases. The results provide a baseline for future benchmarking of other packages and more advanced calculations.

McGaughey, Alan J. H. [Carnegie Mellon Univ., Pitt↗
DOE OSTI2025

Dual-season thermal energy storage-integrated heat pump system for load shifting: experiments, modeling and nationwide analysis

Here, this paper proposes and validates a dual-season thermal energy storage-integrated heat pump (TES-HP) system that shifts building electric load in both cooling and heating modes using a single phase-change TES. The design uses the TES as a heat sink for the condenser in summer and a heat source for the evaporator in winter, thereby reducing on-peak electric demand while ensuring compatible with existing air-distribution systems. A 14-kW TES-HP prototype with six operating modes and novel refrigerant charge management strategies was experimentally tested under varied conditions. Subsequently, data-driven polynomial performance curves were developed and validated against experiments, then coupled with a phase change material (PCM) model and a calibrated building model within a rule-based controller. Single-day case studies demonstrated effective on-peak demand reduction while maintaining thermal comfort. During cooling, hourly power was lowered by up to 1.5 kW; in heating, average hourly power decreased from over 3.5 kW to 1.8 kW, compared to an HP-only system, achieving 64.5 % load shifting, and the need for resistance heating was eliminated. Seasonal simulations showed typical on-peak electricity savings of 25–35 % in cooling and 40–65 % in heating, with the largest benefits on the hottest and coldest days. Extended response-surface analysis and nationwide mapping quantified load shifting as functions of ambient and TES temperatures, revealing a seasonal trade-off in TES phase-change temperature. These results demonstrate the TES-HP system as a practical and scalable solution for grid-interactive buildings that reduces on-peak demand and electricity use while maintaining thermal comfort and enhancing grid flexibility and reliability.

Grid-interactive buildings↗
DOE OSTI2021

Pressure, temperature, and orientation dependent thermal conductivity of pentaerythritol tetranitrate (PETN)

We use reverse non-equilibrium molecular dynamics to determine the thermal conductivity tensor, κ, of tetragonal pentaerythritol tetranitrate (PETN). The most stable form under standard ambient conditions (PETN I, with space group P4¯2 1 c) is considered in the temperature and pressure intervals 200–500 K and 0–8 GPa, which covers approximately the stability range for this polymorph. We compute κ along the [100] and [001] directions, which is enough to construct the full thermal conductivity tensor for a system with tetragonal symmetry. In addition, we also determine κ along the [010] direction and confirm that κ 100 ~ κ 010 , with less than 4% average absolute error between the two quantities. We observe an anisotropic response for κ, with κ 100 > κ 001 across the whole (T,P) interval, and 37% difference at 300 K, 0 GPa. Furthermore, we provide analytical functions to interpolate κ(P,T) within the fitting interval and serve as input for continuum-scale simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Predicting responses to climate change using a joint species, spatially dependent physiologically guided abundance model

Abstract Predicting the effects of warming temperatures on the abundance and distribution of organisms under future climate scenarios often requires extrapolating species–environment correlations to climatic conditions not currently experienced by a species, which can result in unrealistic predictions. For poikilotherms, incorporating species' thermal physiology to inform extrapolations under novel thermal conditions can result in more realistic predictions. Furthermore, models that incorporate species and spatial dependencies may improve predictions by capturing correlations present in ecological data that are not accounted for by predictor variables. Here, we present a joint species, spatially dependent physiologically guided abundance (jsPGA) model for predicting multispecies responses to climate warming. The jsPGA model uses a basis function approach to capture both species and spatial dependencies. We apply the jsPGA model to predict the response of eight fish species to projected climate warming in thousands of lakes in Minnesota, USA. By the end of the century, the cold‐adapted species was predicted to have high probabilities of extirpation across its current range—with 10% of lakes currently inhabited by this species having an extirpation probability >0.90. The remaining species had varying levels of predicted changes in abundance, reflecting differences in their thermal physiology. Though the model did not identify many strong species dependencies, the variation in estimated spatial dependence across species suggested that accounting for both dependencies was important for predicting the abundance of these fishes. The jsPGA model provides a new tool for predicting changes in the abundance, distribution, and extirpation probability of poikilotherms under novel thermal conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Gaussian-preserved, non-volatile shape morphing in three-dimensional microstructures for dual-functional electronic devices

Motile plant structures such as Mimosa pudica leaves, Impatiens glandulifera seedpods, and Dionaea muscipula leaves exhibit fast nastic movements in a few seconds or less. This motion is stimuli-independent mechanical movement following theorema egregium rules. Artificial analogs of tropistic motion in plants are exemplified by shape-morphing systems, which are characterized by high functional robustness and resilience for creating 3D structures. However, all shape-morphing systems developed so far rely exclusively on continuous external stimuli and result in slow response. Here, we report a Gaussian-preserved shape-morphing system to realize ultrafast shape morphing and non-volatile reconfiguration. Relying on the Gaussian-preserved rules, the transformation can be triggered by mechanical or thermal stimuli within a microsecond. Moreover, as localized energy minima are encountered during shape morphing, non-volatile configuration is preserved by geometrically enhanced rigidity. Using this system, we demonstrate a suite of electronic devices that are reconfigurable, and therefore, expand functional diversification.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Handbook on the Material Properties of Yttrium Hydride for High Temperature Moderator Applications

Yttrium hydride is an optimal choice for a high-temperature moderator material in advanced thermal neutron spectrum reactors that require small core volumes. A complete database of the thermomechanical properties of yttrium hydride is not yet available, although it is much needed to understand and predict the moderator performance during service in reactors. This milestone report presents the properties of unirradiated bulk yttrium hydride as a function of hydrogen concentration—including density, crystal structure, specific heat capacity, thermal diffusivity, thermal conductivity, hardness, elastic/shear moduli, Poisson’s ratio, fracture strength, microstructure, and thermal stability. This information provides a baseline measurement for the subsequent neutron irradiation response study of yttrium hydride. The recommended empirical treatment of the data is suggested. In addition, other properties (i.e., hydrogen retention, thermal hydrogen migration, and irradiation response) requiring investigation are discussed. Also included are the Preliminary post-irradiation examination (PIE) data of yttrium hydride irradiated at 600 and 900°C to 0.1 dpa in the High Flux Isotope Reactor (HFIR). The thermophysical properties have insignificant change following this low dpa irradiation for both irradiation temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Analytical solution and parameter estimation for heat of wetting and vapor adsorption during spontaneous imbibition in tuff

An analytical expression is derived for the thermal response observed during spontaneous imbibition of water into a dry core of zeolitic tuff. Sample tortuosity, thermal conductivity, and thermal source strength are estimated from fitting an analytical solution to temperature observations during a single laboratory test. Here, the closed-form analytical solution is derived using Green's functions for heat conduction in the limit of “slow” water movement; that is, when advection of thermal energy with the wetting front is negligible. The solution has four free fitting parameters and is efficient for parameter estimation. Laboratory imbibition data used to constrain the model include a time series of the mass of water imbibed, visual location of the wetting front through time, and temperature time series at six locations. The thermal front reached the end of the core hours before the visible wetting front. Thus, the predominant form of heating during imbibition in this zeolitic tuff is due to vapor adsorption in dry zeolitic rock ahead of the wetting front. The separation of the wetting front and thermal front in this zeolitic tuff is significant, compared to wetting front behavior of most materials reported in the literature. This work is the first interpretation of a thermal imbibition response to estimate transport (tortuosity) and thermal properties (including thermal conductivity) from a single laboratory test.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Topological Control of Triply Periodic Minimal Surfaces for Thermal Design and Advanced Manufacturing: A Gyroid Case Study

Recently, there has been a heightened interest in using triply periodic minimal surfaces (TPMSs) in the design of compact process engineering components. The benefits of high surface area per unit volume, modular form, and inherent periodicity provide a holistic self-supporting network and flow-conducive features. Applications of importance include thermal power management, biomimetic scaffolds and structures, and feasibility of advanced manufacturing. This study presents a novel approach to the manipulation of the characteristic Schwarz-G, or gyroid TPMS, for thermal design in the context of advanced manufacturing. The study presents relationships between design parameters and resulting surface area as a target response using the characteristic equation of a gyroid. Through parametric control, the characteristic equation is manipulated to produce a 20-fold increase in achievable area over a baseline design characteristic of 25.4 mm through controlled combinations of design parameters. A second relationship is presented as a function of the maximum area achieved and manipulated design parameters. Through the analysis, the study presents a framework to identify and maximize the achievable area of TPMSs for advanced manufacturing and thermal management applications.

gyroid↗
DOE OSTI2024

An Exact Expression for Multidimensional Spectroscopy of a Spin‐Boson Hamiltonian

Multidimensional coherent spectroscopy is a powerful tool to characterize nonlinear optical response functions. Typically, multidimensional spectra are interpreted via a perturbative framework that straightforwardly provides intuition into the density matrix dynamics that give rise to specific spectral features. When the goal is to characterize system coupling to a thermal bath however, the perturbative formalism becomes unwieldy and yields less intuition. Here, an approach developed by Vagov et al. is extended to provide an exact expression for multidimensional spectra of a spin-boson Hamiltonian up to arbitrary order of electric field interaction. Furthermore, the utility of this expression is demonstrated by modeling polaron formation and coherent exciton-phonon coupling in quantum dots, which strongly agree with experiment.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Precision Loading and Delivery of Molecular Cargo by Size-Controlled Coacervation of Gold Nanoparticles Functionalized with Elastin-like Peptides

Thermoresponsive elastin-like peptides (ELPs) have been extensively investigated in biotechnology and medicine, but little attention has been paid to the process by which coacervation causes ELP-decorated particles to aggregate. Here, using gold nanoparticles (AuNPs) functionalized with a cysteine-terminated 96-repeat of the VPGVG sequence (V96-Cys), we show that the size of the clusters that reversibly form above the ELP transition temperature can be finely controlled in the 250 to 930 nm range by specifying the concentration of free V96-Cys in solution and using AuNPs of different sizes. We further find that the localized surface plasmon resonance peak of the embedded AuNPs progressively red-shifts with cluster size, likely due to an increase in particle–particle contacts. We exploit this fine control over size to homogeneously load precise amounts of the dye Nile Red and the antibiotic Tetracycline into clusters of different hydrodynamic diameters and deliver cargos near-quantitatively by deconstructing the aggregates below the ELP transition temperature. Beyond establishing a key role for free ELPs in the agglomeration of ELP-functionalized particles, our results provide a path for the thermally controlled delivery of precise quantities of molecular cargo. This capability might prove useful in combination photothermal therapies and theranostic applications, and to trigger spatially and temporally uniform responses from biological, electronic, or optical systems.

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

Dielectric Resonator Design for Low Power and Low Temperature Microwave Plasma

Waveguide-based microwave plasmas generally operate at high temperatures (2000 - 6000K)[1], making it difficult to directly interface solid materials with the plasma without significant thermal damage. Dielectric microwave resonators (DMRs), long studied for wave-based manipulation of electromagnetic radiation for telecom and optics, can focus radiation to extremely small mode volumes, creating intense localized fields with low-power input.[2] This phenomenon can be used for applications ranging from efficient plasma electronics to near-ambient plasma-materials interactions. Such DMR-based plasmas have been demonstrated a handful of times in the literature, but the majority of research towards this utilize the lowest frequency resonance mode.[3], [4], [5] By carefully controlling the geometry of cylindrical resonators, a variety of electromagnetic modes can be excited. In this work, COMSOL Multiphysics simulations are used to study the electric field enhancement and absorption properties of CaTiO3 DMRs as a function of geometry and excitation frequency. Whereas previous studies have utilized the HEM111 resonance frequency to drive low power plasma excitation, we find that higher order resonance frequencies are more effective at field enhancement and result in less power loss within the dielectric material, hence less wasted heating. The effectiveness of these modes is also geometry dependent and can be computationally optimized for plasma generation. Complementing these computational efforts, we demonstrate a new closed-system reactor design built in a WR-650 waveguide and experimentally demonstrate the formation of atmospheric argon microwave plasma using < 30 W input power on DMR dimers. We observe a shifting resonance frequency as the DMRs heat in response to microwave excitation and develop a Python-based lock-in mechanism to effectively track the DMR resonance over time, leading to stable plasma operation. We use infrared thermal imaging to monitor the temperature of the DMR dimers and surrounding quartz chamber, demonstrating thermal temperatures < 60 degreesC. Finally, we utilize optical emission spectroscopy (OES) to probe the plasma properties as a function of the resonance mode.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Intrinsic and induced quantum quenches for enhancing qubit-based quantum noise spectroscopy

Quantum sensing protocols that exploit the dephasing of a probe qubit are powerful and ubiquitous methods for interrogating an unknown environment. They have a variety of applications, ranging from noise mitigation in quantum processors, to the study of correlated electron states. Here, we discuss a simple strategy for enhancing these methods, based on the fact that they often give rise to an inadvertent quench of the probed system: there is an effective sudden change in the environmental Hamiltonian at the start of the sensing protocol. These quenches are extremely sensitive to the initial environmental state, and lead to observable changes in the sensor qubit evolution. We show how these new features give access to environmental response properties. This enables methods for direct measurement of bath temperature, and for detecting non-thermal equilibrium states. We also discuss how to deliberately control and modulate this quench physics, which enables reconstruction of the bath spectral function. Extensions to non-Gaussian quantum baths are also discussed, as is the application of our ideas to a range of sensing platforms (e.g., nitrogen-vacancy (NV) centers in diamond, semiconductor quantum dots, and superconducting circuits).

59 BASIC BIOLOGICAL SCIENCES↗