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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Exceptional thermal stability of additively manufactured CoCrFeMnNi high-entropy alloy with cellular dislocation structures

CoCrFeMnNi high-entropy alloy (HEA) was additively manufactured (AM) by laser powder-bed fusion (L-PBF). The AM CoCrFeMnNi has prominent cellular dislocation structures with a small number of Mn-rich oxides. The thermal stability of the AM CoCrFeMnNi was investigated by isochronal annealing treatment at various temperatures from 400 to 1300°C for 1h. Microstructural analysis shows slow dislocation recovery, retarded recrystallization process, and precipitation of additional Cr-Mn based oxides during thermal annealing, resulting in exceptional thermal stability and retained high hardness at elevated temperatures. Further, by thermodynamic calculations, a low stored energy of 1.31 MJ/m 3 and a high activation energy of 353 kJ/mol for recrystallization were estimated for the AM CoCrFeMnNi. The exceptional thermal stability of the AM CoCrFeMnNi HEA is mechanistically attributed to the low crystallographic misorientations across the dislocation cell walls, sluggish atomic diffusion, and the pinning effects of the oxide nanoprecipitates.

36 MATERIALS SCIENCE↗
DOE OSTI2024

First-Principles Studies of Tritium Species Diffusivity Across the Interface of Nickel-Plated Zircaloy-4

α-Zr and its alloys are known as best 3 H getters due to their excellent corrosion resistance under chemically corrosive environment, low thermal neutron absorption cross-section and robust mechanical strength at high temperature. In this work, we conducted a systematic study on constructing and optimizing Ni and Zr surfaces, interfacing the stable Ni and Zr surfaces to create an optimal Ni-Zr interface, and an understanding the diffusion mechanisms for 3 H across Ni(111)-Zr(0001) interface with and without vacancies and impurities such as Sn and O.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Kinetically Controlled Synthesis of Rhodium Nanocrystals with Different Shapes and a Comparison Study of Their Thermal and Catalytic Properties

Here, we report the synthesis of Rh nanocrystals with different shapes by controlling the kinetics involved in the growth of preformed Rh cubic seeds. Specifically, Rh nanocrystals with cubic, cuboctahedral, and octahedral shapes can all be obtained from the same cubic seeds under suitable reduction kinetics for the precursor. The success of such a synthesis also relies on the use of a halide-free precursor to avoid oxidative etching, as well as the involvement of a sufficiently high temperature to remove Br – ions from the seeds while ensuring adequate surface diffusion. The availability of Rh nanocrystals with cubic and octahedral shapes allows for an evaluation of the facet dependences of their thermal and catalytic properties. The data from in situ electron microscopy studies indicate that the cubic and octahedral Rh nanocrystals can keep their original shapes up to 700 and 500 °C, respectively. When tested as catalysts for hydrazine decomposition, the octahedral nanocrystals exhibit almost 4-fold enhancement in terms of H 2 selectivity relative to the cubic counterpart. As for ethanol oxidation, the order is reversed, with the cubic nanocrystals being about three times more active than the octahedral sample.

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

Effects of Hydrogen Bonding on Nuclear Data Development of Liquid Anhydrous HF

Anhydrous Hydrogen Fluoride (HF) at high temperatures and pressures is used to process and manufacture nuclear fuel. As HF is often used directly with uranium, correct neutron thermal scattering cross sections are crucial to criticality safety applications. Classical molecular dynamics (CMD) simulation of the flexible HF system was used to create the thermal scattering law (TSL) and cross sections. The initial 2-site model is used in LAMMPS, and it can not capture the H-bond. To correctly represent the H-bond effects, a second, 3-site model was constructed in GROMACS. The 3-site model handled H-bonds by connecting a massless charge to the molecule. Key model parameters were compared to experimental data to verify the approach and models. To get the normalized VACF, the model was compared using hydrogen and fluorine bond length, density, potential energy, and diffusion coefficient. The phonon DOSs for both models were derived from the normalized VACF. DOSs were used to estimate the TSL ( S ( α, β )) and neutron thermal scattering cross sections for hydrogen in HF. The TSLs were evaluated using the FLASSH code with the Schofield diffusion model. It was observed that the representation of the hydrogen bonding changes the TSL's diffusional contributions. This is represented in the low energy scattering cross section, where intermolecular binding effects shift the cross section.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Dynamics of deuterium retention and desorption from plasma-facing materials in fusion reactor-relevant conditions

Hydrogen isotopes retention and desorption during and after discharges in fusion devices are still not well understood due to the complex device conditions and limitations of in-situ diagnostics and measurements. Here, we simulated well-diagnosed recent experiments at the DIII-D facility to benchmark our ITMC-DYN integrated package of modeling deuterium diffusion, retention, and desorption during and after D discharge irradiation. Modeling results were compared with detail experimental data of D desorption fluxes for various irradiation conditions. We predicted the temporal evolution of free and trapped D distribution in tungsten (W) plasma-facing material (PFM). Effects of key parameters namely diffusion coefficient, recombination rate, trapping energies against different defect types, were examined in these simulations. Existing experimental data of these parameters in literature varies significantly which makes it harder to identify key mechanisms and physics responsible for hydrogen isotope retention and desorption. The purpose of this work is to accurately simulate recent well-diagnosed reactor experiments given the uncertainties in such parameters and identify mechanisms responsible for the retention and desorption. We implemented the best identified diffusion, recombination, and trapping parameters in ITMC-DYN package that integrate both various collisional and thermal processes. We predicted, for example, that sample cooling between discharges in DIII-D operations can significantly affect the spatial distribution of trapped D in W under reactor irradiation conditions. Correct prediction of desorption spectra from samples irradiated during 10 DIII-D discharges showed that up to 35% of D can be retained in high binding energy defects such as vacancy clusters or voids.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Dark matter in X-rays: revised XMM-Newton limits and new constraints from eROSITA

We investigate two classes of dark matter (DM) candidates, sub-GeV particles and primordial black holes (PBHs), that can inject low-energy electrons and positrons into the Milky Way and leave observable signatures in the X-ray sky. In the case of sub-GeV DM, annihilation or decay into e + e - contributes to the diffuse sea of cosmic-ray (CR) leptons, which can generate bremsstrahlung and inverse Compton (IC) emission on Galactic photon fields, producing a broad spectrum from X-rays to γ-rays detectable by instruments such as eROSITA and XMM-Newton. For PBHs with masses below ∼10 17 g, Hawking evaporation similarly yields low-energy e ± , leading to comparable diffuse emission. Using the first data release from eROSITA and incorporating up-to-date CR propagation and diffusion parameters, we derive new constraints on both scenarios. For sub-GeV DM, we exclude thermally averaged annihilation cross sections in the range ∼ 10 -27 –10 -25 cm 3 /s and decay lifetimes of ∼ 10 24 –10 25 s for masses between 1 MeV and 1 GeV, with eROSITA outperforming previous X-ray constraints below ∼ 30 MeV. For asteroid-mass PBHs, we set new bounds on the DM fraction based on their Hawking-induced emission. Finally, we revisit previous constraints in 10.1088/1475-7516/2023/07/026JCAP 07 (2023) 026 and consequently 10.3847/1538-4357/ad41e0Astrophys. J. 968 (2024) 46 derived from XMM-Newton, finding that they were approximately four orders of magnitude too stringent due to the use of the instrument's geometric solid angle rather than its exposure-weighted solid angle. Upon using the exposure-weighted solid angle, we show that the revised XMM-Newton limits are slightly weaker than those from eROSITA.

X-rays↗
DOE OSTI2024

Absence of magnetoelastic deformation in the spin-chain compound CuBr 2

Here, we investigate a spin-$\frac{1}{2}$ antiferromagnet, CuBr 2 , which has quasi-one-dimensional structural motifs. The system has previously been observed to exhibit unusual Raman modes possibly due to a locally deformed crystal structure driven by the low-dimensional magnetism. Using hard x-ray scattering and neutron total scattering, here we aim to verify a specific form of tetramerizing deformation proposed in the previous study. Apart from diffuse scattering signals, which we can reproduce by performing a thorough modeling of the lattice's thermal vibrations, we do not observe evidence for a tetramerized lattice structure within our detection sensitivity. Consequently, we consider it unlikely that the unusual Raman modes in CuBr 2 arise from quantum spin-singlet-driven lattice deformations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Magnetized ablative Rayleigh-Taylor instability in three dimensions

Three-dimensional extended-magnetohydrodynamics simulations of the magnetized ablative Rayleigh-Taylor instability are presented. Previous two-dimensional (2D) simulations claiming perturbation suppression by magnetic tension are shown to be misleading, as they do not include the most unstable dimension. For perturbation modes along the applied field direction, the magnetic field simultaneously reduces ablative stabilization and adds magnetic tension stabilization; the stabilizing term is found to dominate for applied fields > 5 T, with both effects increasing in importance at short wavelengths. Additionally, for modes perpendicular to the applied field, magnetic tension does not directly stabilize the perturbation but can result in moderately slower growth due to the perturbation appearing to be 2D (albeit in a different orientation to 2D inertial confinement fusion simulations). In cases where thermal ablative stabilization is dominant the applied field increases the peak bubble-spike height. Resistive diffusion is shown to be important for short wavelengths and long timescales, reducing the effectiveness of tension stabilization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Numerical Improvements in Magnetohydrodynamic, Pulsed Power Simulations of Near-Target Plasmas

Magnetohydrodynamic (MHD) simulations of pulsed power experiments frequently result in unphysical runaway heating when extended to low values of mass density. Traditionally, this has been addressed by the use of conductivity floors below which the plasma is given an arbitrarily low conductivity value. Here, a low-density treatment is presented that allows for low-density material to carry current while maintaining a stable temperature. This treatment is implemented in the Ares multiphysics code. Numerical modifications to the conventional vacuum treatment include an energy-conserving density floor, a modified averaging procedure to determine the thermal conductivity at the edge of adjacent zones, and a modified averaging procedure for determining the thermal conductivity in mixed zones that contain material from differing regions. Additionally, anomalous resistivity and Bohm diffusion models are implemented as simplified models for microphysics-induced enhancement of collisional transport. Importantly, the advantage of these various improvements are illustrated in a simple 1-D pulsed power target where the combined changes result in stable temperatures within the lower density regions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Compact Diffusion Bonded Heat Exchanger Fatigue Life Simulations (Final Report)

Compact diffusion bonded heat exchangers are essential for high pressure heat exchange, but they are subject to thermal fatigue and ramp rate limitations. Simulation of these geometries is challenging with a large range of length and time scales from thousands of mm-sized microchannels inside a m-sized heat exchanger. Multi-physics simulations including thermal, fluid, and solid mechanics components are being used to predict stress within the heat exchangers under these conditions. These predictions can then be used to understand thermal ramp rate limitations while keeping maximum stresses low as well as fatigue life predictions from well-known empirical models.

42 ENGINEERING↗
DOE OSTI2023

Dynamic Formation of a Transient Jet from Arched Magnetized Laboratory Plasma

Abstract A laboratory plasma experiment was built to explore the eruptive behavior of arched magnetized plasmas with dimensionless parameters relevant to the Sun’s photosphere ( β ≈ 10 −3 , Lundquist number ≈10 4 , plasma radius/ion gyroradius ≈20, ion–neutral collision frequency ≫ion cyclotron frequency). Dynamic formation of a transient plasma jet was observed in the presence of the strapping magnetic field. The eruption leading to the jet is unintuitive because the arched plasma is both kink- and torus-stable. The jet structure erupts within a few Alfvén transit times from the formation of the arched plasma. Extensive measurements of plasma temperature, density, magnetic field, and flows are presented. In its early stages, the jet plasma flows away from the arch with supersonic speeds (Mach 1.5). This high-speed flow persists up to the resistive diffusion time in the arched plasma and is driven by large gradients in the magnetic and thermal pressures near the birthplace of jets. There are two distinct electric current channels within the jet, one consisting of outgoing electrons and another composed of electrons returning to the anode footpoint. Significant current density around the jet is a consequence of the diamagnetic current produced by a large thermal pressure gradient in the jet. Ion–neutral charge-exchange collisions provide an efficient mechanism to produce the cross-field current and control the dynamics of the complex current channels of the jet.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Internally-cooled compressor diaphragm

An internally-cooled diaphragm for an internally-cooled compressor is provided. The internally-cooled diaphragm may include an annular body configured to cool a process fluid flowing through a fluid pathway of the internally-cooled compressor. The annular body may define a return channel of the fluid pathway, and a cooling pathway in thermal communication with the fluid pathway. The return channel may be configured to at least partially diffuse and de-swirl the process fluid flowing therethrough, and the cooling pathway may be configured to receive a coolant to absorb heat from the process fluid flowing through the return channel.

42 ENGINEERING↗
DOE OSTI2022

Three-dimensional microstructural characterization of FBR MOX fuel and the contribution of microstructural features to the thermal conductivity of the fuel

Combination of microstructural characterization, property measurements, and phase field modeling is used to investigate fast breeder reactor (FBR) mixed oxide (MOX) fuel irradiated to burnup of 13.7% fissions per initial metal atom (FIMA). Here, the fuel was characterized at different radial locations, which revealed that grey phase can be present in the central region if it nucleates on five metal precipitates (FMPs). In addition, in the mid-radial region FMPs do not diffuse out of the region once formed and the size of Pd–Te precipitates is dictated by the porosity present in the region. Thermal conductivity measurements were conducted as a function of radial location and the microstructure of the fuel was correlated with the observed trend. Reconstructions of the 3D solid and gaseous fission product structures in different regions of the fuel were used to simulate the effective thermal conductivity (ETC) of the respective regions and determine which microstructural feature has the strongest impact on thermal conductivity. Based on conducted assessment, FMPs and Pd-Te precipitates improve local conductivity of central and mid-radial regions even in the presence of grey phase, but defects are primary contributor to the degradation of thermal conductivity on the periphery of the fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Time-domain theory of transient heat conduction in the local limit

Ultrafast and nanoscale heat conduction demands a unified theoretical framework that rigorously bridges macroscopic transport equations with microscopic material properties derived from statistical physics. Existing empirical generalizations of Fourier's law often lack a solid microscopic foundation, failing to connect observed non-Fourier behavior with underlying atomic-scale mechanisms. In this work, we present a time-domain theory of transient heat conduction rooted in Zwanzig's statistical theory of irreversible processes. Central to this framework is the time-domain transport function $\overleftrightarrow{𝑍}$⁡(𝑡) defined through equilibrium time-correlation functions of heat fluxes. This function generalizes the conventional concept of steady-state thermal conductivity, governing the transition of conduction dynamics from onset second sound type wave propagation at finite speeds to diffusion-dominated behavior across broad temporal and spatial scales. Unlike phonon hydrodynamic models that rely on mesoscopic constructs such as phonon drift velocity, our approach provides a quantitative and microscopic description of intrinsic memory effects in transient heat fluxes and applies universally to bulk materials at any temperature or length scale. By integrating atomistic-scale first-principles calculations with continuum-level macroscopic equations, this framework offers a robust foundation for numerical simulations of transient temperature fields. Furthermore, it facilitates the interpretation and design of transient thermal grating experiments using nanometer-scale heat sources and ultrafast laser systems in the extreme ultraviolet and x-ray wavelength ranges, advancing our understanding of heat dissipation dynamics.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

The role of Si in determining the stability of the θ; precipitate in Al-Cu-Mn-Zr alloys

Precipitate-strengthened Al-Cu-Mn-Zr (ACMZ) alloys demonstrate improved microstructural stability compared to conventional Al-Cu alloys. Si is typically considered an impurity in these alloys, but a systematic analysis of ACMZ alloys with a range of Si levels revealed that there is a Si content range (~0.05 wt. % < Si < ~0.10 wt. %) for achieving optimal high temperature microstructural stability. A combination of hardness testing, scanning electron microscopy, and atom probe tomography has been used to understand the relationship between Si content, aging response, and thermal stability in ACMZ alloys. It is shown that the optimum Si content range corresponds to a reduced as-aged hardness, but a greater hardness value retained by the alloy after thermal exposure to 350 °C for 200 hours. A mechanism is proposed in which optimum Si levels reduce the number density of nucleation sites for θ' precipitates, resulting in larger precipitates that, on average, provide a reduced as-aged strength but are more coarsening resistant. This inherent θ' coarsening resistance allows more time for slow-diffusing Mn and Zr to reach the θ' interface and further stabilize the precipitates. Designing the optimal as-aged microstructure for improved thermal stability of the alloy by controlling impurity levels is a concept that has implications beyond the investigated Al-Cu alloy system.

36 MATERIALS SCIENCE↗
DOE OSTI2025

An Investigation of Thermal Properties of 2D Materials [Dissertation]

Studying the thermal conductivity of 2D materials is important due to the applications of 2D materials in fields such as thermal management, thermoelectricity, renewable energy, and sensors. As such, measurements of the thermal conductivity of these 2D materials become important to measure. Thermal conductivity is often difficult to measure for 2D materials due to their atomically thin nature and many experimental methods for doing so requiring contact with the sample, which can alter the thermal properties. A non-contact method for calculating the thermal conductivity of 2D materials supported on substrates in order to model the thermal conductivity of 2D materials for devices, is proposed and experimentally performed in this dissertation. The optothermal Raman technique is a useful non-contact diagnostic technique useful in determining the thermal conductivity of 2D materials. The optothermal Raman typically does not account for heat losses due to convection or radiation or substrate resistance, which are shown to be important factors to consider when developing an optothermal Raman model. Additionally, the calculation of the interfacial thermal conductance between the bottom surface of the sample and the top surface of the substrate, plays an important role in determining the final value of the thermal conductivity of a supported sample, and will yield differing results based on whether or not the conductance is calculated using an approach such as the Diffuse Mismatch Model (DMM) or calculated directly by varying the laser heating profile (usually done by changing the laser objective). This is shown to be the case for both graphene on Ni, graphene on Cu, and SnSe 2 on Cu. In addition to experimentally calculating the thermal conductivity of a 2D material with the optothermal Raman technique, the thermal conductivity of 2D materials can also be calculated using computational methods. The three-phonon method is a method which can be used to simulate phonon scattering processes and determine the thermal conductivity of semiconductors, wherein phonon scattering is the dominant mechanism which determines the thermal conductivity. The three-phonon method uses relaxation times for phonon scattering with other phonons, electrons, and other material system elements, such as isotopes or material defects, in order to create a single-mode relaxation time approximation (SMRTA), which is used to calculate the final value of the thermal conductivity. An important consideration when determining the thermal conductivity of a 2D material using this method is the device geometry, which is reflected in this work as the phonon-boundary scattering relaxation time. This inclusion is important along with the inclusion of phonon-electron scattering in accurately determining the thermal conductivity of a 2D material. In both the optothermal Raman experiments and the three-phonon method computations, strain is shown to have a demonstrable effect on the thermal conductivity of 2D materials. When a 1.1% strain was applied to the mechanical properties of SnSe, the three-phonon processes yielded a lower thermal conductivity than the no-strain case. For the optothermal Raman experiments, the strain induced in the Cu substrate and transferred to a single-layer graphene (SLG) sample yields a trend where the thermal conductivity of the SLG decreases with respect to strain applied. In the case where the interfacial thermal conductance was calculated directly, the conductance increased with respect to strain applied. This presents strain as a reliable and viable method for tuning the thermal properties of 2D materials for device applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Numerical thermalization in 2D PIC simulations: Practical estimates for low-temperature plasma simulations

The process of numerical thermalization in particle-in-cell (PIC) simulations has been studied extensively. It is analogous to Coulomb collisions in real plasmas, causing particle velocity distributions (VDFs) to evolve toward a Maxwellian as macroparticles experience polarization drag and resonantly interact with the fluctuation spectrum. This paper presents a practical tutorial on the effects of numerical thermalization in 2D PIC applications. Scenarios of interest include simulations, which must be run for many thousands of plasma periods and contain a population of cold electrons that leave the simulation space very slowly. This is particularly relevant to many low-temperature plasma discharges and materials processing applications. We present numerical drag and diffusion coefficients and their associated timescales for a variety of grid resolutions, discussing the circumstances under which the electron VDF is modified by numerical thermalization. Though the effects described here have been known for many decades, direct comparison of analytically derived, velocity-dependent numerical relaxation timescales to those of other relevant processes has not often been applied in practice due to complications that arise in calculating thermalization rates in 1D simulations. Using these comparisons, we estimate the impact of numerical thermalization in several examples of low-temperature plasma applications including capacitively coupled plasma discharges, inductively coupled plasma discharges, beam plasmas, and hollow cathode discharges. Finally, we discuss possible strategies for mitigating numerical relaxation effects in 2D PIC simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Molecular simulation of the separation of toluene and p-xylene with the thermally-robust ionic liquid triphenyl-p-phenyl sulfonyl phenyl phosphonium

In this work, we use molecular dynamics simulations to evaluate the behavior of toluene (TOL) and p-xylene (XYL) within three different solvents: n -hexane (HEX), n -heptane (HEP), and a thermally-robust ionic liquid (triphenyl-p-phenyl sulfonyl phenyl phosphonium + bis(trifluoromethylsulfonyl)imide (TPSP + Tf 2 N)). Several different temperatures are explored, while the solvation structures, energetics, and dynamics of TOL and XYL are evaluated. The solvation free energy (ΔG solv ) and transfer free energy (ΔG transfer ) of TOL and XYL are calculated using thermodynamic integration (TI). The transfer free energy of XYL from HEP/HEX to TPSP + Tf 2 N is more favorable than that of TOL over the entire temperature range, while the transfer of both TOL and XYL becomes less favorable as the temperature increases. The diffusion rates of TOL and XYL are lower in the ionic liquid when compared to HEP and HEX, but the thermal stability of TPSP + Tf 2 N may allow for higher operating temperatures and accelerated mass transport rates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗