Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermal diffusion”

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

Plasma instabilities dominate radioactive transients magnetic fields: the self-confinement of leptons in Type Ia and core-collapse supernovae, and kilonovae

The light curves of radioactive transients, such as supernovae and kilonovae, are powered by the decay of radioisotopes, which release high-energy leptons through $\beta ^+$ and $\beta ^-$ decays. These leptons deposit energy into the expanding ejecta. As the ejecta density decreases during expansion, the plasma becomes collisionless, with particle motion governed by electromagnetic forces. In such environments, strong or turbulent magnetic fields are thought to confine particles, though the origin of these fields and the confinement mechanism have remained unclear. Using fully kinetic particle-in-cell (PIC) simulations, we demonstrate that plasma instabilities can naturally confine high-energy leptons. These leptons generate magnetic fields through plasma streaming instabilities, even in the absence of pre-existing fields. The self-generated magnetic fields slow lepton diffusion, enabling confinement, and transferring energy to thermal electrons and ions. Our results naturally explain the positron trapping inferred from late-time observations of thermonuclear and core-collapse supernovae. Furthermore, they suggest potential implications for electron dynamics in the ejecta of kilonovae. We also estimate synchrotron radio luminosities from positrons for Type Ia supernovae and find that such emission could only be detectable with next-generation radio observatories from a Galactic or local-group supernova in an environment without any circumstellar material.

instabilities↗
DOE OSTI2025

Mechanistic nuclear fuel performance modeling of uranium nitride

Uranium mononitride (UN) is a nuclear fuel candidate for advanced reactor designs and an alternative being considered for light water reactors due to its higher thermal conductivity and uranium density than UO 2 . As with any nuclear fuel, swelling and fission gas release are important factors for safety, while also being some of the hardest phenomena to predict with a high degree of confidence. Getting a grasp on the gas swelling behavior and release is crucial to lower the barrier for UN utilization. An accelerated swelling rate at high temperatures observed experimentally, sometimes referred to as “breakaway swelling,” further complicates the prediction of fuel performance of UN. A mechanistic model has been developed using a multiscale approach to describe the intragranular and intergranular fission gas behavior. Lower-length-scale calculations have been employed to inform models of the gas and self-diffusion behavior, resolution rate, and bubble shape. Leveraging previous work on high burnup UO 2 , two populations of intragranular bubbles are considered; small bulk bubbles and larger bubbles located along dislocations. The dislocation bubbles were found to be crucial to the overall swelling behavior, and the breakaway swelling transition was associated with the transition in the gas atom diffusion mechanism from an irradiation-induced athermal diffusion regime at lower temperatures to an intrinsic thermal equilibrium regime at higher temperatures, accelerating the growth of the dislocation bubbles. Similarly, the threshold for fission gas release was associated with the grain boundary vacancy diffusivity surpassing the gas atom diffusivity at sufficiently high temperatures, allowing the over-pressurized grain boundary bubble to grow in size and interconnect. Using thermo-mechanical models with the fission gas model, two integral fuel pin assessment cases were simulated. Finally, this work demonstrates the ability of a multiscale approach to accelerate the understanding of advanced fuel forms when experimental data is limited.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Dynamics and structure of the B2→B19’ phase transformation in NiTi revealed through in situ 4D-STEM

The structural evolution of NiTi during the B2→B19’ martensitic phase transformation via thermal cycling is investigated using in situ four dimensional scanning transmission electron microscopy (4D-STEM). With 4D-STEM, we can directly visualize and quantify the nanoscale evolution of the martensitic structure on thermal cycling and also investigate the origin of diffuse scattering of NiTi in the pre-transitional state. Mapping of the martensite orientation and strain visualizes the progression of the transformation front and self-accommodation of the B19’ structure. Diffuse streaking and strain are measured in the pre-transitional austenite (B2) phase and demonstrate no localization or preferential directionality hinting that long-range homogeneous instability rather than nanoscale heterogeneities may be the origin of the pre-transitional anomalies in NiTi. Finally, it is revealed that NiTi does not reform the same martensite nanostructure on thermal cycling but does express similar features. This small variation is likely owing to transformation-induced dislocations.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Multi-Objective design of interlocking metasurfaces using conditional diffusion models

Unit cell design remains a major challenge for interlocking metasurfaces, a promising joining technology for dissimilar materials, due to the complex, competing, multivariate design space and the need for rapid adaptation to varying performance requirements. This study explores Conditional Diffusion Models as a design optimization tool for interlocking metasurfaces. Given the complex, competing, multivariate design space for interlocking metasurfaces, unit cell design remains a major challenge for this joining technology. We trained a conditional diffusion model on 25,000 finite element analysis-simulated interlocking metasurface unit cells to generate designs with tailored thermo-mechanical properties (tensile strength, shear strength, and thermal conductivity) based on specified performance criteria. The model demonstrated a success rate of approximately 72 % in producing designs that met specified property bounds. The conditional diffusion model generated both thermally resistive and conductive designs, revealing clear trends in design characteristics: taller, dendritic structures were advantageous for tensile loads, while shorter, robust designs excelled in shear applications. Our findings indicate that the model's performance is more influenced by the breadth of the design space than by the quantity of training data, highlighting the importance of expansive design domains for generating innovative solutions. This work establishes conditional diffusion models as a highly efficient and adaptable tool for rapid interlocking metasurface unit cell design, paving the way for advancements in multi-material joining technologies, as well as highlighting the justification to leverage conditional diffusion models as design tools across complex design domains.

Conditional diffusion models↗
DOE OSTI2020

Physical effects on compact high-velocity clouds in the circumgalactic medium

ABSTRACT We numerically investigate the evolution of compact high-velocity clouds (CHVCs) passing through a hot, tenuous gas representing the highly ionized circumgalactic medium (CGM) by applying the adaptive-mesh refinement code flash. The model clouds start from both hydrostatic and thermal equilibrium and are in pressure balance with the CGM. Here, we present 14 models, divided into two mass categories and two metallicities each and different velocities. We allow for self-gravity and thermal conduction or not. All models experience mass diffusion, radiative cooling, and external heating leading to dissociation and ionization. Our main findings are (1) self-gravity stabilizes clouds against Rayleigh–Taylor instability, which is disrupted within 10 sound-crossing times without; (2) clouds can develop Jeans-instable regions internally even though they are initially below Jeans mass; (3) all clouds lose mass by ram pressure and Kelvin–Helmholtz instability; (4) thermal conduction substantially lowers mass-loss rates, by this, extending the clouds’ lifetimes, particularly, more than doubling the lifetime of low-mass clouds; (5) thermal conduction leads to continuous, filamentary stripping, while the removed gas is heated up quickly and mixes efficiently with the ambient CGM; (6) without thermal conduction the removed gas consists of dense, cool, clumpy fragments; (7) thermal conduction might prevent CHVCs from forming stars; and (8) clouds decelerated by means of drag from the ambient CGM form head-tail shapes and collapse after they reach velocities characteristic for intermediate-velocity clouds. Conclusively, only sophisticated modelling of CHVCs as non-homogeneous and non-isothermal clouds with thermal conduction and self-gravity explains observed morphologies and naturally leads to the suppression of star formation.

Sander, Bastian↗
DOE OSTI2021

Evidence for a thermally driven charge-density-wave transition in 1T-TaS 2 thin-film devices: Prospects for GHz switching speed

We report on the room-temperature switching of 1T-TaS 2 thin-film charge-density-wave devices, using nanosecond-duration electrical pulsing to construct their time-resolved current–voltage characteristics. The switching action is based upon the nearly commensurate to incommensurate charge-density-wave phase transition in this material, which has a characteristic temperature of 350 K at thermal equilibrium. For sufficiently short pulses, with rise times in the nanosecond range, self-heating of the devices is suppressed, and their current–voltage characteristics are weakly nonlinear and free of hysteresis. This changes as the pulse duration is increased to ~200 ns, where the current develops pronounced hysteresis that evolves nonmonotonically with the pulse duration. By combining the results of our experiments with a numerical analysis of transient heat diffusion in these devices, we clearly reveal the thermal origins of their switching. In spite of this thermal character, our modeling suggests that suitable reduction of the size of these devices should allow their operation at GHz frequencies.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Radiation induced athermal diffusivity in uranium mononitride

Uranium mononitride (UN) is one of the ceramic nuclear fuel alternatives to oxide fuel considered for light water reactors and advanced reactor designs. Properties like self- and fission gas diffusivity need to be better understood, given that they influence key fuel performance phenomena such as fission gas swelling and release. In particular, the radiation induced athermal (D 3 ) diffusivity remains challenging to accurately predict and has only been sparsely characterized in UN, despite its importance as it likely governs diffusion at the low temperatures this high-thermal-conductivity fuel form may operate. Molecular Dynamics simulations are used to estimate the mean square displacement induced by a primary knock-on atom (PKA) with a given kinetic energy. These results are combined with the PKA energy distributions obtained from binary collision approximation calculations to obtain the displacement due to a particular fission fragment. Finally, this is combined with experimental fission fragment yields to determine the displacement due to an average fission event and, thus, express the athermal diffusivity as a function of the fission rate density. These results are in excellent agreement with available experimental data. In conclusion, a particular importance is given to the understanding and the quantification of the variability of these results.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Silicon nanocluster anion-argon cation recombination via hybrid continuum-molecular dynamics calculations

Prediction of the charge distribution on particles in an aerosol is critical not only in electrical mobility based characterization methods, but also in understanding the role that charging plays in particle growth in gas phase synthesis reactors. The latter is particularly important in non-thermal plasma synthesis reactors, wherein nanoclusters form and grow from vapor phase precursors in a high electron and high ion density environment. In plasmas, because free electrons are much less massive and much more mobile than positive gas ions, nanocluster charge distributions are biased negative. Furthermore, neutral or even positively charged nanoclusters may exist, depending on the rate of nanocluster-gas ion recombination, and nanocluster-nanocluster collisions may greatly contribute to nanocluster growth if not Coulombically suppressed.

42 ENGINEERING↗
DOE OSTI2024

Mesoscale modeling and semi-analytical approach for the microstructure-aware effective thermal conductivity of porous polygranular materials

Here we established a comprehensive modeling approach for investigating the microstructure-aware effective thermal conductivity ($κ_{eff}$) for porous microstructures containing solid particles and gaseous pores. Our approach combines the mesoscale computational modeling framework and the semi-analytical method, allowing for efficient prediction of $κ_{eff}$ for realistic porous microstructures, while considering complicated microstructural thermal conduction pathways effectively in the prediction. We used the diffuse-interface mesoscale computational model to generate extensive simulated $κ_{eff}$ data for realistic digital representations of microstructures with wide ranges of porosity ($f_p$), thermal conductivity of the gas phase ($κ_g$), and thermal conductivity of the solid phase ($κ_s$). From the simulated data, we identified two property variation regimes for $κ_{eff}$: (1) a slow $κ_{eff}$ increase for $κ_s ~ κ_g$; and (2) a faster $κ_{eff}$ increase for $κ_s \gg κ_g$. To capture the key features of the relationship between the microstructure and $κ_{eff}$, we derived a semi-analytical model by introducing structure and intensification factors. The two new factors incorporate the calibrated effective contribution of the solid volume with $κ_s$ and additional interfacial effects into the prediction of $κ_{eff}$, respectively, allowing for consideration of parallel, serial, and interfacial conduction mechanisms effectively. Using the selected simulation data, we quantified key model parameters within the semi-analytical model and verified that the parameterized model exhibits excellent agreement with simulated $κ_{eff}$ for the entire range of the parameter space.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Diffuse neutrino background from past core collapse supernovae

Core collapse supernovae are among the most powerful explosions in the Universe, which emit thermal neutrinos that carry away most of the gravitational binding energy released. These neutrinos produce a diffuse supernova neutrino background (DSNB), which is one of the largest energy budgets among all radiation backgrounds. Detecting the DSNB is an important goal of modern high-energy astrophysics and particle physics, which provides valuable insights into core collapse modeling, neutrino physics, and cosmic supernova rate history. In this review, the key ingredients of DSNB calculation and what can be learned from future detections, including black hole formation and non-standard neutrino interactions are discussed. Moreover, an overview of the latest updates in neutrino experiments, which could lead to the detection of the DSNB in the next decade, is provided. With the promise of this breakthrough discovery on the horizon, the study of DSNB has great potential to further our understanding of the Universe.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

FSISPH: An SPH formulation for impacts between dissimilar materials

Here, we present an SPH formulation with several new features designed to better model the fully-compressible interaction of dissimilar materials. We developed the new method to simulate the atmospheric entry and break-up of small celestial bodies in planetary atmospheres. The formulation uses a unity-based, density-energy discretization of the hydrodynamic conservation laws with linear-corrected kernel gradients. To account for variations in compressibility, we use an HLLC approximate Riemann solver to adjust the velocity gradient at material interfaces. To handle large transverse velocity discontinuities, we introduce a simple slip interface model that limits the artificial viscosity at material interfaces. Diffusion is optionally applied through the velocity gradient and this allows the density and specific thermal energy to evolve in a manner more consistent with the first law of thermodynamics in comparison to other more direct diffusion schemes. We also introduce a material-local second-order artificial conduction scheme used to smooth the specific thermal energy field. Material damage fits neatly under this framework by treating the damage front as a material interface. The method has been implemented as a solver, FSISPH, within the code, Spheral++, and is publicly available on github. We test our new solver on a number of classic shock, mixing, and multi-material problem. The components we outline can significantly improve accuracy of SPH for problems with sharp contact discontinuities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Anharmonic lattice dynamics and superionic transition in AgCrSe 2

Intrinsically low lattice thermal conductivity ( κ l a t ) in superionic conductors is of great interest for energy conversion applications in thermoelectrics. Yet, the complex atomic dynamics leading to superionicity and ultralow thermal conductivity remain poorly understood. Here, we report a comprehensive study of the lattice dynamics and superionic diffusion in A g C r S e 2 from energy- and momentum-resolved neutron and X-ray scattering techniques, combined with first-principles calculations. Our results settle unresolved questions about the lattice dynamics and thermal conduction mechanism in A g C r S e 2 . We find that the heat-carrying long-wavelength transverse acoustic (TA) phonons coexist with the ultrafast diffusion of Ag ions in the superionic phase, while the short-wavelength nondispersive TA phonons break down. Strong scattering of phonon quasiparticles by anharmonicity and Ag disorder are the origin of intrinsically low κ l a t . The breakdown of short-wavelength TA phonons is directly related to the Ag diffusion, with the vibrational spectral weight associated to Ag oscillations evolving into stochastic decaying fluctuations. Furthermore, the origin of fast ionic diffusion is shown to arise from extended flat basins in the energy landscape and collective hopping behavior facilitated by strong repulsion between Ag ions. These results provide fundamental insights into the complex atomic dynamics of superionic conductors.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Experimental and modeling study of the effect of confinement on the thermal decomposition of organic materials

Simultaneous Thermal Analysis is widely used to study thermal decomposition of a variety of organic materials, and the choice of experimental conditions is very important to obtain reliable mass loss and heat flow curves. Here in this work we investigate the effect of confinement on the thermal decomposition of organic materials. As an example, we explore high-density polyethylene (HDPE) and one of its pyrolysis products, namely eicosane (C 20 H 42 ), through experiments and modeling in open crucibles and crucibles closed with a pierced lid. A new model is developed for evaporation and sublimation from open and pinhole thermal-analysis pans. It considers the kinetic resistance of evaporation at the liquid-vapor interface, inhibition of evaporation by buildup of vapor inside the pan, thermal expansion of vapor inside the pan, and diffusion of vapor from the pan. The model is validated using simultaneous heat-flow and mass-loss measurements of n-eicosane evaporation for various pinhole sizes, down to 50 µm. The importance of product evaporation inhibition for measuring polymer decomposition was demonstrated using HDPE. This study sheds light on the effect of confinement on the mass loss rate of organic materials.

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

Stable and Catalytically Active Shape-Engineered Cerium Oxide Nanorods by Controlled Doping of Aluminum Cations

Shape-engineered nanocrystals (SENs) promise a better selectivity and a higher activity in catalytic reactions than the corresponding non-shape-engineered ones because of their larger specific surface areas and desirable crystal facets. However, often, it is challenging to apply SENs in practical catalytic applications at high reaction temperatures, where SENs deforms into more stable, less active nanoparticles. In this paper, we show that atomic layer deposition (ALD) of Al 2 O 3 at 200 °C can controllably dope Al cations into the shape-engineered CeO 2 nanorods (NRs) to not only increase their shape transition temperature from 400 °C to beyond 700 °C but also greatly increase their specific reversible oxygen storage capacity (srOSC). Furthermore, the substituted Al3+ ions impede the surface diffusion of Ce ions and therefore improve the thermal stability of CeO 2 NRs. These Al3+ dopants form -Al–O–Ce–O– clusters, which are new Ce species and can be reversibly reduced and oxidized at 500–700 °C. This low-temperature chemical doping method decouples the synthesis process of SENs from the doping process and maintains the shape of the SENs during the activation of dopants. This concept could be adopted to enable the applications of other SENs in challenging high-temperature environments.

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

Prediction of Bi 2 Te 3 -Sb 2 Te 3 Interfacial Conductance and Superlattice Thermal Conductivity Using Molecular Dynamics Simulations

Bismuth telluride (Bi 2 Te 3 ) and its alloys with antimony telluride (Sb 2 Te 3 ) have long been considered to be the best room-temperature bulk thermoelectric (TE) materials. In recent decades, proof-of-concept demonstrations on Bi 2 Te 3 -Sb 2 Te 3 nanostructures have shown high TE performance due to reduction in lattice thermal conductivities. Particularly, ultra-low thermal conductivities have been observed in Bi 2 Te 3 -Sb 2 Te 3 1D superlattices, leading to thermoelectric figures of merit (ZT) as high as 2.4. In contrast, very few computational studies have been performed to provide insight into the phonon transport across these nanostructures. In this work, we use non-equilibrium molecular dynamics simulations with previously developed force fields to simulate thermal transport across Bi 2 Te 3 -Sb 2 Te 3 interfaces and superlattices. We first calculate the thermal conductance associated with a Bi 2 Te 3 -Sb 2 Te 3 interface across a temperature range of 200–400 K. Furthermore, the values are also compared with thermal conductances calculated by a modified Landauer transport formalism using phonon transmission coefficients obtained from the diffuse mismatch model. Our results show that inelastic scattering processes contribute to an increase in interfacial thermal conductance at higher temperatures. Finally, we calculate the thermal conductivities of Bi 2 Te 3 -Sb 2 Te 3 superlattices with varying period lengths from 2 to 18 nm. A minimum thermal conductivity of 0.27 W/mK is observed at a period length of 4 nm, which is attributed to the competition between incoherent and coherent phonon transport regimes. In comparison with previous experimental measurements in the literature, our results show good agreement with respect to the range of thermal conductivity values and the period length corresponding to the minimum superlattice thermal conductivity.

42 ENGINEERING↗
DOE OSTI2025

Unveiling the Origin of Morphological Instability in Topologically Complex Electrocatalytic Nanostructures

Coarsening and degradation phenomena in metals have largely focused on thermally driven processes, such as bulk and surface diffusion. However, dramatic coarsening has been reported in high-surface-area, nanometer-sized Pt-based catalysts during potential cycling in an electrolyte at room temperature─a temperature too low for the process to be explained purely by surface mobility values measured in both vacuum and electrolytes (∼10 –22 and ∼10 –18 cm 2 /s, respectively). This morphological evolution must be due to a different mechanism for mass transport that is sensitive to electrochemical conditions (e.g., electrolyte composition, potential limits, and scan rate). However, there have been no notable studies of electrochemically induced coarsening in nanometer-sized electrocatalysts. Here, we unveil the origins of coarsening in an electrolyte through coupled in situ experiments and atomistic kinetic Monte Carlo (kMC) simulations. Our work demonstrates electrochemical coarsening is driven by two concurrent mechanisms that can be explained at the atomistic level: (i) dissolution/redeposition during the reduction of an oxidized species and (ii) rapid surface diffusion of undercoordinated atoms.

gold↗
DOE OSTI2022

Tuning the Co/Fe ratio in BaCo x Fe 0.8– x Zr 0.1 Y 0.1 O 3– δ , a promising triple ionic and electronic conducting oxide, to boost electrolysis and fuel cell performance

The triple conducting oxide BaCo 0.4 Fe 0.4 Zr 0.1 Y0.1O 3–δ (BCFZY4411), which accommodates simultaneous transport of protons, oxygen ions, and p-type electronic carriers, has been intensively investigated in recent years as a high-performance positive electrode material for fuel cell and electrolysis applications. The heavy Co and Fe-based transition metal doping in BCFZY4411 ensures adequate electrical conductivity while the multiple oxidation states of Co and Fe assist the electrocatalytic and redox ability. Despite the considerable role of Co and Fe transition metal doping in controlling electrochemical activity, however, the study of alternative BCFZY compositions with varying Co/Fe ratios has not yet been pursued. Here, we evaluate the electrochemical performance of a series of BaCo x Fe 0.8–x Zr 0.1 Y0.1O 3–δ compositions with varying Co/Fe ratio (x = 0.1, 0.2, 0.4, 0.6, 0.7) and use oxygen ion tracer diffusion and in situ high-temperature X-ray diffraction to investigate the effect of Co/Fe ratio on electrocatalytic activity, electronic conductivity, oxygen ion incorporation and transport kinetics, and thermomechanical behavior. We find that Co-rich BCFZY7111 yields the highest performance due to exceptionally high oxygen vacancy diffusion and shows a lower and more linear thermal expansion behavior compared to Fe-rich compositions. A protonic ceramic button cell incorporating a BCFZY7111 positive electrode yields a peak power density of 695 mW cm –2 under fuel cell mode and an electrolysis current density of 1976 mA cm –2 at 1.4 V at 600 °C, underscoring the promise of this new BCFZY electrode composition.

30 DIRECT ENERGY CONVERSION↗