Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time of arrival estimation”

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 91 records · Page 5

Seasonality and Albedo Dependence of Cloud Radiative Forcing in the Upper Colorado River Basin

Mountains create and enhance their own clouds, which both scatter and absorb shortwave radiation from the sun and absorb and re-emit land surface and atmospheric longwave radiation. However, the impacts of clouds on the surface radiation balance in high elevation snowy mountain terrain are poorly explored. In this study, we use data collected by the SAIL field campaign and partner organizations in the upper elevations (2,880 m.a.s.l) of the Upper Colorado River Basin (UCRB) over a 21-month period from September 2021 to June 2023 to estimate Cloud Radiative Forcing (CRF) in the shortwave, longwave, and the net effect. Longwave warming effects dominate during the winter when snow albedos are high (0.8–0.9) and the background atmospheric precipitable water vapor is low (<0.5 cm), yielding a maximum monthly average net CRF of +34.7 W·m -2 , meaning that clouds increase the net radiation relative to clear skies during this time period. The sign of net CRF switches in the warm season as snow recedes, sun-angles increase, and the North American monsoon arrives, yielding a minimum monthly average net CRF of -47.6 W·m -2 with hourly minima of -600 W·m -2 . The sign of net CRF is typically positive, even at solar noon, when the surface is snow covered, except for a brief period over melting, low-albedo snow (0.5–0.6) impacted by dust impurities. Sensitivity tests elucidate the role of the surface albedo on the net CRF. The results suggest that net CRF will increase in magnitude and lead to a more persistent cooling effect on the surface net radiation budget as the snow cover declines.

54 ENVIRONMENTAL SCIENCES↗

Data Analysis of the 2020 Central Idaho Mainshock-Aftershock Sequence

In an effort to inform the Senior Seismic Hazard Analysis Committee for the Idaho National Laboratory, we provide an improved aftershock catalog related to the March 31, 2020, Mw6.5 Stanley, Idaho earthquake from picks related to a temporary network of two real-time and 15 non-telemetered seismometers within the epicentral area. From the permanent and temporary (XP) real-time network, the USGS cataloged 1,946 aftershocks between April 1, 2020 and October 31, 2020. To improve aftershock location and magnitudes, we manually picked arrival times of P and S waves from off-line stations in the XP temporary network, generated a new crustal velocity model, and independently relocated each event using the HypoDD double-difference earthquake algorithm. We created our new velocity model from existing broadband and active source seismic campaign data that were acquired near the epicentral region prior to the 2020 earthquake. We compare arrival time differences, epicentral locations and depths between aftershocks recorded with the two catalogs. We find the addition of local stations provides tighter aftershock clustering that suggests an improved aftershock locations. To detect lower magnitude events, we employed deep learning. Our method solves common problems associated with detecting many events that have a low signal-to-noise ratio. From the machine learning database, we detected more than 74,000 aftershocks. Based on the number of identified earthquakes and Gutenberg-Richter relationships derived from the USGS catalog, we estimate that we have reduced the completion magnitude for the Stanley earthquake sequence to below M1 using this machine learning approach. We located each aftershock with our new velocity model. Our new velocity model and picks suggests aftershocks occurred mostly at shallower depths than assessed in the USGS catalog. These aftershocks align along two linear trends that suggest the activation of two unnamed primary faults.

58 GEOSCIENCES↗

Instantaneous tracking of earthquake growth with elastogravity signals

Rapid and reliable estimation of large earthquake magnitude (above 8) is key to mitigating the risks associated with strong shaking and tsunamis. Standard early warning systems based on seismic waves fail to rapidly estimate the size of such large earthquakes. Geodesy-based approaches provide better estimations, but are also subject to large uncertainties and latency associated with the slowness of seismic waves. Recently discovered speed-of-light prompt elastogravity signals (PEGS) have raised hopes that these limitations may be overcome, but have not been tested for operational early warning. Here we show that PEGS can be used in real time to track earthquake growth instantaneously after the event reaches a certain magnitude. We develop a deep learning model that leverages the information carried by PEGS recorded by regional broadband seismometers in Japan before the arrival of seismic waves. After training on a database of synthetic waveforms augmented with empirical noise, we show that the algorithm can instantaneously track an earthquake source time function on real data. Our model unlocks ‘true real-time’ access to the rupture evolution of large earthquakes using a portion of seismograms that is routinely treated as noise, and can be immediately transformative for tsunami early warning.

58 GEOSCIENCES↗

Cloud condensation nuclei activity of internally mixed particle populations at a remote marine free troposphere site in the North Atlantic Ocean

Here, this study reports results from research conducted at the Observatory of Mount Pico (OMP), 2225 m above mean sea level on Pico Island in the Azores archipelago in June and July 2017. We investigated the chemical composition, mixing state, and cloud condensation nuclei (CCN) activities of long-range transported free tropospheric (FT) particles. FLEXible PARTicle Lagrangian particle dispersion model (FLEXPART) simulations reveal that most air masses that arrived at the OMP during the sampling period originated in North America and were highly aged (average plume age > 10 days). We probed size-resolved chemical composition, mixing state, and hygroscopicity parameter ($\mathcal{κ}$) of individual particles using computer-controlled scanning electron microscopy with an energy-dispersive X-ray spectrometer (CCSEM-EDX). Based on the estimated individual particle mass from elemental composition, we calculated the mixing state index, $χ$. During our study, FT particle populations were internally mixed ($χ$ of samples are between 53% and 87%), owing to the long atmospheric aging time. We used data from a miniature Cloud Condensation Nucleus Counter (miniCCNC) to derive the hygroscopicity parameter, $\mathcal{κ}$ CCNC . Combining $\mathcal{κ}$ CCNC and FLEXPART, we found that air masses recirculated above the North Atlantic Ocean with lower mean altitude had higher $\mathcal{κ}$ CCNC due to the higher contribution of sea salt particles. We used CCSEM-EDX and phase state measurements to predict single-particle $\mathcal{κ}$ ($\mathcal{κ}$ CCSEM-EDX ) values, which overlap with the lower range of $\mathcal{κ}$ CCNC measured below 0.15 % SS. Therefore, CCSEM-EDX measurements can be useful in predicting the lower bound of $\mathcal{κ}$, which can be used in climate models to predict CCN activities, especially in remote locations where online CCN measurements are unavailable.

54 ENVIRONMENTAL SCIENCES↗

Machine-Learning-Based High-Resolution Earthquake Catalog Reveals How Complex Fault Structures Were Activated during the 2016–2017 Central Italy Sequence

The 2016–2017 central Italy seismic sequence occurred on an 80 km long normal-fault system. The sequence initiated with the Mw 6.0 Amatrice event on 24 August 2016, followed by the Mw 5.9 Visso event on 26 October and the Mw 6.5 Norcia event on 30 October. We analyze continuous data from a dense network of 139 seismic stations to build a high-precision catalog of ~900,000 earthquakes spanning a 1 yr period, based on arrival times derived using a deep-neural-network-based picker. Our catalog contains an order of magnitude more events than the catalog routinely produced by the local earthquake monitoring agency. Aftershock activity reveals the geometry of complex fault structures activated during the earthquake sequence and provides additional insights into the potential factors controlling the development of the largest events. Activated fault structures in the northern and southern regions appear complementary to faults activated during the 1997 Colfiorito and 2009 L’Aquila sequences, suggesting that earthquake triggering primarily occurs on critically stressed faults. Delineated major fault zones are relatively thick compared to estimated earthquake location uncertainties, and a large number of kilometer-long faults and diffuse seismicity were activated during the sequence. These properties might be related to fault age, roughness, and the complexity of inherited structures. The rich details resolvable in this catalog will facilitate continued investigation of this energetic and well-recorded earthquake sequence.

58 GEOSCIENCES↗

An Online Approach to Solve the Dynamic Vehicle Routing Problem with Stochastic Trip Requests for Paratransit Services

Many transit agencies operating paratransit and microtransit services have to respond to trip requests that arrive in real-time, which entails solving hard combinatorial and sequential decision-making problems under uncertainty. To avoid decisions that lead to significant inefficiency in the long term, vehicles should be allocated to requests by optimizing a non-myopic utility function or by batching requests together and optimizing a myopic utility function. While the former approach is typically offline, the latter can be performed online. We point out two major issues with such approaches when applied to paratransit services in practice. First, it is difficult to batch paratransit requests together as they are temporally sparse. Second, the environment in which transit agencies operate changes dynamically (e.g., traffic conditions can change over time), causing the estimates that are learned offline to become stale. To address these challenges, we propose a fully online approach to solve the dynamic vehicle routing problem (DVRP) with time windows and stochastic trip requests that is robust to changing environmental dynamics by construction. We focus on scenarios where requests are relatively sparse—our problem is motivated by applications to paratransit services. We formulate DVRP as a Markov decision process and use Monte Carlo tree search to evaluate actions for any given state. Accounting for stochastic requests while optimizing a non-myopic utility function is computationally challenging; indeed, the action space for such a problem is intractably large in practice. To tackle the large action space, we leverage the structure of the problem to design heuristics that can sample promising actions for the tree search. Our experiments using real-world data from our partner agency show that the proposed approach outperforms existing state-of-the-art approaches both in terms of performance and robustness.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Pulse-by-pulse transient thermal deformation in crystal optics under high-repetition-rate FEL

Time-domain modeling of the thermal deformation of crystal optics can help define acceptable operational ranges across the pulse-energy repetition-rate phase space. In this paper, we have studied the transient thermal deformation of a water-cooled diamond crystal for a cavity-based X-ray free-electron laser (CBXFEL), either an X-ray free-electron laser oscillator (XFELO) or a regenerative amplifier X-ray free-electron laser (RAFEL), by numerical simulations including finite-element analysis and advanced data processing. Pulse-by-pulse transient thermal deformation of a 50 µm-thick diamond crystal has been performed with X-ray pulse repetition rates between 50 kHz and 1 MHz. Results for temperature and thermal deformation have been compared with the results of transient analysis using a continuous wave (CW) power loading. Temperature and thermal deformation results from pulse-by-pulse transient analysis vary with time about the results for the CW case for the same average power. The variation amplitude increases with pulse energy and decreases with repetition rate. When the repetition rate increases to infinity, both temperature and thermal deformation converge to the results for the CW case. Two critical time scales for the operation of crystal optics in a CBXFEL are (1) first-turn time, i.e. the time for the XFEL pulse to complete the first turn around the cavity so that the crystal sees the recirculated XFEL pulse, and (2) period-end time, i.e. the time that the next electron bunch arrives for the amplification, so that the crystal outcouples the amplified FEL power. For the same average power, simulation results show that the crystal thermal deformation seen by the XFEL beam decreases with repetition rate at the first-turn time of a 300 m-long cavity and increases with repetition rate at the period-end time. For the wavefront preservation requirement of the crystal optics, a pulse-energy versus repetition-rate phase space has been established. The upper bounds of the pulse energy at both first-turn and period-end times decreases with repetition rate, especially at the period-end time. The upper bound of the thermal deformation of the crystal at the period-end time for any repetition frequency can be estimated from the CW case. For a water-cooled diamond crystal of dimension 5 mm × 5 mm × 0.05 mm, the time to reach a quasi steady-state is about 50 ms for temperature and 50 µs for thermal deformation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Vadose and Saturated Zone Flow and Transport Calculations for the Active Trenches of the Low-Level Burial Grounds, Hanford Site, Washington

The purpose of the fate and transport modeling described in this environmental calculation file (ECF) is to evaluate the impacts to groundwater associated with waste disposal operations at Mixed Waste Low-Level Burial Ground (LLBG) Trenches 31 and 34 to satisfy requirements in DOE O 435.1, Radioactive Waste Management. The model integrates the flow and transport in the vadose zone beneath the active trenches with the saturated zone downgradient of the trenches to predict the radionuclide concentration at the point of assessment (POA). DOE M 435.1-1, Radioactive Waste Management Manual, defines the POA as the point of highest projected dose or concentration beyond a 100 m (328 ft) buffer zone surrounding the disposed waste. The modeling is conducted in accordance with the DOE G 435.1, Implementation Guide for Use with DOE M 435.1-1, performance assessment (PA) guidelines. The modeling involves evaluation of the groundwater concentrations and radionuclide arrival times during the 1,000-year compliance and 10,000-year sensitivity-uncertainty periods per DOE O 435.1 and DOE M 435.1. This analysis does not consider radionuclide release during facility operations, only the post-closure impacts of the radionuclides to the environment. The evaluation of potential radiological dose to groundwater receptors caused by releases from a closed facility containing radioactive waste typically includes the following: (1) Release of radionuclides from that facility (2) Transport of those radionuclides through the environment, and (3) Exposure to humans to environmental concentration levels of those radionuclides The fate and transport three-dimensional (3D) model analysis involves the post-closure impacts to the environment of the technetium-99, iodone-129, and uranium (all isotopes in the waste). The residual inventory estimates include several radionuclides, but technetium-99 is typically responsible for almost all of the beta-gamma dose equivalent associated with groundwater (water resources) protection per 40 CFR 141, “National Primary Drinking Water Regulations” (e.g., see the results in WCH-520, Performance Assessment of Environmental Restoration Disposal Facility, Hanford Site, Washington; hereinafter referred to as the ERDF PA), and iodone-129 can also be a significant dose contributor for some waste (e.g., RPP-RPT-59958, Performance Assessment for the Integrated Disposal Facility, Hanford Site, Washington; hereinafter referred to as the IDF PA). Uranium does not typically factor significantly into the impacts to groundwater, even during the 10,000-year sensitivity-uncertainty period, but always remains of interest as a contaminant. This ECF does not address vadose and saturated zone modeling for Trench 94 of the 200 East Area LLBG. Current information confirms the validity of the low corrosion rate of the naval reactor plant carbon steel (HY-80), and the even lower corrosion rate of the nickel-iron-chromium alloy reactor vessel (Inconel Alloy 600) presented in DOE/EIS-0259, Final Environmental Impact Statement on the Disposal of Decommissioned, Defueled Cruiser, Ohio Class, and Los Angeles Class Naval Reactor Plants. Based on these low corrosion rates, the time to breach the reactor vessel to allow release of radionuclides from the activated metal of the reactor vessel internal structure is at least 10,000 years. This time to breach precludes the need to evaluate the vadose and saturated zone transport of contaminants released from the reactor compartment disposal packages in the 200 East Area LLBGs PA (CP-63826, Waste Release Model Package Report for the Active Trenches of the Low-Level Burial Grounds, Hanford Site, Washington).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Search for Infrasound Signals in InSight Data Using Coupled Pressure/Ground Deformation Methods

The unprecedented quality and sampling rate of seismometer and pressure sensors of the InSight Mars mission allow us to investigate infrasound through its pressure and ground deformation signals. This study focuses on compliance effects induced by acoustic waves propagating almost horizontally close to the surface. The compliance of acoustic waves is first estimated using the compliance estimates from pressure perturbations moving at wind speed. Then, a marker of compliance events is used to detect events of ground deformation induced by pressure variations, in three frequency bands from 0.4 to 3.2 Hz, from InSight sol 180 to 690. Additional selection criteria are imposed on the detected events to focus on acoustic waves and to remove various noise sources (e.g., wind effects or seismometer artifacts). After an automated selection, the visual inspection of the records allows us to validate two infrasound candidates that cannot be related to pressure perturbations moving at wind speed nor to known noise sources. For our highest quality infrasound candidate, the relation between this event and a convective vortex occuring 10 s later is tested. The azimuth of the vortex position at the time of infrasound detection is not consistent with the arrival azimuth of the suspected infrasound inferred from the polarization of seismometer records, thus the link between these two phenomena cannot be demonstrated. Further investigations would require a better understanding of wind-related noise impacting InSight sensors and of the effects of lateral variations of subsurface mechanical properties on the ground deformations induced by atmospheric pressure variations.

58 GEOSCIENCES↗

WHOLESCALE: Microseismic Event Catalog for San Emidio, Nevada 2022

This submission includes a high-precision seismic event catalog estimated from seismic data collected at San Emidio, Nevada from April to May 2022. The catalog lists the precise time, location, and magnitude of microseismic events recorded during this period. Both the seismic data and microseismic event catalog were produced as part of the Water & Hole Observations Leverage Effective Stress Calculations and Lessen Expenses (WHOLESCALE) project. Attached here are the microseismic event catalog, a link to the raw seismic data, and a detailed description of methods used to create the catalog.

15 GEOTHERMAL ENERGY↗

SULI Technical Report: WIRE (Draft)

Local climate undoubtedly affects facility maintenance requirements; however, knowledge of the link between weather and facility condition over time is grounded in the 1990s and needs to be updated. WIRE aims to quantify the relationship between climate and infrastructure by calculating climate-based service life estimates to be used in BUILDER, a program the Department of Energy (DOE) uses for facility management. Phase One of this project is a case study on window-or-wall-mounted AC units. The independent variables influencing service life are load, maintenance quality, and effective age, with load calculated from three weather variables: temperature, humidity, and solar radiation. Given sufficient data, service life can be calculated as a function of these variables via regression. WIRE is currently collecting data. Preliminary estimated service life projections as a function of load derived from Lawrence Livermore National Laboratory (LLNL) data support the need for updated analysis on the relationship between climate and facility condition. With the groundwork laid, WIRE can generate deliverables immediately upon arrival of data. The methods used for this case study constitute a framework for deriving the relationship between weather variables and service life for other components.

42 ENGINEERING↗

Fog Formation Related to Gravity Currents Interacting with Coastal Topography

An interesting mixing-fog event was identified during IOP7 of the C-FOG project, where a cold front arriving from the north-east collided with the Downs peninsula in Ferryland, Newfoundland, to produce foggy conditions. A comprehensive set of field observations suggests that this collision caused turbulent mixing of nearly saturated ambient air with an almost saturated cold front, creating conditions for mixing fog. Additionally, laboratory experiments were performed to study the interaction of lock-exchange induced gravity currents with a (rectangular) obstacle to delve into physical processes underlying this phenomenon. Instantaneous velocity and density fields were obtained using particle image velocimetry and planar laser-induced fluorescence. The observations suggest that the obstacle starts affecting the approaching gravity-current propagation at an upstream distance of 2 H and, upon collision, the mixing is taking place over a spatial scale of 0.83 H , where H is the depth of the ambient fluid layer. The time for largerscale turbulent stirring to permeate to the smallest scales of turbulence and activate the condensation nuclei was estimated as 3t * , where t* = √H/g' is the intrinsic time scale of the gravity current and g' the reduced gravity. Extrapolation of results to IOP7 showed a good agreement with observations.

54 ENVIRONMENTAL SCIENCES↗

Time-of-flight vs time-of-arrival in neutron spectroscopic measurements for high energy density plasmas

The neutron time-of-flight (nToF) diagnostic technique has a lengthy history in Inertial Confinement Fusion (ICF) and High Energy Density (HED) Science experiments. Its initial utility resulted from the simple relationship between the full width half maximum of the fusion peak signal in a distant detector and the burn averaged conditions of an ideal plasma producing the flux [Lehner and Pohl, Z. Phys. 207, 83–104 (1967)]. More recent precision measurements [Gatu-Johnson et al., Phys. Rev. E 94(8), 021202 (2016)] and theoretical studies [Munro, Nucl. Fusion 56, 035001 (2016)] have shown the spectrum to be more subtle and complicated, driving the desire for an absolute calibration of the spectrum to disambiguate plasma dynamics from the conditions producing thermonuclear reactions. In experiments where the neutron production history is not well measured, but the neutron signal is preceded by a concomitant flux of photons, the spectrum can be in situ calibrated using a set of collinear detectors to obtain a true “time-of-flight” measurement. This article presents the motivation and overview of this technique along with estimates of the experimental precision needed to make useful measurements in existing and future nToF systems such as the pulsed power Z-machine located in Albuquerque, NM, at Sandia National Laboratories.

Fusion experiments↗

Acoustic waves from the 2023 April 20 SpaceX Starship rocket explosion travelling in the elevated ‘AtmoSOFAR’ channel

SUMMARY The ability to detect low-frequency sounds from distant energetic events depends on the temperature and wind structure of the atmosphere. Thus, from time to time surface-based acoustic detectors may not be able to capture sounds arriving from certain directions. However, the temperature minimum at the tropopause may create an acoustic duct called the ‘AtmoSOFAR’ channel that could transmit acoustic waves laterally—but perhaps not to the ground. If true, then elevated sensors such as those borne aloft by balloons would record the signatures even in regions where ground-based sensors cannot. This has been difficult to prove because high altitude acoustic sources are rare and balloon deployments are sporadic. This work describes the detection and characterization of powerful acoustic waves generated during the launch and terminal explosion of the SpaceX Starship rocket on 2023 April 20 using a pair of microbarometers on a stratospheric balloon. The signals travelled through the AtmoSOFAR channel, carrying information about the size and nature of their source. This channel also appears to leak some acoustic energy to the ground, in agreement with previous studies. The acoustic yield of the Starship terminal explosion was on the order of $10^3$ tons trinitrotoluene equivalent, which agrees with the estimated fuel load of the vehicle to about a factor of 2, but is two orders of magnitude larger than optical estimates. These results support an earlier study that claimed lateral transmission of sound from a smaller rocket through the AtmoSOFAR channel. The transmission of source information through the AtmoSOFAR channel motivates its use for monitoring other natural and anthropogenic events using balloon-borne sensors. This may become increasingly important as more and more private and government entities conduct spacecraft launches and reentries. It may also provide a means of monitoring clear air turbulence and other sound-generating atmospheric phenomena at a distance.

58 GEOSCIENCES↗

Examining Infrasound Propagation at High Spatial Resolution Using a Nodal Seismic Array

Infrasound—acoustic waves in the atmosphere below 20 Hz—is a useful monitoring tool. Topography and atmospheric structure strongly control infrasound propagation, and at common source–receiver distances neither of these effects can be ignored when quantitative source constraints are sought. Detailed spatial measurements of the infrasound wavefield would inform propagation models and improve source estimates. However, the “large-N” deployment strategy now well-known in seismology has not yet been realized for infrasound studies. Here, we use the 900-node seismic array from the 2014 Imaging Magma Under St. Helens (iMUSH) experiment as a proxy for a large-N infrasound network, by leveraging acoustic–seismic coupled arrivals. The active-source component of iMUSH consisted of 23 shallowly buried explosions around Mount Saint Helens volcano; these explosions produced epicentral infrasound recorded on the nodes. We find that the bulk presence of ground-coupled infrasound on the nodes is controlled by wind noise and source–receiver distance, with observed arrivals for eight explosions. Explosions with the most extensive coupling produce complex spatial waveform patterns across the array. These patterns are related to both topographic and atmospheric propagation effects, as well as spatially variable site (coupling) effects. We compare our observations to simple topographic diffraction and high-resolution wind advection models, and full-wave numerical simulations. We find strong spatial correlations between (a) coupled arrival strength and modeled topographic obstruction and (b) coupled arrival time and along-path winds. Our seismoacoustic analyses and results are applicable to other existing and future nodal seismic data sets and can expand the utility of such deployments.

58 GEOSCIENCES↗

A Census of the 32 Ori Association with Gaia*

I have used high-precision photometry and astrometry from the third data release of Gaia (DR3) to identify candidate members of the 32 Ori association. Spectral types and radial velocities have been measured for subsets of the candidates using new and archival spectra. For the candidates that have radial velocity measurements, I have used UVW velocities to further constrain their membership, arriving at a final catalog of 169 candidates. I estimate that the completeness of the survey is ~90% for spectral types of ≲M7 (≳0.06 M ⊙ ). The histogram of spectral types for the 32 Ori candidates exhibits a maximum at M5 (~0.15 M ⊙ ), resembling the distributions measured for other young clusters and associations in the solar neighborhood. The available UVW velocities indicate that the association is expanding, but they do not produce a well-defined kinematic age. Based on their sequences of low-mass stars in color–magnitude diagrams, the 32 Ori association and Upper Centaurus-Lupus/Lower Centaurus-Crux (UCL/LCC) are coeval to within ±1.2 Myr, and they are younger than the β Pic moving group by ~3 Myr, which agrees with results from previous analysis based on the second data release of Gaia. Finally, I have used mid-IR photometry from the Wide-field Infrared Survey Explorer to check for excess emission from circumstellar disks among the 32 Ori candidates. Disks are detected for 18 candidates, half of which are reported for the first time in this work. The fraction of candidates at ≤M6 that have full, transitional, or evolved disks is 10/149 = 0.07$^{+0.03}_{-0.02'}$ which is consistent with the value for UCL/LCC.

79 ASTRONOMY AND ASTROPHYSICS↗

Time-of-Flight Estimation for Nondestructive Evaluation

Ultrasonic testing (UT) for nondestructive evaluation (NDE) is a critical entity necessary to resolve both the quality and precision questions of complex parts evolving from the classical approaches to the innovative additive manufacturing (AM) process. This modality provides the essential quantitative information for acceptance and potential flaw detectionof apart under investigation. A primary ingredient in UT besides the required precision robotic hardware for theacquisition of high-qualitymeasurement data is the underlying signal processing. It is here that much of the system performance capability resides. In this report,we discuss the basic steps in UT signal processing along with current capabilities that must be achieved in order to satisfy the criticaldemands of the various LLNL programs.In itsbasic foundations, signal processing is essentially the “extraction of critical information (signals) from noisy, uncertain data.” Clearly, in a perfect world the best signal processing is none---just make a reliable, noise-free, uncertainty-free, measurement. Unfortunately, even the best of systems is still straddled within the confines and limitations of physical instrumentation. With this in mind, we discuss the development of signal processing techniques toextract the desired UT information from noisy measurement data. We start with a discussion of a simple homogeneous representation of a “part” under investigation and its insonification by an ultrasonicexcitation from a physics-based perspective. Once developed, we briefly discussthe various choices of excitation signals and their tradeoffs.Next,we discuss a simulation approach employing simple models from the signal processing perspective and then move on to the development of the basic signal processing approach to ultrasonic signal processing all based on estimating the pulse arrival estimation. Starting with simple peak detection techniques, progressing to the processing workhorse so-called “matched-filter” and finally to the more sophisticated “model-based matched filter.”The underlying pre-and post-processing of noisy measurement data along with the application of these techniques are subsequently demonstrated on experimental ultrasonic data.Finally, we present some results for processing “weld” data when the ultrasonic excitation signal is not directly available.

42 ENGINEERING↗

Evaluation of a generalized least squares algorithm for infrasound beamforming with coherent background noise

Infrasonic signals of interest can occur during periods with persistent, coherent, background noise, which may be natural or anthropogenic. For high signal-to-noise (SNR) ratio transient signals, an ‘overprinting’ of the coherent background may occur, and the signal may still be detected. However, this approach fails for low SNR signals of interest, which may be obscured by coherent noise. An infrasound beamforming method based on generalized least squares (GLS) is investigated for detecting transient signals of interest in the presence of coherent and incoherent background noise. This approach relies on an estimate of the noise covariance, captured in a covariance matrix, to effectively null contributions to the array response from noisy directions of arrival. Synthetic array data is used to investigate the performance of the GLS beamformer compared to the Bartlett beamformer when coherent and incoherent backgrounds are present. Additionally, the effects of array element number and relative strength of the interfering signal on the GLS estimates is investigated. GLS empirical area under the curve estimates suggest that the beamformer can recover coherent power for a signal of interest lower in amplitude than the coherent background, but this effectiveness degrades more quickly with SNR for a four element array compared to a six or eight element infrasound array. Finally, infrasound from the Forensic Surface Experiment, a bolide signal observed at IMS array I37NO, and a volcanic signal recorded at the Alaska Volcano Observatory array ADKI are used to evaluate GLS performance on recorded data. A ten minute window was used to capture the background noise, and the coherent background signal was nulled in all three examples.

58 GEOSCIENCES↗