Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic temperature”

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

Ultrafine Interwoven Dendritic Cu 2 Se/CuFeSe 2 Composites with Enhanced Thermoelectric Performance

As an effort to develop high-performance and stable copper chalcogenide thermoelectrics, several (1-x)Cu 2 Se/(x)CuFeSe 2 composites were successfully synthesized via a one-step solid-state transformation of CuSe 2 template. Uniform ultrafine dendritic structures consisting of interweaving Cu 2 Se and CuFeSe 2 nanofibers were observed in samples with low Fe content (0.05 ≤ x ≤ 0.1). Increasing the Fe content to x = 0.5 led to phase segregation into Cu 2 Se-rich region with embedded CuFeSe 2 fine structures interwoven with CuFeSe 2 -rich region containing Cu 2 Se nanofibers. The formation of such an entwined dendritic structure is believed to arise from the temperature-dependent solubility of CuFeSe 2 in the Cu 2 Se matrix. The dynamic dissolution of CuFeSe 2 into the Cu 2 Se at high temperatures leads to temperature-dependent doping of the Cu 2 Se matrix, enabling drastic enhancements of thermoelectric power factor at high temperatures. As a result, such a strategy is expected to be a powerful tool for properties modulation in thermoelectric materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Numerical Study of Liquid Piston Compression Using Large-Eddy Simulation and Volume-of-Fluid Approach

Efforts to increase the efficiency of residential and commercial air conditioners and heat pumps have demonstrated that the compressor accounts for most of the system’s electrical energy usage. Therefore, the efficiency of this component should be improved to reduce its energy usage. The US Department of Energy’s Oak Ridge National Laboratory developed a near-isothermal liquid piston compressor (LPC) that uses propylene glycol (PG) to compress CO2. This report presents numerical studies of the LPC in which the compression chamber fills with injected PG from the bottom inlet. Numerical simulations were performed using the large-eddy simulation (LES) with the wall-adapting local eddy-viscosity (WALE) subgrid-scale model coupled with the multiphase volume-of-fluid (VOF) model to simulate the transient interface between gas and liquid and to capture the heat and mass transfer within the compression chamber. In this effort, effects of boundary conditions applied to the LES-VOF calculations (e.g., no wall, an adiabatic wall, and a wall with a heat flux subscribed) to the overall pressure and temperature of CO2 gas as well as the transient evolution of flow and heat transfer evolution within the compression chamber are investigated and discussed. It was found that the LES calculation with no wall have shown no dynamical flow patterns and the volume-averaged temperature of CO2 increased from 305 to 392.7 K, while LES calculations with a constant wall temperature or a wall heat flux had similar increases of CO2 temperatures. Results of LES simulation using a wall heat flux showed different stages in the compression process and revealed dynamical formation and interaction of CO2 gas layers and circulation flow patterns within the chamber that contribute to the overall heat transfer between the solid wall, gas, and liquid surface in the compressor.

Nguyen, Thien D.↗
DOE OSTI2020

Ultrafast hot-hole injection modifies hot-electron dynamics in Au/p-GaN heterostructures

A fundamental understanding of hot-carrier dynamics in photo-excited metal nanostructures is needed to unlock their potential for photodetection and photocatalysis. Despite numerous studies on the ultrafast dynamics of hot electrons, so far, the temporal evolution of hot holes in metal-semiconductor heterostructures remains unknown. Here, we report ultrafast (t < 200 fs) hot-hole injection from Au nanoparticles into the valence band of p-type GaN. The removal of hot holes from below the Au Fermi level is observed to substantially alter the thermalization dynamics of hot electrons, reducing the peak electronic temperature and the electron-phonon coupling time of the Au nanoparticles. First-principles calculations reveal that hot-hole injection modifies the relaxation dynamics of hot electrons in Au nanoparticles by modulating the electronic structure of the metal on timescales commensurate with electron-electron scattering. Overall, these results advance our understanding of hot-hole dynamics in metal-semiconductor heterostructures and offer additional strategies for manipulating the dynamics of hot carriers on ultrafast timescales. Photo-excited gold nanoparticles are shown to provide ultrafast and efficient hot-hole injection to the valence band of p-type GaN, substantially altering hot-electron dynamics in the nanoparticles and forming a basis to design hot-hole-based optoelectronics.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Upper bounds for 21st-century surface air temperatures in the Western United States

The last decade has seen a large number of severe heatwaves that were unprecedented in the observational record, highlighting challenges associated with observationally-based statistical quantification of the likelihood and magnitude of future extreme temperatures. An alternative to such probabilistic assessments is identification of upper bounds that quantify the hottest surface air temperatures that can possibly be achieved by the end of the 21st century. Theory, simulations, and observational analyses support the existence of a finite upper bound for surface air temperature; however, estimates for future upper-bound values that are realistic and usable for planning remain unavailable. Here, we combine atmospheric theory with large ensembles of dynamically downscaled projections to estimate historical and end-of-century upper bounds for surface air temperatures. A number of physical mechanisms can influence upper bounds, and at the end of the 21st century, estimates based on mechanisms that yield more moderate upper-bounds produce values around 60∘C for much of the Western United States and in excess of 80∘C for the hottest parts of the domain. Even cooler high-altitude locations have end-of-century upper bounds over 50∘C. Although these upper-bound estimates might seem implausibly large, increases in the upper bounds over the 21st century are similar to increases in dynamically downscaled peak surface temperatures after adjusting those downscaled temperatures to eliminate the possibly biased model trends in surface specific humidity. While upper bound estimates are high relative to historical observations, they nonetheless suggest that heatwave intensity risk is bounded, with uncertainty dominated by projections of surface and upper-level humidity.

Risser, Mark D↗
DOE OSTI2022

A Process-based Model with Temperature, Water, and Lab-derived Data Improves Predictions of Daily Culex pipiens/restuans Mosquito Density

While the number of human cases of mosquito-borne diseases has increased in North America in the last decade, accurate modeling of mosquito population density has remained a challenge. Longitudinal mosquito trap data over the many years needed for model calibration, and validation is relatively rare. In particular, capturing the relative changes in mosquito abundance across seasons is necessary for predicting the risk of disease spread as it varies from year to year. We developed a discrete, semi-stochastic, mechanistic process-based mosquito population model that captures life-cycle egg, larva, pupa, adult stages, and diapause for Culex pipiens (Diptera, Culicidae) and Culex restuans (Diptera, Culicidae) mosquito populations. This model combines known models for development and survival into a fully connected age-structured model that can reproduce mosquito population dynamics. Mosquito development through these stages is a function of time, temperature, daylight hours, and aquatic habitat availability. The time-dependent parameters are informed by both laboratory studies and mosquito trap data from the Greater Toronto Area. The model incorporates city-wide water-body gauge and precipitation data as a proxy for aquatic habitat. This approach accounts for the nonlinear interaction of temperature and aquatic habitat variability on the mosquito life stages. We demonstrate that the full model predicts the yearly variations in mosquito populations better than a statistical model using the same data sources. This improvement in modeling mosquito abundance can help guide interventions for reducing mosquito abundance in mitigating mosquito-borne diseases like West Nile virus.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

A Modeling and Neutron Diffraction Study of the High Temperature Properties of Sub-Stoichiometric Yttrium Hydride for Novel Moderator Applications

Low-enriched-uranium (LEU) reactor systems utilize moderators to improve neutron economy. Solid yttrium hydride is one of the primary moderator candidates for high-temperature (>700 °C) nuclear reactor applications. This is due to its ability to retain hydrogen at elevated temperatures compared to other metal hydrides. For reactor modeling purposes, both neutronic and thermos-mechanical modeling, several high-temperature properties for sub-stoichiometric yttrium hydride (YH2−x) are needed. In this paper, we present an atomistics and a neutron diffraction study of the high-temperature properties of Y and YH2−x. Specifically, we focus on the thermal lattice expansion effects in yttrium metal and yttrium hydride, which also govern bulk thermal expansion. Previously reported physical and mechanical properties for sub-stoichiometric yttrium hydride at ambient conditions are expanded using lattice dynamics to take into account high-temperature effects. Accordingly, an array of newly generated properties is presented that enables high-fidelity neutronics, and thermomechanical modeling. These properties include various elastic moduli, thermal expansion parameters for yttrium and yttrium hydride, and single-phase (YH2−x) and two-phase (Y + YH2−x) density as a function of stoichiometry and density.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Directed Gas-Phase Formation of Azulene (C 10 H 8 ): Unraveling the Bottom-Up Chemistry of Saddle-Shaped Aromatics

The azulene (C 10 H 8 ) molecule, the simplest polycyclic aromatic hydrocarbon (PAH) carrying a fused seven- and five-membered ring, is regarded as a fundamental molecular building block of saddle-shaped carbonaceous nanostructures such as curved nanographenes in the interstellar medium. However, an understanding of the underlying gas-phase formation mechanisms of this nonbenzenoid 10π-Hückel aromatic molecule under low-temperature conditions is in its infancy. Here, by merging crossed molecular beam experiments with electronic structure calculations and molecular dynamics simulations, our investigations unravel an unconventional low-temperature, barrierless route to azulene via the reaction of the simplest organic radical, methylidyne (CH), with indene (C 9 H 8 ) through ring expansion. This reaction might represent the initial step toward to the formation of saddle-shaped PAHs with seven-membered ring moieties in hydrocarbon-rich cold molecular clouds such as the Taurus Molecular Cloud-1 (TMC-1). These findings challenge conventional wisdom that molecular mass growth processes to nonplanar PAHs, especially those containing seven-membered rings, operate only at elevated pressure and high-temperature conditions, thus affording a versatile low-temperature route to contorted aromatics in our galaxy.

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

Importance of Nuclear Quantum Effects on Aqueous Electrolyte Transport under Confinement in Ti 3 C 2 MXenes

Protons display a high chemical activity and strongly affect the charge storage capability in confined interlayer spaces of two-dimensional (2D) materials. As such, an accurate representation of proton dynamics under confinement is important for understanding and predicting charge storage dynamics in these materials. While often ignored in atomistic-scale simulations, nuclear quantum effects (NQEs), e.g., tunneling, can be significant under confinement even at room temperature. Here, using the thermostatted ring polymer molecular dynamics implementation of path integral molecular dynamics (PIMD) in conjunction with the ReaxFF force field, density functional tight binding (DFTB), and NequIP neural network potential simulations, we investigate the role of NQEs on proton and water transport in bulk water and aqueous electrolytes under confinement in Ti 3 C 2 MXenes. Although overall NQEs are relatively small, especially in bulk, we find that they can alter both quantitative values and qualitative trends on both proton transport and water self-diffusion under confinement relative to classical MD predictions. Therefore, our results suggest the need for NQEs to be considered to simulate aqueous systems under confinement for both qualitative and quantitative accuracy.

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

Chain-length-controllable upcycling of polyolefins to sulfate detergents

Escalating global plastic pollution and the depletion of fossil-based resources underscore the urgent need for innovative end-of-life plastic management strategies in the context of a circular economy. Thermolysis is capable of upcycling end-of-life plastics to intermediate molecules suitable for downstream conversion to eventually high-value chemicals, but tuning the molar mass distribution of the products is challenging. Here, in this study, we report a temperature-gradient thermolysis strategy for the conversion of polyethylene and polypropylene into hydrocarbons with tunable molar mass distributions. The whole thermolysis process is catalyst- and hydrogen-free. The thermolysis of polyethylene and polyethylene/polypropylene mixtures with tailored temperature gradients generated oil with an average chain length of ~C 14 . The oil featured a high concentration of synthetically useful α-olefins. Computational fluid dynamics simulations revealed that regulating the reactor wall temperature was the key to tuning the hydrocarbon distributions. Subsequent oxidation of the obtained α-olefins by sulfuric acid and neutralization by potassium hydroxide afforded sulfate detergents with excellent foaming behaviour and emulsifying capacity and low critical micelle concentration. Overall, this work provides a viable approach to producing value-added chemicals from end-of-life plastics, improving the circularity of the anthropogenic carbon cycle.

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

Thermal conductivity tensor of β -1,3,5,7-tetranitro-1,3,5,7-tetrazoctane ( β -HMX) as a function of pressure and temperature

Here, we have used reverse non-equilibrium molecular dynamics (RNEMD) simulations to determine the full thermal conductivity tensor for the monoclinic high explosive crystal β-1,3,5,7-tetranitro-1,3,5,7-tetrazoctane (β-HMX). In order to do so for the monoclinic crystal, four directions for heat propagation are used. Effects of the temperature and pressure are investigated between 200 and 500 K and 0 and 5 GPa, respectively, which approximately covers the range where the β polymorph is stable. Simulations are carried out with the Smith–Bharadwaj non-reactive empirical potential [Smith and Bharadwaj, J. Phys. Chem. B 103, 3570 (1999)], which is known to reproduce well the thermo-elastic properties of HMX. Our results indicate that the thermal conductivity, κ, is highly anisotropic, with 36% difference between the two extreme values at 300 K and 0 GPa. A simple function is used to interpolate κ in the pressure-temperature regime considered in this study, which can be used in continuum models. The results from RNEMD simulations compare well with available experimental results from the literature and allow the determination of κ for any direction and temperature and pressure within or around the fitting interval.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Control-oriented core-SOL-divertor model to address integrated burn and divertor control challenges in ITER

The real-time regulation of a burning plasma’s temperature and density, or burn control, will be necessary to produce high fusion power in future tokamaks like ITER. This is made more challenging due to the plasma’s nonlinear characteristics and the interdependence between the core-plasma and edge-plasma regions. For example, a raising plasma temperature leads to increasing reactivity and therefore to more alpha-particle heating, which further increases temperature. Furthermore, a raise of the fusion power increases the heat flow through the scrape-off-layer (SOL), which can compromise the integrity of the divertor without proper safeguards. For control design, a model-based approach is attractive because it can directly incorporate the nonlinear, coupled, burning-plasma dynamics into the design. To facilitate this design approach, a control-oriented core-SOL-divertor (CSD) model is presented in this work. In this CSD model, a core-plasma model captures the nonlinear dynamics of the core’s density and temperature, and a SOL-divertor model defines the plasma conditions at the separatrix and divertor including the heat load on the target plates. The core-plasma and SOL-divertor models are coupled through the exchange of various variables. In particular, the SOL-divertor model yields the separatrix temperature and the influx of recycled particles into the core-plasma. Further, these variables influence the power and particle balances captured by the core-plasma model. In return, the core-plasma model determines the intensity of the heat and particles fluxes across the separatrix, and this outflow strongly impacts the SOL-divertor model. Therefore, the power and density of the core-plasma, which can be readily modulated through external heating systems and pellet injection, can be viewed as control knobs for the SOL-divertor region in addition to the gas puffing. In simulations of the CSD model, it is demonstrated how external actuation can be utilized to meet burn control and divertor control objectives simultaneously.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Study of the water dynamics near hydrophilic, nanostructured CuO surfaces by quasielastic and inelastic neutron scattering

We have used quasielastic and inelastic neutron scattering to investigate the structure, dynamics, and phase transitions of water interacting with superhydrophilic CuO surfaces that not only possess a strong affinity for water but also a “grass-like” topography that is rough on both micro- and nanoscales. Here, we report quasielastic neutron scattering (QENS) measurements on two samples differing in water content at five temperatures below 280 K. The QENS spectra show water undergoing two different types of diffusive motion near the CuO surfaces: a “slow” translational diffusion occurring on a nanosecond time scale and a faster rotational motion. Further from the surfaces, there is “fast” translational diffusion comparable in rate to that of bulk supercooled water and the rotational motion occurring in the interfacial water. Analysis of the QENS spectra supports wetting of water to the CuO blades as seen in electron microscopy images. In addition, we observe an anomalous temperature dependence of the QENS spectra on cooling from 270 to 230 K with features consistent with a liquid–liquid phase transition. We suggest that the solvent-like properties of the coexisting bulk-like water in our CuO samples are a significant factor in determining the temperature dependence of the interfacial water’s dynamics. Our results are compared with those obtained from two well-studied substrate classes: (1) silicas that contain ordered cylindrical nanopores but have weaker hydrophilicity and (2) nanoparticles of other transition-metal oxides, such as TiO 2 , which share the strong hydrophilicity of our samples but lack their porosity.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Prototype energy models for data centers

Data centers in the United States consume about two percent of the nation’s electricity. Because heat gains from IT equipment drive cooling demand, data centers offer unique opportunities for energy savings. However, no prototype energy model for data centers is available in the suite of existing U.S. Department of Energy’s Commercial Prototype Building Models. Here we present the development of two new data center prototype models and their implementation in OpenStudio and EnergyPlus. The small-size data center model represents a computer room in a building served by computer room air conditioners (CRACs); while the large-sized model represents stand-alone data centers served by computer room air handlers (CRAHs) with a central chiller plant. For each data center model, two levels of IT equipment (ITE) load density were considered, to cover the wide range of IT power density of data centers: 40 and 100 W/ft 2 (430 and 1076 W/m 2 ) for the computer room, and 100 and 500 W/ft 2 (1076 and 5382 W/m 2 ) for the stand-alone data center. All other assumptions, such as building envelope, lighting, HVAC efficiencies and schedules, were based on the minimal requirements of ASHRAE Standard 90.1 at various vintages. We introduced a novel concept of supply and return air approach temperatures to capture the essential effects of non-uniform airflow and temperature distribution in data centers. The approach temperatures were pre-computed by computational fluid dynamics (CFD) simulations for various configurations of ITE loads and airflow containment management in data centers. A new feature was developed in EnergyPlus to implement the approach temperature method. A case study was conducted to demonstrate the use of the data center models. The two data center models cover all U.S. climate zones and can be used to evaluate energy saving measures for data centers, as well as to support development of data center energy efficiency codes and standards.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

A coarse-grained simulation model for colloidal self-assembly via explicit mobile binders

Colloidal particles with mobile binding molecules constitute a powerful platform for probing the physics of self-assembly. Binding molecules are free to diffuse and rearrange on the surface, giving rise to spontaneous control over the number of droplet–droplet bonds, i.e., valence, as a function of the concentration of binders. This type of valence control has been realized experimentally by tuning the interaction strength between DNA-coated emulsion droplets. Optimizing for valence two yields droplet polymer chains, termed ‘colloidomers’, which have recently been used to probe the physics of folding. To understand the underlying self-assembly mechanisms, here we present a coarse-grained molecular dynamics (CGMD) model to study the self-assembly of this class of systems using explicit representations of mobile binding sites. Further, we explore how valence of assembled structures can be tuned through kinetic control in the strong binding limit. More specifically, we optimize experimental control parameters to obtain the highest yield of long linear colloidomer chains. Subsequently tuning the dynamics of binding and unbinding via a temperature-dependent model allows us to observe a heptamer chain collapse into all possible rigid structures, in good agreement with recent folding experiments. Our CGMD platform and dynamic bonding model (implemented as an open-source custom plugin to HOOMD-Blue) reveal the molecular features governing the binding patch size and valence control, and opens the study of pathways in colloidomer folding. This model can therefore guide programmable design in experiments.

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

Dynamically preferred state with strong electronic fluctuations from electrochemical synthesis of sodium manganate

Electrochemical (de)intercalation is a delicate method to precisely control the alkaline ion composition in alkaline transition metal oxides. Because of complicated interactions, metal charge ordering patterns can form spontaneously at special fractional alkaline compositions and orderings. In this study, we show that this elegant electrochemical process can create dynamically preferred structures in an anharmonic energy landscape that conventional syntheses and computations can rarely visit. Specifically, electrochemically prepared Na 1/2 MnO 2 ordering exhibits abnormal structure distortions, charge orderings, and dynamical activities. Strong magnetic fluctuations and lattice dynamics are observed in an unusually wide temperature range in Na 1/2 MnO 2 , which distinguishes it from all other Na x MnO 2 at higher or lower Na compositions. The results emphasize the unique opportunity of using electrochemical processes to design and create novel quantum states with strongly coupled and mutually enhanced electronic and lattice fluctuations, likely through a special dynamic charge flux functional, as suggested by our computational investigations.

25 ENERGY STORAGE↗
DOE OSTI2021

Cation Dynamics in Hybrid Halide Perovskites

Hybrid halide perovskite semiconductors exhibit complex, dynamical disorder while also harboring properties ideal for optoelectronic applications that include photovoltaics. However, these materials are structurally and compositionally distinct from traditional compound semiconductors composed of tetrahedrally coordinated elements with an average valence electron count of silicon. The additional dynamic degrees of freedom of hybrid halide perovskites underlie many of their potentially transformative physical properties. Neutron scattering and spectroscopy studies of the atomic dynamics of these materials have yielded significant insights into their functional properties. Specifically, inelastic neutron scattering has been used to elucidate the phonon band structure, and quasi-elastic neutron scattering has revealed the nature of the uncorrelated dynamics pertaining to molecular reorientations. Understanding the dynamics of these complex semiconductors has elucidated the temperature-dependent phase stability and origins of defect-tolerant electronic transport from the highly polarizable dielectric response. Furthermore, the dynamic degrees of freedom of the hybrid perovskites provide additional opportunities for application engineering and innovation.

Materials Science↗
DOE OSTI2023

Thermofield Theory for Finite-Temperature Electronic Structure

Wave function methods have offered a robust, systematically improvable means to study ground-state properties in quantum many-body systems. Theories like coupled cluster and their derivatives provide highly accurate approximations to the energy landscape at a reasonable computational cost. Analogues of such methods to study thermal properties, though highly desirable, have been lacking because evaluating thermal properties involve a trace over the entire Hilbert space, which is a formidable task. Besides, excited-state theories are generally not as well studied as ground-state ones. In this mini-review, we present an overview of a finite-temperature wave function formalism based on thermofield dynamics to overcome these difficulties. Thermofield dynamics allows us to map the equilibrium thermal density matrix to a pure state, i.e., a single wave function, albeit in an expanded Hilbert space. Ensemble averages become expectation values over this so-called thermal state. Around this thermal state, we have developed a procedure to generalize ground-state wave function theories to finite temperatures. As explicit examples, we highlight formulations of mean-field, configuration interaction, and coupled cluster theories for thermal properties of Fermions in the grand-canonical ensemble. To assess the quality of these approximations, we also show benchmark studies for the one-dimensional Hubbard model, while comparing against exact results. We will see that the thermal methods perform similarly to their ground-state counterparts, while merely adding a prefactor to the asymptotic computational cost. Furthermore, they also inherit all the properties, good or bad, from the ground-state methods, signifying the robustness of our formalism and the scope for future development.

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

Atomic Diffusivities of Yttrium, Titanium and Oxygen Calculated by Ab Initio Molecular Dynamics in Molten 316L Oxide-Dispersion-Strengthened Steel Fabricated via Additive Manufacturing

Oxide-dispersion-strengthened (ODS) steels have long been viewed as a prime solution for harsh environments. However, conventional manufacturing of ODS steels limits the final product geometry, is difficult to scale up to large components, and is expensive due to multiple highly involved, solid-state processing steps required. Additive manufacturing (AM) can directly incorporate dispersion elements (e.g., Y, Ti and O) during component fabrication, thus bypassing the need for an ODS steel supply chain, the scale-up challenges of powder processing routes, the buoyancy challenges associated with casting ODS steels, and the joining issues for net-shape component fabrication. In the AM process, the diffusion of the dispersion elements in the molten steel plays a key role in the precipitation of the oxide particles, thereby influencing the microstructure, thermal stability and high-temperature mechanical properties of the resulting ODS steels. In this work, the atomic diffusivities of Y, Ti, and O in molten 316L stainless steel (SS) as functions of temperature are determined by ab initio molecular dynamics simulations. The latest Vienna Ab initio Simulation Package (VASP) package that incorporates an on-the-fly machine learning force field for accelerated computation is used. At a constant temperature, the time-dependent coordinates of the target atoms in the molten 316L SS were analyzed in the form of mean square displacement in order to obtain diffusivity. The values of the diffusivity at multiple temperatures are then fitted to the Arrhenius form to determine the activation energy and the pre-exponential factor. Given the challenges in experimental measurement of atomic diffusivity at such high temperatures and correspondingly the lack of experimental data, this study provides important physical parameters for future modeling of the oxide precipitation kinetics during AM process.

Chemistry↗