SEARCH · Engineering Papers
Results for “thermal dynamics”
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.
Advancing simulations of coupled electron and phonon nonequilibrium dynamics using adaptive and multirate time integration
Electronic structure calculations in the time domain provide a deeper understanding of nonequilibrium dynamics in materials. The real-time Boltzmann equation (rt-BTE), used in conjunction with accurate interactions computed from first principles, has enabled reliable predictions of coupled electron and lattice dynamics. However, the timescales and system sizes accessible with this approach are still limited, with two main challenges being the different timescales of electron and phonon interactions and the cost of computing collision integrals. As a result, only a few examples of these calculations exist, mainly for two-dimensional (2D) materials. Here we leverage adaptive and multirate time integration methods to achieve a major step forward in solving the coupled rt-BTEs for electrons and phonons. Relative to conventional (non-adaptive) time-stepping, our approach achieves a 10x speedup for a target accuracy, or greater accuracy by 3–6 orders of magnitude for the same computational cost, enabling efficient calculations in both 2D and bulk materials. This efficiency is showcased by computing the coupled electron and lattice dynamics in graphene up to ~100 ps, as well as modeling ultrafast lattice dynamics and thermal diffuse scattering maps in bulk materials (silicon and gallium arsenide). In addition to improved efficiency, our adaptive method can resolve the characteristic rates of different physical processes, thus naturally bridging different timescales. This enables simulations of longer timescales and provides a framework for modeling multiscale dynamics of coupled degrees of freedom in matter. Our work opens new opportunities for quantitative studies of nonequilibrium physics in materials, including driven lattice dynamics with phonons coupled to electrons, spin, and other degrees of freedom.
Predictions of Cell-to-Cell Propagation and Vent Gas Production in the Thermal Runaway of Lithium-Ion Battery Stacks
This work presents the thermal runaway propagation model LIM1TR (Lithium-ion Modeling with 1-D Thermal Runaway) as an efficient tool to predict different cell-to-cell thermal runaway propagation scenarios. Here, we explored the vent gas volume production and reaction duration highlighting the relationship between these parameters and thermal runaway propagation due to convection by the vented gases. Two metrics based on gas production rate and heating rate are utilized as good indicators of the start and end of thermal runaway. LIM1TR results are compared with and validated by experiments from the literature for single-cell and multicell array experiments of 5 Ah and 10 Ah cells. By accounting for intraparticle diffusion of reacting species in the electrodes, we were able to capture the general dynamics of thermal runaway propagation and estimate acceptable reaction durations compared with the experimental values. Simulation results further demonstrated that varying heating modes lead to distinct reaction durations, consistent with experimental observations. Vent gas volume predictions indicate the need to consider both full and partial oxidation of the electrolyte. The outcomes of this work are building blocks for further investigations of module-to-module propagation by vented gases through convective heat transfer.
Numerical simulation of frost formation and heat transfer on fin-and-tube heat exchangers in turbulent cross-flow
Frost formation in fin-and-tube heat exchangers in turbulent cross-flow presents significant challenges in industrial refrigeration applications, affecting heat transfer efficiency and operational reliability. The purpose of this work is to investigate frost deposition and growth on a staggered bank of a fin-and-tube freezer coil under turbulent forced convection conditions. The focus here is on investigating conditions that closely replicate real-world scenarios in large walk-in industrial freezers. Using a direct numerical simulation approach, we examine the flow dynamics and thermal behaviour in the presence of frost, considering turbulent regimes characterized by a Reynolds number in the range 1050 ≤ R e D , avg ≤ 4800 , with the characteristic length being the outer diameter of the tube and the velocity being the bulk fluid velocity between the plates (fins). Computational fluid dynamics simulations are employed to resolve the interactions between turbulent airflow and the frost layer. Our approach incorporates a modified immersed boundary method and a slow-time acceleration technique to address the complex dynamic interface between the continuously evolving frost layer and the flowing air stream. Our findings indicate that frost forms more on the sides of the finned surfaces (plates) and less on the tubes themselves. This article is part of the theme issue ‘Heat and mass transfer in frost and ice’.
Proliferated Resilient Economical Half-Meter Aperture Space Telescopes (PREEMPT)
The PREEMPT LDRD was motivated by a major national need for better space-based imaging systems that are both high performing and affordable for the US Government. Current optical payloads for intelligence, surveillance and reconnaissance (ISR) and space domain awareness (SDA) missions can cost hundreds of millions of dollars and take years to develop, which makes it difficult to build the large constellations required for persistent, world-wide coverage. To address this, the team aimed to advance a different kind of large-aperture (>25 cm) telescope, called a monolithic telescope, in which key optical surfaces are built into a single piece of fused silica. This design greatly reduces payload and spacecraft complexity, the need for precision focus actuators, improves mechanical and thermal robustness, and lowers cost when compared with traditional Cassegrain telescopes that rely on many precisely aligned components. The project focused on five primary thrusts. The first thrust was to advance the concept of a V10, 25 cm, monolithic telescope forward from optical design to flight-ready stage. This was accomplished in partnership with Optimax, who delivered the first test unit in the early stages of the LDRD. The team developed several technologies necessary for this optic to be integrated into a flight demonstration. These include carbon fiber housings, highly detailed structural and thermal models and stress-reducing elastic averaging Hirth groove designs. These technologies resulted in a successful maturation of the optic, which is now slated to fly in late 2026/early 2027 for a demonstration mission. The second thrust was to advance the manufacturability of these optics. In collaboration with NIF’s optical manufacturing shop, we reduced polishing time from 480 hours to 65 hours through the implementation of optimized processes and new tools. The NIF team utilized a conceptual V8 (18 cm) optic to demonstrate this optimization, though it can be applied to the rest of the monolithic optic portfolio. Third, the team developed the first conceptual V20 (50 cm) payload, which is slated to be the next generation of LLNL optical payload systems. A set of structural, dynamic and thermal simulations were performed to identify potential challenges in the future development of this payload. Early-stage simulations suggest the payload is feasible, though thermal management will be the key focus area to maintain optimal performance. Fourth, a non-linear model of Viton was developed, to further enhance the reliability of our structural and dynamic models for future payloads. Viton acts as the primary interface material between the optic and its housing. Lastly, the team focused on successfully displaying the feasibility of using additively manufactured metal composites for optical space payloads. The team successfully demonstrated layer by layer deposition of Al-SiC composites, which have highly tunable structural and coefficient of thermal expansion (CTE) properties. These are crucial for optical payloads because CTE mismatch is one of the causes for degraded optical performance for telescopes in orbit. Overall, the work showed that monolithic telescopes could become a practical, lower-cost path to high-resolution space imaging for both national security and scientific missions.
Nonequilibrium Thermodynamics of Colloidal Gold Nanocrystals Monitored by Ultrafast Electron Diffraction and Optical Scattering Microscopy
Metal nanocrystals exhibit important optoelectronic and photocatalytic functionalities in response to light. These dynamic energy conversion processes have been commonly studied by transient optical probes to date, but an understanding of the atomistic response following photoexcitation has remained elusive. In this paper, we use femtosecond resolution electron diffraction to investigate transient lattice responses in optically excited colloidal gold nanocrystals, revealing the effects of nanocrystal size and surface ligands on the electron–phonon coupling and thermal relaxation dynamics. First, we uncover a strong size effect on the electron–phonon coupling, which arises from reduced dielectric screening at the nanocrystal surfaces and prevails independent of the optical excitation mechanism (i.e., inter- and intraband). Second, we find that surface ligands act as a tuning parameter for hot carrier cooling. Particularly, gold nanocrystals with thiol-based ligands show significantly slower carrier cooling as compared to amine-based ligands under intraband optical excitation due to electronic coupling at the nanocrystal/ligand interfaces. Finally, we spatiotemporally resolve thermal transport and heat dissipation in photoexcited nanocrystal films by combining electron diffraction with stroboscopic elastic scattering microscopy. Taken together, we resolve the distinct thermal relaxation time scales ranging from 1 ps to 100 ns associated with the multiple interfaces through which heat flows at the nanoscale. Our findings provide insights into optimization of gold nanocrystals and their thin films for photocatalysis and thermoelectric applications.
Atomic dynamics in 𝑀Cr𝑋 2 (𝑀=Ag, Cu; 𝑋 = S, Se) across magnetic and superionic transitions
Here, a systematic study of atomic dynamics and thermal properties of the family of layered chalcogenide compounds 𝑀Cr𝑋 2 (𝑀= Ag, Cu; 𝑋 = S, Se) was performed, including neutron and x-ray scattering, thermal characterization, and first-principles simulations. In all compounds, we observe a breakdown of specific phonon modes across the superionic phase transition, for phonons whose eigenvectors exhibit large contributions of mobile ions. In particular, the nondispersive portions of transverse acoustic (TA) branches at short-wavelengths and the low-energy optical phonons with large contributions from Ag + or Cu + become severely damped in the superionic phase. However, well-defined quasiparticles persist in the superionic state for long-wavelength TA phonons. In the case of AgCrS 2 , the coupling of lattice dynamics with its antiferromagnetic transition was also investigated. The magnetic ordering couples with the monoclinic–rhombohedral structural transition, and the Cr 3+ spin arrangement strongly affects the phonon dispersions. We qualitatively reproduce the magnetic and nuclear components of the INS measurement for antiferromagnetic AgCrS 2 by combining models of spin-waves and spin-polarized first-principles phonon simulations. Quasielastic magnetic fluctuations persist in the paramagnetic phase up to high temperature, but are clearly distinguished from the nuclear component through their momentum dependence. Finally, we report measurements of the thermal properties of the selenide compounds and find good agreement with our DFT simulations.
Toward high laser power beam manipulation with nanophotonic materials: evaluating thin film damage performance
Nanophotonic materials enable unprecedented control of light-matter interactions, including the ability to dynamically steer or shape wavefronts. Consequently, nanophotonic systems such as metasurfaces have been touted as promising candidates for free-space optical communications, directed energy and additive manufacturing, which currently rely on slow mechanical scanners or electro-optical components for beam steering and shaping. However, such applications necessitate the ability to support high laser irradiances (> kW/cm 2 ) and systematic studies on the high-power laser damage performance of nanophotonic materials and designs are sparse. Here, we experimentally investigate the pulsed laser-induced damage performance (at λ ∼ 1 µm) of model nanophotonic thin films including gold, indium tin oxide, and refractory materials such as titanium nitride and titanium oxynitride. We also model the spatio-thermal dissipation dynamics upon single-pulse illumination by anchoring experimental laser damage thresholds. Our findings show that gold exhibits the best laser damage resistance, but we argue that alternative materials such as transparent conducting oxides could be optimized to balance the tradeoff between damage resistance and optical tunability, which is critical for the design of thermally robust nanophotonic systems. We also discuss damage mitigation and ruggedization strategies for future device-scale studies and applications requiring high power beam manipulation.
Interface-independent sound speed and thermal conductivity of atomic-layer-deposition-grown amorphous AlN / Al 2 O 3 multilayers with varying oxygen composition
Dielectric amorphous multilayers (AMLs) play a critical role in a wide array of technologies such as optical coatings, nanoelectronics, energy harvesting, and recovery devices. However, despite their wide applications, a robust understanding of the effect of the interplay between chemical and structural disorder on the thermal properties of AMLs is still lacking. Therefore, in this paper, we experimentally and numerically investigate the effects of composition and interface density on the sound speed and thermal conductivity of a series of amorphous aluminum nitride and aluminum oxide multilayers grown via plasma-enhanced atomic layer deposition. To systematically change the composition, the oxygen content of the AMLs is proportionally varied with interface density during growth. We find that the longitudinal sound speed of these AMLs is dictated by the oxygen content instead of the number of interfaces. The thermal conductivity, in contrast, is dictated by both interface density and oxygen content. The interfaces act to decrease the thermal conductivity, whereas the oxygen content increases the thermal conductivity. Due to the competing influence of the interfaces and oxygen content, the thermal conductivity of the AMLs remains nearly constant as a function of interface density. Furthermore, our study provides crucial insights into the effect of the interplay of composition and interfaces on the sound speed and thermal conductivity of AMLs.
Scaling patch analysis of turbulent planar plume
Proper scaling in turbulent planar plumes is investigated in this study using a scaling patch approach. Based on the scaled boundary conditions, a proper velocity scale for the mean axial flow is the plume centerline velocity U ref = U ctr , and a proper temperature scale for the temperature excess is Θ ref = T ctr – T ∞ , where T ctr is the plume centerline temperature and T ∞ is the ambient fluid temperature. By seeking an admissible scaling, a key concept in the scaling patch approach, for the mean continuity, mean momentum, and mean energy equations, respectively, the following is found: (1) a proper scale for the mean transverse flow is V ref = ( d δ / d x ) U ctr , where d δ / d x is the growth rate of the plume width. (2) A proper scale for the Reynolds shear stress is R vu , ref = U ctr V ref = ( d δ / d x ) U ctr 2 , a mix of the scales for the mean axial and transverse flows. (3) A proper scale for the turbulent heat flux is R v θ , ref = V ref Θ ctr , a mix of the scales for the mean transverse flow and mean temperature excess. The mean transverse flow thus plays a critical role in the scaling of turbulent planar plumes. Approximate functions are developed for the scaled mean transverse flow, Reynolds shear stress, and turbulent temperature flux, and are found to agree favorably with experimental and numerical simulation data. The integral analysis of the mean momentum equation yields a Richardson number Ri, which remains invariant in the axial direction. The Richardson number is defined as Ri = def g β Θ ctr δ t / ( U ctr V ref ) ≈ 1 / 2 , where g is the gravitational acceleration, β is the thermal expansion coefficient, and δ t is the plume half-width based on the mean temperature profile. This Richardson number arises directly from the scaling patch analysis of the mean momentum equation, including both the streamwise and transverse velocity scales.
Quantitative analysis of cost savings and occupants’ preferences in grid-interactive smart home operation
Many utility companies in the United States have introduced time-of-use (TOU) rates for homeowners with the goal of regulating electricity consumption during peak hours. The electrical appliances in homes include various thermostatically controlled devices, such as air conditioners (AC) for thermal comfort, and nonthermostatically controlled devices such as clothes washers. As a result, homeowners face the complicated challenge of economically operating multiple electrical appliances in their homes while maintaining comfort and convenience. This is usually due to the lack of an explicit understanding of the correlation between cost saving and the users’ comfort. To understand the correlation, this article is designed to construct a framework by integrating three major components: a multi-objective optimization method accommodating multiple competing goals with different weights, a learning-based system modeling approach describing the dynamics and thermal coupling effects of appliances, and a novel comfort index method differentiating preferred and acceptable thermal comfort. Our proposed framework can allow the indoor air temperature to fall into the "preferred" range with a marginal cost increase. Furthermore, the simulation result shows that an additional 8 h for the preferred thermal comfort can be achieved with a cost increase of only 1.77%.
Chemical structure and curing dynamics of bisphenol S, PEEK™-like, and resveratrol phthalonitrile thermoset resins
This effort provides a multifaceted analysis of the structural changes and material dynamics of thermally driven softening and curing of three distinct phthalonitrile (PN) resins that cross-link into thermally stable and oxidation-resistant thermosets. Although this material system has yielded a large subset of fire-retardant composites that require facile processing and low-temperature curing, to date, insufficient information had been available on the fundamental processes that drive their softening and curing stages. Our approach conducted a complementary analysis the chemistry, monomer mobility, and rheology of three PN polymers in order to correlate the curing processes with corresponding structural and behavioral transformations of thermosets. Here, we focused on PNs with a bisphenol S backbone, a bisphenol A (PEEK™-like) backbone, and a resveratrol backbone. We relied on quasi-elastic neutron scattering (QENS) in order to analyze the in situ dynamics and self-diffusion properties of PN monomers, and to track changes in their mobilities during cross-linking and staging. Our analysis facilitates proper control over the staging and final curing of these resins and enables more efficient processing of these thermosets.
Elucidating the effect of preheating temperature on melt pool morphology variation in Inconel 718 laser powder bed fusion via simulation and experiment
In laser powder bed fusion (L-PBF) additive manufacturing, the mechanical performance, microstructure and defects of fabricated parts are closely associated with the melt pool morphology, e.g., its dimension and shape through the building process. Past studies have largely focused on how the process parameters such as laser power and scan speed affect melt pool characteristics. In this study, the melt pool morphology variation as a function of preheating temperature in the conduction, transition, and keyhole regimes and the underlying mechanisms in each regime are investigated through ex-situ sample characterization and computation thermal fluid dynamics (CtFD) simulation. Single tracks with different combinations of laser power and scan speed are deposited on an Inconel 718 bare plate preheated to a temperature range of 100–500 °C in the experiment. Significant changes are observed in melt pool morphology as a function of preheating temperature from optical measurements of melt track cross sections. The depth of melt pool in the three regimes increases monotonically with preheating temperature, e.g., at 500 °C, the experimental melt pool depth is increased by 49% in conduction regime, 34% in transition regime and 33% in keyhole regime, respectively, while the variation of melt pool width in each regime does not all follow an increasing trend but depends on the melt pool regimes. Melt pool width variation in the conduction and transition regimes is found to depend on the enhanced heat conduction directly related to temperature dependent thermal properties. Through validated CtFD simulations, it is found that in the keyhole regime the evaporation mass, recoil pressure, and laser drilling effect is enhanced with higher preheating temperature, which gives rise to a deeper melt pool. The simulations also reveal that preheating temperature significantly elongates the melt track length due to the increased flow rate and strong recoil pressure that accelerates the backward flow.
Combustion dynamics of crude and upgraded Thermal DeOxygenation oils in a compression ignition engine
Thermal DeOxygenation (TDO) is a robust thermochemical conversion scheme to produce hydrocarbons from biomass feedstock. The process targets the carbohydrate fraction of biomass to yield a broad mixture of primarily aromatic hydrocarbons within a boiling point range of 348–798 K with low oxygen content (<4 wt%). The resulting materials are amenable to traditional hydrotreating and distillation processes with approximately 70% by mass in the distillate fuel range. The simple conversion scheme combined with commercially viable upgrading routes makes TDO oils a feasible alternative fuel for transportation applications. Here this paper is the first systematic treatment of fit-for-purpose testing of TDO oils in a compression ignition engine. Blends of partially upgraded TDO oils (e.g. whole oil, distilled oil, hydrotreated oil and hydrotreated-distilled oil) are prepared with certified ultra-low sulfur diesel at 5%, 10%, 15% and 20% by volume. The resulting fuels are analyzed for fuel characteristics and combustion dynamics in an instrumented single-cylinder compression ignition engine. All fuel blends at 10% blend level by volume exhibited adequate engine performance. Hydrotreating of the oils is a necessary step to meet cetane number and EPA soot emissions requirements at 20% blend volume. Combined distillation and hydrotreatment meet or exceed all fuel specification and engine performance benchmarks. Partially upgraded TDO oils are suitable fuel options for compression ignition applications.
Landau-Lifshitz-Bloch equation for ferrimagnets with higher-order interaction
In this work, we present a micromagnetic formulation for modeling the magnetization dynamics and thermal equilibrium in ferrimagnetic materials at low and elevated temperatures. The formulation is based on a mean-field approximation (MFA). In this formulation, the ferrimagnet is described micromagnetically by two coupled sublattices with corresponding interactions, including inter- and intrasublattice micromagnetic exchange as well as four-spin interactions described as an intersublattice molecular field with a cubic dependence of the magnetization. The MFA is used to derive a Landau-Lifshitz-Bloch type equation for ferrimagnetic material, including cases with a ferromagnetic-antiferromagnetic phase transitions. For validation, the results obtained via the presented model are compared with recent experimental data for phase transitions in FeRh.
Mitigation of Molten Salt Corrosion
The desire to increase the peak operational temperature of Concentrated Solar Power (CSP) plants to lower the levelized cost of CSP-derived electricity, has led to consideration of an earth-abundant, low-cost, thermally-stable, molten MgCl 2 -KCl-NaCl (40-40-20 mole%) salt for high-temperature thermal energy storage and heat transfer. The reliable, long-term containment of such a molten salt is challenging, however, due to the corrosion of cost-effective, high-temperature structural metal alloys in this salt (particularly if this salt is contaminated with oxygen-bearing species). This project has been focused on: i) development of equipment for evaluating corrosion kinetics in molten chlorides under well-controlled thermal, fluid dynamic, and atmospheric conditions at CSP-relevant temperatures, ii) use of such equipment to evaluate the influence of forced convection on corrosion kinetics (to determine the role of liquid phase diffusion on such corrosion), and iii) identifying a corrosion-mitigation strategy that can be industrially applied to achieve a projected annual corrosion ≤30 μm at CSP-relevant temperatures. The overall goals of this project have been: i) to understand the rate-limiting step(s) and associated kinetic laws governing corrosion in molten chlorides at 750°C, ii) to determine key operational parameters affecting the corrosion kinetics in molten salts, and iii) to demonstrate the tailoring of such parameters to appreciably reduce such corrosion.
A unified understanding of minimum lattice thermal conductivity
Here, we propose a first-principles model of minimum lattice thermal conductivity ($κ^{min}_L$) based on a unified theoretical treatment of thermal transport in crystals and glasses. We apply this model to thousands of inorganic compounds and find a universal behavior of $κ^{min}_L$ in crystals in the high-temperature limit: The isotropically averaged $κ^{min}_L$ is independent of structural complexity and bounded within a range from ~0.1 to ~2.6 W/(m K), in striking contrast to the conventional phonon gas model which predicts no lower bound. We unveil the underlying physics by showing that for a given parent compound, $κ^{min}_L$ is bounded from below by a value that is approximately insensitive to disorder, but the relative importance of different heat transport channels (phonon gas versus diffuson) depends strongly on the degree of disorder. Moreover, we propose that the diffuson-dominated $κ^{min}_L$ in complex and disordered compounds might be effectively approximated by the phonon gas model for an ordered compound by averaging out disorder and applying phonon unfolding. With these insights, we further bridge the knowledge gap between our model and the well-known Cahill–Watson–Pohl (CWP) model, rationalizing the successes and limitations of the CWP model in the absence of heat transfer mediated by diffusons. Finally, we construct graph network and random forest machine learning models to extend our predictions to all compounds within the Inorganic Crystal Structure Database (ICSD), which were validated against thermoelectric materials possessing experimentally measured ultralow κ L . Our work offers a unified understanding of $κ^{min}_L$, which can guide the rational engineering of materials to achieve .
Investigation of Low-Energy Lattice Dynamics and Their Role in Superionic Na Diffusion and Ultralow Thermal Conductivity of Na 3 PSe 4 as a Solid-State Electrolyte
The atomic dynamics of Na 3 PSe 4 were investigated using a combination of neutron scattering experiments and ab initio and machine-learned molecular dynamics simulations to probe the interplay of fast ionic diffusion with atomic vibrations (phonons) of the host lattice. Our results reveal the existence of low-energy vibrational modes, simultaneously involving motions of Na + ions and framework polyanion subunits, and show that these modes become strongly overdamped in the superionic regime as they couple with the Na + hopping process. In particular, the Na+ migration energy landscape is strongly impacted by low-energy phonons derived from a soft acoustic branch of the host lattice, which modulates the diameter of the Na + diffusion channel at the bottleneck. We find that an additional factor for the enhanced Na + conductivity in Na 3 PSe 4 is the presence of Na-vacancies, which also affect the low-frequency dynamics and thermal vibration amplitudes, pointing to an interplay between Na + vacancies and host dynamics, jointly enhancing ionic diffusivity. Finally, we investigate the origin of ultralow thermal conductivities in Na 3 PSe 4 and Na 3 PS 4 using Green–Kubo simulations and find that low-energy acoustic phonon modes of the overall crystal framework provide a dominant contribution to the thermal conductivity.