Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “3D image analysis”

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

Surface Roughness Calculations of Additively Manufactured Cylinders via X-Ray Computed Tomography

X-ray Computed Tomography (X-ray CT) is an analytical technique used in materials science to non–destructively characterize features in variety of materials like polymers, metals, composites, and explosives. It also has the capability of imaging additively manufactured, machined and assembled parts and processes. The non-destructive imaging allows for the analysis of features (voids and cracks), which give a fundamental understanding of material characteristics. Additionally, X-ray CT can obtain accurate measurements of dimensional and topographic variations as a result of different stimuli, and assess accuracy of material production. This study will focus on parts manufactured via metal additive manufacturing (AM). Although, AM produces parts faster and easier, the printing process can produce defects (pores and surface roughness) that undermine the part’s mechanical properties. The analysis of 3D printed objects has a strength in analysis in that the material has an STL file from which the item was printed, which is not available in many manufactured materials (i.e. foams) due to stochastic structures. For this study, the print accuracy of several 3D printed cylinders will be assessed via X-ray CT to determine optimal printing parameters for parts with less defects and stronger mechanical properties. This will be accomplished by having the original STL file serve as the baseline surface to compute difference in the CT rendering to calculate surface roughness for each cylinders.

36 MATERIALS SCIENCE↗
DOE OSTI2022

X-ray Computed Tomography as a Metrology Technique for the Analysis of Additively Manufactured Material

X-ray computed tomography (X-ray CT) is an analytical technique used in materials science to non-destructively characterize features in a variety materials like polymer, metals, composites, and explosives. It also has the capability of imaging additively manufacture, machine and assembled parts. The non-destructive imaging allows for the analysis of features (voids and cracks), which give a fundamental understanding of the material characteristics. Additionally, X-ray CT can obtain accurate measurements of dimensional and topographic variations due to different stimuli and assess the accuracy of material production. This study focuses on parts manufactured via metal additive manufacturing (AM). Although AM produces parts faster and easier, the printing process can produce defects (pores and surface roughness) that undermine the part’s mechanical properties and performance. The analysis of 3D printed objects has an asset in that the material has an STL file from which the item was printed, which is not available in many manufactured materials (i.e., foams) due to stochastic structures. For this study, the print accuracy of four additively manufactured cylinders will be assessed via X-ray CTto approximate the surface roughness and visualize any major morphological changes to assess the dimensional accuracy of complex additively manufactured parts. It was concluded that using X-ray CT to measure surface roughness was affective because reasonable surface roughness values were measured. Additionally, itwas determined that small-scale features can be produced via additive manufacturing with strong dimensional accuracy so long as the features are highly complex with sharp grooves.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Automated segmentation of soft X-ray tomography: Native cellular structure with submicron resolution at high-throughput for whole-cell quantitative imaging in yeast

Soft X-ray tomography (SXT) is an invaluable tool for quantitatively analyzing cellular structures at suboptical isotropic resolution. However, it has traditionally depended on manual segmentation, limiting its scalability for large datasets. Here, we leverage a deep learning-based autosegmentation pipeline to segment and label cellular structures in hundreds of cells across three Saccharomyces cerevisiae strains. This task-based pipeline uses manual iterative refinement to improve segmentation accuracy for key structures, including the cell body, nucleus, vacuole, and lipid droplets, enabling high-throughput and precise phenotypic analysis. Using this approach, we quantitatively compared the three-dimensional (3D) whole-cell morphometric characteristics of wild-type, VPH1-GFP, and vac14 strains, uncovering detailed strain-specific cell and organelle size and shape variations. We show the utility of SXT data for precise 3D curvature analysis of entire organelles and cells and detection of fine morphological features using surface meshes. Our approach facilitates comparative analyses with high spatial precision and statistical throughput, uncovering subtle morphological features at the single-cell and population level. This workflow significantly enhances our ability to characterize cell anatomy and supports scalable studies on the mesoscale, with applications in investigating cellular architecture, organelle biology, and genetic research across diverse biological contexts.

Chen, Jianhua [Lawrence Berkeley National Laborato↗
DOE OSTI2025

Thermal and kinematic properties of ejecta in SN1987A revealed by XRISM

We present an analysis of high-resolution spectra from the shock-heated plasmas in SN 1987A, based on an observation using the Resolve instrument onboard the X-Ray Imaging and Spectroscopy Mission (XRISM). The 1.7–10 keV Resolve spectra are accurately represented by a single-component, plane-parallel shock plasma model, with a temperature of $2.84_{-0.08}^{+0.09}$ keV and an ionization parameter of $2.64_{-0.45}^{+0.58}$ × $10^{11}\,\,{\rm s\,\, cm}^{-3}$. The Resolve spectra are also well reproduced by the 3D magneto-hydrodynamic simulation presented by Orlando et al. (2020, A&A, 636, A22) suggesting substantial contribution from the ejecta. The metal abundances obtained with Resolve align with the Large Magellanic Cloud value, indicating that the X-rays in 2024 originate from “non-metal-rich” shock-heated ejecta and the reverse shock has not reached the inner metal-rich region of ejecta. Doppler widths of the atomic lines from Si, S, and Fe correspond to velocities of 1500–1700 km s$^{-1}$, where the thermal broadening effects in this non-metal-rich plasma are negligible. Therefore, the line broadening seen in Resolve spectra is determined by the large bulk motion of ejecta. For reference, we determined a $90\%$ upper limit on non-thermal emission from a pulsar wind nebula at $4.3 \times 10^{-13}$ erg cm$^{-2}$ s$^{-1}$ in the 2–10 keV range, aligning with NuSTAR findings by Greco et al. (2022, ApJ, 931, 132). Additionally, we searched for the $^{44}$Sc K line feature and found a $1\sigma$ upper limit of $1.0 \times 10^{-6}$ photons cm$^{-2}$ s$^{-1}$, which translates to an initial $^{44}$Ti mass of approximately $2 \times 10^{-4}\, M_{\odot }$, consistent with previous X-ray to soft gamma-ray observations (Boggs et al. 2015, Science, 348, 670; Grebenev et al. 2012, Nature, 490, 373; Leising 2006, ApJ, 651, 1019).

ISM: supernova remnants↗
DOE OSTI2023

Steptoe Valley NV Data Compilation: Understanding a Stratigraphic Hydrothermal Resource through Geophysical Imaging

Sandia National Laboratories partnered with a multi-disciplinary group of subject matter experts to evaluate a stratigraphic geothermal resource in Steptoe Valley, Nevada using both established and novel geophysical imaging techniques. Provided here are a compilation of newly acquired data over the area and select modeling efforts. This encompasses a 3D geological model (inclusive of full Leapfrog files, Leapfrog viewer files, and XYZ data for faults and stratigraphy) with embedded geophysical modeling, controlled-source electromagnetic (CSEM) and magnetotelluric (MT) data packages, aqueous spring geochemistry data, seismic reflection interpretations, and a gravity data package. The stratigraphic reservoir in Steptoe Valley was previously discovered during oil and gas exploration. Subsequent studies, such as the Nevada Play Fairway Analysis, added data which further highlighted potential resource targets in the basin. Geophysical surveys, complimented with refined geologic mapping and geochemical sampling, were deployed to further characterize the resource. The resulting 3D geologic interpretation, conceptual model refinements, and reservoir simulations suggest that a power-capable reservoir is economically accessible in the Paleozoic carbonates of the deep/central basin. Additional geophysical characterization and exploration drilling efforts are recommended to calibrate interpretation and determine where/how to potentially develop the Steptoe resource. The geophysical tools, interpretations, lessons learned, and publicly available data generated by this study establish an exploration methodology to inform decisions for successful development of stratigraphic reservoirs.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Exploring for Superhot Geothermal Targets in Magmatic Settings: 2022 Field Campaign at Newberry Volcano

This paper presents preliminary results from a subset of work carried out as part of a multinational research project entitled DErisking Exploration for multiple geothermal Plays in magmatic ENvironments (DEEPEN), supported by the U.S. Department of Energy (DOE) and Geothermica, a joint effort by EU member states and associated countries. The DEEPEN project will develop a customized approach to exploration for supercritical and superhot geothermal plays in magmatic systems, which will be applied to two demonstration sites. This paper summarizes field activities carried out at the U.S. demonstration site, Newberry Volcano in central Oregon. The objective of this work effort is to refine the subsurface model of Newberry Volcano, with special focus on deeper zones including the magmatic plumbing system and other key geologic elements. New data collection included gravity and wideband magnetotelluric (MT) surveys, as well as reinstallation of a seismic network. The National Renewable Energy Laboratory (NREL) and Enthalpion Energy LLC (Enthalpion) worked with the Deschutes National Forest Fort Rock District to use a low ground disturbance method of MT deployment to collect MT data inside the caldera and other restricted areas inside the National Volcanic Monument. This opened these areas to geophysical exploration for the first time in decades. Sites along and adjacent to the south rim of the volcano constituted the primary survey objectives. A team from Lawrence Berkeley National Laboratory (LBNL), the U.S. Geological Survey (USGS), and AltaRock also began the process of reinstalling the seismic network from the AltaRock enhanced geothermal system (EGS) demonstration in anticipation of further development activities at the site. The data ingestion, reduction, and analysis phase of the project is ongoing. We are currently processing the MT and gravity data and are developing a new, highly GPU-accelerated, 3D joint MT and gravity inversion to better localize the south rim/south flank conductive target and better understand its relationship to deep heat, fluid sources, and surface extrusive features. Joint inversions, which have not yet been undertaken at Newberry, will allow us to obtain constraints on the geologic model that cannot be determined from each method in isolation, improving our ability to image key geologic features at depth.

geophysics↗
DOE OSTI2025

Genetic algorithm-based geometry calibration for dynamic compression x-ray diffraction experiments

An important component of dynamic compression x-ray diffraction (XRD) experiment analysis is geometry calibration: proper data interpretation requires knowledge of the precise detector position and orientation and, if the experiment involves a single-crystal sample, knowledge of the lattice orientation. The determination of these parameters in the arbitrary three-dimensional (3D) scattering geometries often present in dynamic compression facilities is challenging, as the associated optimization problem can be highly nonlinear, nonsmooth, and discontinuous. We present a genetic algorithm-based approach for performing dynamic compression XRD calibrations that overcomes these obstacles. We provide details regarding the image processing, algorithm implementation, and open-source software deployment and demonstrate the capability of the approach to calibrate the detector and crystal parameters in 3D geometries. Notably, we demonstrate the solver’s capacity to find the crystal orientation without a priori rotation constraints.

Brown, Nathan P. [Sandia National Laboratories (SN↗
DOE OSTI2022

Analysis of Defects in Metal Additive Manufacturing with Augmented Data Generation

Laser powder bed fusion (LPBF) is a method of additive manufacturing (AM) that selectively melts and fuses together microscopic metallic powder. LPBF offers the benefit of producing custom structures out of high strength metals that can be difficult to fabricate with conventional methods. The challenge of LPBF is that 3D printed structures often have internal pores due to process flaws. Pulsed thermal tomography (PTT) is a method for reconstructing the depth profile of materials, allowing the visualization internal voids in solids. In prior work, we developed a convolutional neural network (CNN) which, having been trained on simulated 2D PTT images of subsurface elliptical defects, was able to classify the semi-major radii, semi-minor radii, and angular orientation of the best-fit ellipses in previously unseen PTT images. The unseen PTT images contained subsurface irregular defect shapes imported from scanning electron microscopy (SEM) images of metallic LPBF-printed specimens. Training the CNN on irregular defect shapes instead of on elliptical shapes would make the resulting classifications more descriptive of actual defect shapes. However, this requires a much higher volume of SEM images of material defects, which are difficult to obtain because of random occurrence of defects in LPBF. To address this challenge, we developed a generative adversarial network (GAN) to augment the existing dataset of SEM defect images. The GAN model is demonstrated to create novel yet realistic defect shapes that can be used as input for simulated PTT images to train CNN.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Structure–Property Linkage in Alloys Using Graph Neural Network and Explainable Artificial Intelligence

Deep learning tools have recently shown significant potential for accelerating the prediction of microstructure–property linkage in materials. While deep neural networks like convolution neural networks (CNNs) can extract physics information from 3D microstructure images, they often require a large network architecture and substantial training time. In this research, we trained a graph neural network (GNN) using phase field generated microstructures of Ni-Al alloys to predict the evolution of mechanical properties. We found that a single GNN is capable of accurately predicting the strengthening of Ni-Al alloys with microstructures of varying sizes and dimensions, which cannot otherwise be done with a CNN. Additionally, GNN requires significantly less GPU utilization than CNN and offers more interpretable explanation of predictions using saliency analysis as features are manually defined in the graph. We also utilize explainable artificial intelligence tool Bayesian Inference to determine the coefficients in the power law equation that governs coarsening of precipitates. Overall, our work demonstrates the ability of the GNN to accurately and efficiently extract relevant information from material microstructures without having restrictions on microstructure size or dimension and offers an interpretable explanation.

Chemistry↗
DOE OSTI2023

Slip localization behavior at triple junctions in nickel-base superalloys

Incipient slip localization in the vicinity of hundreds of grain boundary triple junctions (TJs) in a lightly deformed nickel-base superalloy IN718 is studied using a combination of three-dimensional (3D) crystal plasticity finite element (CPFE) modeling, high resolution digital image correlation (HR-DIC) and 3D electron back-scatter diffraction tomography (3D EBSD). A 3D reconstruction method enables identification of thousands of TJs and correspondence of any observed slip bands with their originating TJ lines below the specimen surface. Here, we present a large-scale CPFE model of the experimental 3D microstructure composed of high-fidelity representation of the TJ lines and the boundaries and interiors of the parent grains and use it to calculate the local micromechanical response and slip activity of all TJs at the onset of macroscopic yielding. Statistical analysis of the calculated quantities reveal TJs develop larger stress concentration and grain-average re-orientation than grain interiors and grain boundaries, however no substantial differences in cumulative slip were found among these microstructural regions. We find that TJs with observed slip bands generate lower grain-average re-orientation, fewer active slip systems, and more localized slip on a single system than those without. The distinctions in the reorientation and slip activity are stronger in TJs that experience more intense slip.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Analysis of core asymmetries in inertial confinement fusion implosions using three-dimensional hot-spot reconstruction

Three-dimensional effects play a crucial role during the hot-spot formation in inertial confinement fusion (ICF) implosions. A data analysis technique for 3D hot-spot reconstruction from experimental observables has been developed to characterize the effects of low modes on 3D hot-spot formations. In nuclear measurements, the effective flow direction, governed by the maximum eigenvalue in the velocity variance of apparent ion temperatures, has been found to agree with the measured hot-spot flows for implosions dominated by mode ℓ = 1. Asymmetries in areal-density (ρR) measurements were found to be characterized by a unique cosine variation along the hot-spot flow axis. In x-ray images, a 3D hot-spot x-ray emission tomography method was developed to reconstruct the 3D hot-spot plasma emissivity using a generalized spherical-harmonic Gaussian function. The gradient-descent algorithm was used to optimize the mapping between the projections from the 3D hot-spot emission model and the measured x-ray images along multiple views. Furthermore, this work establishes a platform to analyze 3D low-mode core asymmetries in ICF.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

High-Energy X-ray Tomographic Analysis of Precursor Metal Powders (Ti-6Al-4V) Used for Additive Manufacturing

We utilized high-energy x-ray tomography to characterize Ti-6Al-4V metal powders (both as-received and recycled) used in powder-bed additive manufacturing process. The image processing workflow was developed to process and analyze large amount of data objectively by computer program. The distribution of size and shape of the metal particles as well as defect (mainly porosity) inside the particles was analyzed with the statistical representation and resolution in micrometer. The result revealed that circular-shaped porosity with various sizes could be embedded in the powder particles. These porosities could potentially be transferred to the 3D printed part and critically affect the mechanical performance of the component. The present study shows the effectiveness of characterizing metal powders using x-ray imaging techniques where sufficient number of particles can be sampled within tens of minutes with a minimum sample preparation and high accuracy. Clear structural differences in the as-received and recycled powders were delineated that helps in determining the feasibility of using the recycled powders.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Utilizing computer vision and artificial intelligence algorithms to predict and design the mechanical compression response of direct ink write 3D printed foam replacement structures

Additive Manufacturing (AM) of porous polymeric materials, such as foams, recently became a topic of intensive research due their unique combination of low density, impressive mechanical properties, and stress dissipation capabilities. Conventional methods for fabricating foams rely on complex and stochastic processes, making it challenging to achieve precise architectural control of structured porosity. In contrast, AM provides access to a wide range of printable materials, where precise spatial control over structured porosity can be modulated during the fabrication process enabling the production of foam replacement structures (FRS). Current approaches for designing FRS are based on intuitive understanding of their properties or an extensive number of finite element method (FEM) simulations. These approaches, however, are computationally expensive and time consuming. As such, in this work, we present a novel methodology for determining the mechanical compression response of direct ink write (DIW) 3D printed FRS using a simple cross-sectional image. By obtaining measurement data for a relatively small number of samples, an artificial neural network (ANN) was trained, and a computer vision algorithm was used to make inferences about foam compression characteristics from a single cross-sectional image. Finally, a genetic algorithm (GA) was used to solve the inverse design problem, generating the AM printing parameters that an engineer should use to achieve a desired compression response from a DIW printed FRS. The methods developed herein present an avenue for entirely autonomous design and analysis of additively manufactured structures using artificial intelligence.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Synchrotron‐source micro‐x‐ray computed tomography for examining butterfly eyes

Comparative anatomy is an important tool for investigating evolutionary relationships among species, but the lack of scalable imaging tools and stains for rapidly mapping the microscale anatomies of related species poses a major impediment to using comparative anatomy approaches for identifying evolutionary adaptations. We describe a method using synchrotron source micro-x-ray computed tomography (syn-μXCT) combined with machine learning algorithms for high-throughput imaging of Lepidoptera (i.e., butterfly and moth) eyes. Our pipeline allows for imaging at rates of ~15 min/mm 3 at 600 nm 3 resolution. Image contrast is generated using standard electron microscopy labeling approaches (e.g., osmium tetroxide) that unbiasedly labels all cellular membranes in a species-independent manner thus removing any barrier to imaging any species of interest. To demonstrate the power of the method, we analyzed the 3D morphologies of butterfly crystalline cones, a part of the visual system associated with acuity and sensitivity and found significant variation within six butterfly individuals. Despite this variation, a classic measure of optimization, the ratio of interommatidial angle to resolving power of ommatidia, largely agrees with early work on eye geometry across species. We show that this method can successfully be used to determine compound eye organization and crystalline cone morphology. Our novel pipeline provides for fast, scalable visualization and analysis of eye anatomies that can be applied to any arthropod species, enabling new questions about evolutionary adaptations of compound eyes and beyond.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Defining the scope for altering rice leaf anatomy to improve photosynthesis: a modelling approach

Summary Leaf structure plays an important role in photosynthesis. However, the causal relationship and the quantitative importance of any single structural parameter to the overall photosynthetic performance of a leaf remains open to debate. In this paper, we report on a mechanistic model, eLeaf, which successfully captures rice leaf photosynthetic performance under varying environmental conditions of light and CO 2 . We developed a 3D reaction‐diffusion model for leaf photosynthesis parameterised using a range of imaging data and biochemical measurements from plants grown under ambient and elevated CO 2 and then interrogated the model to quantify the importance of these elements. The model successfully captured leaf‐level photosynthetic performance in rice. Photosynthetic metabolism underpinned the majority of the increased carbon assimilation rate observed under elevated CO 2 levels, with a range of structural elements making positive and negative contributions. Mesophyll porosity could be varied without any major outcome on photosynthetic performance, providing a theoretical underpinning for experimental data. eLeaf allows quantitative analysis of the influence of morphological and biochemical properties on leaf photosynthesis. The analysis highlights a degree of leaf structural plasticity with respect to photosynthesis of significance in the context of attempts to improve crop photosynthesis.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

A Novel Method for Preoperative Positioning of Total Ankle Replacement Using 3D Digital Model

Abstract Objective To establish a digital model of the ankle joint through 3D imaging technology and explore the preoperative placement of ankle replacement prostheses. Methods Computed tomography images of intact ankle joints from 54 cases in the outpatient and inpatient departments of our hospital were collected; according to the INBONE® total ankle system surgery process, the surgery model and surgical osteotomy were finished using MIMICS based on computer simulation method. The shortest distance was measured between the center point and the anterior, posterior, medial, and lateral, respectively, to ensure the precise position of the ankle replacement prosthesis by digital simulation surgery. The relationship between the two variables was analyzed by bivariate correlation analysis. Results The dataset of this study included 48 cases of the sub‐data set (26 males and 22 females) and included 27 cases of left ankle and 21 cases of right ankle. The average medial malleolar angle was 18.67°± 2.87°, the average amount of bone resection was 12.13 ± 1.86 cm 3 , the mid‐anterior distance was 1.72 ± 0.19 cm, the mid‐posterior distance was 2.00 ± 0.19 cm, the ratio of mid‐anterior to mid‐posterior was 0.87, the mid‐medial distance was 1.26 ± 0.17 cm, the mid‐lateral distance was 1.19 ± 0.16 cm, and the ratio of mid‐medial to mid‐lateral was 1.06. After osteotomy, the anteroposterior diameter was 3.73 ± 0.32 cm, the transverse diameter was 2.46 ± 0.27 cm, and the ratio of anteroposterior diameter to transverse diameter was 1.53. In the bottom view, the shape after osteotomy is rectangular. The mid‐anterior distance was strongly negatively correlated with age, the mid‐anterior distance and the amount of bone resection, the mid‐medial distance and the amount of bone resection, the mid‐lateral distance and the amount of bone resection, the mid‐lateral distance and the anteroposterior diameter, the anteroposterior diameter and the transverse diameter were all strongly positively correlated. Conclusion The projection point of the lower tibia centerline on the tibial horizontal osteotomy surface is located at a position slightly anterior to the midpoint of the transverse diameter after ankle arthroplasty. The rational positioning of the total ankle replacement is located at both a position slightly anterior to the midpoint of the transverse diameter and midpoint of the anteroposterior diameter, which can be used as a reference method before total ankle arthroplasty surgery.

Wu, Shi‐xun↗
DOE OSTI2024

Cellular response of keratinocytes to the entry and accumulation of nanoplastic particles

Plastic accumulation in the environment is rapidly increasing, and nanoplastics (NP), byproducts of environmental weathering of bulk plastic waste, pose a significant public health risk. Particles may enter the human body through many possible routes such as ingestion, inhalation, and skin absorption. However, studies on NP penetration and accumulation in human skin are limited. Loss or reduction of the keratinized skin barrier may enhance the skin penetration of NPs. The present study investigated the entry of NPs into a human skin system modeling skin with compromised barrier functions and cellular responses to the intracellular accumulations of NPs. Two in vitro models were employed to simulate human skin lacking keratinized barriers. The first model was an ex vivo human skin culture with the keratinized dermal layer (stratum corneum) removed. The second model was a 3D keratinocyte/dermal fibroblast cell co-culture model with stratified keratinocytes on the top and a monolayer of skin fibroblast cells co-cultured at the bottom. The penetration and accumulation of the NPs in different cell types were observed using fluorescent microscopy, confocal microscopy, and cryogenic electron microscopy (cryo-EM). The cellular responses of keratinocytes and dermal fibroblast cells to stress induced by NPs stress were measured. The genetic regulatory pathway of keratinocytes to the intracellular NPs was identified using transcript analyses and KEGG pathway analysis. The cellular uptake of NPs by skin cells was confirmed by imaging analyses. Transepidermal transport and penetration of NPs through the skin epidermis were observed. According to the gene expression and pathway analyses, an IL-17 signaling pathway was identified as the trigger for cellular responses to internal NP accumulation in the keratinocytes. The transepidermal NPs were also found in co-cultured dermal fibroblast cells and resulted in a large-scale transition from fibroblast cells to myofibroblast cells with enhanced production of α-smooth muscle actin and pro-Collagen Ia. The upregulation of inflammatory factors and cell activation may result in skin inflammation and ultimately trigger immune responses.

36 MATERIALS SCIENCE↗
DOE OSTI2025

STILGAR End-of-Project Report

The Subsurface Tunnel Imaging LeveraGed by Analysis of Rayleigh wave ellipticity (STILGAR) project demonstrated an integrated geophysical approach for detecting, locating, and characterizing underground structural changes using dense seismic arrays and advanced inversion techniques. Field campaigns were conducted at two operational mines—the Redmond salt mine (Utah) and Graymont Pleasant Gap limestone mine (Pennsylvania)—providing real-world testbeds for monitoring anthropogenic subsurface activity. At the Redmond salt mine, seismic interferometry combined with back-projection inversion successfully identified continuous, low-amplitude signals from mining operations. The approach differentiated stationary from migrating anthropogenic sources, captured daily operational cycles, and validated the potential of passive seismic monitoring for remote detection of underground activity. At the Graymont Pleasant Gap mine, two dense seismic deployments in the spring and fall of 2023 generated over 4 TB of high-resolution data. Key outcomes included the relocation of 199 underground and 8 surface explosions with accuracies within tens of meters and the development of a 3D P-wave velocity model using the triple-difference tomography algorithm (tomoTD) that resolved major structural features such as the mine entrance, low-velocity tunnels, and roof-collapse areas. Ambient noise cross-correlation and back-projection analyses revealed persistent sources linked to ongoing mining activity, whereas horizontal-to-vertical spectral ratio (HVSR) and ellipticity studies confirmed stable site responses across seasons and identified soil thickness trends consistent with regional erosional and depositional processes. Checkerboard and sensitivity tests further validated the robustness of the tomographic results. Overall, the findings emphasize that although significant progress has been made in subsurface imaging, further work is needed to enhance the detection and localization of underground structures. Accurate imaging requires higher frequencies, yet anthropogenic sources tend to dominate the seismic record at those frequencies, and high-frequency surface waves are affected by higher modes that complicate interpretation. The improved detection and localization of human-induced signals enabled detailed temporal and spatial mapping of daily mine operations, demonstrating the feasibility of continuous anthropogenic source monitoring. Sensitivity to signals from nontraditional sources, such as fan operations, highlights the broader applicability of this approach to other industrial environments in which continuous and impulsive signals are present. The field campaigns produced a substantial volume of high-quality seismic data, supporting the development and testing of new methods for seismic source characterization and subsurface imaging. Future deployments should include sensors capable of recording lower frequencies to probe deeper structures, increase bandwidth to enhance resolution and sensitivity to both shallow and deep targets, and collect additional large-scale datasets to refine imaging and source characterization techniques. Moreover, conducting 3D modeling studies of seismic wavefields at higher frequencies will provide a better understanding of wave scattering and cavity–wavefield interactions in complex underground environments. In conclusion, the STILGAR project demonstrated that integrated seismic monitoring can effectively characterize underground operations, capturing both natural and anthropogenic signals. The approaches developed provide a foundation for improved detection, localization, and imaging of subsurface structures and are directly transferable to broader industrial monitoring applications.

58 GEOSCIENCES↗