Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “acoustical waves”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Monitoring of corrosion effects in pipes with multi-mode acoustic signals

Current commercial long-range guided wave inspection and monitoring technologies for pipes usually require establishment of certain specific acoustic modes for wave propagation and use a large number of sensors in a certain geometric arrangement (for example, uniformly-spaced arrays of sensors on a collar, placed around a pipe). This paper, in contrast, presents a long-range corrosion detection method in pipes based on multi-mode guided acoustic waves, generated and received using simple, single-element sensors. In this technique, available guided wave modes in any interrogated zone of the structure are automatically down-selected by the medium, and changes in energy distribution over the selected modes are analyzed through signal processing to monitor material loss or addition. Here, to investigate the feasibility of the proposed method, we have performed both experimental and numerical simulation studies on 50.8-mm (2-in)-diameter schedule-40 carbon steel pipes. In these studies, a linear chirp signal is injected in compression mode, perpendicular to the surface of the pipe, generating multi-mode, multi-directional guided waves that travel within the pipe wall. Corrosion effects in the pipes are simulated numerically and experimentally. To reduce errors in the measurement technique due to changes in the local or global temperature of the measured media, a universal temperature compensation algorithm is introduced. The experimental results show good agreement with simulation results, indicating that the proposed multi-mode wave generation and detection in a pipe could be used for long-range detection and monitoring of corrosion effects, such as pipe wall losses or solid material deposits.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Determination of acoustic nonlinearity parameters using thermal modulation of ultrasonic waves

This study presents a test method and its theoretical framework to determine the acoustic nonlinearity parameters (α, β, δ) of material using thermal modulation of ultrasonic waves. Temperature change-induced thermal strain excites the nonlinear response of the material and modulates the ultrasonic wave propagating in it. Experimental results showed a strong correlation between the relative wave velocity change and the temperature change. With a quadratic polynomial model, the acoustic nonlinearity parameters were obtained from the polynomial coefficients by curve fitting the experimental data. Their effects on thermal-induced velocity change were discussed. The parameters α, β, and δ govern the hysteretic gap, average slope, and curvature of the correlation curve, respectively. The proposed theory was validated on aluminum, steel, intact and damaged concrete samples. The obtained nonlinear parameters show reasonable agreement with values reported in the literature. Compared to other nonlinear acoustic methods using vibration or acoustic excitation, the thermal modulation method generates more uniform, slow changing, and larger strain field in the test sample. Employing the thermal effect as the driving force for nonlinearity instead of an undesired influencing factor, this method can measure the absolute values of α, β, and δ with good accuracy using a simple ultrasonic test setup.

Physics↗
DOE OSTI2021

A Self-Synchronizing Underwater Acoustic Network for Mooring Load Monitoring of a Wave Energy Converter

This paper reports on the development of a self-synchronizing underwater acoustic network developed for remote monitoring of mooring loads in Wave Energy Converters (WECs). This network uses Time Division Multiple Access and operates self-contained with the ability for users to remotely transmit commands to the network as needed. Each node is a self-contained unit, consisting of a protocol adaptor board, an FAU-DPAM underwater acoustic modem and a battery pack. A node can be connected to a load cell, to a topside user or to the WEC. Every node is swapable. The protocol adaptor board, named Protocol Adaptor for Digital LOad Cell (PADLOC) supports a variety of digital load cell message formats (CAN, MODBUS, custom ASCII) and underwater acoustic modem serial formats. PADLOC enables topside users to connect to separate load cells through a user-specific command. This is especially important if the user is monitoring multiple load cells during deployment or maintenance, when the primary data system may be offline. Each PADLOC board handles formatting, buffering and has a one-on-one serial connection with each pair (node) of a digital load cell and acoustic modem. In addition, each PADLOC board handles the timekeeping and power saving features for each node. The only limitation is the data bit rate and delay limitations associated with the underwater acoustic modem. A four node self-synchronizing network has been developed to demonstrate the load cell monitoring capability using the PADLOC technology on the CalWave WEC.

acoustic↗
DOE OSTI2022

Feature Extraction for Pipeline Defects Inspection Based Upon Distributed Acoustic Fiber Optic Sensing Data

Fiber-optic distributed acoustic sensing (DAS) is becoming an increasingly important tool for real-time monitoring of energy and civil infrastructure structural health such as pipelines. We present a systematic theoretical study of the potential for DAS to be directly coupled with guided ultrasonic waves typically used in conventional acoustic non-destructive evaluation (NDE) methods for real-time pipeline health monitoring. We are referring to this innovative new NDE technique as ultrasonic guided wave and optical fiber sensor fusion. In the practical application of DAS coupled with guided ultrasonic waves, the structural design of (1) the specific guided waves excited, (2) the physical installation of the acoustic transducers and the fiber optic sensors, and (3) the functional performance specifications (gauge length, sensitivity, Etc.) of fiber optic DAS have an important influence on overall capabilities of the monitoring system. Meanwhile, physics-based analysis of acoustic waves is still a challenge due to the complex nature of the Lamb wave when it propagates, scatters, and disperses in the presence of structural defects. In this work, we simulate carbon steel pipes relevant for oil and gas pipeline applications with diameters of approximately 6-12” and wall thickness of 0.5” as the objects to be monitored. By establishing and implementing these capabilities, we seek to pursue an in-depth study on structural parameter optimization of DAS network, measurement range, and signal processing with an ultimate goal of increasing the sensitivity and efficacy of DAS to defect identification for various modes of corrosion expected in practice. To study the characteristics of scattered acoustic waves and performance of DAS for defect identification, we simulated the response of DAS for multiple pipe structures, defect types, and DAS sensor network configuration using finite element software Ansys, then the properties of signal response are extracted to construct defect-sensitive features. The raw data simulated, and the associated features extracted can ultimately be utilized as annotated training data to benchmark various designs for DAS applications, guided acoustic excitation sources, and learning model parameters to enhance early detection of potentially problematic defects.

Pipeline Defects Inspection, Fiber-optic sensors, ↗
DOE OSTI2020

Validation of finite element analysis strategy to investigate acoustic levitation in a two-axis acoustic levitator

A two-axis acoustic levitator can be used to generate a standing pressure wave capable of levitating solid and liquid particles at appropriate input conditions. This work proposes a simulation framework to investigate the two-axis levitation particle stability using a commercial, computational fluid dynamics software based on the harmonic solution to the acoustic wave equation. The simulation produced predictions of the standing wave that include a strong "+" shaped pattern of nodes and anti-nodes that are aligned with the levitator axes. To verify the simulation, a levitator was built and used to generate the standing wave. The field was probed with a microphone and a motorized-scanning system. After scaling the simulated pressure to the measured pressure, the magnitudes of the sound pressure level at corresponding high-pressure locations were different by no more than 5%. This is the first time a measurement of a two-axis levitator standing pressure wave has been presented and shown to verify simulations. Furthermore, as an additional verification, the authors consulted high speed camera measurements of a reference-levitator transducer, which was found to have a maximum peak-to-peak displacement of 50 +/- 5 mu m. The reference-levitator is known to levitate water at 160 dB. The system for this work was simulated to match the operation of the reference-levitator so that it produced sound pressure levels of 160 dB. This pressure was achieved when the transducer maximum peak-to-peak displacement was 50.8 mu m. The agreement between the two levitators' displacements provides good justification that the modeling approach presented here produces reliable results.

2 axis levitator↗
DOE OSTI2020

Infrasound and Gravity Waves Over the Andes Observed by a Pressure Sensor on Board a Stratospheric Balloon

The study of infrasound (acoustic) and gravity waves sources and propagation in the atmosphere of a planet gives us precious insight on atmosphere dynamics, climate, and even internal structure. The implementation of modern pressure sensors with high rate sampling on stratospheric balloons is improving their study. In our work we analyzed the data from the National Aeronautics and Space Administration Ultra Long Duration Balloon mission (16 May to 30 June 2016). Here, we focus on the balloon's transit of the Andes Mountains. We detected gravity waves that are associated to troposphere convective activity and mountain waves. An increase of the horizontal wavelengths from 50 to 70 km with increasing distance to the mountains is favoring the presence of mountain waves. We also report on the detection of infrasounds generated by the mountains in the 0.01–0.1 Hz range with a pressure amplitude increase by a factor 2 relative background signal. Besides, we characterized the decrease of microbaroms power when the balloon was flying away from the ocean coast. These observations suggest, in a way similar to microseisms for seismometers, that microbaroms are the main background noise sources recorded in the stratosphere even far from the ocean sources. Finally, we observed a broadband signal above the Andes, between 0.45 and 2 Hz, probably associated with a thunderstorm. The diversity of geophysical phenomena captured in less than a day of observation stresses the interest of high rate pressure sensors on board long-duration balloon missions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Acoustic Codes in 2D Spherical Coordinate

Finite-difference methods are widely used to simulate infrasound propagation in the atmosphere. Flexibility of finite-difference scheme allows implementation of highly heterogeneous media for sound propagation as well as complex source models for sound generation. While full 3-D finite-difference methods have been utilized for local infrasound propagation with pronounced topography, 2-D modeling approach has been preferred for regional and global propagation as full 3-D methods generally require enormous computational resources. Infrasound propagation is often simulated with a second-order finite difference scheme. This lowest-order finite-difference scheme is robust and straightforward to implement complex boundary conditions, but the solution includes large error with numerical dispersion and dissipation. This large numerical error may make the second-order finite-difference not optimal for long range infrasound propagation modeling as the numerical dispersion degrades the accuracy of the solution unacceptably. Here, we developed a high-order finite-difference solver for long-range infrasound simulation. The high-order scheme is particularly popular for linear wave modeling in aeroacoustics owing to its low-dispersive and low-dissipative behavior. We develop and evaluate a high-order finite difference scheme in 2-D axisymmetric coordinates. The axisymmetry allows to approximate 3-D spherical sound propagation and amplitude attenuation by a 2-D method. AC2Dr is developed to simulate infrasound propagation in realistic atmosphere, but can be used for linear acoustic waves in general materials with background flow. AC2Dr in an axisymmetric coordinates allows for spherical radiation of acoustic waves from compact sources.

Sjogreen, Bjorn↗
DOE OSTI2020

Significance of the Direct Excitation Mechanism for High-Frequency Response of Premixed Flames to Flow Oscillations

Premixed flames are sensitive to flow disturbances, which can arise from acoustic or vortical fluctuations. For transverse instabilities, it is known that a dominant mechanism for flame response is “injector coupling,” whereby pressure oscillations associated with transverse waves excite axial flow disturbances. These axial flow disturbances then excite heat release oscillations. The objective of this paper is to consider another mechanism—the direct sensitivity of the unsteady heat release to transverse acoustic waves—and to compare its significance relative to the induced axial disturbances, in a linear framework. The rate at which the flame adds energy to the disturbance field is quantified using the Rayleigh criterion and evaluated over a range of control parameters, such as flame length and swirl number. The results show that radial modes induce heat release fluctuations that always add energy to the acoustic field, whereas heat release fluctuations induced by mixed radial-azimuthal modes can add or remove energy. These amplification rates are then compared to the flame response from induced axial fluctuations. For combustor-centered flames, these results show that the direct excitation mechanism has negligible amplification rates relative to the induced axial mechanism for radial modes. For transverse modes, the fact that the nozzle is located at a pressure node indicates that negligible induced axial velocity disturbances are excited; as such, the direct mechanism dominates. As a result, for flames that are not centered on pressure nodes, the direct mechanism for mixed modes dominates for certain nozzle locations and flame angles.

42 ENGINEERING↗
DOE OSTI2023

Chemical Blast Standard (1 kg)

Chemical explosions create blast waves with large overpressure disturbances. It is important to develop a standard blast model based on data to accurately predict acoustic blast-wave amplitudes near detonations and invert for explosion energy from distant observations of blast-wave signals. However, open data from large, controlled chemical explosions with reliable ground truth can be challenging to find. The lack of access to such data could limit the number of contributions to related research and potentially stifle the rate of discoveries or validation of existing models. Here, to address these data scarcity problem, we have curated and compiled a standardized set of 817 blast-wave waveforms from 19 distinct high-explosive events. The blast-wave waveforms are standardized to a 1 kg trinitrotoluene explosion using scaling laws and corrections for location effects. A brief overview of the dataset is presented along with explosion feature models as well as recommendations for extracting explosion features. The resulting dataset is distributed to an open repository in both Seismic Analysis Code and pandas DataFrame formats containing the waveforms, the scaled distances, and the sample rates.

58 GEOSCIENCES↗
DOE OSTI2021

Characterization of damage processes in Montney siltstone under triaxial compression using acoustic emission and diagnostic imaging

SUMMARY Crack nucleation and rock failure processes in a fine-grained siltstone (Montney Formation) under triaxial compression are investigated using combined diagnostic techniques, including ultrasonic-wave measurement, acoustic-emission (AE) monitoring, computed tomography (CT) scanning, and thin-section imaging. The sample displays a weak-to-moderate inherent seismic anisotropy and noticeable stress-induced anisotropy prior to failure. No AE event was detected until the applied axial stress reached 95 per cent of the peak value. The signal-to-noise ratio is relatively low, however, and detectable AE events are more diffuse than those observed in highly brittle rocks. The AE locations correlate with a shear fracture zone imaged by CT scanning. AE moment-tensor analysis reveals that events with larger relative magnitudes are characterized by high volumetric (tensile or compressive) components, and the initiation of the failure zone is dominated by combined shear-tensile (ST) failure. Stress inversion of the AE events with high tensile components is in good agreement with the known applied stress. Microscopic imaging of thin sections from the failed sample shows that the failure zone is an en echelon structure consisting of a major fracture with branching micro and minor cracks. This failure mechanism is consistent with a ST source mechanism and is interpreted to be associated with the fine granular structure and mineral composition of Montney siltstone.

Jia, Suzie Q. (ORCID:0000000326496190)↗
DOE OSTI2021

Failure prediction and estimation of failure parameters

Machine-learning methods and apparatus are disclosed to determine frictional state or other parameters in an earthquake zone or other failing medium, using acoustic emission, seismic waves, or other detectable indicators of microscopic processes. Predictions of future failures are demonstrated in different regimes. A classifier is trained using time series of acoustic emission data along with historic data of frictional state or failure events. In disclosed examples, random forests and gradient boost trees are used, and grid-search or EGO procedures are used for hyperparameter tuning. Once trained, the classifier can be applied to testing or live data in order to assess a frictional state, assess seismic hazard, or make predictions regarding a future failure event. The technology has been developed in a double direct shear apparatus, but can be widely applied to seismic faults, other terrestrial failures, or failures in man-made structures. Variations are disclosed.

Johnson, Paul Allan↗
DOE OSTI2023

Seismoacoustic Signatures Observed During a Long-Term Deployment of Infrasound Sensors at the Nevada National Security Site

Earthquakes have repeatedly been shown to produce inaudible acoustic signals (<20 Hz), otherwise known as infrasound. These signals can propagate hundreds to thousands of kilometers and still be detected by ground-based infrasound arrays depending on the source strength, distance between source and receiver, and atmospheric conditions. Another type of signal arrival at infrasound arrays is the seismic induced motion of the sensor itself, or ground-motion-induced sensor noise. Measured acoustic and seismic waves produced by earthquakes can provide insight into properties of the earthquake such as magnitude, depth, and focal mechanism, as well as information about the local lithology and atmospheric conditions. Large earthquakes that produce strong acoustic signals detected at distances greater than 100 km are the most commonly studied; however, more recent studies have found that smaller magnitude earthquakes ($M$ w <2.0) can be detected at short ranges. In that vein, this study will investigate the ability for a long-term deployment of infrasound sensors (deployed as part of the Source Physics Experiments [SPE] from 2014 to 2020) to detect both seismic and infrasonic signals from earthquakes at local ranges (<50 km). Methods used include a combination of spectral analysis and automated array processing, supported by U.S. Geological Survey earthquake bulletins. Further, this investigation revealed no clear acoustic detections for short range earthquakes. However, secondary infrasound from an $M$ w 7.1 earthquake over 200 km away was detected. Important insights were also made regarding the performance of the SPE networks including detections of other acoustic sources such as bolides and rocket launches. Finally, evaluation of the infrasound arrays is performed to provide insight into optimal deployments for targeting earthquake infrasound.

58 GEOSCIENCES↗
DOE OSTI2024

DAISY Benchmark Performance Data

This repository contains the underlying data from benchmark experiments for Drifting Acoustic Instrumentation SYstems (DAISYs) in waves and currents described in "Performance of a Drifting Acoustic Instrumentation SYstem (DAISY) for Characterizing Radiated Noise from Marine Energy Converters" (https://link.springer.com/article/10.1007/s40722-024-00358-6). DAISYs consist of a surface expression connected to a hydrophone recording package by a tether. Both elements are instrumented to provide metadata (e.g., position, orientation, and depth). Information about how to build DAISYs is available at https://www.pmec.us/research-projects/daisy. The repository's primary content is three compressed archives (.zip format), each containing multiple MATLAB binary data files (.mat format). A table relating individual data files to figures in the paper, as well as the structure of each file, is included in the repository as a Word document (Data Description MHK-DR.docx). Most of the files contain time series information for a single DAISY deployment (file naming convention: [site]_DAISY_[Drift #].mat) consisting of processed hydrophone data and associated metadata. For a limited number of DAISY deployments, the hydrophone package was replaced with an acoustic Doppler velocimeter (file naming convention: [site]_DAISY_[Drift #]_ADV.mat). Data were collected over several years at three locations: (1) Sequim Bay at Pacific Northwest National Laboratory's Marine & Coastal Research Laboratory (MCRL) in Sequim, WA, the energetic tidal channel in Admiralty Inlet, WA (Admiralty Inlet), and the U.S. Navy's Wave Energy Test Site (WETS) in Kaneohe, HI. Brief descriptions of data files at each location follow. - MCRL - (1) Drift #4 and #16 contrast the performance of a DAISY and a reference hydrophone (icListen HF Reson), respectively, in the quiescent interior of Sequim Bay (September 2020). (2) Drift #152 and #153 are velocity measurements for a drifting acoustic Doppler velocimeter in in the tidally-energetic entrance channel inside a flow shield and exposed to the flow, respectively (January 2018). (3) Two non-standard files are also included: DAISY_data.mat corresponds to a subset of a DAISY drift over an Adaptable Monitoring Package (AMP) and AMP_data.mat corresponds to approximately co-temporal data for a stationary hydrophone on the AMP (February 2019). - Admiralty Inlet - (1) Drift #1-12 correspond to tests with flow shielded DAISYs, unshielded DAISYs, a reference hydrophone, and drifting acoustic Doppler velocimeter with 5, 10, and 15 m tether lengths between surface expression and hydrophone recording package (July 2022). (2) Drift #13-20 correspond to tests of flow shielded DAISYs with three different tether materials (rubber cord, nylon line, and faired nylon line) in lengths of 5, 10, and 15 m (July 2022). - WETS - (1) Drift #30-32 correspond to tests with a heave plate incorporated into the tether (standard configuration for wave sites), rubber cord only, and rubber cord, but with a flow shielded hydrophone (November 2022). (2) Drift #49-58 and Drift #65-68 correspond to measurements around mooring infrastructure at the 60 m berth where time-delay-of-arrival localization was demonstrated for different DAISY arrangements and hydrophone depths (November 2022).

16 TIDAL AND WAVE POWER↗
DOE OSTI2024

2-D seismic wave propagation using the distributional finite-difference method: further developments and potential for global seismology

SUMMARY We present a time-domain distributional finite-difference scheme based on the Lebedev staggered grid for the numerical simulation of wave propagation in acoustic and elastic media. The central aspect of the proposed method is the representation of the stresses and displacements with different sets of B-splines functions organized according to the staggered grid. The distributional finite-difference approach allows domain-decomposition, heterogeneity of the medium, curvilinear mesh, anisotropy, non-conformal interfaces, discontinuous grid and fluid–solid interfaces. Numerical examples show that the proposed scheme is suitable to model wave propagation through the Earth, where sharp interfaces separate large, relatively homogeneous layers. A few domains or elements are sufficient to represent the Earth’s internal structure without relying on advanced meshing techniques. We compare seismograms obtained with the proposed scheme and the spectral element method, and we show that our approach offers superior accuracy, reduced memory usage, and comparable efficiency.

Geochemistry & Geophysics↗
DOE OSTI2025

Revealing Defect-Seeded and Interfacial Generation Mechanisms of Photoinduced Coherent Phonons with 4D Ultrafast Electron Microscopy

Ultrafast photoexcitation of coherent phonons is driven by an impulsive, collective displacement of constituent atoms from their average equilibrium lattice positions [1]. Models describing the generation of coherent acoustic modes typically invoke the creation of an anisotropic strain profile arising from the relatively instantaneous absorption of an ultrafast laser pulse [2]. If the skin depth is shallow relative to the specimen thickness, a steep tensile strain gradient perpendicular to the surface ($\frac{∂ε}{∂z}$) results. Initial relaxation occurs via rapid contraction of the surface layers followed by subsequent coherent oscillations of the lattice and launch of a train of coherent elastic strain waves [i.e., coherent acoustic phonons (CAPs)]. Macroscopically, responses are generally well-described by constitutive relations as gleaned from data gathered using ultrasonic methods or ultrafast spectroscopies. Furthermore, at the atomic to nanoscale level, individual lattice discontinuities and their impact on CAP behaviors can be modeled using multiscale methods [3,4]. Further, average unit-cell level responses on ultrafast timescales can be probed using femtosecond electron and X-ray scattering [5,6].

Flannigan, David J. [University of Minnesota, Minn↗
DOE OSTI2024

Tailored ultrasound propagation in microscale metamaterials via inertia design

The quasi-static properties of micro-architected (meta)materials have been extensively studied over the past decade, but their dynamic responses, especially in acoustic metamaterials with engineered wave propagation behavior, represent a new frontier. However, challenges in miniaturizing and characterizing acoustic metamaterials in high-frequency (megahertz) regimes have hindered progress toward experimentally implementing ultrasonic-wave control. Here, we present an inertia design framework based on positioning microspheres to tune responses of 3D microscale metamaterials. We demonstrate tunable quasi-static stiffness by up to 75% and dynamic longitudinal-wave velocities by up to 25% while maintaining identical material density. Using noncontact laser-based dynamic experiments of tunable elastodynamic properties and numerical demonstrations of spatio-temporal ultrasound wave propagation, we explore the tunable static and elastodynamic property relation. This design framework expands the quasi-static and dynamic metamaterial property space through simple geometric changes, enabling facile design and fabrication of metamaterials for applications in medical ultrasound and analog computing.

42 ENGINEERING↗
DOE OSTI2025

Multi-frequency progressive refinement for learned inverse scattering

Interpreting scattered acoustic and electromagnetic wave patterns is a computational task that enables remote imaging in a number of important applications, including medical imaging, geophysical exploration, sonar and radar detection, and nondestructive testing of materials. However, accurately and stably recovering an inhomogeneous medium from far-field scattered wave measurements is a computationally difficult problem, due to the nonlinear and non-local nature of the forward scattering process. We design a neural network, called Multi-Frequency Inverse Scattering Network (MFISNet), and a training method to approximate the inverse map from far-field scattered wave measurements at multiple frequencies. We consider three variants of MFISNet, with the strongest performing variant inspired by the recursive linearization method — a commonly used technique for stably inverting scattered wavefield data — that progressively refines the estimate with higher frequency content. MFISNet outperforms past methods in regimes with high-contrast, heterogeneous large objects, and inhomogeneous unknown backgrounds.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Joint Bayesian Inference for Near-Surface Explosion Yield and Height-of-Burst

Forensic capabilities to understand chemical and nuclear explosions are greatly aided by an accurate estimate of explosive yield with uncertainty. The relationship between explosive size and geophysical observations of seismic, acoustic, and optical waves can be exploited to provide an estimate of yield. Any near-surface yield estimate is complicated by the surface interaction, so an estimate for the explosion height-of-burst is necessarily included in the relationship. Additionally, the relationship dictates a trade-off between estimates of yield and height-of-burst. Fortunately, the surface interaction for each type of observation is different, which breaks the trade-off, and the inclusion of height-of-burst with multiple data types improves yield estimation. We define simple parametric forward models to relate seismoäcoustoöptic observations from a data set of known explosive yields and height-of-bursts. The parameters of the models and a prediction for the yield and height-of-burst of a new event can then be estimated given new observations via Bayesian inference. We report posterior distribution estimates of the parametric models using a Markov chain Monte Carlo sampling technique. These models are then used to predict the yield and height-of-burst of SUGAR, a historical near-surface nuclear explosion, using its reported historical observations. The reported yield of 1.2 ktonne Trinitrotoluene (TNT)-equivalent (Department of Energy, 2015) is within the estimated posterior. Yield uncertainty can be estimated from the spread of the posterior, which is between 0.9 and 2.1 ktonne TNT-equivalent. The posterior for height-of-burst has a wider range between 10 m below and 8 m above ground that includes the true height-of-burst of 1 m.

58 GEOSCIENCES↗