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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Hidden phonon highways promote photoinduced interlayer energy transfer in twisted transition metal dichalcogenide heterostructures

Vertically stacked van der Waals (vdW) heterostructures exhibit unique electronic, optical, and thermal properties that can be manipulated by twist-angle engineering. However, the weak phononic coupling at a bilayer interface imposes a fundamental thermal bottleneck for future two-dimensional devices. Using ultrafast electron diffraction, we directly investigated photoinduced nonequilibrium phonon dynamics in MoS 2 /WS 2 at 4° twist angle and WSe 2 /MoSe 2 heterobilayers with twist angles of 7°, 16°, and 25°. We identified an interlayer heat transfer channel with a characteristic timescale of ~20 picoseconds, about one order of magnitude faster than molecular dynamics simulations assuming initial intralayer thermalization. Atomistic calculations involving phonon-phonon scattering suggest that this process originates from the nonthermal phonon population following the initial interlayer charge transfer and scattering. Our findings present an avenue for thermal management in vdW heterostructures by tailoring nonequilibrium phonon populations.

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

Multiphysics and Multiscale Simulation Methods for Electromagnetic Energy Assisted Fossil Fuel to Hydrogen Conversion (Final Scientific/Technical Report)

This report summarizes the technical accomplishments of the four-year research project “Multiphysics and Multiscale Simulation Methods for Electromagnetic Energy Assisted Fossil Fuel to Hydrogen Conversion” (Award No. DE-FE0032092), conducted at Howard University and the University of Houston (subawardee) from September 2021 to August 2025. The project successfully achieved all four major objectives: 1. 3D Structural Characterization – Developed 3D optical imaging and mechanical sectioning methods to characterize catalyst distribution and support morphology in nickel foam substrates. Successfully reconstructed 3D geometries and imported them into COMSOL Multiphysics for electromagnetic simulations. 2. EM Hotspot Simulation – Created all-frequency stable electromagnetic formulations and 3D nodal discontinuous Galerkin (NDG) methods for coupled electromagnetic-thermal-fluid problems in multiscale catalytic media. Demonstrated stable solutions from DC to microwave frequencies. 3. Multiphysics Coupling – Developed multiscale simulation methods coupling FEM electromagnetic solvers with thermal transport equations. Reactive molecular dynamics (ReaxFF MD) simulations were performed to investigate catalytic reaction mechanisms at the atomistic level. Demonstrated electromagnetic-thermal co-simulation capabilities for porous catalyst structures. 4. System Optimization – Designed and optimized EM-assisted catalytic systems using nickel foam and carbon foam structures, demonstrating significant temperature increases due to microwave heating. Observed and characterized plasma generation in carbon fiber catalysts. Investigated multiple reaction chamber geometries for improved microwave energy deposition. The project produced significant scientific contributions including 15+ peer-reviewed publications, trained multiple Ph.D. students and undergraduate researchers, and advanced the understanding of microwave-assisted hydrogen production from fossil fuels.

08 HYDROGEN↗
DOE OSTI2023

How thermal fluctuations influence the function of the FeMo cofactor in nitrogenase enzymes

The catalytic mechanism of N 2 fixation by nitrogenase remains unresolved in how the strong N≡N bond is activated and why the reductive elimination of H 2 is required. Here, we use density functional theory and physiologically relevant thermal simulations to elucidate the mechanism of the complete nitrogenase catalytic cycle. Over the accumulation of four reducing equivalents, we find that protons and electrons transfer to the FeMo cofactor to weaken and break its bridge Fe–S bond, leading to temporary H 2 S formation that exposes the Fe sites to weakly bind N 2 . Remarkably, we find that subsequent H 2 formation is responsible for chemical activation to an N=N double bond accompanied by a low barrier for H 2 release. We emphasize that finite temperature effects smooth out mechanistic differences between DFT functionals observed at 0 K, thus leading to a consistent understanding as to why H formation is an obligatory step in N 2 adsorption and activation.

DFT↗
DOE OSTI2026

Impact of Forest Canopy Structure on Buoyant Plume Dynamics During Wildland Fires

Heterogeneous forest canopies can generate complex turbulent structures, but in the presence of a fire plume, these interactions are not fully understood. This study investigates the influence of forest canopy heterogeneity on buoyant plume dynamics resulting from surface thermal anomalies representing wildland fires, utilizing Large Eddy Simulation (LES). The Parallelized Large-Eddy Simulation Model (PALM) was employed to simulate six canopy configurations: no canopy, homogeneous canopy, external plume-edge canopy, internal plume-edge canopy, 100 m gap canopy, and 200 m gap canopy. Each configuration was analyzed with and without a static surface heat flux patch of 5000 W ∙ m -2 , resulting in a resting buoyant plume. Simulations were conducted under three crosswind speeds: 0, 5, and 10 m ∙ s -1 . Results show that canopy structure significantly modifies plume behavior, mean flow, and turbulent kinetic energy (TKE) budgets. Plume updraft speed and tilt varied with canopy configuration and crosswind speed. Horizontal pressure gradients associated with plume-atmosphere interaction were modified based on the canopy configuration, resulting in varying crosswind speed reductions at the plume region. Strong momentum absorption was observed above the canopy for the crosswind cases, with the greatest enhancement in the gap canopies. Momentum injection from below the canopy due to the heat source was also observed, resulting in plume structure modulation based on canopy configuration. TKE was found to be the largest in the gap canopy configurations. TKE budget analysis revealed that buoyant production dominated over shear production. At the center of the heat patch, the gap canopy configurations showed enhanced buoyancy within the gap. These results improve our knowledge of fire-canopy-atmosphere interactions that can inform fire models on the impacts of canopy heterogeneity on plume dynamics and ember ejections.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Modeling cell venting and gas-phase reactions in 18650 lithium ion batteries during thermal runaway

A numerical model is developed to study cell venting, internal pressure, and gas-phase dynamics behavior of 18650 Li-ion cells undergoing thermal runaway. A k-e Reynolds-Averaged Navier-Stokes (RANS) model is adopted to describe the turbulent flow out of the cells, while the fluid dynamics inside the cells is described by Darcy-Forchheimer's equation. Thermal abuse reactions and gas generation kinetics are described by a single-step lumped reaction model. Then, a series of computational fluid dynamics (CFD) simulations are conducted on a single 18650 cell at various states-of-charge (100%, 50%, 25%) to study detailed flow and thermal behavior as a function of quantity of gas generated during cell venting. Venting events are categorized into two stages: i) breaching of the cell container and ii) thermal runaway reactions. It is found that the cell response is dominated by the second stage since most of the gases are generated during thermal runaway. Also, the propensity for propagation is highly affected by state-of-charge (SOC). Cells at higher SOCs produce more heat and gas during the venting event, owing to higher mass and concentrations of reacting gases, and consequently reach higher internal cell pressures which increase the risk of side-wall breaching.

25 ENERGY STORAGE↗
DOE OSTI2024

Efficient and Robust Dynamic Crosslinking for Compatibilizing Immiscible Mixed Plastics through In Situ Generated Singlet Nitrenes

Abstract Creating a sustainable economy for plastics demands the exploration of new strategies for efficient management of mixed plastic waste. The inherent incompatibility of different plastics poses a major challenge in plastic mechanical recycling, resulting in phase‐separated materials with inferior mechanical properties. Here, this study presents a robust and efficient dynamic crosslinking chemistry that effectively compatibilizes mixed plastics. Composed of aromatic sulfonyl azides, the dynamic crosslinker shows high thermal stability and generates singlet nitrene species in situ during solvent‐free melt‐extrusion, effectively promoting C─H insertion across diverse plastics. This new method demonstrates successful compatibilization of binary polymer blends and model mixed plastics, enhancing mechanical performance and improving phase morphology. It holds promise for managing mixed plastic waste, supporting a more sustainable lifecycle for plastics.

Chemistry↗
DOE OSTI2021

The benefits of diligence: how precise are predicted gravitational wave spectra in models with phase transitions?

Models of particle physics that feature phase transitions typically provide predictions for stochastic gravitational wave signals at future detectors and such predictions are used to delineate portions of the model parameter space that can be constrained. The question is: how precise are such predictions? Uncertainties enter in the calculation of the macroscopic thermal parameters and the dynamics of the phase transition itself. We calculate such uncertainties with increasing levels of sophistication in treating the phase transition dynamics. Currently, the highest level of diligence corresponds to careful treatments of the source lifetime; mean bubble separation; going beyond the bag model approximation in solving the hydrodynamics equations and explicitly calculating the fraction of energy in the fluid from these equations rather than using a fit; and including fits for the energy lost to vorticity modes and reheating effects. The lowest level of diligence incorporates none of these effects. We compute the percolation and nucleation temperatures, the mean bubble separation, the fluid velocity, and ultimately the gravitational wave spectrum corresponding to the level of highest diligence for three explicit examples: SMEFT, a dark sector Higgs model, and the real singlet-extended Standard Model (xSM). In each model, we contrast different levels of diligence in the calculation and find that the difference in the final predicted signal can be several orders of magnitude. Our results indicate that calculating the gravitational wave spectrum for particle physics models and deducing precise constraints on the parameter space of such models continues to remain very much a work in progress and warrants care.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Advancing stream temperature prediction with a generalizable large-sample framework across CONUS river reaches

Accurately predicting stream temperature in ungauged basins remains a critical challenge for water resource management, thermoelectric power plant cooling, and ecosystem conservation. Large-sample machine learning models trained on hundreds of well-monitored river basins have shown remarkable performance; however, such models have yet to be developed solely using forcing data that can be readily extracted to simulate stream temperatures anywhere in the contiguous United States (CONUS). In this study, we present a scalable, large-sample deep learning framework using Long Short-Term Memory (LSTM) networks to simulate daily stream temperatures in ungauged basins across the CONUS. The framework leverages both modeled reanalysis of meteorological and streamflow inputs as well as static attributes available for all 2.7 million CONUS river reaches in the National Hydrography Dataset Plus (NHDPlusV2). By generating dynamical inputs from predefined thermally relevant upstream contributing areas, rather than the entire upstream basin, the model also offers improvements in very large basins where full-basin averaging can dilute the most important influences on stream temperature. Evaluated across 300 basins, the model achieves a median Mean Absolute Error (MAE) of 1.1 °C and a Nash-Sutcliffe Efficiency (NSE) of 0.95 on temporally and spatially distinct test folds—comparable to models trained exclusively using meteorological and streamflow observational data. The flexible, high-performing framework generalizes to any unmonitored river reach without significant regulation or unnatural thermal input immediately upstream, substantially expanding predictive capabilities in data-scarce regions.

Hydrology↗
DOE OSTI2022

Temperature-dependent formation of gradient structures with anomalous hardening in an Al–Si alloy

The temperature effect on forming gradient structure in an Al-Si alloy during surface severe deformation is studied in this report. The intermediate temperature (473 K) produces the steepest gradient structure by an anomalous hardening on the top layer compared to lower (300 K) and higher temperature (673 K) counterparts. Our analysis shows profuse aluminum oxide particles in the top layer of the gradient structure under 473 K and thus lead to anomalous hardening by oxide-dispersion strengthening. Here, the counterintuitive enhancement of strengthening at the intermediate temperature is explained by the dynamic interplay between thermal-driven, mechanical-aided oxidation and the wear-induced loss of surface materials as a function of temperature, which yields a critical processing temperature to achieve the steepest gradient structure.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Finely tunable dynamical coloration using bicontinuous micrometer-domains

Nanostructures similar to those found in the vividly blue wings of Morpho butterflies and colorful photonic crystals enable structural color through constructive interference of light waves. Different from commonly studied structure-colored materials using periodic structures to manipulate optical properties, we report a previously unrecognized approach to precisely control the structural color and light transmission via a novel photonic colloidal gel without long-range order. Nanoparticles in this gel form micrometer-sized bicontinuous domains driven by the microphase separation of binary solvents. This approach enables dynamic coloration with a precise wavelength selectivity over a broad range of wavelengths extended well beyond the visible light that is not achievable with traditional methods. The dynamic wavelength selectivity is thermally tunable, reversible, and the material fabrication is easily scalable.

42 ENGINEERING↗
DOE OSTI2021

The effect of substrate temperature on the critical velocity in microparticle impact bonding

A number of nascent coating and micromanufacturing processes involve particles sprayed at a high velocity, which impact a substrate, deform, and adhere. Successful bonding between the particles and the substrate requires impact velocities higher than the so-called critical adhesion velocity. This critical velocity is influenced by the temperature of the substrate, a variable that we isolate in this work by conducting single-particle impacts on a variable-temperature substrate for three systems (Al–Al, Sn–Sn, and Ti–Ti) with particles individually selected within a narrow size distribution. Our results quantitatively connect the increase in substrate temperature to a significant lowering of the critical velocity, which we attribute to the lower dynamic strength of the thermally softened substrate. The data are generally consistent with expectations for bonding being controlled by a hydrodynamic process of jetting upon impact.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Towards modelling AR Sco: calibration – reproducing high-energy pulsar emission and testing convergence to Aristotelian electrodynamics

In recent years, kinetic simulations have been crucial to further our understanding of pulsar electrodynamics. Yet, due to the large-scale separation between the gyro-period and the stellar rotation period, resolving the particle gyration has been computationally unfeasible for realistic pulsar parameters. The main aim of this work is comparing our gyro-phase-resolved model with a gyro-centric pulsar model, where our model solves the general equations of motion with included radiation reaction using a higher order numerical solver with adaptive time-steps. Specifically, we aim to (i) reproduce a pulsar’s high-energy emission maps, namely one with 10 per cent of the surface B-field strength of Vela, and the spectra produced by an independent gyro-centric pulsar emission model; and (ii) test convergence of these results to the radiation-reaction limit of Aristotelian electrodynamics. (iii) Additionally, we identify the effect that a large $E_{\parallel }$-field has on the trajectories and radiation calculations. We find that we can reproduce the curvature radiation emission maps and spectra well, using 10 per cent field strengths of the Vela pulsar and injecting our particles at a higher altitude in the magnetosphere. Using sufficiently large $E_{\parallel }$-fields, our numeric results converge to the analytic radiation-reaction limit trajectories. Additionally, we illustrate the importance of accounting for the $\mathbf {E}\times \mathbf {B}$-drift in the particle trajectories and radiation calculations, validating the Harding and collaborators’ model approach. Lastly, we found that our model deals very well with the high-radiation-reaction and high-field regimes present in pulsars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Large spin Hall effect in Si at room temperature

Silicon's weak intrinsic spin-orbit coupling and centrosymmetric crystal structure are a critical bottleneck to the development of Si spintronics, because they lead to an insignificant spin Hall effect (spin current generation) and inverse spin Hall effect (spin current detection). In this work, we undertake current, magnetic field, crystallography dependent magnetoresistance, and magnetothermal transport measurements to study the spin transport behavior in freestanding Si thin films. We observe a large spin Hall magnetoresistance in both $\mathcal{p}$ – Si and $\mathcal{n}$ – Si at room temperature and it is an order of magnitude larger than that of Pt. One explanation of the unexpectedly large and efficient spin Hall effect is spin-phonon coupling instead of spin-orbit coupling. The macroscopic origin of the spin-phonon coupling can be large strain gradients that can exist in the freestanding Si films. This discovery in a light, earth abundant and centrosymmetric material opens a new path of strain engineering to achieve spin dependent properties in technologically highly developed materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Molecular Recognition with Resolution Below 0.2 Angstroms Through Thermoregulatory Oscillations in Covalent Organic Frameworks

Crystalline materials with uniform molecular-sized pores are desirable for a broad range of applications, such as sensors, catalysis, and separations. However, it is challenging to tune the pore size of a single material continuously and to reversibly distinguish small molecules (below 4 angstroms). We synthesized a series of ionic covalent organic frameworks using a tetraphenoxyborate linkage that maintains meticulous synergy between structural rigidity and local flexibility to achieve continuous and reversible (100 thermal cycles) tunability of “dynamic pores” between 2.9 and 4.0 angstroms, with resolution below 0.2 angstroms. This results from temperature-regulated, gradual amplitude change of high-frequency linker oscillations. These thermoelastic apertures selectively block larger molecules over marginally smaller ones, demonstrating size-based molecular recognition and the potential for separating challenging gas mixtures such as oxygen/nitrogen and nitrogen/methane.

crystalline materials↗
DOE OSTI2026

$\overline{\Sigma }^{\pm }$ production in $\text {pp}$ and $\text {p}{-}\text{Pb}$ collisions at $\sqrt{s_{\textrm{NN}}} = 5.02$ TeV with ALICE

The transverse momentum spectra and integrated yields of anti-$Σ$ hyperons ($\overline{\Sigma}^{\pm}$) have been measured in and collisions at $\sqrt{s_{\textrm{NN}}} = 5.02$ TeV with the ALICE experiment. Measurements are performed via the newly accessed decay channel $\overline{\Sigma}^{\pm}$ → $\bar{\textrm{n}}π^±$. A new method of antineutron reconstruction with the PHOS electromagnetic spectrometer is developed and applied to this analysis. The p T spectra of $\overline{\Sigma}^{\pm}$ are measured in the range 0.5 < p T < 3 GeV/c and compared to predictions of the PYTHIA 8, DPMJET, PHOJET, EPOS LHC and EPOS4 models. The EPOS LHC and EPOS4 models provide the best descriptions of the measured spectra both in pp and p-Pb collisions, while models which do not account for multiparton interactions provide a considerably worse description at high p T . The total yields of $\overline{\Sigma }^{\pm }$ in both pp and p-Pb collisions are compared to predictions of the Thermal-FIST model and dynamical models PYTHIA 8, DPMJET, PHOJET, EPOS LHC and EPOS4. All models reproduce the total yields in both colliding systems within uncertainties. The nuclear modification factors R pPb for both $\overline{\Sigma}^{+}$ and $\overline{\Sigma}^{-}$ are evaluated and compared to those of protons, $Λ$ and $Ξ$ hyperons, and predictions of EPOS LHC and EPOS4 models. No deviations of R pPb for $\overline{\Sigma}^{\pm}$ from the model predictions or measurements for other hadrons are found within uncertainties.

Abualrob, I. J. [University of Houston] (ORCID:000↗
DOE OSTI2024

Capture Cavities for the CW Polarized Positron Source Ce+BAF

The initial design of the capture cavities for a continuous wave (CW) polarized positron beam for the Continuous Electron Beam Accelerator Facility (CEBAF) up-grade at Jefferson Lab is presented. A chain of standing wave multi-cell copper cavities inside a solenoid channel are selected to capture positrons in CW mode. The cavity shunt impedance is surveyed by tuning the cavity geometry while considering accommodating large phase space distribution positron beams with large beam pipe radius while ensuring a large enough passband mode separation. The RF field wall loss power and maximum wall loss power density are considered in cavity and waveguide design. A range of design parameters are given for larger system optimization when the capture cavities are considered together with thermal calculation and beam dynamics in next phase of work.

Wang, S.↗