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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Thermal transport exceeding bulk heat conduction due to nonthermal micro/nanoscale phonon populations

While classical size effects usually lead to a reduced effective thermal conductivity, we report here that nonthermal phonon populations produced by a micro/nanoscale heat source can lead to enhanced heat conduction, exceeding the prediction from Fourier's law. We study nondiffusive thermal transport by phonons at small distances within the framework of the Boltzmann transport equation (BTE) and demonstrate that the transport is significantly affected by the distribution of phonons emitted by the source. We discuss analytical solutions of the steady-state BTE for a source with a sinusoidal spatial profile, as well as for a three-dimensional Gaussian “hot spot,” and provide numerical results for single crystal silicon at room temperature. If a micro/nanoscale heat source produces a thermal phonon distribution, it gets hotter than that predicted by the heat diffusion equation; however, if the source predominantly produces low-frequency acoustic phonons with long mean free paths, it may get significantly cooler than that predicted by the heat equation, yielding an enhanced heat transport beyond bulk heat conduction.

42 ENGINEERING↗
DOE OSTI2021

Thermal stability of filler‐based ethylene/vinyl acetate/vinyl alcohol terpolymer in nitroplasticizer oxidative environment

Abstract To study the thermal stability of ethylene/vinyl acetate/vinyl alcohol (EVA‐OH) composite with filler particles, systematic thermogravimetric analysis (TGA) investigation was conducted in both nonisothermal and isothermal modes. The effect of polymer concentration on the aging behavior of EVA‐OH was investigated in the nitroplasticizer (NP) environment. Fourier transform infrared spectroscopy was used to probe chemical structural changes in the EVA‐OH polymer before and after thermal aging. The results suggest that filler addition accelerates the rate of NP exudation out of the EVA‐OH composites and deteriorates the thermal stability of NP at moderate temperatures (up to 70°C). The degradation of NP, in turn, accelerates the degradation of EVA‐OH polymer in its composite form through oxidation and hydrolysis indicating the importance of antioxidants in such phase blends.

36 MATERIALS SCIENCE↗
DOE OSTI2025

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.

25 ENERGY STORAGE↗
DOE OSTI2020

The effect of submicron grain size on thermal stability and mechanical properties of high-entropy carbide ceramics

(Hf 0.2 Zr 0.2 Ta 0.2 Nb 0.2 Ti 0.2 )C high-entropy ceramics (HEC) with a submicron grain size of 400 to 600 nm were fabricated by spark plasma sintering using a two-step sintering process. Both X-ray and neutron diffractions confirmed the formation of single-phase with rock salt structure in the as-fabricated (Hf 0.2 Zr 0.2 Ta 0.2 Nb 0.2 Ti 0.2 )C samples. The effect of submicron grain size on the thermal stability and mechanical properties of HEC was investigated. The grain growth kinetics in the fine-grained HEC was small at 1300 and 1600°C, suggesting high thermal stability that was possibly related to the compositional complexity and sluggish diffusion in HEC. Compared to the coarse-grain HEC with a grain size of 16.5 µm, the bending strength and fracture toughness of fine-grained HEC were 25% and 20% higher respectively. The improvement of mechanical properties in fine-grained HEC may be attributed to micromechanistic mechanisms such as crack deflection.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Atomistic simulations of the thermal conductivity of liquids

We present a method based on sinusoidal approach to equilibrium molecular dynamics (SAEMD) to compute the thermal conductivity of liquids Similar to nonequilibrium molecular dynamics, and unlike equilibrium simulations based on the Green-Kubo formalism, the method only requires the calculation of forces and total energies. The evaluation of heat fluxes and energy densities is not necessary, thus offering the promise of efficiently implementing first principles simulations based on density functional theory or deep molecular dynamics. Our approach is a generalization of SAEMD for solids, where the thermal conductivity is computed in the steady state, instead of a transient regime, thus properly taking into account diffusive terms in the heat equation. We present results for liquid water at ambient conditions and under pressure and discuss simulation requirements to obtain converged values of the thermal conductivity as a function of size and simulation time.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Ultra-thin passivation layers in Cu(In,Ga)Se 2 thin-film solar cells: full-area passivated front contacts and their impact on bulk doping

In the search for highly transparent and non-toxic alternative front layers replacing state-of-the-art CdS in Cu(In,Ga)Se 2 thin-film solar cells, alternatives rarely exceed reference devices in terms of efficiency. Full-area ultra-thin aluminium oxide tunnelling layers do not require any contact patterning and thus overcome the main drawback of insulating passivation layers. Even a few monolayers of aluminium oxide can be deposited in a controlled manner by atomic layer deposition, they show excellent interface passivation properties, low absorption, and suitable current transport characteristics on test devices. Depositing a ZnO-based transparent front contact, however, results in extremely poor solar cell performance. The issue is not necessarily a low quality of the alternative front layer, but rather the intricate relation between front layer processing and electronic bulk properties in the absorber layer. We identify three challenges critical for the development of novel front passivation approaches:(i) both Cd and Zn impurities beneficially reduce the high native net dopant concentration in the space charge region, (ii) sputter deposition of ZnO damages the passivation layer resulting in increased interface recombination, (iii) thermal treatments of devices with ZnO layer result in substantial Zn diffusion, which can penetrate the full absorber thickness already at moderate temperatures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Discrete ion stochastic continuum overdamped solvent algorithm for modeling electrolytes

In this paper we develop a methodology for the mesoscale simulation of strong electrolytes. The methodology is an extension of the fluctuating immersed-boundary approach that treats a solute as discrete Lagrangian particles that interact with Eulerian hydrodynamic and electrostatic fields. In both algorithms the immersed-boundary method of Peskin is used for particle-field coupling. Hydrodynamic interactions are taken to be overdamped, with thermal noise incorporated using the fluctuating Stokes equation, including a "dry diffusion" Brownian motion to account for scales not resolved by the coarse-grained model of the solvent. Long-range electrostatic interactions are computed by solving the Poisson equation, with short-range corrections included using an immersed-boundary variant of the classical particle-particle particle-mesh technique. Also included is a short-range repulsive force based on the Weeks-Chandler-Andersen potential. This methodology is validated by comparison to Debye-Hückel theory for ion-ion pair correlation functions, and Debye-Hückel-Onsager theory for conductivity, including the Wien effect for strong electric fields. In each case, good agreement is observed, provided that hydrodynamic interactions at the typical ion-ion separation are resolved by the fluid grid.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Optical pulse-induced ultrafast antiferrodistortive transition in SrTiO 3

The ultrafast dynamics of the antiferrodistortive phase transition in perovskite SrTiO 3 is monitored via time-domain Brillouin scattering. Using femtosecond optical pulses, we initiate a thermally driven tetragonal-to-cubic structural transformation and detect the crystal phase through changes in the frequency of Brillouin oscillations (BO) induced by propagating acoustic phonons. Coupling the measured BO frequency with a spatiotemporal heat diffusion model, we demonstrate that, for a sample kept in the tetragonal phase, deposition of sufficient thermal energy induces a rapid transformation of the heat-affected region to the cubic phase. The initial phase change is followed by a slower reverse cubic-to-tetragonal phase transformation occurring on a timescale of hundreds of picoseconds. Here, we attribute this ultrafast phase transformation in the perovskite to a structural resemblance between atomic displacements of the R-point soft optic mode of the cubic phase and the tetragonal phase, both characterized by anti-phase rotation of oxygen octahedra. The structural relaxation time exhibits a strong temperature dependence consistent with the prediction of the equation of motion describing collective oxygen octahedra rotation based on the energy landscape of the phenomenological Landau theory of phase transitions. Evidence of such a fast structural transition in perovskites can open up new avenues in information processing and energy storage sectors.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Advanced Moderation Module for High-Temperature Micro-Reactor Applications

Development and deployment of micro-reactors that provide competitive efficiency, prevailing compactness, and inherent safety is an immediate focus of the U.S. nuclear industry. For thermal neutron micro-reactors operating at elevated temperature for optimized efficiency, such as micro molten-salt reactors (MSRs), heat-pipe reactors, and very-high temperature reactors (VHTRs), high-performance moderator based on metal hydride can enhance the neutron economy and therefore achieve reduced weight and enhanced portability. As unclad metal hydride inevitably decomposes at elevated temperature, an enclosure is required for hydride moderator to deliver desired performance at elevated temperatures. Conventional moderator enclosure solutions based on high-temperature alloys introduce extraneous neutron penalty into the micro-reactor, affecting the neutronic benefits provided by the hydride moderator. Additionally, the compatibility between the high-temperature alloy enclosure and the micro-reactor matrix is also a potential issue. In this report, we disclose an Advanced Moderator Module (AMM) concept enabled by an innovative enclosure solution combining the advantages of refractory metals, ceramic matrix composites (CMCs), and advanced coating technology to serve as hydrogen permeation barrier up to very-high temperatures. A schematic description of the AMM structure is illustrated in Figure A1. The AMM contains a moderating material core made of metal hydride with high thermal stability, such as YH 2-x . The hydride core is enclosed by an H 2 barrier layer coated on a ductile refractory metal liner to minimize hydrogen loss during high-temperature operation. A ceramic matrix composite (CMC) cladding is adopted to provide further structural strength, especially during power transients. Between the CMC cladding and metal liner, an extra diffusion barrier coating is inserted to suppress the chemical interaction at elevated temperatures. Hence, based on a series of innovative material solutions, the AMM is capable of containing the metal hydride core at elevated temperature (>900oC) inside coated and lined CMC envelop with negligible hydrogen loss throughout the microreactor lifetime. The benefit of the AMM technology was assessed based on a comprehensive multi-stage reactor-physics analysis. Advanced moderation based on hydride metals (such as YH2-x) enables reaching optimum moderation with higher fuel content than traditional VHTR technology, which is required to design compact micro-reactor cores targeting long-life operation. The AMM enclosure technology provides lower thermal neutron absorption rates than traditional solutions (based on stainless steel or Mo-based alloy such as TZM), which reduces the fissile enrichment requirements by 6-8% on a TRISO-fueled design based on the EMPIRE core. Finally, combining the hydride moderator with neutron transparent enclosure solutions provides significant potential to boost neutronics performance of micro-reactors in terms of increased core lifetime or reduced size and weight by 30-50% on a micro-reactor based on the Holos Quad technology. The progress and plans for on-going development and demonstration efforts are also discussed in this report. The current demonstration of AMM is focused on the hydrogen diffusion barrier demonstration under thermal cycling, while radiation tolerance demonstration is planned. The next step of the demonstration plan will be the assembly of miniature AMM for high-temperature testing, which can be used as a prototype for future pilot scale demonstration.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Chemical Mixing Induced by Internal Gravity Waves in Intermediate-mass Stars

Internal gravity waves can cause mixing in the radiative interiors of stars. We study this mixing by introducing tracer particles into 2D hydrodynamic simulations. Following the work of Rogers & McElwaine, we extend our study to different masses (3, 7, and 20 M ⊙ ) and ages (ZAMS, midMS, and TAMS). The diffusion profiles of these models are influenced by various parameters such as the Brunt–Väisälä frequency, density, thermal damping, the geometric effect, and the frequencies of waves contributing to these mixing profiles. We find that the mixing profile changes dramatically across age. In younger stars, we noted that the diffusion coefficient increases toward the surface, whereas in older stars the initial increase in the diffusion profile is followed by a decreasing trend. We also find that mixing is stronger in more massive stars. Hence, future stellar evolution models should include this variation. In order to aid the inclusion of this mixing in 1D stellar evolution models, we determine the dominant waves contributing to these mixing profiles and present a prescription that can be included in 1D models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Developing the Science Basis for Understanding Polymer Encapsulant Degradation Mechanisms: DuraMAT 2.0 Final Project Report

Polymeric encapsulants are essential materials in photovoltaic modules, protecting sensitive electronics from the environment while providing mechanical integrity to the multilayered assembly. However, these polymeric materials are susceptible to degradation processes driven by the ingress of environmental species, ultraviolet radiation, thermal stresses, and mechanical loading. In this study, we employ a combined atomistic simulation and accelerated aging experimental approach to study the molecular-scale mechanisms of encapsulant degradation. Classical molecular dynamics simulations quantify the diffusion of environmental and degradation species through the polymer matrix, producing composition-specific diffusion coefficients. Reactive simulations characterize activation energy barriers and reaction rate constants for key chemical pathways. In parallel, thermal-desorption analyses coupled with mass spectrometry monitor the emergence and concentration profiles of degradation products under controlled stressor conditions. By integrating simulation and experiment, we establish quantitative correlations between polymer composition, species diffusivity, and chemical reactivity. We anticipate that these relations and quantitative values could serve as high-fidelity inputs to reaction-diffusion models, enabling physics-informed lifetime predictions and guiding the design of more durable encapsulant materials for solar energy applications.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Digital quantum magnetism on a trapped-ion quantum computer

Digital quantum matter—realized when discrete quantum gates approximate continuous time evolution—is susceptible to heating into chaotic, structureless states. If digitization errors are adequately suppressed, a long-lived transient regime of approximately energy-conserving dynamics can be observed on gate-based quantum computers. Conservation of energy, in turn, enables the exploration of a wide variety of complex behaviours observed in equilibrium systems, ranging from the non-trivial microscopic origins of thermalization itself to the stabilization of effective models hosting exotic emergent properties. Here we use Quantinuum’s H2 quantum computer to simulate digitized dynamics of the quantum Ising model, suppressing digitization errors well enough to observe thermalization on timescales that severely challenge classical simulation methods. Relaxation of an inhomogeneous state reveals an emergent hydrodynamics owing to approximate energy conservation and we compute the associated diffusion constant. By reprogramming our simulations to take place on a triangular lattice with periodic boundary conditions, we observe thermalization consistent with emergent gauge and topological constraints resulting from lattice frustration. Furthermore, our results were enabled by continued advances in two-qubit gate quality (native partial entangler fidelities of 99.94(1)%) and establish digital quantum computers as powerful tools for studying (effectively) continuous-time dynamics.

Information theory and computation↗
DOE OSTI2020

Universal and versatile morphology engineering via hot fluorous solvent soaking for organic bulk heterojunction

After explosive growth of efficiency in organic solar cells (OSCs), achieving ideal morphology of bulk heterojunction remains crucial and challenging for advancing OSCs into consumer market. Herein, by utilizing the amphiphobic nature and temperature-dependent miscibility of fluorous solvent, hot fluorous solvent soaking method is developed to optimize the morphology with various donor/acceptor combinations including polymer/small-molecule, all-polymer and all-small-molecule systems. By immersing blend film into hot fluorous solvent which is utilized as liquid medium with better thermal conductivity, the molecular reorganization is accelerated. Furthermore, fluorous solvent can be miscible with the residue of chloroform and chloronaphthalene above upper critical solution temperature. This mixed solvent diffuses around inside the active layer and selectively promotes molecular reorganization, leading to optimized morphology. Compared to widely-used thermal annealing, this approach processed under mild conditions achieves superior photovoltaic performance, indicating the practicality and universality for morphological optimization in OSCs as well as other optoelectronic devices.

14 SOLAR ENERGY↗
DOE OSTI2025

Quantum Ornstein-Zernike theory for two-temperature two-component plasmas

Laboratory plasma production almost always preferentially heats either the ions or electrons, leading to a two-temperature state. In this state, density functional theory molecular dynamic simulation is the state of the art for modeling bulk material properties. We construct a statistical mechanics model for the two temperature limit that is theoretically consistent with the molecular dynamics method. We proceed to derive the electron-ion multi-temperature quantum Ornstein-Zernike equations for the first time. This allows the construction of a two-temperature two-component plasma model using the average atom from which we can compute bulk material properties at a fraction of the computation time of the two-temperature density functional theory simulation. The accuracy of the model is benchmarked against ion pair correlation and self-diffusion results from ab initio simulation. Here, we proceed to compute the viscosity and ion thermal conductivity as a function of both ion and electron temperature.

Ab initio molecular dynamics↗
DOE OSTI2023

A Data-driven, Physics-based Transport Model of Solar Energetic Particles Accelerated by Coronal Mass Ejection Shocks Propagating through the Solar Coronal and Heliospheric Magnetic Fields

Abstract In an effort to develop computational tools for predicting radiation hazards from solar energetic particles (SEPs), we have created a data-driven physics-based particle transport model to calculate the injection, acceleration, and propagation of SEPs from coronal mass ejection (CME) shocks traversing through the solar corona and interplanetary magnetic fields. The model runs on an input of corona and heliospheric plasma and magnetic field configuration from a magnetohydrodynamic model driven by solar photospheric magnetic field measurements superposed with observed CME shocks determined from coronagraph images. SEP source particles are injected at the shock using the result of diffusive shock acceleration formulation from a characteristic obliquity-dependent injection from a heated solar wind thermal tail population. With several advanced computation techniques involving stochastic simulation and integration, the model obtains the particle intensity at any location in interplanetary space through the rigorous solution to the time-dependent 5D focus transport equation in the phase space that includes perpendicular diffusion. We apply the model to the 2011 November 3 CME event. The calculation results reproduce multispacecraft SEP observations at Earth and STEREO-B reasonably well without normalization of particle flux. The observations at STEREO-A can be reproduced by rescaling particle energy or modified energy dependence of particle diffusion coefficients. This circumsolar SEP event seen by spacecraft at Earth, STEREO-A, and STEREO-B at widely separated longitudes can be explained by diffusive shock acceleration by a single CME shock with a moderate speed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Nonlinear Elimination Applied to Radiation Diffusion

We apply a nonlinearly preconditioned, quasi-Newton framework to accelerate the numerical solution of the thermal radiative transfer (TRT) equations. This framework was inspired by the unpublished method that has existed for years in Teton, Lawrence Livermore National Laboratory’s deterministic TRT code. In this paper, we cast this iteration scheme within a formal nonlinear preconditioning framework and compare its performance against other iteration schemes in the framework. With proper choices of iteration controls for the various levels of the solver, we can recover the standard linearized one-step method, a full nonlinear Newton scheme, as well as the method in Teton. In brief, the nonlinear preconditioning TRT scheme formally eliminates the material temperature equation from the nonlinear system in a nonlinear analog of a Schur complement. This nonlinear elimination step involves solving a decoupled nonlinear equation for each spatial degree of freedom and is therefore inexpensive. By applying a quasi-Newton iteration scheme on the new system, we obtain a three-level iteration scheme that is at least as efficient as commonly used TRT schemes. The new method allows full convergence to the nonlinear backward Euler time-discretized system, increasing accuracy and robustness, while using a similar number of linear iterations as the more common linearized one-step methods Eq. (4).

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2026

Mechanistic transition of stainless steel 304L corrosion in MgCl 2 phase change material under thermal cycling

Magnesium chloride (MgCl 2 )-based molten salt phase change materials (PCMs) are promising candidates for thermal energy storage (TES) owing to their advantageous thermophysical properties, however, their high corrosivity in the molten state limits material compatibility. In this study, the corrosion behavior of 304?L stainless steel (SS304L) exposed to MgCl2 PCM was investigated as a function of thermal cycling (up to 9000 cycles) relevant to PCM-based TES systems. Microstructural changes and elemental redistribution at the alloy-salt interface under thermal cycling conditions were characterized using synchrotron-based X-ray fluorescence (XRF) and scanning electron microscopy coupled with energy dispersive X-ray spectroscopy (SEM-EDS). Results show that corrosion is governed by selective Cr-dissolution and elemental diffusion near the interface and grain boundaries. A clear transition in the dominant corrosion mechanism is observed with increasing thermal cycles. Surface corrosion is the dominant mechanism up to 1000 cycles, whereas intergranular corrosion becomes the primary degradation mode beyond 2000 cycles. These mechanistic changes are interpreted in terms of evolving interfacial chemistry, microstructure and the diffusion-driven redistribution of alloying elements.

36 MATERIALS SCIENCE↗
DOE OSTI2025

In situ investigation of high-pressure hydrogen-induced swelling in elastomers and its correlation with material properties

The resistance of elastomeric materials to high-pressure hydrogen-induced damage is essential for ensuring the reliability of hydrogen infrastructure. Here, in this study, we systematically investigated the swelling behavior and hydrogen transport properties of four elastomer types – EPDM, NBR, FKM, and HNBR – using a custom in-situ view cell system capable of real-time monitoring during decompression from pressures up to 96.5 MPa. Each elastomer was formulated with and without fillers and plasticizers to assess the effects of formulation on swelling response. Thermal desorption analysis (TDA) was employed to determine equilibrium hydrogen content and diffusion coefficients, providing insight into gas uptake and mobility within each material. Correlation analyses using Pearson and Spearman coefficients revealed that the diffusion coefficient showed a stronger relationship with swelling behavior than hydrogen content, highlighting the dominant role of hydrogen mobility. Filled elastomers, particularly those with carbon black, consistently showed reduced swelling due to enhanced stiffness and reduced diffusivity. These results deepen our understanding of diffuso-mechanical interactions in elastomers and support the rational design of sealing materials for high-pressure hydrogen systems.

EPDM↗