Engineering Papers⌕ Search

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.

At least 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

High-throughput screening of hypothetical metal-organic frameworks for thermal conductivity

Thermal energy management in metal-organic frameworks (MOFs) is an important, yet often neglected, challenge for many adsorption-based applications such as gas storage and separations. Despite its importance, there is insufficient understanding of the structure-property relationships governing thermal transport in MOFs. To provide a data-driven perspective into these relationships, here we perform large-scale computational screening of thermal conductivity k in MOFs, leveraging classical molecular dynamics simulations and 10,194 hypothetical MOFs created using the ToBaCCo 3.0 code. We found that high thermal conductivity in MOFs is favored by high densities (> 1.0 g cm -3 ), small pores (< 10 Å), and four-connected metal nodes. We also found that 36 MOFs exhibit ultra-low thermal conductivity (< 0.02 W m -1 K -1 ), which is primarily due to having extremely large pores (~65 Å). Furthermore, we discovered six hypothetical MOFs with very high thermal conductivity (>10 Wm -1 K -1 ), the structures of which we describe in additional detail.

36 MATERIALS SCIENCE↗
DOE OSTI2026

From dynamics to kinetics for the dissociation of •QOOH radicals derived from isopentane

Thermal dissociation rates for hydroperoxyalkyl (•QOOH) radicals are challenging to measure due to their inherent instability. Recent dynamics experiments have demonstrated that in situ synthesis of •QOOH radicals in a supersonic expansion in combination with IR action spectroscopy provides an effective route to accurately map their energy-resolved dissociation rates. Collaborative theoretical analyses provide a route for converting from those microcanonical rates to accurate temperature- and pressure-dependent thermal rate constants. Here, we present the results of this conversion for the dissociation of two isopentyl β-QOOH isomers of relevance to isopentane oxidation. The results are compared with existing literature values, as well as with results obtained in a similar fashion for related •QOOH systems. Extrapolations from benchmark calculations of barrier properties for smaller systems to accurate predictions for larger systems and multidimensional hindered rotor treatments are both shown to be important components of the theoretical analysis. For completeness, we extend the theoretical analysis to a treatment of the full RO 2 /QOOH system for isopentyl (1,1-dimethyl-propyl) radical.

Ab initio kinetics↗
DOE OSTI2025

Stress evolution and creep deformation in solid-oxide electrolysis cell systems – Dynamic modeling and multi-objective optimization to maximize stack life and efficiency

Here, this study develops a thermal stress model of solid-oxide electrolysis cells (SOECs) including a model for creep strain and failure probability that is integrated with a dynamic plant-wide model of a hydrogen production process. Uncertainties in key material properties of the cell are quantified to assess their impact on stress profile variability. The oxygen electrode is found to have about 10 times higher failure probability compared to the fuel electrode. The study shows that if the stack operation is not optimized, cycling operation would lead to stress build-up eventually leading to catastrophic failure. A dynamic optimization problem is set up for obtaining the optimal operational profile considering a variable hydrogen production rate. Due to the tradeoff between the efficiency and stress build-up, the dynamic optimization problem is multi-objective. It is observed that the optimizer can considerably reduce the stress build-up (i.e., can increase the stack life) albeit at the cost of a lower efficiency thus exhibiting strong tradeoffs between capital and operating costs. For example, if the stack would be replaced in 0.5 yr, specific energy requirement would be 48.5 kWh/kg H 2 while for a stack replacement time of about 6 yr, the specific energy requirement rises by about 4.2 %.

SOEC↗
DOE OSTI2024

Inelastic Deformation in Methylammonium Lead Iodide Perovskite and Mitigation by Additives during Thermal Cycling

Metal halide perovskite thin films are promising materials for next-generation photovoltaic applications, but the thermomechanical instabilities of these materials are critical barriers to device longevity. In this study, we measure the evolution of stresses in methylammonium lead iodide (MAPbI3) films thermally cycled from 25 to 105 °C using a multibeam optical stress sensor (MOSS) system. We demonstrate that these films undergo residual stress buildup due to inelastic deformation as well as a simple mitigation strategy by incorporating low concentrations of 5-aminovaleric acid hydrochloride (5-AVACl) as an additive. Adding 5-AVACl increases the tensile stresses measured at room temperature but beneficially decreases the biaxial modulus from 9.60 ± 0.44 GPa to 7.95 ± 0.55 GPa, resulting in greater accommodation of thermal stresses. The hysteresis loop is ~63% smaller with the inclusion of 5-AVACl, and the stress relaxation at 50 °C decreases from 86.46 to 51.49 MPa. The larger stress relaxation in pristine MAPbI 3 is correlated with grain boundary opening after five cycles, while the inclusion of 5-AVACl mitigates this degradation. Furthermore, our findings underscore the importance of studying the dynamic response of perovskite films during thermal cycling and lay the foundation for further exploration into the mechanisms governing the thermomechanical behavior of perovskite thin films.

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

Persistent dynamic magnetic state in artificial honeycomb spin ice

Abstract Topological magnetic charges, arising due to the non-vanishing magnetic flux on spin ice vertices, serve as the origin of magnetic monopoles that traverse the underlying lattice effortlessly. Unlike spin ice materials of atomic origin, the dynamic state in artificial honeycomb spin ice is conventionally described in terms of finite size domain wall kinetics that require magnetic field or current application. Contrary to this common understanding, here we show that a thermally tunable artificial permalloy honeycomb lattice exhibits a perpetual dynamic state due to self-propelled magnetic charge defect relaxation in the absence of any external tuning agent. Quantitative investigation of magnetic charge defect dynamics using neutron spin echo spectroscopy reveals sub-ns relaxation times that are comparable to the relaxation of monopoles in bulk spin ices. Most importantly, the kinetic process remains unabated at low temperature where thermal fluctuation is negligible. This suggests that dynamic phenomena in honeycomb spin ice are mediated by quasi-particle type entities, also confirmed by dynamic Monte-Carlo simulations that replicate the kinetic behavior. Our research unveils a macroscopic magnetic particle that shares many known traits of quantum particles, namely magnetic monopole and magnon.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Material Resilience in Harsh Service Conditions

Resilience describes the attributes of a material that allow it to withstand or resist detrimental environmental effects degrading properties and performance. In service, materials may experience harsh or extreme conditions, but even modest thermal or load conditions experienced over a long period can degrade performance. Thus, the National Nuclear Security Administration mission requires predictive understanding of materials performance in harsh and extreme conditions over long periods. This performance is particularly relevant for applications in which replacement is impractical, impossible, or costly. This area of leadership addresses the evolution of material properties in environments that include static and dynamic stress, radiation, and chemical or thermal extremes. A particular focus is on situations when environments coexist or for which collection of experimental data is challenging or impossible. The capability to predict and control the nature and evolution of properties to allow designing resilience is a crucial aspect of mission success in national nuclear, global, and energy security.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Dynamical ejecta synchrotron emission as a possible contributor to the changing behaviour of GRB170817A afterglow

ABSTRACT Over the past 3 yr, the fading non-thermal emission from the GW170817 remained generally consistent with the afterglow powered by synchrotron radiation produced by the interaction of the structured jet with the ambient medium. Recent observations by Hajela et al. indicate the change in temporal and spectral behaviour in the X-ray band. We show that the new observations are compatible with the emergence of a new component due to non-thermal emission from the fast tail of the dynamical ejecta of ab-initio binary neutron star merger simulations. This provides a new avenue to constrain binary parameters. Specifically, we find that equal mass models with soft equations of state (EOSs) and high-mass ratio models with stiff EOSs are disfavoured as they typically predict afterglows that peak too early to explain the recent observations. Moderate stiffness and mass ratio models, instead, tend to be in good overall agreement with the data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

In situ melt pool measurements for laser powder bed fusion using multi sensing and correlation analysis

Laser powder bed fusion is a promising technology for local deposition and microstructure control, but it suffers from defects such as delamination and porosity due to the lack of understanding of melt pool dynamics. To study the fundamental behavior of the melt pool, both geometric and thermal sensing with high spatial and temporal resolutions are necessary. This work applies and integrates three advanced sensing technologies: synchrotron X-ray imaging, high-speed IR camera, and high-spatial-resolution IR camera to characterize the evolution of the melt pool shape, keyhole, vapor plume, and thermal evolution in Ti–6Al–4V and 410 stainless steel spot melt cases. Aside from presenting the sensing capability, this paper develops an effective algorithm for high-speed X-ray imaging data to identify melt pool geometries accurately. Preprocessing methods are also implemented for the IR data to estimate the emissivity value and extrapolate the saturated pixels. Quantifications on boundary velocities, melt pool dimensions, thermal gradients, and cooling rates are performed, enabling future comprehensive melt pool dynamics and microstructure analysis. The study discovers a strong correlation between the thermal and X-ray data, demonstrating the feasibility of using relatively cheap IR cameras to predict features that currently can only be captured using costly synchrotron X-ray imaging. Such correlation can be used for future thermal-based melt pool control and model validation.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Occupational disorder as the origin of flattening of the acoustic phonon branches in the clathrate Ba 8 ⁢Ga 16 ⁢Ge 30

In the search for high-performance thermoelectrics, materials such as clathrates have drawn attention due to having both glasslike low phonon thermal conductivity and crystal-like high electrical conductivity. Ba 8 ⁢Ga 16 ⁢Ge 30 (BGG) has a loosely bound guest Ba atom trapped inside rigid Ga-Ge cage structures. Avoided crossings between acoustic phonons and the flat guest atom branches have been proposed to be the source of the low lattice thermal conductivity of BGG. Ga-Ge site disorder with Ga and Ge exchanging places in different unit cells has also been reported. We used time-of-flight neutron scattering to measure the complete phonon spectrum in a large single crystal of BGG and compared these results with predictions of density functional theory to elucidate the effect of the disorder on heat-carrying phonons. Experimental results agreed much better with the calculation assuming the disorder than with the calculation assuming the ordered configuration. Although the atomic masses of Ga and Ge are nearly identical, we found that disorder strongly reduces phonon group velocities, which significantly reduces thermal conductivity. Finally, our work points to a path towards optimizing thermoelectrics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Proton Collective Quantum Tunneling Induces Anomalous Thermal Conductivity of Ice under Pressure

Proton tunneling is believed to be nonlocal in ice, but its range has been shown to be limited to only a few molecules. Here, we measured the thermal conductivity of ice under pressure up to 50 GPa and found it increases with pressure until 20 GPa but decreases at higher pressures. We attribute this nonmonotonic thermal conductivity to the collective tunneling of protons at high pressures, supported by large-scale quantum molecular dynamics simulations. The collective tunneling loops span several picoseconds in time and are as large as nanometers in space, which match the phonon periods and wavelengths, leading to strong phonon scattering at high pressures. Finally, our results show direct evidence of global quantum motion existing in high-pressure ice and provide a new perspective to understanding the coupling between phonon propagation and atomic tunneling.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Impact of Thermal Stress on Attachment and Stability of High Temperature Strain Sensors

Static and dynamic strain sensing is required in high-temperature (HT), harsh-environments (HE) for industrial, aerospace, and energy sector applications to ensure equipment and process safety, reduce the cost of operation and maintenance, and increase process efficiency. Challenges that arise in HT HE sensing applications include device mounting, packaging, integrity, and stability to HT HE conditions. In previously reported work, static and dynamic strain surface acoustic wave resonator (SAWR) sensors were fabricated on langasite (LGS) and mounted on Inconel 625 strain beams for wireless testing up to 400°C. In this work, it has been identified that after subjecting the mounted SAWR strain sensor to thermal cycling between 100°C and 425°C, the measured sensitivity to dynamic strain decreased by 73% due to cracking at the adhesive/LGS interface, further deteriorating after additional thermal cycles. Strain modeling of the mounted LGS sensor chip up to 400°C revealed the existence of concentrated strain at the borders of the LGS chip. Microcracks caused by dicing make the chip boarders the most susceptible location for cracks to initiate when the sensor is subjected to thermal stress. In an attempt to mitigate the high strain at the LGS chip borders due to heating, adhesive shaping is proposed in this work. Simulations indicate a strain reduction of 50% at the border is achieved using both circular and triangular adhesive shapes, while also reducing the maximum strain over the entire adhesive/LGS interface by around 30%. Furthermore, the technique is thus promising for improving the integrity, reliability, and stability of static and dynamic strain sensors, particularly while operating under HT HE with several hundreddegree Celsius temperature excursions.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Computational fluid dynamic analysis of a novel particle-to-air fluidized-bed heat exchanger for particle-based thermal energy storage applications

Long-duration energy storage technologies are being targeted to enable cost-effective, decarbonized energy systems. Particle-based thermal energy storage systems are one promising technology by storing excess electricity or heat as sensible thermal energy in inexpensive, solid, inert particles. These systems are only possible if an effective and economical particle-to-working fluid heat exchanger exists. This study predicts the performance of a proposed, direct-contact, particle-to-air, pressurized fluidized-bed heat exchanger using computational fluid dynamics. The common Eulerian-Eulerian framework for modeling fluidized beds is first benchmarked to experimental results at a previously untested operating condition and application. Then, the benchmarked model evaluates the performance of a proposed design for a commercial-scale version of the novel particle-to-air heat exchanger. The results show pressure drop and gas-phase approach temperatures are advantageous compared to other proposed designs for particle-to-air heat exchangers in the literature; approach temperatures were less than 5 °C and gas-phase pressure drop across the fluidized bed was 32 kPa. The model also highlights the importance of gas distributor design and representation in the Eulerian-Eulerian framework to control fluidization behavior. In conclusion, the model built and benchmarked in this study can be leveraged to advance the design and analysis of these heat exchangers critical to the deployment of a promising long-duration energy storage technology.

25 ENERGY STORAGE↗
DOE OSTI2021

Flattening the Radial Temperature Profile across the Transformational Challenge Reactor Core

The Transformational Challenge Reactor (TCR) program is demonstrating an agile development approach to advanced nuclear reactor design, which has traditionally utilized a linear design process. In leveraging artificial intelligence, additive manufacturing, advanced materials, and cutting-edge modeling and simulation, the TCR program aims to minimize the high cost and lengthy deployment timelines now standard in the nuclear industry. Within a relatively short period of time, a robust and mature advanced gas-cooled reactor was iteratively designed under the TCR program, using these cutting-edge technologies. The TCR is a 3 MWt gas-cooled microreactor fueled with uranium nitride (UN) tristructural isotropic (TRISO) fuel particles. Though manufactured via traditional means, these UN TRISO particles are loaded into additively manufactured silicon carbide (SiC) cans [4]. Once loaded with TRISO particles, the SiC cans are densified using a chemical vapor infiltration process. The additively manufactured SiC enables significantly more freedom in the design of the fuel form than could ever be achieved using traditionally manufactured SiC. The helium coolant, pressurized to 5 MPa, enters the core at 300°C and nominally exits it at 500°C. Typically, the most thermally limiting components in any reactor design are the fuel assemblies in the core center. To provide a wide thermal margin in these central fuel assemblies, the flow may be biased toward the center of the core to more effectively cool these fuel assemblies with more power deposition and flatten the core’s radial temperature distribution. An analytical fluid model of the TCR core was developed to explore methods for biasing the flow away from the cooler outer fuel assemblies and towards the hotter inner ones. Higher-fidelity models developed in STAR-CCM+ 2020.3.1, a computational fluid dynamics code, were then utilized to verify the analytical model’s findings.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Microtearing instabilities and electron thermal transport in low and high collisionality NSTX discharges

Microtearing mode (MTM) real frequency, growth rate, magnetic fluctuation amplitude and resulting electron thermal transport are studied in systematic NSTX scans of relevant plasma parameters. Here, the dependency of the MTM real frequency and growth rate on plasma parameters, suitable for low and high collision NSTX discharges, is obtained by using the reduced MTM transport model. The plasma parameter dependencies are compared and found to be consistent with the results obtained from MTM using the Gyrokinetic GYRO code. The scaling trend of collision frequency and plasma beta is found to be consistent with the global energy confinement trend observed in the NSTX experiment. The strength of the magnetic fluctuation is found to be consistent with the gyrokinetic estimate. In earlier studies, it was found that the version of the Multi-Mode (MM) anomalous transport model, which did not contain the effect of MTMs, provided an appropriate description of the electron temperature profiles in standard tokamak discharges and not in spherical tokamaks. When the MM model, which involves transport associated with MTMs, is incorporated in the TRANSP code and is used in the study of electron thermal transport in NSTX discharges, it is observed that the agreement with the experimental electron temperature profile is substantially improved.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Giant isotope effect on phonon dispersion and thermal conductivity in methylammonium lead iodide

Lead halide perovskites are strong candidates for high-performance low-cost photovoltaics, light emission, and detection applications. A hot-phonon bottleneck effect significantly extends the cooling time of hot charge carriers, which thermalize through carrier–optic phonon scattering, followed by optic phonon decay to acoustic phonons and finally thermal conduction. To understand these processes, we adjust the lattice dynamics independently of electronics by changing isotopes. We show that doubling the mass of hydrogen in methylammonium lead iodide by replacing protons with deuterons causes a large 20 to 50% softening of the longitudinal acoustic phonons near zone boundaries, reduces thermal conductivity by ~50%, and slows carrier relaxation kinetics. Phonon softening is attributed to anticrossing with the slowed libration modes of the deuterated molecules and the reduced thermal conductivity to lowered phonon velocities. Our results reveal how tuning the organic molecule dynamics enables control of phonons important to thermal conductivity and the hot-phonon bottleneck.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Machine learning maximized Anderson localization of phonons in aperiodic superlattices

Nanostructuring materials to achieve ultra-low lattice thermal conductivity has proven to be extremely attractive for numerous applications such as thermoelectric energy conversion. Anderson localization of phonons due to aperiodicity can reduce thermal conductivity in superlattices, but the lower limit of thermal conductivity remains elusive due to the prohibitively large design space. In this work, we demonstrate that an intuition-based manual search for aperiodic superlattice structures (random multilayers or RMLs) with the lowest thermal conductivity yields only a local minimum, while a genetic algorithm (GA) based approach can efficiently identify the globally minimum thermal conductivity by only exploring a small fraction of the design space. Our results show that this minimum value occurs at an average RML period that is, surprisingly, smaller than the period corresponding to the minimum SL thermal conductivity. Interestingly, above this critical period, scattering of incoherent phonons at interfaces is less, whereas below this period, the room for randomization becomes less, thus putting more coherent phonons out of Anderson localization and causing increased thermal conductivity. Moreover, the lower limit of the thermal conductivity occurs at a moderate rather than maximum randomness of the layer thickness. Our machine learning approach demonstrates a general process of exploring an otherwise prohibitively large design space to gain non-intuitive physical insights.

42 ENGINEERING↗
DOE OSTI2023

Real-time spin systems from lattice field theory

We construct a lattice field theory method for computing the real-time dynamics of spin systems in a thermal bath. This is done by building on previous work of Takano with Schwinger-Keldysh and functional differentiation techniques. We derive a Schwinger-Keldysh path integral for generic spin Hamiltonians, then demonstrate the method on a simple system. Our path integral has a sign problem, which generally requires exponential run time in the system size, but requires only linear storage. The latter may place this method at an advantage over exact diagonalization, which is exponential in both. Our path integral is amenable to contour deformations, a technique for reducing sign problems.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Large-Eddy Simulation Study of Flow and Combustion Dynamics in a Full-Scale Hydrogen–Air Rotating Detonation Combustor-Stator Integrated System

In the present work, a first-of-its-kind three-dimensional (3D) large-eddy simulation (LES) study is conducted to numerically investigate the combustion dynamics as well as aero-thermal phenomena in a full-scale nonpremixed hydrogen–air rotating detonation engine (RDE) (with a diverging-shaped lower-end wall), when integrated with nozzle guide vanes (NGV) acting as the turbine stator. The wall-modeled LES framework incorporates hydrogen–air detailed chemical kinetics and adaptive mesh refinement (AMR). A comparative analysis is carried out for two operating conditions with different fuel/air mass flow rates but global equivalence ratio of unity, and considering RDE configurations without and with stator. The LES model is validated against available experimental data for the low mass flux condition with respect to detonation wave speed/height, wave dynamics, and axial static pressure distribution. Numerical results indicate significant deflagrative combustion occurring in the fill region near the inner wall due to formation of recirculation zones in the injection near-field driven by the backward facing step. The leading detonation wave is found to be trailed by an azimuthal reflected-shock combustion (ARSC) wave, consistent with experimental observations, which consumes unburned vitiated reactants that leak through the main detonation wave. The main detonation wave characteristics, such as detonation wave speed/height and combustion efficiency, do not change appreciably with the presence of NGV. A novel combustion diagnostic technique based on chemical explosive mode analysis (CEMA) is employed to quantify the fraction of heat release occurring in the detonative mode versus deflagrative mode for the simulated conditions. The exit flow is found to be nearly fully subsonic and supersonic for the low and high mass flux conditions, respectively. Further analysis of the exit flow profiles shows that the presence of NGV renders the flow more axial and significantly impacts the exit Mach number and total pressure, while the total temperature shows negligible change. In addition, the low mass flux operating point, despite exhibiting more deflagrative losses within the combustor, yields overall lower pressure drop from plenum to exhaust, which is mainly attributed to lower pressure drop across the injectors. Lastly, the rotating detonation engine-nozzle guide vanes (RDE-NGV) configuration exhibits higher total pressure loss compared to rotating detonation engine (RDE) without stator across both the mass flux conditions. In conclusion, this study extends the state-of-the-art in numerical modeling of pressure gain combustion (PGC) systems by demonstrating high-fidelity 3D reacting LES of full-scale RDE-NGV systems relevant to RDE-turbine integration for stationary power generation.

Combustion↗