Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “waveform 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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Integration of Waveform Simulation Methods

The generation of synthetic seismograms through simulation is a fundamental tool of seismology required to run quantitative hypothesis tests. A variety of approaches have been developed throughout the seismological community and each has their own specific user interface based on their implementation. This causes a challenge to researchers who will need to learn new interfaces with each new software they wish to use and create substantial challenges when attempting to compare results from different tools. Here we provide a unified interface that facilitates interoperability amongst several simulation tools through a modern containerized Python package. Further, this package includes post-processing analysis modules designed to facilitate end-to-end analysis of synthetic seismograms. In this report we present the conceptual guidance and an example implementation of the new Waveform Simulation Framework.

58 GEOSCIENCES↗
DOE OSTI2022

Yield Estimation of the August 2020 Beirut Explosion by Using Physics‐Based Propagation Simulations of Regional Infrasound

Abstract The August 2020 Beirut explosion is the largest single‐fired ammonium nitrate explosion documented in history. The massive explosion excited loud infrasound in the atmosphere, and clear waveforms were recorded by a regional infrasound array at an epicentral distance of 100 km, allowing for accurate measurements of explosion energy. We estimate the explosion size based on the infrasound waveform inversion. Unlike conventional inversions using empirical models, we perform full 3‐D finite‐difference simulations to obtain a physics‐based propagation model for the inversion. Accurate numerical modeling of infrasound is challenging as the propagation is substantially affected by the turbulent atmosphere. Instead of a single deterministic prediction, we provide a range of waveform predictions by running multiple simulations with stochastic weather forecast models, which allows for comprehensive uncertainty analysis of numerical modeling and estimated yields. Finally, we expand the yield estimation technique for seismoacoustic analysis and demonstrate the substantial advantage of the joint approach.

58 GEOSCIENCES↗
DOE OSTI2025

Catalytic resonance theory for parametric uncertainty of programmable catalysis

Microkinetic models are useful tools for screening catalytic materials; however, errors in their input parameters can lead to significant uncertainty in model predictions of catalyst performance. Here, in this work, we investigate the impact of linear scaling and Brønsted-Evans-Polanyi relation parametric uncertainty on microkinetic predictions of programmable-catalyst performance. Two case studies are considered: a generic A-to-B prototype reaction and the oxygen evolution reaction (OER). The results show that error-unaware models can accurately predict trends and, for the prototype reaction, values of optimal waveform parameters. The specific model parameters driving output uncertainty are identified via variance-based global sensitivity analysis. However, predictions of dynamic rate enhancement can decrease when uncertainty is propagated into the models. In both cases, we identify operating conditions where the programmable catalyst achieves a rate enhancement of at least one order of magnitude despite parametric uncertainty in the model, supporting programmable catalysis as a viable strategy for exceeding the Sabatier limit.

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

Uncovering extreme nonlinear dynamics in solids through time-domain field analysis

Time-domain analysis of harmonic fields with sub-cycle resolution is now experimentally viable due to the emergence of sensitive, on-chip techniques for petahertz-scale optical-field sampling. We demonstrate how such a time-domain, field-resolved analysis uncovers the extreme nonlinear electron dynamics responsible for high-harmonic generation within solids. Time-dependent density functional theory was used to simulate harmonic generation from a solid-state band-gap system driven by near- to mid-infrared waveforms. Particular attention was paid to regimes where both intraband and interband emission mechanisms play a critical role in shaping the nonlinear response. We show that a time-domain analysis of the harmonic radiation fields identifies the interplay between intra- and interband dynamical processes underlying the nonlinear light generation. With further analysis, we show that changes to the dominant emission regime can occur after only slight changes to the peak driving intensity and central driving wavelength. Furthermore, time-domain analysis of harmonic fields also reveals, for the first time, the possibility of rapid changes in the dominant emission mechanism within the temporal window of the driving pulse envelope. Finally, we examine the experimental viability of performing time-domain analysis of harmonic fields with sub-cycle resolution using realistic parameters.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Unsupervised Clustering of Microseismic Events and Focal Mechanism Analysis at the CO 2 Injection Site in Decatur, Illinois

Characterization of induced microseismicity at a carbon dioxide (CO 2 ) storage site is critical for preserving reservoir integrity and mitigating seismic hazards. We apply a multilevel machine learning (ML) approach that combines the nonnegative matrix factorization and hidden Markov model to extract spectral representations of microseismic events and cluster them to identify seismic patterns at the Illinois Basin-Decatur Project. Unlike traditional waveform correlation methods, this approach leverages spectral characteristics of first arrivals to improve event classification and detect previously undetected planes of weakness. By integrating ML-based clustering with focal mechanism analysis, we resolve small-scale fault structures that are below the detection limits of conventional seismic imaging. Our findings reveal temporal bursts of microseismicity associated with brittle failure, providing insights into the spatio-temporal evolution of fault reactivation during CO 2 injection. This approach enhances seismic monitoring capabilities at CO 2 injection sites by improving fault characterization beyond the resolution of standard geophysical surveys.

Willis, Rachel Marie [Sandia National Laboratories↗
DOE OSTI2026

Waveform Simulation Framework: User Manual with Tutorials

This manuscript describes the Waveform Simulation Framework (WSF), a Python-based framework that provides a unified, programmable interface for generating synthetic seismograms for applications such as seismic array design, method development, and special event analysis. WSF standardizes how users define sources, receivers, and velocity models while abstracting simulator-specific configuration details, enabling workflows that are largely independent of the underlying numerical engine. The document provides installation guidance and tutorial-driven examples for three WSF simulator wrappers—WSF PyFK, WSF SW4, and WSF SPECFEM2D—illustrating end-to-end workflows from forward waveform simulation to common post-processing tasks (e.g., visualization and backprojection) using consistent data products (e.g., ObsPy Stream objects and SAC files).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Elastic Stochastic Full Waveform Inversion (eSFWI)

This collaboration between Lawrence Livermore National Security, LLC (LLNS) as manager and operator of Lawrence Livermore National Laboratory (LLNL) and Chevron USA Inc., acting through its Chevron Technical Center division, aimed at developing next-generation computational methods for the Elastic Stochastic Full Waveform Inversion (eSFWI). Seismic imaging is heavily used in the oil and gas industry for identifying and operating subsurface reservoirs. Improved seismic imaging methods can improve productivity, lower costs, and improve operational and environmental safety. This CRADA demonstrated that new high-performance computing (HPC) architectures being rolled out over the next five years can enable unprecedented seismic imaging resolution when using eSFWI techniques to process active seismic data. An open-source computational mini-application was developed, capable of demonstrating near-peak performance for eSFWI algorithms on CPU and GPU enabled HPC platforms. Performance was demonstrated on LLNL HPC systems such as Lassen, as well as on Chevron systems. This project benefited Chevron USA Inc. by demonstrating the potential computational efficiency of their full waveform inversion capabilities used to characterize oil/gas reservoirs, which in turn benefits the public through potential increases in capabilities to perform analysis of leasing sites.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2024

Comparative Analysis of Explosion Signals on Smartphones and Seismometers [Poster]

Seismoacoustic signatures from a chemical explosion was recorded on two seismometers and eight smartphone accelerometers at Idaho National Laboratory (INL). The captured waveform’s propagation and its time frequency representations are compared. The smartphone accelerometers were able to pick up a clear seismoacoustic explosion signal at within 1 km (150 m/kg 1/3 ) and a noisy signal up to 5.3 km (790 m/kg 1/3 ).

98 - NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL↗
DOE OSTI2025

Data from: 'Abiotic influences on continuous conifer forest structure across a subalpine watershed'

This package archives the core data used for analysis and inference in 'Abiotic influences on continuous conifer forest structure across a subalpine watershed' (Worsham et al., 2025). All data were collected in the East River, Washington Gulch, Slate River, and Coal Creek watersheds of Colorado. In the paper, we quantified the relative influence of climate, topographic, edaphic, and geologic factors on conifer stand structure and composition, and their functional relationships, at the watershed scale. We used waveform LiDAR data to derive spatially continuous stand structure metrics. We fused these with a species-level classification map to estimate tree species abundance. We applied generalized additive and generalized boosted models to evaluate the covariability of structural and compositional metrics with abiotic variables. The package contains the essential products required for reproducing our analysis and the tables and figures reported in the publication. The products comprise four classes: (1) geospatial data, (2) tabular data used for inferential analysis, (3) tabular data describing analytical results and performance statistics, and (4) a data user guide. (1) includes discretized waveform LiDAR data, locations and attributes of individual tree crowns, sampling locations and domain boundaries, a canopy height model, and raster files of estimated forest structural and compositional metrics at 100 m grid scale. (2) includes all response and explanatory variable values applied in inferential models. Response variables include conifer forest stand density, basal area, 95th percentile height, quadratic mean diameter, and others. Explanatory variables include climatic water deficit, actual evapotranspiration, elevation, heat load, soil available water content, and others. (3) includes results of training and testing several individual tree detection (ITD) algorithms, as well as inferential modeling results. (4) is a PDF user guide for this data package, including detailed descriptions and data dictionaries for all files. The data package root contains 17 assets: 8 compressed tape archive (.tar.gz) files, 5 comma-separated values (.csv) files, 3 Geographic Tagged Image File Format (GeoTIFF) (.tif) files, and 1 Portable Document Format (.pdf) file. The compressed .tar.gz archives contain ESRI shapefiles (.shp) .tif, compressed LASer (.laz), and .csv files. The archives must first be decompressed using the widely distributed command-line software utility TAR. All other files, including constituent files within the .tar.gz archives, can be opened in the open-source R statistical computing environment. Alternatively, .csv files may also be read in any simple text editor software or Microsoft Excel. Geospatial files including .shp and .tif files can also be opened in GIS software, such as QGIS (open-source) or ESRI ArcGIS (proprietary). The .pdf Data User Guide can be read with Adobe Acrobat Reader or other compatible readers.

2018 NEON and 2025 CHESS Campaigns↗
DOE OSTI2023

Acoustic signal based analysis of batteries

Systems and methods for acoustic signal based analysis, include obtaining acoustic response signal data of at least a portion of a battery cell, the acoustic response signal data comprising waveforms generated by transmitting one or more acoustic excitation signals into at least the portion of the battery cell and recording response vibration signals to the one or more acoustic excitation signals. One or more metrics are determined from at least the acoustic response signal data, the one or more metrics being determined based on correlation of the one or more metrics to one or more characteristics of battery cells and a reference model is generated from the one or more metrics. A test battery can be evaluated using the reference model. Actionable insights or recommendations can be generated based on the evaluation. The reference model can also be updated based on the evaluation.

Van Tassell, Barry J.↗
DOE OSTI2024

Machine Learning Decoding of Full Duplex Signals

A full duplex signal is when two endpoints (server one and server two) transmit on a single conductor pair simultaneously and with the same frequency. This results in the waveforms created from each server to be merged with one another when observed at any point along the transmission line making physical analysis of the wave unobtainable. This project was orchestrated to find the means to separate the merged signal into two separate signals which represent the signals originally sent from each server without an active tap.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Developing new analysis tools for near surface radio-based neutrino detectors

The ARIANNA experiment is an Askaryan radio detector designed to measure high-energy neutrino induced cascades within the Antarctic ice. Ultra-high-energy neutrinos above 10 16 eV have an extremely low flux, so experimental data captured at trigger level need to be classified correctly to retain as much neutrino signal as possible. We first describe two new physics-based neutrino selection methods, or "cuts", (the updown and dipole cut) that extend the previously published analysis to a specialized ARIANNA station with 8 antenna channels, which is double the number used in the prior analysis. For a standard trigger with a threshold signal to noise ratio at 4.4, the new cuts produce a neutrino efficiency of > 95% per station-year of operation, while rejecting 99.93% of the background (corresponding to 53 remaining experimental background events). When the new cuts are combined with a previously developed cut using neutrino waveform templates, all background is removed at no change of efficiency. In addition, the neutrino efficiency is extrapolated to 1,000 station-years of operation, obtaining 91%. This work then introduces a new selection method (the deep learning cut) to augment the identification of neutrino events by using deep learning methods and compares the efficiency to the physics-based analysis. The deep learning cut gives 99% signal efficiency per station-year of operation while rejecting 99.997% of the background (corresponding to 2 remaining experimental background events), which are subsequently removed by the waveform template cut at no significant change in efficiency. Here the results of the deep learning cut were verified using measured cosmic rays which shows that the simulations do not introduce artifacts with respect to experimental data. The paper demonstrates that the background rejection and signal efficiency of near surface antennas meets the requirements of a large scale future array, as considered in baseline design of the radio component of IceCube-Gen2.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

User Guide: A Curated Dataset of Regional Meteor Events with Simultaneous Optical and Infrasound Observations

This user guide supports a curated dataset of 71 meteor events recorded between 2006 and 2011 in Southwestern Ontario, Canada. Each event was simultaneously observed by ground-based optical cameras and an infrasound array, providing a rare opportunity to examine meteor trajectories and acoustic signals from the same atmospheric entry events. The dataset includes raw and processed optical data, meteor trajectories, photometric light curves, infrasound waveforms, and atmospheric specifications relevant for acoustic modeling. The archive is structured to support reproducible research in meteor physics, atmospheric acoustics, and shock wave analysis. It is organized following transparent file naming conventions and structured folders to facilitate scientific reuse, comparison, and integration across research domains. The dataset is freely available on Zenodo, doi: 10.5281/zenodo.15868512.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Time-Frequency Warped Waveforms for Well-Contained Massive Machine Type Communications

Herein this paper proposes a novel time-frequency warped waveform for short symbols, massive machine-type communication (mMTC), and internet of things (IoT) applications. The waveform is composed of asymmetric raised cosine (RC) pulses to increase the signal containment in time and frequency domains. The waveform has low power tails in the time domain, hence better performance in the presence of delay spread and time offsets. The time-axis warping unitary transform is applied to control the waveform occupancy in time-frequency space and to compensate for the usage of high roll-off factor pulses at the symbol edges. The paper explains a step-by-step analysis for determining the roll-off factors profile and the warping functions. Gains are presented over the conventional Zero-tail Discrete Fourier Transform-spread-Orthogonal Frequency Division Multiplexing (ZT-DFT-s-OFDM), and Cyclic prefix (CP) DFT-s-OFDM schemes in the simulations section.

42 ENGINEERING↗
DOE OSTI2024

A comparison of smartphone and infrasound microphone data from a fuel air explosive and a high explosive

For prompt detection of large (>1 kt) above-ground explosions, infrasound microphone networks and arrays are deployed at surveyed locations across the world. Denser regional and local networks are deployed for smaller explosions, however, they are limited in number and are often deployed temporarily for experiments. With the expanded interest in smaller yield explosions targeted at vulnerable areas such as population centers and key infrastructures, the need for more dense microphone networks has increased. An “attritable” (affordable, reusable, and replaceable) and flexible alternative can be provided by smartphone networks. Explosion signals from a fuel air explosive (thermobaric bomb) and a high explosive with trinitrotoluene equivalent yields of 6.35 and 3.63 kg, respectively, were captured on both an infrasound microphone and a network of smartphones. The resulting waveforms were compared in time, frequency, and time-frequency domains. The acoustic waveforms collected on smartphones produced a filtered explosion pulse due to the smartphone's diminishing frequency response at infrasound frequencies (<20 Hz) and was found difficult to be used with explosion characterization methods utilizing waveform features (peak overpressure, impulse, etc.). However, the similarities in time frequency representations and additional sensor inputs are promising for other explosion signal identification and analysis. As an example, a method utilizing the relative acoustic amplitudes for source localization using the smartphone sensor network is presented.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

A Comprehensive Analysis of Real-World Accelerometer Data Quality in a Global Smartphone-based Seismic Network

The proliferation of low-cost sensors in smartphones has facilitated numerous applications; however, large-scale deployments often encounter performance issues. Sensing heterogeneity, which refers to varying data quality due to factors such as device differences and user behaviors, presents a significant challenge. In this research, we perform an extensive analysis of 3-axis accelerometer data from the MyShake system, a global seismic network utilizing smartphones. We systematically evaluate the quality of approximately 22 million 3-axis acceleration waveforms from over 81 thousand smartphone devices worldwide, using metrics that represent sampling rate and noise level. We explore a broad range of factors influencing accelerometer data quality, including smartphone and accelerometer manufacturers, phone specifications (release year, RAM, battery), geolocation, and time. Our findings indicate that multiple factors affect data quality, with accelerometer model and smartphone specifications being the most critical. In addition, we examine the influence of data quality on earthquake parameter estimation and show that removing low-quality accelerometer data enhances the accuracy of earthquake magnitude estimation.

58 GEOSCIENCES↗
DOE OSTI2023

Electrorefining Bismuth Using AC Superimposed DC Waveforms

Bismuth has been successfully electrorefined using both direct current (DC) and alternating current (AC) superimposed DC waveforms on a kilogram scale at Brigham Young University (BYU) in Provo, Utah. Finding a suitable surrogate for Pu Electrorefining has always been a technical challenge. This work allows the possibility to improve Pu electrorefining with nonradioactive material at universities, which greatly reduces the cost at Lawrence Livermore National Laboratory (LLNL) for improving Pu Electrorefining. Several possible surrogates (Ce, In, Sn, Zn, and Bi) were selected based on their ability to mimic the Pu electrorefining process. Ce and In electrorefining experiments were conducted at the same temperature as Pu electrorefining in CaCl 2 while Sn, Zn, and Bi, electrorefining experiments were conducted at a lower temperature in a eutectic molten salt composed of LiCl, KCl, and CaCl 2 . The electrorefining experiments using Ce with a Ga impurity did not produce a cathode ring. While the remaining Ce anode after electrorefining was significantly less than what was initially added, Scanning Electron Microscopy with Energy Dispersive X-Ray (SEM-EDX) analysis showed that the Ce had been oxidized and formed a colloid with the surrounding molten salt. This is supported in the literature by analysis of Ce-rich CaCl 2 mixtures. All the electrorefining experiments except for the Ce electrorefining experiments did not introduce an impurity to the system. Electrorefining experiments with In showed that the InCl 3 volatilized out of the molten salt to an extent that only a low current (<<1 A) could be supported without decomposing the molten salt. Zn electrorefining experiments had marginal success, however an easily separable product was not formed. The Zn had to be rinsed from the salt after the experiment. Of the identified potential surrogates. Investigated, only Bi and Sn yielded an easily separable product ring, which enables quantitative analysis of yields and coulombic efficiencies. Bismuth was first identified as a surrogate at BYU, while Sn was identified as a surrogate at LLNL later. Thus, the electrorefining experiments conducted at BYU using AC superimposed DC waveforms used Bi as the metal to be electrorefined.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Joint Inversion of Regional Waveform, First-Motion Polarity, and Synthetic Aperture Radar Surface Displacement for the Fourth and Sixth North Korean Declared Nuclear Explosions

Here, this study analyzed the Democratic People’s Republic of Korea’s (DPRK) fourth (DPRK4, 6 January 2016 M w 4.49) and sixth (DPRK6, 7 September 2017 M w 5.2) declared nuclear tests, employing a joint seismic and Interferometric Synthetic Aperture Radar (InSAR) inversion to improve understanding of these events and enhance moment tensor (MT) inversion capabilities. The recent efforts have focused on employing seismic waveform and InSAR geodetic deformation data separately to analyze these and the previous nuclear tests (e.g., Chiang et al., 2018; Myers et al., 2018; Wang et al., 2018). Building upon our previous work (Chi-Durán et al., 2021), we performed a joint regional waveform, first-motion (FM) polarity, and surface displacement inversion, which demonstrated improved source-type discrimination, a revised MT solution with reduced scalar moment uncertainty, and an independently constrained location. In this article, we build on the previous results for DPRK6 by including an analysis using a four-layered velocity model with free-surface topography to compute the near-source static deformation Green’s functions. The model consists of a 50 m basalt layer (⁠V P = 2.07 km/s, V S = 1.2 km/s⁠), a 250 m stratified volcanic deposit layer (⁠V P = 1.73 km/s, V S = 1.0 km/s⁠), a 700 m weathered granodiorite layer (⁠V P = 2.5 km/s, V S = 1.3 km/s⁠⁠), and a granodiorite half-space (⁠V P = 5.35 km/s, V S = 3.09 km/s⁠⁠). The half-space shares the velocity of the regional MDJ2 velocity model (Ford et al., 2010), which has proven effective for waveform inversion in the region. This model considers the range of reported values for various lithologies and weathering effects. Our findings show that using the layered velocity model enhances the recovery of source location and depth for both the explosions by improving fits and reducing uncertainties. The joint inversion also improves source-type discrimination and better constrains the scalar seismic moment necessary for downstream yield estimation.

58 GEOSCIENCES↗