Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Technology Trends”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Surface ionization waves propagating over non-planar substrates: wavy surfaces, cut-pores and droplets

Atmospheric pressure plasmas intersecting with dielectric surfaces will often transition into surface ionization waves (SIWs). Several applications of these discharges are purposely configured to be SIWs. During propagation of an SIW over a dielectric surface, the plasma charges the surface while responding to changes in geometrical and electrical material properties. This is particularly important for non-planar surfaces where polarization of the dielectric results in local electric field enhancement. In this paper, we discuss results from computational investigations of negative and positive SIWs propagating over nonplanar dielectrics in three configurations—wavy surfaces, cuts through porous materials and water droplets on flat surfaces. Here we found that negative SIWs are particularly sensitive to the electric field enhancement that occurs at the crests of non-planar surfaces. The local increase in ionization rates by the electric field enhancement can result in the SIW detaching from the surface, which produces non-uniform plasma exposure of the surface. Positive SIWs tend to adhere to the surface to a greater degree. These trends indicate that treatment of pathogen containing droplets on surfaces may be best performed by positive SIWs. The same principles apply to the surfaces cut through pores. Buried pores with small openings to the SIW may be filled by plasma by either flow of plasma into the pore (large opening) or initiated by photoionization (small opening), depending on the size of the opening compared to the Debye length.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Dependence of the boundary heat flux width on core and edge profiles in Alcator C-Mod

This work presents new evidence that the heat flux width, λ q , in the Alcator C-Mod tokamak scales with the edge electron pressure, as observed in the ASDEX Upgrade (AUG) tokamak, but the scaling with volume-averaged pressure, $\bar{p}$, from the plasma stored energy, is a better predictor of λ q in Alcator C-Mod than the edge electron pressure. These previous studies, which find that λ q decreases with increasing plasma pressure, imply that a high performance core at high pressure will lead to challenging heat and particle exhaust due to very small λ q . This concern has led to our significant enlargement of the C-Mod database with the electron density, temperature, and pressure profile data from the Thomson scattering and electron cyclotron emission diagnostics. Using the C-Mod database augmented with new profile data, we find that λ q decreases with increasing edge electron pressure as ${\lambda }_{q}\propto {p}_{\mathrm{e},95}^{-0.26}$, similar to results from AUG, and showing the strength of cross-machine comparisons. We also find that ${\lambda }_{q}\propto {p}_{\mathrm{e},\mathrm{c}\mathrm{o}\mathrm{r}\mathrm{e}}^{-0.56}$, consistent with the original finding from C-Mod that the heat flux width scales as ${\bar{p}}^{-0.48}$. The scalings of λ q with separatrix pressure and gradient scale length are found to match the AUG results qualitatively. The C-Mod scalings with edge plasma quantities have more scatter than the $\bar{p}$ scaling, and, importantly, show different trends for H-modes relative to L- and I-mode. Furthermore, investigating the source of this discrepancy presents an opportunity for further study that may improve our ability to predict the heat flux width in different confinement scenarios in the pursuit of optimizing core-edge performance in future reactors.

Alcator C-Mod↗
DOE OSTI2025

Influence of counterion substitution on the properties of imidazolium-based ionic liquid clusters

Due to their unique physiochemical properties that may be tailored for specific purposes, ionic liquids (ILs) have been investigated for various applications, including chemical separations, catalysis, energy storage, and space propulsion. The different cations and anions comprising ILs may be selected to optimize a range of desired properties, such as thermal stability, ionic conductivity, and volatility, leading to the designation of certain ILs as designer “green” solvents. The effect of counterions on the properties of ILs is of both fundamental scientific interest and technological importance. Herein, we report a systematic experimental and theoretical investigation of the size, charge, stability toward dissociation, and geometric/electronic structure of 1-ethyl-3-methyl imidazolium (EMIM)-based IL clusters containing two different atomic counterions (i.e., bromide [Br − ] and iodide [I − ]). This work extends our studies of EMIM + cations with atomic chloride (Cl − ) and molecular tetrafluoroborate (BF 4 − ) anions reported previously by Baxter et al. [Chem. Mater. 34, 2612 (2022)] and Zhang et al . [J. Phys. Chem. Lett. 11, 6844 (2020)], respectively. Distributions of anionic IL clusters were generated in the gas phase using electrospray ionization and characterized by high mass resolution mass spectrometry, energy-resolved collision-induced dissociation, and negative ion photoelectron spectroscopy experiments. The experimental results reveal anion-dependent trends in the size distribution, relative abundance, ionic charge state, stability toward dissociation, and electron binding energies of the IL clusters. Complementary global optimization theory provides molecular-level insights into the bonding and electronic structure of a selected subset of clusters, including their low energy structures and electrostatic potential maps, and how these fundamental characteristics are influenced by anion substitution. Collectively, our findings demonstrate how the fundamental properties of ILs, which determine their suitability for many applications, may be tuned by substituting counterions. These observations are critical in the sub-nanometer cluster size regime where phenomena do not scale predictably to the bulk phase, and each atom counts toward determining behavior.

cluster↗
DOE OSTI2025

Assessment of metadynamic recrystallization in single copper particle impacts by focused ion beam tomography

We study single Cu-on-Cu impacts relevant to cold spray deposition and quantitatively analyze the metadynamic recrystallization (mDRX) that takes place after the impact by virtue of lingering impact adiabatic heat. Unlike prior studies, the current full 3D tomographic analysis of the mDRX volume shows that mDRX is extremely common in such impacts, although it is often missed when examining 2D sections. We also report an unexpected trend: there is a “sour spot” for mDRX at velocities about 20–40 % above the velocity for particle adhesion. This non-monotonic trend is contrary to the expectations based on increasing adiabatic heating with velocity. With a schematic model, we show that the trend can be explained on the basis of heat transfer: cooling of the heat-affected region is limited by transport through the bonded regions at the particle-substrate interface. Thus, bonding has a prominent role in the heat dissipation process and the best bonded particles most rapidly bulk quench, avoiding mDRX. Here, the developed semi-empirical model aligns with the experimental findings and may help inform microstructural evolution during cold spray and post-spray processing.

FIB-SEM tomography↗
DOE OSTI2025

Shock compression of additively manufactured high solids loaded polymer composites with aligned porosity

The effects of aligned porosity on the shock-compression response of additively manufactured (AM) composite structures with collinear and log-cabin geometries are investigated. The composites contain two inorganic and two organic particle populations and a UV-curable polymer binder. X-ray phase contrast imaging (X-PCI) is used as an in-material diagnostic in plate-impact experiments performed at the Dynamic Compression Sector of the Advanced Photon Source at the Argonne National Laboratory. Macroscale shock and particle velocities and microscale pore collapse velocities are measured from X-PCI images obtained from samples impacted in various orientations relative to the AM build directions. The linear shock and particle velocity relation for the structures with collinear build geometry shows directional differences, with the build direction exhibiting the highest slope and the filament print direction having the shallowest slope. Composites with log-cabin geometry, due to their slightly lower porosity and smaller sized pores, show higher shock velocities than the collinear structures at the same particle velocity. Microscale directionality effects are also observed, with the pore collapse velocity dependent on the pore configuration. Specifically, the pore collapse velocity is highest for pores elongated along the impact direction which is the case for the colinear structures impacted in the filament print direction, which is consistent with the directionality trend observed in the data for the shock and particle velocity equation of state. The pore collapse velocity approaches ∼80 % of the shock velocity and is several times the average bulk particle velocity, suggesting a highly hydrodynamic mechanism of pore collapse. Overall, both the macroscale response as measured by the shock and particle velocity relationship, and the microscale response as measured by the pore collapse velocity, are found to be anisotropic and dependent on the AM composite structure.

Wagner, K. B. [Georgia Institute of Technology, At↗
DOE OSTI2022

Electron temperature gradient driven transport model for tokamak plasmas

A new model for electron temperature gradient (ETG) modes is developed as a component of the multi-mode anomalous transport module [Rafiq et al., Phys Plasmas 20, 032506 (2013)] to predict a time-dependent electron temperature profile in conventional and low aspect ratio tokamaks. This model is based on two-fluid equations that govern the dynamics of low-frequency short- and long-wavelength electromagnetic toroidal ETG driven drift modes. A low collisionality NSTX discharge is used to scan the plasma parameter dependence on the ETG real frequency, growth rate, and electron thermal diffusivity. Electron thermal transport is discovered in the deep core region where modes are more electromagnetic in nature. Several previously reported gyrokinetic trends are reproduced, including the dependencies of density gradients, magnetic shear, β and gradient of β (β'), collisionality, safety factor, and toroidicity, where β is the ratio of the plasma pressure to the magnetic pressure. Furthermore, the electron heat diffusivity associated with the ETG mode is discovered to be on a scale consistent with the experimental diffusivity determined by power balance analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Tuning plasma parameters to control reactive species fluxes to substrates in the context of plasma catalysis

The key reactive plasma-produced gas phase species responsible for the enhanced conversion of chemicals in plasma catalysis compared to thermal catalysis have to date not been identified. This outstanding question is mainly due to the inherent large variety of plasma-produced species and the challenge of controlling and measuring the flux of each constituent of the cocktail of reactive species to a (catalytic) substrate. In this paper, we explore the possibility to control the dominant reactive species fluxes, relevant for plasma-catalysis, to a substrate in the effluent of an RF driven Ar-O 2 plasma jet. The absolute species densities of the major reactive species (O, O 2 (a 1 Δ g ), O 3 and ions) were quantified by molecular beam mass spectrometry (MBMS) to assess the possibility of using treatment distance, O 2 admixture concentration, plasma dissipated power, RF modulation frequency and duty cycle as well as the feed gas flow rate to alter the dominant species densities. Selected experimental results were also compared with a pseudo-1D plug flow model. The short-lived and long-lived species can be effectively separated by changing the treatment distance and the RF modulation frequency. Furthermore, adjusting the O 2 admixture concentration enables to change the ratio of the O 2 (a 1 Δ g ) and O 3 density. The changes in the trend of ion and O flux were found to be very similar for nearly all investigated parameters. Nonetheless the gas flow rate was able to significantly change the ratio of the O and ion density in the plasma jet effluent. Here, the impact of the surface-dependent loss probability and boundary layer reactions on the species flux to a substrate and how this qualitatively relates to the MBMS density measurements is further addressed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Sensor enabled data-driven predictive analytics for modeling and control with high penetration of DERs in distribution systems

The electric power grid is undergoing a tremendous transformation due to the increasing penetration of renewable energy resources beginning with wind and more recently with the distributed energy resources (DERs) such as solar and battery storage. DERs have dramatically changed the role of the distribution systems in the overall power grid, and they are expected to contribute a significant portion of power generation in the future. If current trends for DERs continue, system operation and control will need to change dramatically for improved grid reliability and resiliency. As renewable resources increase in penetration, new and challenging operational, planning, and design problems are expected to emerge. Some of the key challenges that arise in the planning and operation of the future grid are: 1) Quantifying the impact of high DER penetration in distribution systems on bulk grid behavior over multiple time scales. 2) Identifying whether a particular DER configuration/settings have a large impact on the overall grid behavior. These challenges can be addressed in an offline manner using detailed T&D grid models and they can also be addressed in an online manner using sensor measurements. In particular, the advancement and planned growth in sensor technology in power grid over various voltage levels provide us with a unique opportunity to tackle these challenges from a data analytic perspective without needing detailed T&D grid models. A few questions that naturally arise when addressing the challenges from DERs using sensor data are: 1) How can we use limited sensor measurements to monitor & control voltage stability and small signal stability of the bulk system? 2) How can we ensure that the developed data analytic methods are robust to data availability and quality issues? 3) How can we compute the developed analytics in a scalable manner using streaming measurements? In this project, we addressed the aforementioned challenges arising from DERs and answered the questions raised above on how to effectively use the sensor measurements to enhance the reliability and performance of the electric grid. Thus, the overarching goal of this project is to develop effective reduced/representative system models from data that make the computational complexity sufficiently manageable so as to be useful to simulate, analyze, and even control complex non-linear power systems dynamics with large penetrations of DERs. In order to achieve the objective, the project team established a four-fold technical approach 1) Formulated a combined transmission-distribution co-simulation framework for data generation and validation, 2) Derived reduced/representative models of power systems based on data-driven methods for efficient computation and appropriate representation of system behavior, 3) Developed data driven characterization of power system behavior based on transfer operator theory, machine learning and optimization for model estimation, 4) Incorporated a scalable data management and processing architecture using distributed Kafka streaming applications that coordinate input data streams to the developed data analytics. The key accomplishments of the project are: 1) Development of a scalable multi-timescale T&D co-simulation framework (both for steady state and for dynamic co-simulation) using commercial solvers (PSSE and GridLAB-D). The steady-state T&D co-simulation interface is shared with our industry partner (PJM). 2) A structured reduced order dynamic model of distribution systems that can represent partial motor stalling along with a systematic procedure to derive the model parameters. 3) A PMU based online method to monitor, localize and mitigate fault-induced delayed voltage recovery using DER reactive support and load control in distribution systems. 4) Development of linear operator based robust methodologies for dynamic state estimation, uncertainty quantification, system identification and trajectory prediction for power system dynamics. 5) An adaptive damping control for utilizing wind energy resources to provide oscillation damping and system stability. 6) Implementation of Kafka-based framework for efficient processing of streaming data using Linux-based local virtual environment.

DER integration↗
DOE OSTI2025

Excitation of whistler and slow-X waves by runaway electrons in a collisional plasma

Runaway electrons are known to provide robust ideal or collisionless kinetic drive for plasma wave instabilities in both the whistler and slow-X branches, via the anomalous Doppler-shifted cyclotron resonances. In a cold and dense post-thermal-quench plasma, collisional damping of the plasma waves can compete with the collisionless drive. Previous studies have found that, due to their higher wavelength and frequency, slow-X waves suffer stronger collisional damping than the whistlers, while the ideal growth rate of slow-X modes is higher. Here, we study runaway avalanche distributions that maintain the same eigen distribution and increase only in magnitude over time. The distributions are computed from the relativistic Fokker–Planck–Boltzmann solver, upon which a linear dispersion analysis is performed to search for the most unstable or least damped slow-X and whistler modes. Taking into account the effect of plasma density, plasma temperature, and effective charge number, we find that the slow-X modes tend to be excited before the whistlers in a runaway current ramp-up. Furthermore, even when the runaway current density is sufficiently high that both branches are excited, the most unstable slow-X mode has a much higher growth rate than the most unstable whistler mode. The qualitative and quantitative trends uncovered in the current study indicate that even though past experiments and modeling efforts have concentrated on whistler modes, there is a compelling case that slow-X modes should also be a key area of focus in the runaway self-mediation through wave instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

The Kimberlina synthetic multiphysics dataset for CO 2 monitoring investigations

Abstract We present a synthetic multi‐scale, multi‐physics dataset constructed from the Kimberlina 1.2 CO 2 reservoir model based on a potential CO 2 storage site in the Southern San Joaquin Basin of California. Among 300 models, one selected reservoir‐simulation scenario produces hydrologic‐state models at the onset and after 20 years of CO 2 injection. Subsequently, these models were transformed into geophysical properties, including P‐ and S‐wave seismic velocities, saturated density where the saturating fluid can be a combination of brine and supercritical CO 2 , and electrical resistivity using established empirical petrophysical relationships. From these 3D distributions of geophysical properties, we have generated synthetic time‐lapse seismic, gravity and electromagnetic responses with acquisition geometries that mimic realistic monitoring surveys and are achievable in actual field situations. We have also created a series of synthetic well logs of CO 2 saturation, acoustic velocity, density and induction resistivity in the injection well and three monitoring wells. These were constructed by combining the low‐frequency trend of the geophysical models with the high‐frequency variations of actual well logs collected at the potential storage site. In addition, to better calibrate our datasets, measurements of permeability and pore connectivity have been made on cores of Vedder Sandstone, which forms the primary reservoir unit. These measurements provide the range of scales in the otherwise synthetic dataset to be as close to a real‐world situation as possible. This dataset consisting of the reservoir models, geophysical models, simulated time‐lapse geophysical responses and well logs forms a multi‐scale, multi‐physics testbed for designing and testing geophysical CO 2 monitoring systems as well as for imaging and characterization algorithms. The suite of numerical models and data have been made publicly available for downloading on the National Energy Technology Laboratory's (NETL) Energy Data Exchange (EDX) website.

58 GEOSCIENCES↗
DOE OSTI2023

Continuum-Level Modeling of Li-Ion Battery SEI by Upscaling Atomistically Informed Reaction Mechanisms

Understanding and controlling solid-electrolyte interphase (SEI) formation to stabilize cell performance is a significant challenge for next-generation Li-ion battery technologies. In recent years, computational modeling has become an essential tool in providing fundamental insights into SEI properties and dynamics. However, neither atomistic nor continuum-level approaches alone can capture the complexities of SEI chemistry across all relevant length and time scales. In this work, a continuum-level model is developed that is informed by reaction mechanisms obtained from first-principle calculations. The atomistically informed continuum-level model is used to understand electrolyte degradation, including the decomposition of ethylene carbonate (EC), ethyl methyl carbonate (EMC), and fluoroethylene carbonate (FEC). The model presented here is the most chemically complex continuum-level SEI model in the literature to date. The SEI model is calibrated against experimental irreversible leakage currents and shows qualitative agreement with expected SEI growth trends. The model framework is expected to accelerate fundamental understanding of SEI formation, facilitate mechanism development feedback, and dynamically interact with experimental insights.

continuum-level model↗
DOE OSTI2024

How turbulence spreading improves power handling in quiescent high confinement fusion plasmas

Viable magnetic fusion devices necessitate combining good confinement with effective power flux handling. A major concern for ITER, and devices beyond, is the divertor heat load width, which sets peak boundary heat loads on the plasma-facing materials. Current estimates of the heat flux width are narrow for future reactors. Here, we demonstrate how pedestal turbulence can expand into, or entrain, the stable scrape-off-layer and so broaden the heat flux width beyond these neoclassical predictions. Employing combined theoretical, computational, and experimental approaches, we focus on quiescent high confinement discharges on the DIII-D tokamak, but the results are of broader significance. Our findings uncover common trends in the edge turbulence intensity flux, the pressure perturbation skewness, and the turbulence mixing length, which together determine the heat flux width. This research demonstrates the physics of scrape-off-layer broadening by turbulence and highlights the promise of a turbulent pedestal for successful core-edge integration in ITER and future fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Impact of toroidal rotation on the resistive ballooning modes in ASDEX Upgrade tokamak

Abstract In this work, we investigate the behavior of instabilities appearing between type-I edge localized modes (ELMs), with increasing neutral beam injection (NBI) power concomitant increase in toroidal rotation, and compare it to the modeling result of the linear magneto-hydrodynamic (MHD) code CASTOR3D. An injection of one NBI beam, increasing toroidal rotation, results in the mode slowing down from 12 kHz to 7 kHz, and its associated radial displacement decreases from 5 mm to 3.5 mm. In addition, modes shift radially outwards towards higher q , decreasing their poloidal mode numbers. The mode velocity is measured to be close to the E × B velocity with significant uncertainties. Through a set of CASTOR3D simulations with varying profiles, resistivity has been identified as the primary contributor to the growth rates. Only a small stabilizing effect due to toroidal rotation has been observed. While experimental results show a decrease of mode frequency with rotation, the opposite trend is observed in modeling. Reasons for discrepancies between modeling and experiment are discussed. Nevertheless, a main contributor to the mode frequency has been identified to be rotation velocity. CASTOR3D classifies modes as resistive ballooning modes as they do not appear unstable in ideal MHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

High-Throughput Statistical Interrogation of Mechanical Properties with Build Plate Location and Powder Reuse in AlSi10Mg

Abstract Additive manufacturing (AM) allows agile, rapid manufacturing of geometrically complex components that would otherwise be impossible through traditional manufacturing methods. With this maturing manufacturing technology comes the need to adopt testing methods that are commensurate with the speed of additive manufacturing and take advantage of its geometric flexibility. High-throughput tensile testing (HTT) is a technique that allows a large number of tensile bars to be tested in a short amount of time. In the present study, HTT is used to evaluate AM AlSi10Mg produced using powder bed fusion with a Renishaw AM250 machine. Three parameters were varied in this study: (1) powder reuse history, (2) location on the build plate, and (3) size of the tensile specimen. For all parameter combinations, at least 22 specimens were tested; in several cases, over 40 were tested. This large dataset, consisting of over 500 tensile tests, permits Weibull statistical analysis and provides sufficient fidelity to isolate subtle trends that would have likely been missed in smaller, traditional datasets. The observed trends are rationalized in terms of the role of porosity and surface crust on mechanical response.

Carroll, Jay D. (ORCID:0000000258184709)↗
DOE OSTI2021

A multi-agent approach to distribution system fault section estimation in smart grid environment

We report that Multi-Agent Systems (MAS) are seen from different areas as one of the paramount trends for the next generation of power systems. Numerous published studies about MAS discuss its utilization in power distribution networks but none focuses on the prior step to restoration and self-healing that is fault section estimation. This paper aims to show how MAS can improve the utilities’ reliability indexes and consumer satisfaction by overcoming the multiple fault section estimation problem. In order to do this, the authors considered using MAS as a means of communication between smart meters. The purpose of smart meters usage is to employ devices that are already present in smart grids, mainly because of their reading and saving data capacity. The proposed method was tested on a radial feeder generated by the authors. The network was built on HYPERSIM, a software platform of OPAL-RT Technologies. The simulation results show that this MAS provides speed, efficiency, and automation for the process of fault section estimation.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Effect of outer divertor leg detachment on the high field side scrape-off layer in DIII-D and ASDEX Upgrade

In this work, evidence is presented that detachment of the outer divertor leg leads to a reduction of electron density and neutral pressure in the high-field side (HFS) scrape-off layer (SOL) of the ASDEX Upgrade (AUG) and DIII-D tokamaks with ion $B$ x $∇B$ drift directed toward the X-point (favorable configuration). These results are observed across multiple diagnostics and without the use of impurity seeding to reach detachment. In AUG, outer divertor leg detachment correlates with a decrease in electron density near the separatrix at the inner midplane, measured with HFS reflectometry. A concurrent reduction in inner divertor density and neutral pressure at the inner target is observed using divertor Thomson scattering and neutral pressure gauges. These effects are present in both L- and H-mode plasmas. In DIII-D, a similar reduction is detected through analysis of line-integrated hydrogenic emission measured by multiple diagnostics in the HFS SOL close to the separatrix. The consistent trends in both devices indicate that high electron density and strong hydrogenic emissivity in the HFS SOL are common features of H-mode plasmas in the favorable configuration, independent of wall material and divertor geometry. In L-mode plasmas, the reduction in electron density and neutral pressure is not observed in DIII-D, possibly due to differences in wall material. These results emphasize the importance of the divertor state in determining the two-dimensional neutral distribution and edge fueling.

detachment↗
DOE OSTI2022

Aero-Optics of Hypersonic Turbulent Boundary Layers

Aero-optics refers to optical distortions due to index-of-refraction gradients that are induced by aerodynamic density gradients. At hypersonic flow conditions, the bulk velocity is many times the speed of sound and density gradients may originate from shock waves, compressible turbulent structures, acoustic waves, thermal variations, etc. Due to the combination of these factors, aero-optic distortions are expected to differ from those common to sub-sonic and lower super-sonic speeds. This report summarizes the results from a 2019-2022 Laboratory Directed Research and Development (LDRD) project led by Sandia National Laboratories in collaboration with the University of Notre Dame, New Mexico State University, and the Georgia Institute of Technology. Efforts extended experimental and simulation methodologies for the study of turbulent hypersonic boundary layers. Notable experimental advancements include development of spectral de-aliasing techniques for highspeed wavefront measurements, a Spatially Selective Wavefront Sensor (SSWFS) technique, new experimental data at Mach 8 and 14, a Quadrature Fringe Imaging Interferometer (QFII) technique for time-resolved index-of-refraction measures, and application of QFII to shock-heated air. At the same time, model advancements include aero-optic analysis of several Direct Numerical Simulation (DNS) datasets from Mach 0.5 to 14 and development of wall-modeled Large Eddy Simulation (LES) techniques for aero-optic predictions. At Mach 8 measured and predicted root mean square Optical Path Differences agree within confidence bounds but are higher than semi-empirical trends extrapolated from lower Mach conditions. Overall, results show that aero-optic effects in the hypersonic flow regime are not simple extensions from prior knowledge at lower speeds and instead reflect the added complexity of compressible hypersonic flow physics.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Machine Learning Vacancy Formation Energy in Nickel-Based Superalloys

Thermal vacancies play a critical role in high-temperature Ni-based superalloys and influence various properties such as creep resistance, oxidation, etc. This study systematically investigates the impact of commonly used transition metals (Cr, Co, Fe), refractory metals (Nb, Ta, Mo, W) and other elements (Al, Cu, Ti, Mn) on the thermodynamic stability of 36 binary, 20 ternary, 11 quaternary, 9 quinary, and 3 senary FCC Ni-based alloys covering various elemental combinations. Density functional theory-based studies on Ni-X binary alloys show that higher concentrations of Cr, Nb, Ta, Al, and Ti introduce significant lattice distortions and broaden the distribution of vacancy formation energies (standard deviation up to 0.15 eV). These elements partially donate electrons, reducing their self-consistent chemical potentials relative to single-element reference values and lowering vacancy formation energies, while Co, Fe, Mo, and W show lower charge localization. These trends extend from 3-6 element alloys, where Cr, Nb, and Ta-rich compositions have low-energy states (~0.5 eV) that increase vacancy concentrations. Finally, graph neural network models are developed to screen over 5000 virtual alloys. Eleven leading compositions are identified with mean vacancy formation energy higher than 1.75 eV and vacancy concentration ~2 orders of magnitude lower than pure Ni at 1000 K. These results provide valuable guidelines to achieve controlled defect engineering in structural alloys.

DFT↗