Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “statistical change detection”

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

Detecting Thermally Induced Spinodal Decomposition with Picosecond Ultrasonics in Cast Austenitic Stainless Steels

Given the existential climate crisis faced by mankind and the world, the lifetime and sustainability of nuclear reactors as a carbon-free source of renewable energy depend on the susceptibility of their structural components to environmental degradation. In particular, critical components for light water reactors (LWRs) evolve over decades in service, losing ductility and toughness due to thermal and irradiation aging. Techniques to monitor their health cannot be easily applied in the field due to their destructive, expensive, or immobile nature. Thus, non-destructive evaluation (NDE) methods are sought to monitor and evaluate the health of major LWR components such as core barrels, steam generator tubes, or primary coolant pipes and are often required by policy, such as NRC policy #10-CFR-50.65. In this work we demonstrate the use of gigahertz, non-contact ultrasonics to gauge the state of cast austenitic stainless steels (CASS), used in some of the largest components in LWR primary systems. We do so by linking changes in their surface acoustic wave (SAW) characteristics using transient grating spectroscopy (TGS) to transmission electron microscopy (TEM)-verified evidence of spinodal decomposition and G-phase precipitation. In this thesis, thermal aging is shown to induce SAW peak splitting in spinodally decomposed CASS alloys, correlated strongly with lowered toughness and decreased ductility. Furthermore, statistical testing on the number of SAW peak splits observed show that the second SAW peak significantly appears more frequently and is significantly different in frequency in comparison to counts and frequencies measured in unaged specimens. The ability of this technique to non-destructively detect microstructural degradation at a distance in a predictive manner in the case of CASS motivates extending gigahertz ultrasonics to detect other LWR material degradation modes as an in-vessel inspection technique, such as reactor pressure vessel (RPV) embrittlement. This allows for the greater use of NDE techniques for confident monitoring of LWR structural material health to 80 years and beyond, saving costs by minimizing structural replacements until needed and maximizing energy production by preventing early decommission until necessary.

36 MATERIALS SCIENCE↗
DOE OSTI2022

An Ensemble Approach to Computationally Efficient Radiological Anomaly Detection and Isotope Identification

Radiological source search is a challenging task involving detection and identification of weak sources in a constantly changing radiological background. As of now, many radiological source detection algorithms have been proposed; however, their computational complexity, and hence reliance on power intensive processing units inhibit low-power applications of radiological source search systems. In this work, we introduce the anomaly filter (AF) algorithm; a computationally light, yet effective time-series source detection algorithm based on exponential weighted moving average (EWMA) and Poisson deviance statistics. Then, we demonstrate that the proposed algorithm can be used in ensemble with other more computationally intensive source detection and identification algorithms to achieve both increased detection performance and reduced power consumption. The proposed AF algorithm and the ensemble algorithms were thoroughly benchmarked against several existing source detection and identification algorithms. The results show that the AF algorithm outperforms existing conventional source detection algorithms, and the ensemble approach improves the overall performance of existing source detection and isotope identification algorithms. Furthermore, the AF algorithm and the Non-negative Matrix Factorization approach based source identification (NMF-ID) algorithm were combined and implemented on a singleboard microcontroller and the power consumption was measured. This ensemble algorithm reduced the power consumption of the NMF-ID algorithm almost by a factor of 100, while improving the detection performance of the overall system.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Self-diagnosis of model suitability for continuous measurements of stream-dissolved organic carbon derived from in situ UV–visible spectroscopy

Application of high-frequency monitoring of dissolved organic carbon (DOC) is difficult in instances where training datasets are challenging to develop (e.g., remote locations) and the relationship between optical features and DOC concentration changes due to environmental or landscape shifts (e.g., climate or land-use change). We developed and compared three partial least squares (PLS) models using in situ water level measurements, conductivity, and UV–Vis spectral attenuation to predict DOC. Two site-specific models were developed using data from a hillslope-dominated forest or a low-relief wetland-pond-dominated stream catchment. The third model, using data from both sites, exhibited the best performance (DOC range = 4–15.5 mg C L −1 , mean = 8.38 mg C L −1 , training RMSE = 0.34 mg C L −1 , internal validation RMSE = 0.50 mg C L −1 , external validation RMSE = 2.43 mg C L −1 ). We further demonstrate using PLS model statistics to monitor performance and elucidate when and how models should be updated. These statistics, Hotelling's T 2 and squared prediction errors, are useful consistency checks for the predictions made and detect underlying inconsistencies that, if undetected, can reduce the robustness of DOC prediction. For example, via the T 2 statistic, we identified the summer–autumn transition as a period when DOC composition differed from what was represented in the training dataset. We also determined that elevated SUVA 254 values contributed to the overall bias observed in predictions made during the subsequent year as part of the external validation. This enabled the application of a bias correction that reduced the RMSE from 2.43 to 0.89 mg C L −1 . The method presented here could be applied to future monitoring programs enabling model updates to monitor DOC fluxes accurately from optical datasets (e.g., attenuance or fluorescence) in the face of developing datasets in remote locations or environmental change. In conclusion, implementation of this approach may also identify possible regime shifts or landscape and hydrologic change associated with climate and other environmental changes relevant to terrestrial to aquatic fluxes.

UV-VIS spectroscopy↗
DOE OSTI2021

Mass Spectrometry Adduct Calculator

We describe the Mass Spectrometry Adduct Calculator (MSAC), an automated Python tool to calculate the adduct ion masses of a parent molecule. Here, adduct refers to a version of a parent molecule [M] that is charged due to addition or loss of atoms and electrons resulting in a charged ion, e.g. [M+H] + . MSAC includes a database of 2,341 potential ions and their mass-to-charge ratios (m/z) as extracted from the NIST/EPA/NIH Mass Spectral Library (NIST17), the Global Natural Products Social Molecular Networking Public Spectral Libraries (GNPS), and MassBank of North America (MoNA). The calculator relies on user-selected subsets of the combined database to calculate expected m/z for adducts of molecules supplied as formulas This tool is intended to help researchers create identification libraries to collect evidence for the presence of molecules in mass spectrometry data. While the included adduct database focuses on adducts typically detected during liquid chromatography-mass spectrometry analyses, users may supply their own lists of adducts and charge states for calculating expected m/z. We also analyzed statistics on adducts from spectra contained in the three selected mass spectral libraries. MSAC is freely available at https://github.com/pnnl/MSAC.

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

CAHS: Context-Aware Homology Search

Protein homology search is foundational to bioinformatics: it supports annotation transfer, structure/function inference, and evolutionary analysis over rapidly expanding sequence repositories (e.g., UniProtKB). Profile hidden Markov models (pHMMs), as implemented in HMMER, remain the most widely trusted approach because they provide statistically calibrated E-values; however, their gap behavior is fixed once a profile is trained, despite biological evidence that insertion/deletion tolerance varies across flexible loops and intrinsically disordered regions. We present CAHS (Context-Aware Homology Search), a lightweight query-time adapter for pHMM search that incorporates learned and biologically motivated signals without changing HMMER's downstream search pipeline or its calibrated E-value reporting. Given a query sequence, CAHS computes per-residue representations from a protein language model and a disorder predictor, maps these to profile coordinates, and modulates only match-state transition rows (gap-open and gap-extension probabilities) while preserving Plan7 constraints. We comprehensively evaluate CAHS across six structurally diverse protein families and multi-domain architectures against a 570k-sequence target corpus. CAHS expands detection capability, retrieving thousands of additional remote homologs at relaxed thresholds by maintaining alignment quality through flexible regions. For multi-domain proteins, context-aware modulation resolves 94% of fragmented alignments. Crucially, CAHS preserves hit-set invariance at stringent operating points (E<10-10), demonstrating increased statistical confidence without inflating false positives. Furthermore, sharper statistical distinction between homologs and background noise during early filter stages yields up to a 3.87× acceleration in end-to-end wall-clock time on high-performance computing clusters. Overall, CAHS illustrates a practical AI-for-science design pattern: augmenting a trusted probabilistic model with query-specific learned signals to improve interpretable, reproducible inference in data-rich biology.

Bhattaram, Swethasree [Georgia Institute of Techno↗
DOE OSTI2022

Structural Characterization of a Toppling Rock Slab From Array-Based Ambient Vibration Measurements and Numerical Modal Analysis

Accurate assessments of the internal structure and boundary conditions of unstable rock slopes are imperative for evaluating landslide hazard scenarios. However, instability characterization at depth remains challenging and is often limited by costly or invasive subsurface investigations. Here, we develop a new approach coupling array-based ambient vibration modal analysis and numerical modeling to improve structural characterization of rock slope instabilities at depth. We used ambient noise cross-correlation on 4 hr of seismic data recorded by an array of 30 nodal geophones at a 500-m-long toppling rock slab in Utah, USA to identify modal frequencies between 0.8 and 3.5 Hz and derive modal displacements. We show that transverse and longitudinal bending modes span the length of the instability, indicating an interconnected slab. Further, statistical comparison of field results with outputs from >1,000 finite element models with varying boundary conditions showed that the instability depth varies between 40–70 and 10–20 m in the middle and lateral regions, respectively. Our approach yields new information on the structural conditions of rock cliff and column instabilities at depth, which is not easily obtained by other means but is imperative for change detection monitoring and improved hazard assessments.

58 GEOSCIENCES↗
DOE OSTI2021

Analysis of Long-Term Quality Control Data for a 137 Cs Dosimetry Calibration Source

Strict quality assurance programs are required for many radiological applications, but these seldom exist for verifying dosimetry calibration sources. After initial characterization of a dosimetry calibration facility, quality control procedures are recommended to ensure the early detection of any changes or malfunctions. These also result in refined knowledge about average dose rate and experimental variations in dose delivery. This paper describes the implementation of a phase I quality control protocol for a 137 Cs dosimetry calibration source and includes an analysis of the resulting data collected over a 24-mo period. During this time, substantial data was collected to establish trial control limits. Air kerma rate measurements were obtained using an ion chamber and were adjusted for decay, corrected for ambient temperature, pressure and humidity, and then analyzed using quality control charts. Three variations of rational subgrouping methods were used in order to find assignable causes of error, and Nelson's Rules were followed to detect any non-random statistical variations. Measurements were subgrouped according to same-day measurements in order to detect positional errors as well as atmospheric correction errors. Additionally, measurements were subgrouped according to analogous experimental setups in order to detect failure in equipment or incorrect settings. Both were analyzed using the X-bar and R chart method. Similarly, individuals and moving ranges charts were used to carefully examine each position in order to observe any situational errors that may occur which include timing, positional, or interference errors. Each method was successful in identifying unique out-of-control data points that occurred during the phase I application of forming control limits. Furthermore, over the 24-mo period, enough data points were deemed in-control to establish reliable trial limits. Future experiments will include the phase II application of gaining more reliable measurements in order to fine-tune the limits, as well as performing a designed experiment, where variables are purposefully changed in order to test the variation of the data.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2022

Cross-Layered Distributed Data-Driven Framework for Enhanced Smart Grid Cyber-Physical Security

Smart Grid (SG) research and development has drawn much attention from academia, industry and government due to the great impact it will have on society, economics and the environment. Securing the SG is a considerably significant challenge due the increased dependency on communication networks to assist in physical process control, exposing them to various cyber-threats. In addition to attacks that change measurement values using False Data Injection (FDI) techniques, attacks on the communication network may disrupt the power system's real-time operation by intercepting messages, or by flooding the communication channels with unnecessary data. Addressing these attacks requires a cross-layer approach. In this paper a cross-layered strategy is presented, called Cross-Layer Ensemble CorrDet with Adaptive Statistics(CECD-AS), which integrates the detection of faulty SG measurement data as well as inconsistent network inter-arrival times and transmission delays for more reliable and accurate anomaly detection and attack interpretation. Numerical results show that CECD-AS can detect multiple False Data Injections, Denial of Service (DoS) and Man In The Middle (MITM) attacks with a high F1-score compared to current approaches that only use SG measurement data for detection such as the traditional physics-based State Estimation, Ensemble CorrDet with Adaptive Statistics strategy and other machine learning classification-based detection schemes.

cyber-physical security↗
DOE OSTI2023

Small-Magnitude Seismic Swarms in Central Utah (US): Interactions of Regional Tectonics, Local Structures and Hydrothermal Systems

Swarms in Central Utah are situated in the complex transition between the Basin and Range (BR) province and the Colorado Plateau. Transecting transverse structures, volcanic deposits, and hydrothermal systems complicate the extensional BR horst and graben structures and provide a multitude of plausible triggering mechanisms. Revisiting the catalog of the University of Utah Seismograph Stations (1981–2022), we analyze spatio-temporal patterns and characteristic features of seismic sequences. Swarms with alternating seismicity rates, bursts, and longer swarms with persistent moment release exhibit a remarkable diversity in temporal evolution. Swarm durations do not scale with cumulative seismic moment: swarms lasting less than 1 day can have similar cumulative seismic moments as month-long swarms. Here, we observe stationary swarms re-occurring for years (e.g., Mineral Mountains), as well as singular swarms in low-seismicity areas (e.g., activating a local structure). The swarms show a pronounced heterogeneity in triggering and driving mechanisms, observed in the detailed analysis of exemplary sequences (detections, relocations, moment tensors, waveform-based clustering, and repeater analysis). The 2022 Sevier Valley sequence activated a BR-related normal fault, the first resolved fault plane in the valley since 1983. The 2011 Circleville sequence is interpreted as a swarm triggered by mainshock-aftershock activity characterized by increasing magnitudes, changing rupture mechanisms, and a concentration of highly similar events in the second part of the sequence. By jointly discussing exemplary sequences and catalog statistics, we draw a comprehensive picture of swarm activity and its relation to geothermal and tectonic activity.

58 GEOSCIENCES↗
DOE OSTI2020

Planck 2018 results

We analyse the Planck full-mission cosmic microwave background (CMB) temperature and E -mode polarization maps to obtain constraints on primordial non-Gaussianity (NG). We compare estimates obtained from separable template-fitting, binned, and optimal modal bispectrum estimators, finding consistent values for the local, equilateral, and orthogonal bispectrum amplitudes. Our combined temperature and polarization analysis produces the following final results: f NL local = -0.9 ± 5.1; f NL equil = -26 ± 47; and f NL ortho = -38 ± 24 (68% CL, statistical). These results include low-multipole (4 ≤ ℓ < 40) polarization data that are not included in our previous analysis. The results also pass an extensive battery of tests (with additional tests regarding foreground residuals compared to 2015), and they are stable with respect to our 2015 measurements (with small fluctuations, at the level of a fraction of a standard deviation, which is consistent with changes in data processing). Polarization-only bispectra display a significant improvement in robustness; they can now be used independently to set primordial NG constraints with a sensitivity comparable to WMAP temperature-based results and they give excellent agreement. In addition to the analysis of the standard local, equilateral, and orthogonal bispectrum shapes, we consider a large number of additional cases, such as scale-dependent feature and resonance bispectra, isocurvature primordial NG, and parity-breaking models, where we also place tight constraints but do not detect any signal. The non-primordial lensing bispectrum is, however, detected with an improved significance compared to 2015, excluding the null hypothesis at 3.5 σ . Beyond estimates of individual shape amplitudes, we also present model-independent reconstructions and analyses of the Planck CMB bispectrum. Our final constraint on the local primordial trispectrum shape is g NL local = (-5.8 ± 6.5) × 10 4 (68% CL, statistical), while constraints for other trispectrum shapes are also determined. Exploiting the tight limits on various bispectrum and trispectrum shapes, we constrain the parameter space of different early-Universe scenarios that generate primordial NG, including general single-field models of inflation, multi-field models (e.g. curvaton models), models of inflation with axion fields producing parity-violation bispectra in the tensor sector, and inflationary models involving vector-like fields with directionally-dependent bispectra. Our results provide a high-precision test for structure-formation scenarios, showing complete agreement with the basic picture of the ΛCDM cosmology regarding the statistics of the initial conditions, with cosmic structures arising from adiabatic, passive, Gaussian, and primordial seed perturbations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Resolution dependence of the turbulent atmospheric boundary layer in global storm-resolving climate simulations

The current generation of state-of-the-art global climate models are being run at increasingly higher horizontal resolutions, with the goal of resolving organised deep convection explicitly. How a kilometre-scale resolution impacts the representation of the atmospheric boundary layer is, however, not well known. Using statistical analysis on global fields as well as high-frequency data at selected locations, produced with the Integrated Forecasting System (IFS) model for the Next Generation Earth-system Models (nextGEMS) project, we investigate the horizontal resolution dependence of some boundary-layer processes. We find that a change in resolution from 9 to 2.8 km causes no substantial changes to boundary-layer properties and processes at most of the locations studied, although some global changes are detected that indicate circulation changes. Small changes to the boundary-layer depth and structure are found in the Tropics. The short simulation length and lack of data for optimal boundary-layer analysis limits the conclusions, especially in relation to the connection between the boundary layer and the atmosphere general circulation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Can ERA5 Be Used to Study Mesoscale Convective System Climatological Characteristics?

Mesoscale convective systems (MCSs) produce more than half of tropical rainfall and are central to the global hydrologic cycle. As the climate warms, environments favorable for MCSs may become more common; however, limited observational records hamper understanding of how MCSs respond to variations and changes in their environments. Here, we evaluate how well MCSs are represented in ERA5, a widely used global high‐resolution reanalysis product. Using PyFLEXTRKR, which jointly tracks top‐of‐atmosphere infrared brightness temperature and surface precipitation, we identified MCSs in ERA5 and compared them with those identified in satellite observations using the same detection algorithm. This comparison analysis spans 2007–2020 using hourly data at 0.25° horizontal resolution focusing over the tropics. ERA5 reproduces observed brightness‐temperature statistics and captures the geographic distribution and seasonal and diurnal cycles of MCS cold cloud shields. However, ERA5 precipitation exhibits an intensity bias—too much light rain and too little heavy rain—which shifts the rain‐rate distribution and reduces the frequency of MCSs relative to observations. Within MCSs, ERA5 precipitation exhibits the same pattern of bias, yielding a systematic underestimation of MCS precipitation intensity. Consistent with these biases, ERA5 underestimates the contribution of MCS to tropical rainfall by 25%–34% in key regions. Overall, ERA5 is suitable for studying MCS cold cloud‐shield climatology and evolution, but precipitation‐based MCS characteristics (including event‐level precipitation features and the geospatial distribution of MCS precipitation) should be interpreted with caution. These findings clarify which aspects of MCS behavior are robustly represented in ERA5 for climatological applications.

mesoscale convection↗
DOE OSTI2025

First observations of solar halo gamma rays over a full solar cycle

We analyze 15 years of Fermi-LAT data and produce a detailed model of the Sun’s inverse-Compton scattering emission (solar halo), which is powered by interactions between ambient cosmic-ray electrons and positrons with sunlight. By developing a novel analysis method to analyze moving sources, we robustly detect the solar halo at energies between 31.6 MeV and 100 GeV, and angular extensions up to 45° from the Sun, providing new insight into spatial regions where there are no direct measurements of the Galactic cosmic-ray flux. The large statistical significance of our signal allows us to subdivide the data and provide the first 𝛾-ray probes into the time variation and azimuthal asymmetry of the solar modulation potential, finding time-dependent changes in solar modulation both parallel and perpendicular to the ecliptic plane. Our results are consistent with (but with independent uncertainties from) local cosmic-ray measurements, unlocking new probes into astrophysical processes near the solar surface.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Longitudinal analysis of host protein serum signatures of treatment and recovery in pulmonary tuberculosis

A better understanding of treatment progression and recovery in pulmonary tuberculosis (TB) infectious disease is crucial. This study analyzed longitudinal serum samples from pulmonary TB patients undergoing interventional treatment to identify surrogate markers for TB-related outcomes. Serum that was collected at baseline and 8, 17, 26, and 52 weeks from 30 TB patients experiencing durable cure were evaluated and compared using a sensitive LC-MS/MS proteomic platform for the detection and quantification of differential host protein signatures relative to timepoint. The global proteome signature was analyzed for statistical differences across the time course and between disease severity and treatment groups. A total of 676 proteins showed differential expression in the serum over these timepoints relative to baseline. Comparisons to understand serum protein dynamics at 8 weeks, treatment endpoints at 17 and 26 weeks, and post-treatment at 52 weeks were performed. The largest protein abundance changes were observed at 8 weeks as the initial effects of antibiotic treatment strongly impacted inflammatory and immune modulated responses. However, the largest number of proteome changes was observed at the end of treatment time points 17 and 26 weeks respectively. Post-treatment 52-week results showed an abatement of differential proteome signatures from end of treatment, though interestingly those proteins uniquely significant at post-treatment were almost exclusively downregulated. Patients were additionally stratified based upon disease severity and compared across all timepoints, identifying 461 discriminating proteome signatures. These proteome signatures collapsed into discrete expression profiles with distinct pathways across immune activation and signaling, hemostasis, and metabolism annotations. Insulin-like growth factor (IGF) and Integrin signaling maintained a severity signature through 52 weeks, implying an intrinsic disease severity signature well into the post-treatment timeframe. Previous proteome studies have primarily focused on the 8-week timepoint in relation to culture conversion status. While this study confirms previous observations, it also highlights some differences. The inclusion of additional end of treatment and post-treatment time points offers a more comprehensive assessment of treatment progression within the serum proteome. Examining the expression dynamics at these later time periods will help in the investigation of relapse patients and has provided indicative markers of response and recovery.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Preventing Failures By Dataset Shift Detection in Safety-Critical Graph Applications

Dataset shift refers to the problem where the input data distribution may change over time (e.g., between training and test stages). Since this can be a critical bottleneck in several safety-critical applications such as healthcare, drug-discovery, etc., dataset shift detection has become an important research issue in machine learning. Though several existing efforts have focused on image/video data, applications with graph-structured data have not received sufficient attention. Therefore, in this paper, we investigate the problem of detecting shifts in graph structured data through the lens of statistical hypothesis testing. Specifically, we propose a practical two-sample test based approach for shift detection in large-scale graph structured data. Our approach is very flexible in that it is suitable for both undirected and directed graphs, and eliminates the need for equal sample sizes. Using empirical studies, we demonstrate the effectiveness of the proposed test in detecting dataset shifts. We also corroborate these findings using real-world datasets, characterized by directed graphs and a large number of nodes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Exploring variability in seasonal average and extreme precipitation using unsupervised machine learning.

Focal Area(s): We will use unsupervised machine learning methods to identify and quantify the influence of large scale natural modes of climate variability to gain insight into the observed and simulated seasonal average and extreme precipitation changes. Science Challenge: A recent paper, led by co-PI Mark Risser, finds that although much of the variability in seasonal average and extreme precipitation over CONUS is unforced, the effect of large-scale modes of circulation variability (such as ENSO, AMO, PNA, etc.) can be detected and attributed. However, it is unclear whether or not unsupervised learning methods can (a) replicate this finding or (b) yield insight into possible nonlinear behavior that was not captured in the initial statistical analysis. Further work would entail extending this framework to other global land areas.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A detection of the sea level fingerprint of Greenland Ice Sheet melt

Rapid melting of ice sheets and glaciers drives a unique geometry, or fingerprint, of sea level change. However, the detection of individual fingerprints has been challenging because of sparse observations at high latitudes and the difficulty of disentangling ocean dynamic variability from the signal. Here we predict the fingerprint of Greenland Ice Sheet (GrIS) melt using recent ice mass loss estimates from radar altimetry data and model reconstructions of nearby glaciers and compare this prediction to an independent, altimetry-derived sea surface height trend corrected for ocean dynamic variability in the region adjacent to the ice sheet. A statistically significant correlation between the two fields (P < 0.001) provides an unambiguous observational detection of the near-field sea level fingerprint of recent GrIS melting in our warming world.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Quantifying Microstructure Variability in Laser Powder Bed Fusion 316 L Stainless Steel Microstructures with Spatial Statistics

Here, we have explored data-driven methods for material microstructure quantification that improve sensitivity to microstructural changes compared to traditional approaches. The methods integrate multiple microstructural properties, including grain morphology, crystallographic orientation, and material phase information. The simpler method employs maps of the Euclidean distance transformation metric to evaluate the morphology of grain boundary networks. The more intensive approach employs generalized spherical harmonic mapping for crystallographic orientations, per-pixel phase information, and a variational auto-encoder for dimensionality reduction and results in a multidimensional clustering of by microstructure similarity. Applied to an experimental dataset of additively manufactured steel, both methods detected slight variations in samples produced under nominally identical processing conditions. Both methods were able to distinguish between samples from multiple (nominally identical) builds, while the generalized spherical harmonics-based method could additionally cluster data samples rotated at two orientations on the build plate. The improved sensitivity of the methods, demonstrated through comparison with traditional microstructure characterization techniques, offers advantages for microstructure quantification and comparisons in advanced manufacturing applications.

SS316L↗